Archaean Plate Tectonics in the North Atlantic Craton of West Greenland Revealed by Well-Exposed Horizontal Crustal Tectonics, Island Arcs and Tonalite-Trondhjemite-Granodiorite Complexes

The 700 km-long North Atlantic Craton (NAC) in West Greenland is arguably the best exposed and most continuous section of Eo-to Neoarchaean crust on Earth. This allows a close and essential correlation between geochemical and isotopic data and primary, well-defined and well-studied geological relationships. The NAC is therefore an excellent and unsurpassed stage for the ongoing controversial discussion about uniformitarian versus non-uniformitarian crustal evolution in the Archaean. The latest research on the geochemistry, structural style, and Hf isotope geochemistry of tonalite-trondhjemite-granodiorite (TTG) complexes and their intercalated mafic to intermediate volcanic belts strongly supports previous conclusions that the NAC formed by modern-style plate tectonic processes with slab melting of wet basaltic oceanic crust in island arcs and active continental margins. New studies of the lateral tectonic convergence and collision between juvenile belts in the NAC corroborate this interpretation. Nevertheless, it has repeatedly been hypothesised that the Earth’s crust did not develop by modern-style, subhorizontal plate tectonics before 3.0 Ga, but by vertical processes such as crustal sinking and sagduction, and granitic diapirism with associated dome-and-keel structures. Many of these models are based on supposed inverted crustal density relations, with upper Archaean crust dominated by heavy mafic ridge-lavas and island arcs, and lower Archaean crust mostly consisting of felsic, supposedly buoyant TTGs. Some of them stem from older investigations of upper-crustal Archaean greenstone belts particularly in the Dharwar craton, the Slave and Superior provinces and the Barberton belt. These interpreted interactions between these upper and lower crustal rocks are based on the apparent down-dragged greenstone belts that wrap around diapiric granites. However, in the lower crustal section of the NAC, there is no evidence of any low-density granitic diapirs or heavy, downsagged or sagducted greenstone belts. Instead, the NAC contains well-exposed belts of upper crustal, arc-dominant greenstone belts imbricated and intercalated by well-defined thrusts with the protoliths of the now high-grade TTG gneisses, followed by crustal shortening mainly by folding. This shows us that the upper and lower Archaean crustal components did not interact by vertical diapirism, but by subhorizontal inter-thrusting and folding in an ambient, mainly convergent plate tectonic regime.


INTRODUCTION
The Archaean craton in West Greenland is a key part of the North Atlantic Craton (NAC) (Bridgwater et al., 1973a, Bridgwater et al., 1973bNutman, 1997). The excellent three-dimensional exposures along this N-S, 700 km-long craton makes this one of the most informative transects of Archaean crust on Earth, especially because it contains many intercalated sections of both upper-and lowercrustal rocks and their connected structures that are rarely available in granite-greenstone terrains. It was first realised by Windley and Bridgwater (1971) that high-grade Precambrian gneiss belts primarily consist of strongly deformed, mostly TTG (tonalitetrondhjemite-granodiorite) orthogneisses, and rarely include rocks of sedimentary origin as normally seen in granite-greenstone belts. The distinction between Eo-and Mesoarchaean orthogneisses in the Nuuk region (McGregor, 1973) led to the recognition of large-scale, subhorizontal, tectonic intercalation of older and younger crustal segments (Bridgwater et al., 1974), and later to the concept of Archaean tectono-stratigraphy (e.g., Friend et al., 1988;Nutman et al., 1989;Friend and Nutman, 2019) where individual terranes with specific lithologies, ages and tectonic histories were recognised within the central NAC, and later found to be applicable to the Canadian Superior Province (e.g., Polat et al., 2009).
Subsequently, the Greenlandic NAC continued to yield abundant field, geochemical and geochronological evidence of Archaean crustal accretion/evolution by plate tectonic processes, including identification of andesitic calc-alkaline metavolcanic rocks in supra-subduction zone island arcs intercalated with TTG gneisses. Following the terrane accretion models of Friend et al. (1987), Friend et al. (1988), McGregor et al. (1991) made the first specific Archaean plate tectonic model for the Godthåbsfjord region. These and other papers with relevant data were summarised by Windley and Garde (2009) who reported strong evidence for the operation of plate tectonic processes in the Meso-to Neoarchaean of West Greenland. In the following we particularly examine more recent research, except for the NAC in East Greenland which still lacks many detailed studies. We finally address if non-uniformitarian models of Archaean crustal evolution might be viable alternatives to the accumulated field and geochemical data from the large, wellexposed, representative rocks in West Greenland.

Components of the Eoarchaean Itsaq Gneiss Complex
The largely Mesoarchaean craton of West Greenland contains the extensively studied, NE-trending Eoarchaean Itsaq Gneiss Complex. McGregor (1973) made his outstanding breakthrough in Archaean history at the southern end of this belt, around the town of Nuuk (formerly Godthåb). The 3.8-3.6 Ga Eoarchaean belt, which is about 180 km long (Nutman et al., 1996;Nutman et al., 2013), consists primarily of TTG orthogneisses, within which there are many relict layers, lenses and large bodies of a variety of rocks, particularly in the Isua supracrustal belt in the north and on Akilia island and nearby areas in the south (e.g., Nutman et al., 2020;van de Löcht et al., 2020). Eoarchaean orthogneisses also occur in the Aasivik terrane north of the Itsaq Gneiss Complex (Figure 1; Rosing et al., 2001).
The existence of surface water at 3.8-3.7 Ga is indicated in Isua by the presence of pillow lavas (Maruyama et al., 1991;Appel et al., 1998;Maruyama and Komiya, 2011) and of bedded sediments such as pelagic cherts and banded iron formation (BIF) (Moorbath et al., 1973;Nutman et al., 1997b). According to Maruyama and Komiya (2011) the total thickness of the largest unit of pillow basalts in Isua is up to 1.0 km, suggesting that the ambient ocean depth was at least 1.0 km; from the presence of vesicular basalts de Wit and Furnes (2013) suggested a depth of up to 2-4 km. By using the source mantle potential temperature of 1,450°C and the ocean surface temperature of 100°C in the Early Archaean by Komiya (2004), Komiya (2007) and Komiya et al. (1999), which compare with the corresponding temperatures of the modern Earth of 1,330 and 0°C, respectively, Maruyama and Komiya (2011) calculated the oceanic geotherm of the Early Archaean lithosphere and concluded that the lithosphere was sufficiently rigid to operate modern-style plate tectonics in Isua times. Frei and Polat (2007) used isotopic relations in the Isua BIF to demonstrate that two rare-earth element (REE) sources were derived from sea-floor vented hydrothermal fluids and from ambient surface seawater that were responsible for the iron-and silica-rich layers, implying that the silica was derived from a now-unexposed mafic landmass. Friend et al. (2008) demonstrated that the Isua BIF has modern seawater-like trace element signatures consistent with chemical sedimentation from seawater with little clastic input at c. 3.7 Ga. In a more detailed geochemical study of the Isua BIF Aoki et al. (2018) showed that with stratigraphic height the ratio of hydrothermal fluid to seawater component decreased throughout the period of deposition with increasing distance from a ridge to a trench in common with the changes in modern sediments on a spreading oceanic plate, confirming the role of ocean plate stratigraphy (OPS, see below) in the plate-tectonic history of Isua.

Lithological, Structural, Geochronological and Geochemical Models for the Isua Supracrustal Belt
The tectonic significance of the Isua belt revolves around three main conceptual models of how to interpret the complex interrelationships of an orogenic belt. The lithology/structure-based, the chronology/structure-based, and the geochemicallyconstrained model. Each model takes account of field-based structures, albeit in different ways, and they all interpret the Isua belt in terms of modern-style plate tectonics (for a different view see Webb et al., 2020).

The Lithology/Structure-Based Model
The first major model to explain the Isua geology came from the Maruyama-Komiya team, who with their experience of the Mesozoic OPS in Japan discovered essentially identical OPS in NE Isua (Maruyama et al., 1991;Komiya et al., 1999) followed by Maruyama and Komiya (2011) and Arai et al. (2015). Ocean plate stratigraphy (OPS) is a bulk term that describes the travel history of an opening-closing oceanic plate from basalts, gabbros and ultramafic rocks that form at a mid-oceanic ridge to pelagic cherts and/or carbonates that accumulate during the movement of the ocean toward its trench. When the outer trench slope is reached hemipelagic mudstones are deposited on the cherts, and succeeded by voluminous clastic turbidites, greywackes and conglomerates that may be derived by erosion of the accretionary wedge, the adjacent volcanic arc, or an eroded old craton behind the arc (Kusky et al., 2013). The key point relevant for this discussion is that the sequentially developed package of basalt-chert-turbidite can be recognised in many accretionary orogens from the Eoarchaean to the Cenozoic worldwide. As all the maps of any origin show, this tripartite package is clearly present in NE Isua (but see Geochemically-Constrained Models regarding basalt compositions). It represents the core of the argument that Isua developed as an accretionary orogen little different from those of the Phanerozoic.

The Chronology/Lithology-Based Model
The second major model to explain Isua geology comes from the Nutman-Friend team (Nutman et al., 1996;Nutman et al., 1997a;Nutman et al., 1997b;Nutman et al., 2004;Nutman et al., 2010;Nutman et al., 2013;Nutman 2015a;Nutman et al., 2015b;Friend and Nutman, 2019), who used abundant zircon ages mostly of (meta-) tonalites and diorites to work out the sequential development of two main terranes in NE Isua; a southern c. 3.8 Ga and a northern c. 3.7 Ga terrane. The central boundary between the terranes is a mylonitised belt of chert, BIF and carbonate rocks that contain 3,940-3,750 Ma detrital zircons ). The northern, younger terrane is an imbricate package with no preserved stratigraphy, composed of arc components that young from c. 3,720-3,715 Ma boninites (Polat et al., 2002), tholeiites and picrites (Polat et al., 2002), to c. 3,710-3,700 Ma andesites, lesser dacites and sediments derived from them, to c. 3,700-3,695 Ma BIF and minor zircon-bearing siliceous detritus from a distal source ). The older southern terrane contains a similar arc assemblage and is similar in origin to the northern terrane, but it contains minor >3,850 Ma crust. This (now-thin) central boundary belt probably acted as a décollement when the northern younger, c. 3.7 Ga, terrane was thrust southwards over the southern older, c. 3.8 Ga, terrane.
Although Nutman and considered that Isua developed by accretionary processes at a convergent plate margin similar to those operating in the modern Earth, their tectonic development by upward younging of the accretionary belt is the opposite of the younging downward polarity recognised in the accretionary wedges of modern subduction zones, and in accretionary orogens of any age (e.g., Leggett et al., 1979;Platt 1986;Kawai et al., 2007;Isozaki et al., 2010;Kusky et al., 2013). underlain by massive and pillow basalts, representing the whole OPS history of Isua. Because the early traverses of Komiya et al. (1999) crossed some of the later isotopically-defined units and boundaries of , and because Polat et al. (2002) and  found that all the basalts have island arc chemistry, inevitably  do not believe the OPS units and structural set-up of Komiya et al. (1999), but consider that their three zones developed separately and at different times and were juxtaposed by accretionary collisions. Kamber et al. (2001) made a critical assessment of the isotopic constraints on the evolution of Isua rocks. However, we will not comment further, as it is beyond the scope of this paper to present a comparative detailed analysis of these opposing geochemical arguments.

The Mesoarchaean Sections of the Greenland Craton
A plate-tectonic model for the central part of the West Greenland NAC was first proposed by McGregor et al. (1991). In the entire NAC of West Greenland Windley and Garde (2009) proposed a repetitive pattern of crustal tilting into six large blocks with tectonic boundaries (Ivittuut, Kvanefjord, Bjørnesund, Sermilik, Fiskefjord and Maniitsoq, Figure 1). Each block except Maniitsoq comprises a transition from granulite facies, lower zones in the north to unretrogressed amphibolite facies upper zones in the south, with intervening rocks partially retrogressed to upper amphibolite facies. The block boundaries were originally drawn to bring out this general pattern. Some have been revised in Figure 1, and further revisions may be required as new information from less intensely studied regions becomes available.
The lower, deep-crustal zones of the blocks are dominated by highly ductile, double/triple fold interference patterns that overprint earlier thrusts and recumbent isoclinal folds. They comprise TTG gneisses with subordinate, kilometre-sized belts, small rafts and enclaves of homogeneous to indistinctly layered amphibolites that rarely preserve primary structures. Like mafic counterparts in the upper zones they have MORB-like compositions that contain minor but important layers of anorthositic complexes, like Fiskenaesset ( Figure 1), and layers and enclaves of amphibolites, which like their mafic counterparts in the upper zones have MORB-like compositions and flat chondrite-normalised REE patterns.
The upper zones display simpler, less ductile structures and contain better-preserved intrusive contacts and thrustintercalated volcanic and plutonic components. They contain well-defined remnants of volcanic arcs with recognisable pillow basalts, lapilli tuffs and andesitic volcaniclastic rocks. The arcs are bimodal in composition, with tholeiitic mafic and calc-alkaline andesitic components (Table 1; Testing Mesoarchaean Plate Tectonics in West Greenland: Recent Field and Laboratory Results; Szilas et al., 2018). Many of these metavolcanic belts are associated with contemporaneous, intrusive, layered anorthositic complexes believed to have crystallised from hydrous tholeiitic magmas. The Fiskenaesset complex in the Bjørnesund block (Figure 1) is the most extensive, best exposed and most well-studied Mesoarchaean layered anorthositic complex worldwide (see Fiskenaesset Complex).
Most parts of the craton also contain I-type granitic rocks, which are contemporaneous with, or younger than the TTG complexes, and generally interpreted as melts of mafic crust (Rapp et al., 2010;Friend and Nutman, 2019).

Tonalite-Trondhjemite-Granodiorite Orthogneisses
Mesoarchaean TTG orthogneisses make up about 80% of the West Greenland NAC (Kalsbeek and Myers, 1973;Wedepohl et al., 1991;Friend and Nutman, 2019). Much earlier work on the TTG orthogneisses addressed general magmatic and tectonic aspects, demonstrating that they are invariably younger than and intrusive into the metavolcanic belts, and/or they are mutually thrust-intercalated (Myers, 1976;Kisters et al., 2012;Keulen et al., 2014;Polat et al., 2016;Szilas et al, 2018). Early thrusts and recumbent isoclinal folds were refolded by several phases of ductile folds in the lower crust and by simpler, upright compressional folds in the middle-upper crust (e.g., Berthelsen, 1960;Garde, 1997;Windley and Garde, 2009). The TTG protoliths were intruded as subhorizontal sheets synchronous with, but mainly following, the emplacement of the layered anorthositic complexes (Myers, 1976), which enabled the widespread breakup of the complexes and their host metavolcanic rocks into inclusions in the TTG gneisses, most spectacularly seen in the Ivittuut region (Berthelsen and Henriksen, 1975;Windley and Garde, 2009). We note that exactly the same process occurs in the upper crust of modern accretionary orogens, like Japan where peraluminous tonalites were intruded along thrust planes during subhorizontal tectonics and granulite facies metamorphism in the Tertiary Hidaka arc-trench system (Komatsu et al., 1989).
The major and trace element compositions of the Greenlandic TTG orthogneisses repeatedly indicate a characteristic depletion in Nb, Ta and heavy rare-earth elements, which is consistent with slab melting in a subduction zone with local mantle wedge interaction (Steenfelt et al., 2005). Experimental studies have commonly demonstrated that the magmatic generation of such TTG rocks requires a component of slab melting from hydrated mafic crust and a melt component from previously depleted mantle (Arth and Hanson, 1972;Winther and Newton, 1991;Rapp et al., 2010;Moyen and Martin, 2012). This important observation supports Mesoarchaean plate subduction, which might indicate a high potential mantle temperature in the early Earth (Komiya, 2007;Herzberg et al., 2010), or alternatively it might mean that warmer, ridge-subducted, crust was being subducted, as in the modern Earth (Arndt, 2013). Gardiner et al. (2019) presented major, trace element, isotope and geochronological data for a large number of orthogneisses in the northern part of the Finnefjeld block, which are similar to those in other parts of the NAC. They also compared the average composition of the TTG gneisses surrounding the Finnefjeld domain with that of the rocks within this domain, which had previously been shown to consist of cataclastic rocks comprising both TTG and metavolcanic rocks (Garde et al., 2014; see The Maniitsoq Structure). The Finnefjeld rocks have ∼3,000 Ma zircon U-Pb ages which are identical to those in the surrounding orthogneisses, but lower SiO 2 and higher TiO 2 , FeOt, MgO and CaO, as well as significantly higher HREE.
Abundant additional age data from the TTG gneisses in the different blocks have been published in the last decades, too numerous to be reviewed here; the most important age groups in the Nuuk region can be gleaned from Examination of Continental Growth Using Trace Element and Hafnium Isotope Compositions, see also Friend and Nutman (2019). Significant in the present context are the relative ages of adjacent metavolcanic belts (older) and TTG orthogneisses (younger).

Examination of Continental Growth Using Trace Element and Hafnium Isotope Compositions
The development of more precise trace element, geochemical and isotope analytical tools and decreasing funding for systematic field work changed the nature of recent studies in the Earth sciences and led to new tests of older crustal accretion models.
A Nb/Yb vs. Th/Yb diagram ( Figure 2) of the metabasaltic rocks from each of the Meso-to Neoarchaean supracrustal belts in the West Greenland NAC ( Table 1) illustrates potential relationships between their basaltic source magmas and Archaean continental crust as explained in the figure caption. See also the generalised Figure 10 in Szilas (2018). The element ratios straddle the present-day MORB-OIB array, forming diagonal trends that extend toward an approximate Archaean continental crustal composition. We interpret this to reflect interaction between uncontaminated magmas along the oceanic basalt array and crustal compositions above it. The ranges of the measured ratios indicate wide variations in the crustal input. There is also a remarkable age and geographical coupling: the Graedefjord, Ikkattup Nunaa and Nigerlikasik belts in the south generally have younger magmatic ages (see Andesitic Metavolcanic Arcs and Gold Mineralisation) and higher Th/Nb ratios than the northern belts, suggesting more pronounced crustal material in the former belts; the southern belts may comprise inboard, continental arc systems, with outboard, more oceanic arc systems in the northern belts.
In the last decade several studies of zircon Hf isotope data from the West Greenland NAC have been published. Hiess et al. (2009) showed that metaplutonic and metavolcanic rocks from the 3.85-3.69 Ga Itsaq Gneiss Complex (see Main Components of the North Atlantic Craton in West Greenland), all have chondritic εHf values, i.e., close to zero within analytical uncertainty, and argued that processes of crustal growth must have operated incrementally during the Hadean and Eoarchaean. Naeraa et al. (2012) presented complementary Hf isotope data from the Nuuk region including zircon in TTG gneisses and detrital zircons in rocks of Eo-to Neoarchaean age, thus providing a more complete view of the TTG melt source evolution ( Figure 3). Their Hf isotope dataset replicated the chondritic Eoarchaean signature of Hiess et al. (2009), but also revealed a broad trend of declining ε Hf values from 0 at ∼3.85 Ga to −8 at 3.25 Ga, with a distinct lack of depleted Hf isotope signatures (i.e., high ε Hf values) between c. 3.6-3.25 Ga. Naeraa et al. (2012) noted there are no known exposures of felsic crust and only very few detrital zircons in this region within the latter time interval, and thus coupled the lack of depleted Hf isotope signatures to a paucity in juvenile mantle input to the crust during this time. Although some of the detrital grains plotting along the declining ε Hf trend might reflect post-crystallisation disturbance (see, e.g., Guitreau et al., 2019), the enriched anchor point, which is represented by Fe-rich gneisses from the inner part of Ameralik Fjord, is unquestionably robust. The isotopically enriched Fe-rich gneisses clearly demonstrate that a ∼3.85 Ga old mafic protolith was available at 3.25 Ga as a source for TTG magmatism  and much later for the 2.55 Ga Qôrqut granite complex. In contrast, the zircons in TTG gneisses that formed after 3,250 Ma have highly variable ε Hf compositions ranging intermittently from supra-to highly sub-chondritic values, which indicates renewed influx of juvenile material derived from a depleted mantle source along with contributions from evolved Eoarchaean crust with negative ε Hf compositions. Naeraa et al. (2012) pointed out a clear coupling between the existence of volcanic supracrustal belts and the formation of depleted TTGs and noted that this pattern of periodic juvenile arc accretion and TTG flare-ups is more consistent with Mesoarchaean continental crust evolving in accretionary orogens than by amalgamation of truly exotic terranes. Overall, Naeraa et al. (2012) interpreted their data to suggest that continental crust grew from 3.2 Ga ago in a geodynamic regime analogous to that of modern plate tectonics.
In another study of zircon ε Hf compositions in TTG gneisses from the northern Fiskefjord block with a narrower age range, Gardiner et al. (2019) made similar arguments. As shown in Figure 3, their zircon ε Hf data between 3.25 and 2.95 Ga overlap with those of Naeraa et al. (2012), whereas younger, remelted granitic rocks have more negative ε Hf compositions consistent with reworking of old crust to form these rocks.
Several other studies in West Greenland have constrained the sources of metavolcanic belts using bulk-rock Lu-Hf and Sm-Nd isotope data combined with best estimates of their zircon U-Pb magmatic ages (e.g., Hoffmann et al., 2011;Szilas et al., 2015;Szilas et al., 2016a;Szilas et al., 2017). In general, the initial ε Hf compositions of the mafic and andesitic rocks corroborate the interpretation of Naeraa et al. (2012) that the Mesoarchaean crust in West Greenland was primarily formed by melt extraction from a depleted mantle source rather than by recycling of older crust, and that the underlying magmatic processes were both variable and complex (see also Friend et al., 2009). Some whole-rock isotopic data even suggest extraction from sources with ε Hf compositions above that of depleted mantle ( Figure 3). Studying mafic and andesitic rocks of the Qussuk-Bjørneøen metavolcanic belt in the south-eastern Fiskefjord block, Szilas et al. (2017) combined whole-rock Lu-Hf isotope data with trace element compositions and showed that the andesitic rocks have  Table 1, suggesting input of continental crust (CC) into the precursor magmas in agreement with field evidence of arc magmatism (see main text). In this diagram (originally developed to investigate crustal influence on basaltic rocks), uncontaminated modern MORB and OIB samples define an oceanic basalt array that reflects variable amounts of recycled oceanic crust in their mantle sources. Crustal rocks, as well as basaltic magmas affected by partial crustal melting or crustal contamination, have higher Th/Nb ratios and plot above this array. Their relative Th/Nb enrichment reflects a stronger Nb than Th retention in the subducting slab. All modern arc-related basalts plot above the oceanic basalt array, forming a high-Th/Nb volcanic arc array which univocally signifies a crustal influence on these subduction-related basalts, either through recycling detrital sediment into the mantle source region via subduction or subsequent assimilation in the crust. The diagram can also portray relative depletion and enrichment of basalts, as more depleted basalts have lower Nb/Yb and Th/Yb ratios. Notably, a crustal input through a subduction component will affect the Nb/Yb and Th/Yb ratios of depleted magmas more strongly than in enriched magmas, due to their much lower incompatible element content. This means that for any given amount of crustal addition to a mantle source, a depleted magma will be shifted to a point significantly higher above the oceanic basalt array than an enriched magma.
Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 540997 enriched incompatible trace element compositions but nearchondritic initial ε Hf compositions, whereas the mafic rocks have more depleted ε Hf compositions about +4 and lessenriched trace element compositions. They concluded that the andesitic metavolcanic rocks could only have formed by mixing of contemporaneous mafic and felsic magmas. This is consistent with the most widely accepted general models for the generation of modern andesitic arcs (e.g., Isozaki et al., 2010;Stern, 2010).
In summary, the large variations in whole-rock and zircon ε Hf compositions of juvenile Mesoarchaean volcanic belts and TTG gneisses clearly show that the generation of new Mesoarchaean crust was not simply controlled by direct extraction of melts from an increasingly depleted mantle source, but also involved incorporation of older crustal material, pointing to contamination and/or mixing of melts from different sources. The observation that the periods with the most abundant TTG  Gardiner et al. (2019). Best estimates of the magmatic ages of the respective supracrustal belts are shown with green arrows. The combined measurement of Lu-Hf and U-Pb isotope compositions in zircon have proven a valuable tool for extracting information about the timing of mineral growth and the radiogenic enrichment of the source from which the zircon grew. While the U-Pb system dates the timing of crystallisation or recrystallisation of the zircon Hf isotope composition expresses the degree to which the host melt was derived from a radiogenic reservoir (i.e., positive ε Hf values; e.g., depleted mantle) or an unradiogenic reservoir (i.e., negative ε Hf values; e.g., ancient continental crust), or some mixture of these sources (Spencer et al., 2020). (A) The ε Hf zircon data from TTG rocks suggest episodic growth of the continental crust (TTG formation) from a depleted mantle source at specific times, marked by the yellow bars. Between ∼3.65 and 3.25 Ga there are no known exposures of felsic crust in the NAC and very few detrital zircons. The shift toward higher (depleted) ε Hf values after ∼3.25 Ga indicates a change toward more pristine sources for TTG magmatism, which in turn suggests a change in the crustal evolution toward a more modern style of plate tectonics. (B) The whole-rock ε Hf values in the metavolcanic belts show a surprisingly large range, indicating different, variously depleted mantle sources and, in some rocks, admixture of crustal components with negative ε Hf values. Note that the highest ε Hf values require a more extreme depleted mantle evolution line than that shown in diagram (A). magmatism are associated with negative excursions in ε Hf suggests increasing degrees of contemporaneous crustal reworking .

Andesitic Metavolcanic Arcs and Gold Mineralisation
As noted in the introduction, andesitic metavolcanic arcs occur in the upper part of all of the crustal blocks identified by Windley and Garde (2009);Szilas et al. (2018) recently compiled an overview of their geochemical characteristics. Properties of the best-known arcs are shown in Table 1. The Greenlandic arcs also contain some of the oldest known epithermal mineralisation, including gold mineralisation and hydrothermal alteration in the Ivisaartoq metavolcanic belt (Polat et al., 2011), the Qussuk-Bjørneøen metavolcanic belt (Garde et al., 2012a), the Bjørnesund Supracrustal belt, the Godthåbsfjord-Ameralik belt, and the Tartoq Group, all described below. Epithermal gold mineralisation is an important, diagnostic constituent of modern island arcs. In high-grade Archaean terrains it can be difficult to distinguish between metamorphosed, deformed epithermal mineralisation as at Qussuk-Bjørneøen, and orogenic gold. For instance, in the Tartoq Group, a recent study of the initial gold mineralisation associated with hydrothermal arsenopyrite dated at 3.18-3.13 Ga is not "orogenic gold" but part of the arc-related magmatism (Saintilan et al., 2020). Unspecified gold mineralisation also occurs in the Eoarchaean Isua supracrustal belt, particularly in the areas most strongly affected by later tectono-metamorphic events (Olsen and Grahl-Madsen, 1994).

Qussuk-Bjørneøen Metavolcanic Belt
The arc around 64°30′N in Godthåbsfjord ( Figure 1) hosts a large epithermal gold occurrence (Garde et al., 2012a). The Qussuk arc contains a folded and metamorphosed, shallow-crustal argillic to advanced argillic alteration zone with a minimum strike length of 20 km, embedded within the core of a less altered, bimodal andesitic-mafic metavolcanic belt. The leached and metamorphosed varieties of the andesite precursors contain abundant biotite, garnet and sillimanite, but very little or no feldspar. These rocks were originally mapped as pelitic metasediment (see also Schumacher et al., 2011). However, anomalously high contents of immobile trace elements including REEs, negative Eu anomalies (related to leaching of plagioclase), and anomalous Ga/Al ratios show that up to 80% of the precursor rock was lost by leaching. The geochemical signature is best explained by near-surface boiling and dissociation of hydrothermal fluids into a barren, highly acidic phase and an enriched brine, as found in many modern circum-Pacific arcs (Simmons et al., 2005). Garde et al. (2012a) recorded in these altered rocks, zoned zircon grains, which contain igneous cores with high Th/U ratios typical of primary volcanic zircons, yielding ages up to c. 3,070 Ma ( Figure 4); the oldest nearby TTG rocks are c. 3,060 Ma in age (Garde et al., 2000). The surrounding zircon shells are strongly enriched in U and Th, with REE patterns typical of hydrothermal zircons. Thin metamorphic rims yield ages of c. 3,000 Ma, contemporaneous with a regional thermal event recorded in most of the Akia terrane (see below). Schlatter and Christensen (2010) had previously shown that the intensity of alteration in the gold-mineralised rocks is highly variable, as is commonly the case in Phanerozoic epithermal gold deposits. The latter authors preferred an interpretation as "orogenic gold" related to the Palaeoproterozoic Ataneq fault 10 km away (Thøgersen et al., 2006;Garde et al., 2012a), although this is barren of gold.

Storø Supracrustal Belt
The Storø supracrustal belt is part of a complex zone of thrustimbricated and folded slices of supra-and infracrustal rocks of Eo-to Neoarchaean age (Friend and Nutman 2019 and older references therein) at the northern margin of the Godthåbsfjord-Ameralik belt at around 64°30′N ( Figure 1). It includes TTG orthogneiss, gabbro and anorthosite, mafic amphibolite, metasedimentary rocks, and altered garnet-rich lithologies. Dispersed gold mineralisation associated with quartz veins is hosted by supracrustal amphibolite in central Storø, in the vicinity of chemically altered rocks of volcanic origin Szilas and Garde, 2013;Szilas et al., 2014, Szilas, 2016b. In spite of the geological complexity, it is widely agreed that the alteration of the aluminous and garnet-rich rocks was pre-metamorphic and that the gold was mobilised into its present location in the core of a late antiform (Szilas et al., 2016b). Arsenopyrite Re-Os and previously published zircon U-Pb age data led Scherstén et al. (2012) to propose that the gold mineralisation formed at 2.71-2.64 Ga, i.e. somewhat younger than its host rocks. Szilas et al. (2016b) argued on geochemical grounds that the altered and subsequently metamorphosed rocks represent transported and weathered sediment derived from mafic volcanic rocks and that these rocks were unrelated to the gold mineralisation, although field observations suggest that the alteration took place in situ in volcanic rocks (Szilas and Garde, 2013). Irrespective of the source, mode and age of the gold mineralisation, there is consensus that the Storø supracrustal belt formed in a convergent plate tectonic setting, where continentally derived sediments were deposited along with arc-related mafic and intermediate volcanic rocks.

Ameralik Supracrustal Belt
The Ameralik supracrustal belt at 64°05′N in the Sermilik block ( Figure 1) was mapped at a scale 1:100,000 (Rehnström, 2011), but it has not been described in detail. Szilas et al. (2015) made a geochronological and geochemical investigation of the belt. It is not a coherent unit but consists of numerous inclusions of amphibolite, ultramafic rocks, garnet-mica schists, sillimaniteand cordierite-bearing schists, quartzites and remnants of layered anorthosite complexes that are sheared and boudinaged with the adjacent TTG gneisses. The sillimanite-cordierite-bearing rocks and the quartzites were interpreted as metasedimentary rocks, but might also represent hydrothermally altered (leached) rocks of volcanic origin as documented at Qussuk (see Qussuk-Bjørneøen Metavolcanic Belt); notably Szilas et al. (2015) showed that several of the rocks recognised as of volcanic origin have been hydrothermally altered. The least-altered metavolcanic rocks have a narrow major-element compositional range with SiO 2 contents between 48.0-51.2 wt%, but Szilas et al. (2015) identified two different geochemical groups based FIGURE 4 | Epithermal (pre-metamorphic) leaching of andesitic volcaniclastic rocks in the Qussuk-Bjørneøen metavolcanic belt north of Qussuk. (A) Unaltered meta-tuffaceous rock with partial melt seams. (B) Garnet-sillimanite-rich rock almost devoid of plagioclase, interpreted as a moderately leached volcaniclastic rock originally similar to that in (A). (C) Quartz-sillimanite rock interpreted as an altered andesitic volcaniclastic rock subjected to pre-metamorphic advanced argillic alteration. (D) REE diagrams of an unaltered volcaniclastic rock with a normal, LREE-enriched pattern (green) and a quartz-sillimanite rock strongly enriched in all REE except Eu, caused by leaching of up to 80% of the original rock and removal of all plagioclase, hence the negative Eu anomaly. (E,F) Zircon with three different age and geochemical components, comprising 1) original ∼3,075 Ma 207 Pb/ 206 Pb igneous cores with moderate U and Th contents and high Th/U ratios, 2) pre-metamorphic, hydrothermal zircon with abundant quartz inclusions and very high U and Th contents and younger apparent 207 Pb/ 206 Pb ages (partially metamict and probably partially affected by early lead loss), and (3) rims with lower U and Th contents, low Th/U ranges and yielding ∼3,000 Ma 207 Pb/ 206 Pb ages known from the regional metamorphism. See Garde et al. (2012a) for further details.
Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 540997 on trace elements. Like many other belts described here, as pointed out in the geochemical review by Szilas et al. (2018), the two groups do not represent a single magmatic suite, because they contain rocks with different crustal and mantle source components (Figure 2), as also indicated by the large range of initial εNd 2970Ma values from 0.0 to +5.7 and the initial εHf (2,970 Ma) values from +0.7 to +10.4 in the least-altered amphibolites (Szilas et al., 2015).

Nunatak 1390
As its name implies, Nunatak 1390 is a nunatak at 63°40′N in the Greenland ice cap in the Sermilik block, east-north-east of Sermilik- Figure 1 (Stendal and Scherstén, 2007;Scherstén et al., 2008;Szilas et al., 2012b). The nunatak contains a wellpreserved bimodal volcanic succession of mafic to ultramafic rocks, as well as associated acid volcanic rocks and granite intrusions. The nunatak was first described by Escher and Pidgeon (1976) and later mapped in detail by Stendal and Scherstén (2007), who established a relative stratigraphy; the thickness of the entire volcanic package is about 2 km. Being only slightly to moderately deformed, the rocks show wellpreserved primary pillow lavas, finely laminated tuffites, ignimbrites ( Figure 5), and relict vesicles in pillows and flows.
The lower half of the supracrustal package is dominated by deformed mafic pillow lavas and pillow breccias that are extensively altered and contain calc-silicate minerals, and cut by an E-W-trending swarm of 1-5 m-thick, gabbroic to noritic dykes. In the upper pillow lavas ultramafic greenstones and soapstones have been interpreted as intrusive sills. Above the upper mafic pillow lavas there is a ca 80 m-thick unit of felsic volcanic and pyroclastic rocks that include ignimbrites with recognizable fiammé textures; these rocks are cut by felsic dykes (<1 m) interpreted as feeders to higher felsic lavas in the pile. The felsic rocks are overlain by mafic flows and ash beds that are intercalated with ultramafic sills and mafic rusty layers containing sulfides and tourmalinites. The top 700-800 m of the sequence consist of finely bedded tuffs, which pass upwards increasingly to subvolcanic granite intrusions.
Geochemically the lower mafic-ultramafic rocks have whole rock compositions that plot along a well-defined tholeiitic trend (high Fe/Mg series, Figure 2), whereas the felsic rocks, which are intercalated with tholeiites, form a loosely defined calc-alkaline group (low Fe/Mg) with an andesitic to dacitic composition. The tholeiites have Nb/La ratios of 0.6 ± 0.1 (1σ; n 9), whereas the felsic volcanics have a substantially lower ratio of 0.14 ± 0.03 (1σ; n 5); the tholeiitic ratios are thus reminiscent of those of lower crust, whereas the low ratios of the felsic volcanics are typical of arc-related volcanic rocks. Based on the close spatial relationship between the felsic and tholeiitic rocks, Stendal and Scherstén (2007) considered that an arc setting was most plausible. They contended that the intercalation of the mafic and felsic rocks could by analogous to that of a modern arc system, reflecting input from volcanoes in the trench side (tholeiites) and calcalkaline rocks in an arc to back-arc.

Graedefjord Supracrustal Belt
The Graedefjord supracrustal belt, situated at 63°20′N on the southern side of Graedefjord in the Sermilik block (Figure 1), is bimodal comprising about 50% amphibolite (metabasalts) and 40-50% leuco-amphibolite (mostly meta-andesite) together with minor ultramafic rocks and late granitoids (see later). The Graedefjord area is characterised by intense cataclastic and brittle-ductile mylonites (Dziggel et al., 2017). The degree of deformation varies throughout the supracrustal belt, being strongest toward the sheared contacts with the surrounding TTG gneisses. Toward the interior of the belt deformation is less pronounced and primary pyroclastic lithologies such as agglomerates, tuffs and volcanic breccias are preserved.
The supracrustal belt has a poorly constrained Sm-Nd and Lu-Hf isochron formation age of ∼2,970 Ma, similar to that of the Fiskenaesset Complex (Szilas et al., 2013a). Whole rock geochemical data, which confirm the bimodal compositions inferred from field observations, define the two main groups Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 540997 that belong to different chemical differentiation series. The amphibolites have a tholeiitic (high Fe/Mg) affinity coupled with a moderately depleted incompatible trace element signature of (La/Sm) N ∼1) and weak negative Nb-Ta anomalies of (Nb/Nb* 0.6-1.1), which plot above the mantle array in a Th/Yb vs. Nb/Yb Pearce diagram (Figure 2). In contrast, the leuco-amphibolites are calc-alkaline with andesitic compositions that show moderately enriched incompatible trace element signatures with (La/Sm) N 2.7-5.7 and negative Nb-Ta anomalies of Nb/Nb* 0.01-0.45) that plot well within the arc-field on a Th/Yb vs. Nb/Yb Pearce diagram. Geochemical modelling suggests that the trace element variation of the leuco-amphibolites is due to mixing between basaltic mantle-derived magma and crust-derived felsic magma with at least 50% contribution from the latter. Hf-isotope constraints point to contamination with pre-existing continental crust that has an age of c. 3,250 Ma, a model that is underpinned by the preservation of coeval inherited zircons in local TTG gneiss. Ultramafic rocks of the Graedefjord supracrustal belt have platinum-group element (PGE) patterns that are similar to those of magmas generated from a high-degree melting of mantle, but they have relatively enriched trace element patterns. The rocks are interpreted to represent arc-related picrites or alternatively high degree partial melts of metasomatised subcontinental lithospheric mantle. Szilas et al. (2013a) concluded that the proposed mixing model for the Graedefjord supracrustal belt is comparable to that in modern continental subduction zone environments, where mixing and crustal contamination processes are univocally important. A possible difference may be the extent of crustal melting, which was higher in the Archaean due to a hotter Archaean sub-continental mantle, compared to that of modern arcs where assimilation of pre-existing continental crust is relatively more important.

Bjørnesund and Ravns Storø (Ikkattup Nunaa) Supracrustal Belts
The Bjørnesund and Ravns Storø supracrustal belts that crop out at 63°N just north of Frederikshåb Isblink (Figure 1) consist of upper amphibolite facies, andesitic metavolcanic rocks with prominent pyroclastic lapilli tuffs, agglomerate bombs and pillows (Andersen and Friend, 1973;Keulen et al., 2014, including a summary of previous work). The two belts belong to the same mafic volcanic sheet that was subsequently folded. Leuco-amphibolites in the Ravns Storø belt have a whole-rock Sm-Nd isochron age of 3,020 ± 78 Ma and juvenile values of εNd (2,970 Ma) of +2.1 to +3.5 and εHf (2,970 Ma) of +3.5 to +4.3 (Szilas et al., 2012a), while Nutman et al. (2004) reported a SHRIMP zircon Pb-Pb age of 2,908 ± 13 Ma for a volcanic sediment in the same supracrustal belt. The supracrustal belts were intruded by the anorthositic-gabbroic-ultramafic rocks of the Fiskenaesset complex (see below), by subhorizontal sheets of diorite, and later protoliths of the granodiorite-trondhjemite orthogneisses (Keulen et al., 2014) that make up about 75-85% of the Archaean crust in West Greenland (Kalsbeek and Myers, 1973;Steenfelt et al., 2005). The main regional deformation at Bjørnesund is associated with major folds (F2 in most literature), thrusts and shear zones.
Large-scale regional F2 folds throughout the entire area between Fiskenaesset and Frederikshåb Isblink (e.g., Myers, 1985) belong to a deformation event that moved from south (Paamiut area), north-eastwards (Bjørnesund region), to north-west (north of Fiskenaesset) between 2.86 and 2.83 Ga . The folding was associated with intrusion of granites in fold axial planes.
The Bjørnesund supracrustal belt is mineralised with gold, which represents an original arc-related mineralisation that was metamorphosed and variably remobilised. According to Kolb et al. (2013) thrusts and shear zones that were active during the main deformation event acted as pathways for fluids, which created interactions with the host amphibolites and transported the gold into the amphibolites. The gold accumulation is thought to have occurred during the intrusion of granites into the Bjørnesund supracrustal belt and during deformation associated with coeval thrusts and strike slip faults (Keulen et al., 2014). Metasomatic aluminous gneisses containing cordierite, garnet and sillimanite were generated during hightemperature hydrothermal alteration of the volcanic rocks Schlatter and Stensgaard, 2014); these are difficult to distinguish from metamorphosed rocks that had already been altered in a surficial environment.

Nigerlikasik Belt
The Nigerlikasik metavolcanic belt at 62°05′N is the easternmost of three subparallel supracrustal belts in the Paamiut region of the Kvanefjord block (Figure 1). The entire sequence was metamorphosed at upper amphibolite facies and subjected to at least two stages of folding (Klausen et al., 2017). An E-W, c. 550 m-profile shown in Figure 6 demonstrates that the belt consists of a mafic to ultramafic lower part and an intermediate to felsic upper part. From base to top the belt comprises: 1. Ultramafic serpentinites (hornblende ± phlogopite ± diopside ± serpentinite/talc). 2. Mafic amphibolites, generally fine-grained, layered hornblende ± biotite ± plagioclase rocks. 3. Intermediate, andesitic-dacitic, fine-grained, layered biotite ± hornblende ± plagioclase ± quartz schists and gneisses. The surrounding felsic TTG gneisses are clearly intrusive into the supracrustal sequence. The Nigerlikasik belt has a minimum age of 2,930 Ma, constrained by zircon U-Pb ages of intrusive TTG-like sheets within the metavolcanic sequence (Klausen et al., 2017). Geochemically the Nigerlikasik belt records a distinct geochemical change in magmatism from tholeiitic (lower part) to calc-alkaline (upper part) (Figure 2). This indicates a fundamental change in the process of derivative magma generation, because it is not possible to link the two rock suites simply by magmatic differentiation. Klausen et al. (2017) interpreted the entire metavolcanic sequence in the framework of a single maturing island arc, with an initial tholeiitic stage characterised by more juvenile (mantlesourced) magmatism, followed by a differentiated calc-alkaline stage when crustal sources were involved in the generation of the andesite-dacite suite. The final stage of TTG intrusion was interpreted as a separate stage (unrelated to the andesite-dacite suite) of melt formation in the lower crust, which by that time Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 540997 would have been substantially thickened by magmatic accretion and might constitute a proto-craton.

Tartoq Group
The Au-mineralised Tartoq Group at 61°45′N (Higgins and Bondesen, 1966;Higgins, 1968;Appel and Secher, 1984;Evans and King, 1993) comprises greenschist-to granulite-grade, mainly mafic metavolcanic rocks with well-preserved subaqueous lava flows, basaltic pillows, andesitic lapilli tuffs, shallow crustal mafic sills and dykes, gabbroic rocks, and serpentinites that may have been derived from harzburgitic cumulates or from the mantle, as well as sedimentary banded iron formation; this combined rock assemblage closely resembles the lithological units that typically occur in oceanic crust (Szilas et al., 2013b). The mafic lavas have chondrite-normalised REE patterns and Th/Y-Nb/Yb ratios indicative of an arc affinity when compared with modern island arcs (Figure 2). Nutman and Kalsbeek (1994) reported a SHRIMP zircon U-Pb age of 2,944 ± 7 Ma for an intrusive TTG, which provides a minimum crystallisation age, while Lu-Hf regression gives a model age of c. 3,189 ± 65 Ma for the Tartoq Group (Szilas et al., 2013b). Saintilan et al. (2020) obtained Re-Os model ages of c. 3.18-3.13 Ma from arsenopyrite in gold-sulphide-mineralised shear zones and interpreted the mineralisation as associated with arc volcanism (and thus likely originally epithermal), whereas Nutman et al. (2004) found detrital zircons as young as 2,842 ± 6 Ma. The various latter data may indicate that the Tartoq Group contains significantly different age components. Kisters et al. (2012) and Polat et al. (2016) described the postdepositional geological history of the Tartoq Group. It is interpreted as beginning with burial and subduction of supracrustal rocks, followed by exhumation and expulsion, and finally overprinting by transcurrent shearing. Most of the observed features comply with those occurring during hightemperature subduction and retrogresssion during exhumation at a modern convergent margin (see also van Hinsberg et al., 2018). Different from current subduction settings are the absence of high-P metamorphic rocks and the close spatial and temporal association between granitoid plutonism and the structural burial and return flow of the Tartoq Group, both of which are explained by a higher geothermal gradient than today in this Archaean subduction zone. Slab melting and the formation of granitoids at Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 540997 a relatively shallow depth in the subduction channel move the magmatic front to the forearc region of the convergent margin. The high geothermal gradient would have caused a higher volume of felsic melts and a weaker crust, enhancing the disintegration of the downgoing slab and thus producing a short-lived subduction event alternating with slab break-off and backstepping of the subduction front in a wide accretionary complex. This scenario was replicated by Polat et al. (2016) who made a convincing comparison between the accretionary Mesoarchaean Tartoq Group and the Upper Cretaceous Chugach-Prince William accretionary complex in Alaska, having formed by subductionaccretion processes.

Magma Chambers and Layered Complexes in Island Arcs
Both the lower and upper crustal sections of many Archaean cratons contain minor, but prominent and distinctive layered intrusions that consist predominantly of megacrystic, calcic anorthosites (often chromite-layered) and anorthositic leucogabbros that are typically accompanied by gabbros, melanogabbros, dunites and harzburgites. These complexes are often amazingly well-preserved, and contain important diagnostic evidence for their mode of origin and tectonic environment. There is at least 500 km strike-length of such anorthositic complexes in the West Greenland NAC. Yet, they are surprisingly poorly referenced in many papers on Archaean crustal evolution. Overall, it is widely considered that these intrusions were generated in oceanic island arcs and derived from sub-arc hydrous mantle sources.

Fiskenaesset Complex
West Greenland contains the most extensive and best-known examples of these layered intrusions, of which the Fiskenaesset complex at 63°20′N (Figure 1) is the largest and most complete. Its original maximum stratigraphic thickness is close to 540 m. The current tectonic thickness reaches 2 km and some folded layers have an exposed length of at least 30 km (Myers 1985). Although isoclinally and tightly folded three times, metamorphosed up to granulite facies in some areas and commonly retrogressed to high amphibolite facies, the complex has a remarkably well-preserved magmatic stratigraphy, which from top to bottom is: upper gabbro, chromite-layered anorthosite with megacrystic plagioclase (An 75-95 ), upper leucogabbro, middle gabbro, lower leucogabbro, ultramafic unit (peridotite, dunite, hornblendite), lower gabbro; the weighted lithological average is leucogabbro (Myers, 1985). Excellently preserved primary structures include hundreds of mineral/size-graded layers, cumulate textures, scour channels and trough layering. The leucogabbros are characterised by a cumulate plagioclasehornblende assemblage in which plagioclase megacrysts may reach 30 cm across, and undeformed chromitite layers in anorthosites are up to 20 m thick. The Fiskenaesset complex has a Sm-Nd age of 2,973 ± 28 Ma, and a Pb-Pb isotope regression age of 2,945 ± 36 Ma (Polat et al., 2010). Best-preserved anorthosites and leucogabbros contain zircons with a weighted mean U-Pb age of 2,936 ± 13 Ma, which provides a minimum crystallisation age of the complex (Souders et al., 2013), or 2,973 ± 28 Ma (Polat et al., 2010), or 2,985 Ma . Huang et al. (2012) demonstrated that the complex is characterised by depletions of high field strength elements (HFSE; Nb, Ta, Zr, and Hf) relative to Th, U, LREE and LILE (Rb, Ba), as well as enrichments of Sr and Pb relative to MREE, which are consistent with generation in an oceanic island arc. They also showed that abundant magmatic hornblende, the high anorthite content of plagioclase, and negative anomalies of Nb in hornblende are collectively consistent with a sub-arc, hydrous, parental liquid of high-alumina basalt derived from a mantle source (see also, e.g., Hoffmann et al., 2012;Souders et al., 2013;Polat et al., 2015). The Fiskenaesset complex was emplaced into island-type volcanic rocks (now amphibolites) at 2.97-2.95 Ga (Polat et al., 2010;Keulen et al., 2010), intruded by voluminous tonalites at 2,863 ± 10 and c. 2,802 Ma (Friend and Nutman, 2001;Steenfelt et al., 2005;Hoffmann et al., 2012), and likely underwent collision with other arcs and continental blocks in the late Neoarchaean. A characteristic and diagnostic feature of the Fiskenaesset complex is the presence over hundreds of kilometres of spectacular zones with a remarkable metasomatic assemblage with inter alia sapphirine, ruby, kornerupine, phlogopite, pargasite, enstatite, and red spinel-see Figure 7. The ruby and other aluminous phases were formed by metasomatic reactions at the contacts between the anorthosite and ultramafic rocks (Herd et al., 1969;Appel and Ghisler, 2014;Keulen et al., 2014;Smith et al., 2016;Keulen et al., 2020). The deformation and crustal thickening associated with early thrusting, subduction, and collision in the Fiskenaesset region gave rise to large-scale double and triple fold interference patterns and broadly synchronous amphibolite-granulite facies metamorphism (Windley and Garde, 2009).

Other Anorthositic Complexes
Many smaller anorthosite-leucogabbro-gabbro-ultramafic complexes lack the diagnostic stratigraphy, chromitites and associated distinctive sapphirine-bearing rocks of the Fiskenaesset complex, and thus are considered to be derived from separate island arcs. Several were summarised by Windley and Garde (2009). These include.
A. The Naajat Kuuat complex in the eastern Tasiusarsuaq terrane (Figure 1), which has been subjected to the most detailed geochemical and isotopic study of all the anorthositic complexes in Greenland . B. The Buksefjorden and Tre Brødre complexes in the Tasiusarsuaq terrane (Friend et al., 1987). C. Ivittuut, South Greenland. This up to 15 km wide and 55 km long anorthosite-gabbro complex was defined entirely by mapping innumerable anorthosite-gabbro inclusions in the host TTG orthogneisses (Berthelsen and Henriksen, 1975). The original igneous complex was invaded passively by voluminous protoliths of the gneisses, leaving the inclusions in their original magmatic stratigraphic positions; this illustrates the remarkable degree of preservation of original material that is possible in the Earth's deep crust, and emphasises the importance of not neglecting the role of inclusions in TTG gneisses.

E. A conspicuous isolated anorthosite mountain called
Fuglefjeldet occurs at the south-eastern boundary of the Fiskefjord block on the eastern side of the Qussuk peninsula.

Dunitic Bodies in the Fiskefjord Block
The south-western part of the Fiskefjord block contains several kilometre-sized lenses of ultrabasic rocks comprising peridotite, dunite and occasional layered chromitite. They are older than the regional TTG orthogneisses and mostly associated with layered metagabbroic and noritic rocks; a few are isolated enclaves in orthogneiss. Some contact zones are jumbled mixtures of ultrabasic rocks, amphibolite and orthogneiss. The dunitic component of the Seqi body (Figures 1 and 8) is highly refractory (Fo 91-93 ; Szilas et al., 2018). The dunitic bodies may represent residual upper mantle like the nearby kimberlite xenoliths (Bernstein et al., 1998;Bernstein et al., 2007), but their common association with metagabbros suggests they are more likely basal parts of layered complexes as indicated by their chromite compositions and platinum-group element patterns (Garde 1997;Guotana et al., 2018;Szilas et al., 2018). The Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 540997 dunitic bodies are clearly not metamorphosed komatiites and do not yield support for plume or sagduction models.

The Norite Belt in the Fiskefjord Block
About 75 noritic intrusions ∼3.0 Ga old with hybrid margins to their hosts occur in the northern Fiskefjord block, and related post-kinematic diorites farther south (Secher 1983;Garde 1991;Garde et al., 2013b;Waterton et al., 2020). The very variable geochemistry of these highly unusual rocks with high Ni and Cr contents suggests derivation from ultramafic magmas contaminated with 30% TTG material. They are unrelated to the older mafic-ultramafic complexes and difficult to explain, since the crustal evolution of the Fiskefjord block up to 3.0 Ga was normal for the NAC. Origins by crustal thinning and extraterrestrial impacting have been suggested (see below).

Collision Between the Godthåbsfjord-Ameralik Belt and the Tasiusarsuaq Terrane (Sermilik and Bjørnesund Blocks), and High-Pressure Metamorphism
New evidence of extensive high-pressure metamorphism (Dziggel et al., 2014(Dziggel et al., , 2017Friend and Nutman 2019;Nutman et al., 2020) has reinforced the importance of collisional tectonics in the boundary region between the Godthåbsfjord-Ameralik belt and the Sermilik block (McGregor et al., 1991), supplementing lowpressure metamorphism interpreted by Wells (1979) as a result of injection of voluminous tonalitic magma into the middle crust. Prior to the final assembly between ca. 2,720-2,700 Ma (Nutman and Friend 2007;Dziggel et al., 2014;Dziggel et al., 2017;Dziggel, 2019), the structurally upper Tasiusarsuaq terrane underwent a prolonged period of medium-and high-pressure granulite facies metamorphism from ca. 2,800 Ma, which terminated in nearisobaric cooling and amphibolite facies reworking (Crowley 2002;Kolb et al., 2012;Dziggel et al., 2012;Dziggel et al., 2014;Dziggel et al., 2017). An excellent explanation of the history of the terrane concept in the Godthåbsfjord region and how it solved controversies in the 1980s between field observations and geochronological data is by Friend and Nutman (2019). The Tasiusarsuaq terrane (Sermilik block, Figure 1) extends from the southern Nuuk region at Ameralik fjord all the way southward to Frederikshåb Isblink ( Figure 1) and comprises TTG gneisses with intrusion ages of 2,880-2,860 Ma (Friend and Nutman 2005a;Naeraa and Scherstén 2008). Geochronological data suggest that a significant earlier crust-forming event took place at around 2,920 Ma (Schiøtte et al., 1989;Kokfelt et al., 2011). McGregor and Friend (1992) recognised crustal tilting in the present Bjørnesund block and reported a rare prograde amphibolite to granulite facies boundary in inner Bjørnesund. Dziggel et al. (2017) described deep-crustal granulite-facies rocks in the central Tasiusarsuaq terrane as ductile fold nappes like those in the orogens of the Grenville and Himalayas. The interpretation by Dziggel et al. (2017) was based on the estimated thickness of crust in the Neoarchaean (≥60 km), an apparent extrusion of deep-crustal nappes into the middle crust, and the evidence for contemporaneous underplating of cool and older continental crust, as suggested by McGregor et al. (1991) for the Godthåbsfjord region. Dziggel et al. (2019) argued that deformation of the Tasiusarsuaq terrane evolved from initial accretion to terminal continent-continent collision.

Evidence of Collision in the Kvanefjord Block
The Kvanefjord block was studied by McGregor and Friend (1997), Friend and Nutman (2001) and Nutman et al. (2004) who proposed several new tectono-metamorphic terranes (the Sioraq, Paamiut, Neria and Sermiligaarsuk blocks in their terminology).
Using published maps and subsequent observations, Windley and Garde (2009) re-examined the regional distribution of amphibolite facies, granulite facies and amphibolite facies rocks retrogressed from granulite facies within the Kvanefjord block and showed that the region can be viewed as one contiguous crustal segment of prograde amphibolite to granulite facies rocks, overlain by a major fold nappe, the Neria nappe, which was thoroughly retrogressed from the granulite facies and refolded by younger, upright folds. In contrast to early, recumbent isoclinal folds in other areas, which document horizontal tectonics during the crustal accretion, the proposed nappe affected rocks that had already been metamorphosed at granulite-facies conditions. It may well constitute evidence for collisional tectonics, comparable to that described from the Sermilik block farther north.
New Crustal Accretion Models for the Fiskefjord Block Kirkland et al. (2018) proposed a new Meso-/Neoarchaean orogenic cycle in the north-western Fiskefjord block based on zircon geochronology of thin intercalations of homogeneous quartzo-feldspathic sheets in metavolcanic rocks. The model inferred complete denudation, a new regional unconformity, formation of a younger series of metavolcanic rocks than that established by previous workers and a second phase of granulite facies metamorphism after ∼2.8 Ga, see Discussion. Yakymchuk et al. (2020) examined static partial melting of mafic rocks at four localities in the same region. They dated the partial melts to 3,010-2,990 Ma and modelled partial melt compositions in these rocks and compared them to the compositions of the regional TTG gneisses. Their modelling was consistent with an apparent, unusually high geothermal gradient of >1,050°C/GPa at around 3.0 Ga, implying thin crust. Along with the static metamorphism this prompted them to propose a non-uniformitarian, "stagnant lid" model (Sizova et al., 2015) with generation of TTGs from melting of coeval lower-crustal basic rocks without subduction; they also noted widespread regional evidence of early zircon lead loss. They expressed caution about their interpretation due to limitations in their modelling and a rather imperfect geochemical match of their modelled melt compositions with the regional TTG gneisses.

The Maniitsoq Structure
The most important but also the most controversial discovery in the northern Fiskefjord block is the 3.0 Ga Maniitsoq impact structure, which occupies most of its northern part including the study area addressed in New Crustal Accretion Models for the Fiskefjord Block. (Figures 1 and 9; Garde et al., 2012b;Garde et al., 2013a, Garde et al. 2013bScherstén and Garde, 2013;Garde et al., 2014;Keulen et al., 2015). This structure comprises several anomalous, crustal-scale features ( Table 2) which have recently given rise to different interpretations. The Finnefjeld domain 35 by 50 km large in the centre of the structure, consists of mechanically mixed cataclastic material, first described by Berthelsen (1962). The norite belt, which is not cataclastic, is located just east of this domain. The melt zone surrounding the Finnefjeld domain ( Figure 1) contains evidence of direct melting of K-feldspar and biotite and in places localised bulk melting of felsic orthogneiss, besides further cataclasis and widespread, intense hydrothermal alteration in some parts. A minimum age from intensely altered zircon in the eastern part of the structure gave a hydrothermally reset U-Pb age of 3,000.9 ± 1.9 Ma (Scherstén and Garde, 2013).
The Maniitsoq structure is interpreted as the lower-crustal expression of a giant extraterrestrial impact structure with a surface diameter of ∼500 km, estimated from iSALE impact modelling ( Figures 9D-F; Trowbridge et al., 2017). Direct (shock) melting of K-feldspar and biotite requires an extremely rapid and short-lived temperature excursion to temperatures not attainable by endogenic crustal processes (>1,300 and >1,400°C, respectively) and signifies an extraterrestrial impact (French and Koeberl, 2010). Whereas unequivocal shock lamellae in quartz have not been preserved, shock lamellae in zircon were described by Garde et al. (2019). Glikson (2013) recorded the Maniitsoq structure as the oldest known impact structure on Earth. Because of its large size and deep depth of preservation it constitutes an important target for studies of impact processes approaching those that shaped the outer Earth during the ∼3.85 Ga Late Heavy Bombardment.
A Relict Suture Between the Fiskefjord and Maniitsoq Blocks Windley and Garde (2009) placed their tectonic boundary between the Fiskefjord and Maniitsoq blocks where a metamorphic facies change and previous age determinations suggest a major break, namely a 2.5 ± 0.2 Ga metasedimentary 87 Sr/ 86 Sr whole-rock age on Hamborgerland (Kalsbeek and Hansen, 1989), a metamorphic zircon age of 2,738 ± 6 Ma near Maniitsoq (Friend and Nutman, 1994), and detrital and metamorphic zircon ages of 3,180-2,700 and 2,550 Ma respectively, of a metasedimentary gneiss at Kangerluarsuk (Garde et al., 2000). With new zircon ages of ∼2,560 Ma from migmatite leucosomes about 5 and 30 km north-west of Kangerluarsuk, respectively, Dyck et al. (2015) proposed a S-dipping suture (a palaeo-subduction zone, Figure 10) along Kangerluarsuk fjord near the proposed block boundary (Figure 1). This zone includes medium-to high-pressure, kyanite-bearing rocks, and disrupted lenses of dunite surrounded by metasomatic corundum-bearing rocks, Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 540997 suggesting tectonic extrusion along the proposed suture. Adjacent anthophyllite-gedrite-bearing amphibolites most likely represent hydrothermally altered tholeiitic lithologies similar to those containing semi-precious ortho-amphibole gemstones near Nuuk (Appel and Jensen, 1987). Dyck et al. (2015) also described a lithological change and different stream sediment compositions across the tectonic boundary, and a gravity high, from which they proposed a new tectonostratigraphic unit, the Majorqaq belt, between the Fiskefjord and Maniitsoq blocks (Figure 1). This belt supposedly constitutes the last phase of amalgamation between the Fiskefjord and Maniitsoq blocks. The belt probably contains several unrelated imbricated components, given metamorphic zircon ages of ∼2,815 Ma near Kangerluarsuk and ∼3,010 Ma just north of Maniitsoq (Garde et al., 2000;Yakymchuk et al., 2020). The metamorphic history of the Majorqaq belt indicates a clockwise P-T evolution, interpreted as evidence for rifting and subduction, rather than over-accretion that would produce an anti-clockwise evolution (Dyck et al., 2015). The inferred ∼2.56 Ga subduction is contemporaneous with the intrusion of the crust-derived Qôrqut granite complex in the Nuuk region (Nutman et al., 2010;Naeraa et al., 2014); Dyck et al. (2015) suggested that dehydration fluids from the subducted oceanic plate triggered lower-crustal remelting that produced the granite. In summary, the southern boundary of the Majorqaq belt constitutes a strong candidate for Neoarchaean plate tectonic subduction and collision, essentially undisturbed by later crustal events.
Xenocrystic zircons in a 557 Ma kimberlite on the southern side of the Majorqaq belt yielded a wide range of Eo-to Neoarchaean ages that were interpreted as lithospheric growth by tectonic stacking after the assembly of the NAC .   Berthelsen (1962); Garde et al. (2012b, Garde et al. 2014 Direct melts of K-feldspar and biotite Observed in granitic gneisses and with micro-cataclasites. Note: direct melting is a shockinduced, instantaneous phase transformation from mineral to melt without chemical reaction with other phases. Scherstén and Garde (2013, Figure 2); Garde et al. (2014): Keulen et al. (2015) Micro-cataclasites Widespread micro-cataclastic zones akin to pseudotachylytes, formed in situ. Interpreted as reflecting numerous pulses of intense, deviatoric pure shear stress. First described as "grey dykes" in 1976 in NW Godthåbsfjord. See also Garde and Klausen (2016). Garde et al. (2012b, Figures 4a, 5);Bridgwater et al. (1976) K-feldspar veins lining microcataclasites Veins of K-feldspar lining micro-cataclasites, crystallised from K-feldspar melts and undeformed. Interpreted as direct melts formed during cataclasis and separated from the solid, micro-cataclastic constituents.  Planar microstructures in zircon External, very closely spaced, planar microstructures in several directions, resembling those in, e.g., the Vredefort and Sudbury impact structures and indicative of shock metamorphism. Note: quartz microstructures described by Garde et al. (2012) are poorly preserved and not diagnostic. Garde et al. (2019) Hydrothermal alteration Very widespread, very intense and commonly associated with brecciation of mafic rock components. Not studied in detail. Prior to 2012 interpreted as endogenic metamorphic retrogression from granulite facies. Scherstén and Garde (2013); Garde et al. (2014, Figure 11) Coincidence of radiometric ages close to 3.0 Ga A highly unusual regional preponderance of ∼3.0 Ga and slightly younger zircon U-Pb ages in TTG orthogneisses and crustal melts, including reset ages. Note: this regional pattern of radiometric ages has not yet been fully understood. Garde et al. (2014, Garde (1991); Garde et al. (2012b) Aeromagnetic anomaly First described in 2012 as a potential primary impact-related feature, but possibly related to the subsequent hydrothermal alteration. Garde et al. (2012b) Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 540997

DISCUSSION
This brief review of the geology of, and recent research in, the Greenlandic section of the NAC provides us with the opportunity to evaluate that geology in terms of the current, polemic debate about the controversial uniformitarian versus nonuniformitarian models.

Contrasting Interpretations in the Northern Fiskefjord Block
We begin with an examination of contrasting interpretations of the northern Fiskefjord block, where Yakymchuk et al. (2020) recently proposed a non-uniformitarian, stagnant lid model of 3.0 Ga crustal accretion within the area of the Maniitsoq structure. Non-uniformitarian crustal accretion models for the Meso-to Neoarchaean NAC are rare, and the new model is interesting. Yakymchuk et al. (2020) examined static partial melting in mafic lithologies and argued that their observations are compatible with production of voluminous TTG partial melts at around 3.0 Ga from mafic lower crust in a static or extensional scenario with an unusually steep geothermal gradient and relatively thin crust. They argued that the effect of supposed crustal thinning was effective simultaneously over a large area and can therefore not only be explained by accretionary processes. They did not, however, supply observational support for crustal extension or test their model against the Maniitsoq structure. In the Maniitsoq structure, the biotite melt seams in cataclastic plagioclase and disequilibrium microstructures of three coexisting feldspars in granitic gneisses document a very large, but short-lived temperature excursion, to at least 1,300-1,400°C (Scherstén and Garde, 2013;Keulen et al., 2015), succeeded by long-term lower-crustal thermal perturbation under tectonically static conditions. Whereas the concept of large-scale mechanical mixing of cataclastic rocks in the lower crust is clearly beyond all normal endogenic crustal paradigms, it is familiar to scientists studying collisional planetary processes and impacts (e.g. Morgan et al., 2016). As argued by the latter authors and Garde et al. (2012b) a meteorite impact is the only way to achieve these relations. In regard to the stagnant lid model of Yakymchuk et al. (2020) this implies that at 3.0 Ga heat was supplied to the lower crust from above and not only from below, which makes the calculation of a normal crustal gradient redundant. On the other hand, their field observations of rock fragmentation (their Figure 2A), partial melting in the mafic lithologies, the ages of the partial melts, and the local P-T modelling are all fully compatible with a giant meteorite impact: Keulen et al. (2015) showed what happened to lower-crustal granitic gneisses in this scenario, whereas Yakymchuk et al. (2020) have demonstrated some high-temperature effects in equivalent mafic lithologies. In conclusion, a stagnant lid model for the northern Fiskefjord block may be speculatively possible, but it does not provide a satisfactory solution to the available salient information. The proposition by Kirkland et al. (2018) of a new orogenic cycle in the northern Fiskefjord block relies heavily on a highly unusual and controversial interpretation of homogeneous felsic layers as sedimentary in origin, and on complex zircon age data with lead loss, must be viewed with caution. Citing Yakymchuk et al. (2020), "Evidence for radiogenic-Pb loss, in U-Pb age data and zircon morphologies, is very common in zircon from Mesoarchaean orthogneiss of the Akia Terrane (Friend and Nutman, 1994;Garde et al., 2012b;Kirkland et al., 2018;Gardiner et al., 2019)." Further discussion of the shortcomings of this proposition is not relevant for a general readership and is not pursued here.
First, we emphasise that for the last 50 years or so there has been a near-total consensus by a huge number of workers in all Earth Science disciplines, who studied the Archaean geology of West Greenland, that the geological relationships are without doubt compatible with equivalent accretionary or collisional processes in the Mesozoic-Cenozoic Earth. As pointed out by Friend and Nutman (2019), the evidence of partial melting of the upper mantle by a hydrous fluid flux, and the absence of komatiites in West Greenland strongly indicate a plate-tectonic regime rather than plume-driven vertical tectonics. A recent boron isotope study by Smit et al. (2019) is a possible exception to the general consensus; the authors suggested that the boron isotopic composition of TTGs seems to be incompatible with derivation primarily from subducted, hydrated oceanic crust. There are several ways in which non-uniformitarian models may impinge upon interpretation of the geology of West Greenland. In principle, there are two types of orogens: accretionary and collisional. Accretionary orogens (Cawood et al., 2009), like the Phanerozoic Japanese Islands (Isozaki et al., 1990;Isozaki et al., 2010) or the Altaids of Central Asia (Wilhem et al., 2012), form by the accretion of small fragments of mid-ocean ridges, ophiolites, ocean plate stratigraphy, seamounts, oceanic plateaus and island arcs. These are fragmental orogens because subduction erosion during offscraping accretion caused most ocean floor rocks to be subducted, and only a few to be accreted (Kimura and Ludden, 1995). Therefore, accretionary orogens do not contain or belong to a global network of plate boundaries (the putative requirements for plate tectonics of Cawood et al., 2018 andBrown et al., 2020), but rather only short boundaries (which arguably cannot even be called plate boundaries) that surround small relicts of island arcs, ophirags, and other rocks. Significantly, accretionary orogens were predominant in the Archaean, when the Earth's lithospheric plates were forming by oceanic tectonics, and that was before the development of large, stable and rigid continental plates in the Palaeoproterozoic, which evolved via the Wilson Cycle (Windley et al., 2020).
Collisional orogens, on the other hand, form by the collision of two continental blocks, like the Himalayas between India and Eurasia, during the final stage in the evolution of a Wilson Cycle (Windley 1995). Accordingly, collisional orogens typically contain long sutures along convergent plate boundaries, which have been maintained for considerable periods. The Himalayan example demonstrates that such a definable boundary can extend for thousands of kilometres from the Pyrenees to Malaysia with a global network of linked plate boundaries.
Returning to West Greenland, the Archaean rocks formed in two main periods, the Eo-and Mesoarchaean. In the Eoarchaean from ca. 3.9-3.6 Ga ago the Itsaq Gneiss Complex contains both the Isua belt and younger plutonic TTG gneisses (Nutman et al., 1996;Komiya et al., 1999;Nutman et al., 2004;Friend and Nutman 2005b;Friend and Nutman, 2013;Friend and Nutman, 2015a;Friend and Nutman, 2019), as well as a few other Eoarchaean fragments, the largest of which is the Aasivik terrane. Although disputed by Friend and Nutman (2019 and references therein) different parts of the Isua belt contain an upward succession from basal low-K tholeiitic basalts with or without pillows, thinly bedded cherts and magnetite-BIF, to uppermost mixed turbidites and conglomerates, which Komiya et al. (1999Komiya et al. ( , 2004 pointed out is essentially identical to modern Ocean Plate Stratigraphy (Kusky et al., 2013), which is transported horizontally to a trench. Bolhar et al. (2004) reported that the REE, yttrium patterns, and diagnostic anomalies (Ce/Ce* and Pr/Pr*) of the magnetite-quartz banded iron formation (BIF), magnetitecarbonate BIF, and banded magnetite-rich quartz rocks are entirely comparable with the same diagnostic trace-element patterns of modern sedimentary equivalents and seawater proxies. These relations strongly suggest that these Eoarchaean rocks in Greenland formed in a horizontal tectonic regime in one or more subduction-accretion complexes (Friend and Nutman, 2019). Furthermore, from their detailed structural studies Hanmer and Greene (2002) reported that the Isua belt is underlain by a major 3.64 Ga mylonitic thrust-nappe stack. The whole geologicalstructural regime demonstrates that the deformational behaviour, rheology, and overall strength of the Palaeoarchaean and modern continental crust were fundamentally similar.
The bulk of the Greenland craton consists largely of Mesoarchaean TTG gneisses, the protoliths of which were intruded into and intercalated with metavolcanic arc-type belts (such as Qussuk-Bjørneøen, Tartoq, Ivisaartoq, Ravns Storø and Nigerlikasik, described above) and layered anorthositic intrusions (such as the Fiskenaesset and Ivittuut complexes) that formed in island arcs. All these upper and lower crustal rocks were thrustimbricated by subhorizontal shortening, as evidenced by abundant early thrusts and associated recumbent isoclinal folds, followed by upright folding by crustal shortening. As shown by Grocott et al. (1999) and Grocott and McCaffrey (2017) this is also the case in the two Palaeoproterozoic orogens bounding the NAC in West Greenland, where large, dome-like compressional folds had previously been interpreted as crustal diapirs. All the wellpublished geological, structural, geochemical and isotopic data of these Mesoarchaean rocks in Greenland point to a subhorizontal tectonic regime. Confirmation of this tectonic setting comes from 2.7 ± 0.3 Ga eclogite xenoliths in a West Greenland kimberlite, which are coeval with the regional TTG magmatism (Tappe et al., 2011 were derived from seafloor-altered oceanic crust and grew by melting of a basaltic slab in a subduction zone. Now we return to the question of the relevance of the West Greenland rocks to the current non-uniformitarian models for Archaean crustal evolution, most of which are based on vertical, gravity-driven tectonics. It has been argued that the start or early operation of plate tectonics on Earth can "only be defined if the rocks can be shown to belong to a global system of plate boundaries" (Cawood et al., 2018;Brown et al., 2020). We fundamentally disagree with this conclusion, because all the papers and observations used to construct this model considered only the application of long-boundary collisional orogens and tectonics, and noticeably failed to consider the fact that the bulk of early to mid-Archaean orogens did not form by collisional tectonics because they are accretionary orogens, and accretionary orogens do not belong to "a globally linked system of plates," as described above and by Windley et al. (2020).
The North Atlantic craton as defined by Bridgwater et al. (1973 a, b) included not only the Archaean craton of Greenland, but also the Saglek block of Labrador in Canada. The Nulliak supracrustal belt at Saglek has been well described and comprehensively analysed by inter alia Komiya et al. (2015), who showed that it consists of an Eoarchaean accreted package of (in ascending primary stratigraphic order) ultramafic rocks, basalts, bedded white and green cherts, pelites, conglomerates and psammites; this is a typical ridge-ocean floor pelagic-trench assemblage that evolves on a moving oceanic plate. This Nulliak pile of volcanicsedimentary rocks contains imbricate thrust duplexes with roof, floor and link thrusts; this type of structure and stratigraphy can only be explained tectonically by the accretion of horizontaltravelled ocean floor rocks in a subduction trench setting. Accordingly, the Labrador segment of the NAC provides indubitable evidence of an Archaean horizontal tectonic regime, which is entirely consistent with that in West Greenland.
(2) Other accretionary orogens that formed in the early Archaean (e.g. Nulliak, Barberton, Nuvvuagittuq, and Yilgarn) are not significantly different from those in the Neoarchaean belts of the Superior and Slave Cratons in Canada. There are some minor differences, but they are not greater than the differences between different Phanerozoic accretionary orogens (3) The period of 3.2-3.0 Ga did not mark the onset of modern-style plate tectonics (e.g., van Kranendonk, 2011;Dhuime et al., 2015;Tang et al., 2016), nor the start of the Wilson Cycle (e.g., Shirey and Richardson, 2011), see Windley et al. (2020). The evidence that the Rb-Sr ratio and silica content of juvenile crust increased dramatically at 3.2-3.0 Ga (Dhuime et al., 2015) and that the crust before 3.0 Ga had a highly mafic bulk composition, and a more felsic bulk composition after that mirrors exactly what takes place in the modern Earth, when juvenile oceanic crust and primitive oceanic island arcs change to a more Si-and Rb-rich felsic crust in Andean-type tonalitic magmatic arcs in active continental margins. The period 3.2-3.0 Ga marks the time in all Archaean cratons of the progressive development of a rigid lithosphere; from a juvenile oceanic crustal regime to a more mature continental-influenced crust at a convergent plate boundary; the period 3.2-3.0 Ga does not mark the onset of plate tectonics, and certainly not the start of the Wilson Cycle, as so presciently and well reported by de Wit and Hart (1993) see also Windley et al. (2020). The major change at 3.2-3.0 Ga is well recorded in West Greenland when the more juvenile Eoarchaean oceanic crust in the Isua supracrustal belt changed to the TTGdominant protoliths of the widespread Mesoarchaean gneisses. (4) Eighty per cent of the Archaean craton consists of 3.25-2.82 Ga lower-crustal TTG gneisses, which have been inter-thrusted with upper crustal arc-derived volcanic rocks and anorthositic complexes. (5) Neoarchaean thrusting and folding of previously accreted terranes took place in the Nuuk region and the Sermilik block at 2.85-2.65 Ga, and probably also in the Neria nappe in the Kvanefjord block (Windley and Garde, 2009), most likely as a result of collision tectonics (Dziggel et al., 2019;Friend and Nutman, 2019). (6) With regard to the current controversy of whether plate tectonics started at c. 3.9 Ga or 3.2-3.0 Ga, like Windley et al. (2020) we reiterate and emphasise the well-documented fact that the evolution of the Archaean craton of West Greenland took place in two main periods: A. In the Eoarchaean Itsaq Gneiss Complex, which is extremely well documented by many publications over at least 20 years by Allen Nutman, Clark Friend, their colleagues and many others, which led to the conclusion that juvenile crust formation took place between 3,900 and 3,600 Ma in a proto-arc by episodic growth at a convergent, incipient plate boundary (Friend and Nutman, 2019). B. At 3.25-2.82 Ga when the bulk of the Archaean craton formed when upper crustal island arc-generated volcanic rocks and anorthositic complexes were inter-thrusted with lower crustal TTG gneisses, the protoliths of which were generated in active continental magmatic margins.
Thus, the NAC contains abundant evidence of semicontinuous crustal growth by modern-style plate tectonic processes throughout the bulk of the Archaean. The geological relations in West Greenland and Labrador within the North Atlantic Craton were fortunately mapped and interpreted from more assured, basic principles in the field constrained by appropriate geochemical and isotopic age data, without recourse to theoretical and speculative laboratory-based models.

AUTHOR CONTRIBUTIONS
All four authors participated in writing the manuscript, with AG and BW being responsible for the major part. AG prepared and compiled the illustrations. All four authors contributed with their intimate knowledge of various parts of the North Atlantic craton based on previous and ongoing field and laboratory work.

FUNDING
There was no specific funding for this work, whereas previous field and laboratory work was mainly supported by the Geological Survey of Denmark and Greenland, with additional contributions for field work from the Ministry of Mineral Resources of the Greenland Government.