Anticlockwise metamorphic pressure–temperature paths and nappe stacking in the Reisa Nappe Complex in the Scandinavian Caledonides, northern Norway: evidence for weakening of lower continental crust before and during continental collision

Abstract. This study investigates the tectonostratigraphy and metamorphic and tectonic
evolution of the Caledonian Reisa Nappe Complex (RNC; from bottom to top:
Vaddas, Kåfjord, and Nordmannvik nappes) in northern Troms, Norway.
Structural data, phase equilibrium modelling, and U-Pb zircon and titanite
geochronology are used to constrain the timing and pressure–temperature
(P–T) conditions of deformation and metamorphism during nappe stacking
that facilitated crustal thickening during continental collision. Five
samples taken from different parts of the RNC reveal an anticlockwise
P–T path attributed to the effects of early Silurian heating (D1)
followed by thrusting (D2). At ca. 439 Ma during D1 the
Nordmannvik Nappe reached the highest metamorphic conditions at
ca. 780 ∘C and ∼9–11 kbar inducing kyanite-grade partial
melting. At the same time the Kåfjord Nappe was at higher, colder, levels
of the crust ca. 600 ∘C, 6–7 kbar and the Vaddas Nappe was
intruded by gabbro at > 650 ∘C and ca. 6–9 kbar. The
subsequent D2 shearing occurred at increasing pressure and decreasing
temperatures ca. 700 ∘C and 9–11 kbar in the partially molten
Nordmannvik Nappe, ca. 600 ∘C and 9–10 kbar in the Kåfjord
Nappe, and ca. 640 ∘C and 12–13 kbar in the Vaddas Nappe.
Multistage titanite growth in the Nordmannvik Nappe records this evolution
through D1 and D2 between ca. 440 and 427 Ma, while titanite
growth along the lower RNC boundary records D2 shearing at 432±6 Ma. It emerges that early Silurian heating (ca. 440 Ma) probably
resulted from large-scale magma underplating and initiated partial melting
that weakened the lower crust, which facilitated dismembering of the crust
into individual thrust slices (nappe units). This tectonic style contrasts
with subduction of mechanically strong continental crust to great depths as
seen in, for example, the Western Gneiss Region further south.



Introduction
Large-scale thrusting and nappe stacking are the main processes responsible for crustal shortening during continental collision. Crustal rocks can be subducted to depths of ∼ 200 km (e.g. Chopin, 2003;Spengler et al., 2009;Hacker et al., 2010), or be included in large-scale nappe stacks (e.g. Escher et al., 1993;Escher and Beaumont, 1997). The crucial factor that controls the style of large-scale deformation and whether rocks are subducted or are stacked in mid-to lower-crustal ductile nappe piles is the rheology, which is dependent on temperature, availability of fluids, composition of the rocks involved, and the presence or absence of melt (e.g. Hollister and Crawford, 1986;Beaumont et al., 2006;Gerya and Meilick, 2010;Labrousse et al., 2010;Jeřábek et al., 2012). The immediate pre-collisional history of rocks involved in nappe stacking plays a significant role in determining these factors. In active continental collision zones such as the Himalaya, mid-and lower-crustal processes are generally inferred from surface deformation, geophysical information, and the geochemistry of erupted volcanic rocks (e.g. Schulte-Pelkum et al., 2005;Zhao et al., 2009). In contrast, older and deeply eroded continental collision zones such as the Caledonides allow for direct study of nappe thrusting at mid-and lower-crustal levels, allowing for insight into the high-temperature processes of large-scale continental subduction and crustal shortening. The nappe concept of largescale thrust units was developed in the Alps (Bertrand, 1884;Schardt, 1893). The nappe units were defined on the basis of stratigraphy (in the sediments), or, in higher-grade metamorphic "basement" units, by nappe dividers of cover sediments between these (e.g. Pfiffner, 2016). When the metamorphic temperatures during thrusting are high in all units, the distinction between basement and cover units becomes virtually impossible. Under such high-grade conditions, the thrust units typically are thin parallel rock slices of variable extent, and the distinction of individual nappe units can only be made on the basis of metamorphic grade or age of metamorphism and deformation. This is the situation in the western part of the northern Caledonides, where only high-grade rocks of the deep part of the orogen form the nappe stack.
The Caledonides were formed by convergence and collision between Baltica and Laurentia in Silurian to Devonian times (Fig. 1a). The resulting large-scale nappe stacks provide access to the study of mid-to lower-crustal processes during a continental collision of Himalayan style and extent (e.g. Streule et al., 2010;Labrousse et al., 2010). In northern Norway, a well preserved section, from the autochthonous Baltica basement to exotic terranes and ophiolites, is exposed ( Fig. 1b), displaying large gradients in metamorphic grade and deformational style . The Reisa Nappe Complex (RNC) represents a large viscously deformed nappe stack, metamorphosed at amphibolite-to granulite-facies conditions, displaying pervasive deformation and possible Caledonian partial melting (Roberts and Sturt, 1980;Andresen, 1988). Based on its similarity to rocks at the same tectonostratigraphic level (e.g. the Magerøy Nappe, located in northern Finnmark; Andersen et al., 1982;Corfu et al., 2006), the RNC is considered equivalent to other Iapetus-derived or outer Baltica margin (upper allochthon) units in the Caledonides. It may also preserve an early Silurian history, recording events immediately prior to or during early continental collision between Baltica and Laurentia (e.g. Andréasson et al., 2003;Slagstad and Kirkland, 2018). In addition, its tectonostratigraphic position directly below the ophiolitic rocks of the Lyngsfjellet Nappe places it in a distinctive palaeogeographic position. Understanding the structure and composition of the RNC together with its metamorphic, deformation, and magmatic history will help to establish how lower-crustal nappe stacking during continental collision takes effect. The pre-Caledonian and Caledonian evolution of the rocks affected the deformation behaviour during nappe stacking in this particular part of the Caledonian orogen.

Geological framework
The Scandinavian Caledonides consist of a series of allochthonous nappes, with, from bottom to top, Baltican, Iapetus, and Laurentian affinities (Roberts and Gee, 1985;Stephens and Gee, 1989;Fig. 1a).

The north Norwegian Caledonides
In northern Norway, the Baltican basement with its metasedimentary autochthonous cover is overlain by the parautochthonous metasedimentary rocks of the Gaissa Nappe Complex ( Fig. 1b; e.g. Ramsay et al., 1985). Above these, the allochthonous rocks form a nappe stack including the following units from bottom to top: (1) the metasedimentary and metaigneous rocks of the Kalak Nappe Complex (KNC), (2) the metasedimentary and metaigneous rocks of the RNC, (3) the ophiolitic and metasedimentary rocks of the Lyngsfjellet Nappe, (4) gneisses and migmatites of the Nakkedalen Nappe, and (5) metasedimentary rocks of the Tromsø Nappe ( Fig. 1b ;Andresen, 1988;Kvassnes et al., 2004;Kirkland et al., 2006Kirkland et al., , 2007aAugland et al., 2014;Rice, 2014;Gee et al., 2017). Caledonian metamorphic conditions associated with nappe thrusting generally increase upwards through the KNC and RNC, from low greenschist facies in the parautochthonous metasediments to upper amphibolite or lower granulite facies at the top of the RNC, with local granulite facies metamorphism generally interpreted to be pre-Caledonian (Andresen, 1988;Elvevold et al., 1993Elvevold et al., , 1994). The Lyngsfjellet Nappe displays greenschist facies metamorphism at its base with higher-grade metamorphism in the overlying metasediments, and the overlying Nakkedalen and Tromsø nappes show amphibolite to eclogite facies metamorphism (Andresen and Bergh, 1985;Andresen and Steltenpohl, 1994;Corfu et al., 2003). The predominant Caledonian deformation, associated with continental collision (often referred to as Scandian; e.g. Corfu et al., 2014) in the allochthons in northern Norway is associated with top-to-SE and top-to-E shearing (Rice, 1998).
The palaeogeographic origin of several of the Caledonian nappes in northern Norway is debated. The Gaissa Nappe Complex is interpreted to represent telescoped Baltica margin cover, and the KNC has traditionally also been interpreted to represent more outboard Baltica basement rocks and its late Precambrian to Paleozoic cover (e.g. Stephens and Gee, 1989;Gee et al., 2017). However, recent evidence Solid Earth, 10, 117-148, 2019 www.solid-earth.net/10/117/2019/  Gee et al. (2010). LN is the Lyngsfjellet Nappe, WGR is the Western Gneiss Region, NNC is the Narvik Nappe Complex, KNC is the Kalak Nappe Complex, RNC is the Reisa Nappe Complex, and TNC is the Trondheim Nappe Complex. (b) Map showing the extent of the Reisa Nappe Complex in northern Norway, based on Zwaan (1988) and own correlations. The Vaddas, Kåfjord, and Nordmannvik nappes are shown in colour. The study area is denoted by the black box (Fig. 2).
shows that the sedimentary cover in the KNC had already been deposited between 1000 and 900 Ma, and the KNC underwent metamorphism, deformation, and magmatism during several events at about 970-950, 870-830, 700, 600, and 570 Ma (Daly et al., 1991;Kirkland et al., 2006;Corfu et al., , 2011. These Neoproterozoic events are not known in the Baltica basement elsewhere, implying that either the KNC has an exotic origin or that a still unknown type of suture (typical suture rocks are lacking beneath the KNC) separated it from the Archean-Paleoproterozoic northern Baltic basement before the final collision (e.g. Kirkland et al., 2006Kirkland et al., , 2007aKirkland et al., , b, 2008Corfu et al., 2007Corfu et al., , 2011. The RNC has previously been interpreted to represent either Iapetusderived or outermost Baltican rocks, whereas the Lyngsfjellet Nappe is traditionally considered as Iapetus-derived, marking the transition towards the Laurentia-derived Nakkedalen and Tromsø nappes (e.g. Stephens and Gee, 1989;Corfu et al., 2003).

Field results
We investigated the RNC east and north of Lyngen in coastal areas, where it is exposed around a window of KNC rocks (Figs. 1, 2). Detailed structural sections through the RNC were investigated on Uløya, supplemented by structural mapping along the east coast of Lyngen, on Kågen and around Straumfjord, and through the nappe stack on Arnøya (Figs. 2, 3). The nappes on Arnøya were previously assigned to the KNC (Roberts, 1973;Zwaan, 1988), but later reconsidered as being part of the Vaddas Nappe (Andresen, 1988). More recent maps have classified northern Arnøya as part of the Magerøy or Vaddas Nappe and southern Arnøya as KNC (e.g. Corfu et al., 2007;Augland et al., 2014;Gasser et al., 2015). Our own field observations indicate that all three RNC nappes (Vaddas, Kåfjord, and Nordmannvik) are present on Arnøya (Figs. 1b, 2, 3).

Lithologies
The lower boundary of the RNC towards the KNC was investigated on Uløya and Arnøya (Fig. 2). The underlying KNC comprises either quartzofeldspathic rocks or metapelitic paragneisses (Figs. 3,4a), whereas the boundary itself consists of a ∼ 40 m thick strongly mylonitized zone comprising intercalated amphibolite facies quartzofeldspathic rocks, gneisses, and schists (Fig. 4b). Marbles are locally present together with the amphibolite and schist. On Uløya, a migmatitic paragneiss that forms part of the KNC below the boundary is subsequently mylonitized in the solid state, with centimetre-sized porphyroclasts of garnet and feldspar (Fig. 4b, c). Kinematic indicators consistently show top-to-SE shearing. Above the mylonitic layer, sheared metaconglomerate ( Fig. 4d) occurs together with a marble unit, which marks the base of the Vaddas Nappe on Arnøya, Kågen, and Uløya (Fig. 3). This marble layer is overlain by amphibolite and calc-schist (Fig. 4e). Two main marble units are found in the Vaddas Nappe at Straumfjord: one in the lower Vaddas and one in the upper Vaddas (Fig. 3). The marble at the www.solid-earth.net/10/117/2019/ Solid Earth, 10, 117-148, 2019 Figure 3. Representative tectonostratigraphy of the RNC from Lyngseidet (west) to Straumfjord and Kvaenangen (east) showing the location of major structures and relative thickness and spatial variation in tectonic units. Based on Zwaan (1988), Lindahl et al. (2005), and own work. The tectonostratigraphic position of the investigated samples is indicated.
base of the lower Vaddas is variably sheared and overlain by intercalated metasediments and metapsammites similar to those in the KNC. The upper marble layer marks the base of the upper Vaddas and is associated with sheared conglomerates. It is overlain by amphibolite and calc-schist, similar to the Vaddas succession on Arnøya, Kågen, and Uløya, and we consider that only the upper Vaddas is found on these three islands (Fig. 3 A well developed strongly mylonitic foliation in muscovite-rich garnet-mica-schist marks the Vaddas-Kåfjord nappe boundary on Uløya and Arnøya (Figs. 2, 3). The boundary cuts a gabbro body in the Vaddas Nappe on Arnøya (Fig. 3). The Kåfjord Nappe is comprised of a lower unit of homogenous garnet-mica-zoisite-schist with calcsilicate lenses and some amphibolite layers ( Fig. 4g-i), and a homogenous upper unit comprised of garnet-mica-schist and gneiss. This upper unit often displays strongly sheared layers rich in quartz-feldspar sigma clasts, possibly indicating the presence of leucosome prior to mylonitization. The lower Kåfjord is similar to the calc-schists in the underlying Vaddas Nappe, whereas the upper Kåfjord is similar to strongly sheared rocks in the Nordmannvik Nappe (e.g. Figs. 3, 5a). Both units are pervasively and strongly sheared at amphibolite facies conditions. The Kåfjord-Nordmannvik boundary consists of mylonitic gneisses and garnet-micaschist at least 50 m thick. The Nordmannvik Nappe is mainly comprised of garnet-biotite-gneiss with layers and lenses of amphibolite and local calc-silicates, generally dominated by a mylonitic foliation (Fig. 5a, b). Further away from the nappe boundary rare relict lenses (< 50 m) display a higher-grade migmatitic foliation (Fig. 5c). On Arnøya, the frequency of relict lenses increases upwards away from the nappe boundary. The migmatite comprises felsic leucosome and biotite, garnet, and kyanite in the restite. The amount of leucosome varies spatially between 5 % and 25 % (Fig. 5d, e). Layers of amphibolite in the nappe also contain minor tonalitic leucosome suggesting that the mafic rocks have also been migmatized.

Structures
The earliest observed structural element is the synmigmatitic foliation observed in the low-strain lenses in the Nordmannvik Nappe (S 1 ; Figs. 2a, e, 5c-e). In most cases it is sheared and overprinted by the solid-state amphibolite-  Table 1. facies S 2 mylonitic foliation, which is pervasive throughout all the nappes of the RNC and underlying KNC ( Fig. 2a-f). In the lower strain lenses the S 1 migmatitic foliation is found in three structural orientations: (1) as a steep foliation that shows a variable trend (rare; Fig. 5c); (2) folded in open to closed folds with axial planes (containing leucosome) parallel to the S 2 foliation (Fig. 5d); and (3) most commonly as a shallowly dipping foliation, parallel to the S 2 foliation (Fig. 2a,d,e). Fold hinges in the migmatites (F 2 ) plunge towards the SE or NW, parallel to the L 2 stretching lineation. An L 2 intersection lineation between S 1 and S 2 is observed in the Nordmannvik and upper Kåfjord nappes. On Arnøya it is mostly parallel to the L 2 stretching lineation, whereas on Uløya it also plunges shallowly towards the SW and W as well, consistent with the presence of a variable S 1 foliation in the upper Kåfjord Nappe (Fig. 2a, d, e).
The amphibolite-facies solid-state S 2 foliation is generally flat-lying, dipping slightly to the west or north-west in the northern and western parts of the area, and toward the east and south-east in the eastern part of the field area. The main variation is around the larger gabbro bodies, where the S 2 foliation is deflected around them (Fig. 2). Gabbro interiors are nearly unaffected by S 2 shearing, with most of the deformation confined along their boundaries. The S 2 foliation is associated with a generally NW-SE-trending stretching lineation (L 2 ; Fig. 2) and always top-to-SE shear sense indicators ( Fig. 4b, f, g). In the Vaddas and lower Kåfjord metasediments local micro-to meso-scale isoclinal folds transpose a pre-S 2 foliation. Dismembered fold hinges are common, with fold axes normally parallel to the L 2 stretching lineation and axial planes parallel to S 2 ( Fig. 2a, b, d, e). The KNC-Vaddas, Vaddas-Kåfjord, and Kåfjord-Nordmannvik nappe boundaries generally display a relatively stronger S 2 mylonitic foliation than within the nappes, and the L 2 stretching lineation at the nappe boundaries is often more pronounced than within nappe interiors. We associate S 1 in the RNC with www.solid-earth.net/10/117/2019/ Solid Earth, 10, 117-148, 2019 D 1 and S 2 , and L 2 with D 2 . The S 1 in the KNC may be older than the S 1 in the RNC (pre-D 1 ).

Tectonostratigraphic interpretation
Based on our field observations we propose a correlation of the different parts of the RNC across the study area (Fig. 3).
Lithological similarities indicate that the basal marble and conglomerate layer of the Vaddas Nappe on Arnøya, Kågen, and Uløya correlates with the uppermost marble and conglomerate layer at the base of the upper Vaddas Nappe around Straumfjord at a similar tectonostratigraphic level to where Binns and Gayer (1980) identified early Silurian fossils (Fig. 3). This indicates that the lower Vaddas rocks, which are relatively thick in the eastern study area, thin out towards the north and west and are not present elsewhere in the field area (Fig. 3). The lower Vaddas unit includes the Rappesvarre granite, and, based on the structural and litho-logical similarity to the underlying KNC, we favour the interpretation that it is either a composite nappe of KNC and upper Vaddas rocks, or reworked KNC rocks, as suggested by Corfu et al. (2007) and Gee et al. (2017). More work needs to be done to establish whether it has more RNC or KNC affinity (perhaps Corrovarre Nappe) to constrain its extent south of the field area. Gabbros are characteristic of the upper Vaddas Nappe, and the continuation of the Kågen gabbro on southern Arnøya indicates that the upper Vaddas Nappe is present there as well (Figs. 2,3). Our investigation on Arnøya shows an identical nappe stack as present on Uløya, indicating that the north-western half of Arnøya is part of the RNC, including the upper Vaddas, Kåfjord, and Nordmannvik nappes (Figs. 2, 3). The kyanite-bearing migmatites and gneisses on the northern part of Arnøya are assigned to the Nordmannvik Nappe based on their similarity to migmatites at Lyngseidet and on southern Uløya, and their tectonostrati-Solid Earth, 10, 117-148, 2019 www.solid-earth.net/10/117/2019/ graphic position (Fig. 3). We separate the Kåfjord Nappe, based on lithological differences, into upper and lower parts. The upper part is similar to the strongly sheared rocks of the Nordmannvik Nappe, and records a Rb-Sr whole-rock age of ∼ 440 Ma age, possibly for anatexis (Dangla et al., 1978). It could alternatively be considered as part of the Nordmannvik Nappe. The lower part consists of calcareous metasediments that show no evidence of prior anatexis and are more similar to the underlying metasediments of the Vaddas Nappe.
Although the lower part of the Kåfjord Nappe in the field area is comprised of similar metasediments to the Vaddas Nappe (suggesting they may be related), it is defined as its own nappe because its base cuts the upper part of the gabbro in the Vaddas Nappe on Arnøya. All nappes show significant thickness variations, which are partly due to the undeformed Kågen and Kvaenangen gabbros that are large boudins. Erosion has removed all units above them, e.g. the Kåfjord and Nordmannvik nappes (Fig. 3). The entire Vaddas Nappe, as it is currently defined, shows a thickness of 100-1500 m with its lower part wedging out towards the west and north. The Kåfjord Nappe is thickest in its southern part (> 1000 m) and thins towards Arnøya (< 1000 m), whereas the Nordmannvik Nappe is much thicker (> 2000 m) than the Vaddas or Kåfjord nappes and thickest on Arnøya (Figs. 2, 3). This results in flat-lying laterally extensive but lensoid thrust sheets, which are separated by mylonitic shear zones but which are internally pervasively sheared (Fig. 2).

Metamorphism in the RNC
Metamorphic conditions were investigated throughout the entire RNC in order to resolve variations in pressure and temperature (P -T ), and deformation conditions throughout the nappe stack. Sample sites and details are shown in Figs. 2, 3, and Table 1. Table 2 gives a summary of the mineralogy in the samples used for P -T modelling. Mineral abbreviations are according to Whitney and Evans (2010).

Methods
P -T conditions for all samples were estimated using phase equilibrium modelling with the Perple_X software (Connolly, 2005: version 6.6.6) using the internally consistent thermodynamic dataset of Holland and Powell (1998). The calculations were performed in the MnNCKF-MASHTi system using X-ray fluorescence (XRF) wholerock compositions (Table S1 in the Supplement). The following solution mixing models were used: garnet, staurolite, chloritoid (Holland and Powell, 1998), biotite (Tajčmanová et al., 2009), ternary feldspar (Fuhrman and Lindsley, 1988), ilmenite (ideal mixing of ilmenite, geikielite, and pyrophanite endmembers), melt (Holland and Powell, 2001), and white mica (Coggon and Holland, 2002). Iron was assumed to be Fe 2+ as the Fe 3+ content of the miner-als considered is negligible and Fe 3+ oxides occur in negligible amounts. Apatite was observed in all samples and therefore the corresponding amount of CaO bonded to P 2 O 5 observed in the whole-rock analyses was subtracted from the bulk compositions. Measured chemical compositions of the relevant minerals (Tables 3 and 4) were compared with model isopleths in the calculated pseudosections. Molar percent of grossular (Grs) and spessartine (Sps) endmembers and the X Mg (Mg/Mg + Fe tot ) value in garnet are used to constrain P -T conditions, and are shown on the P -T sections. Estimated P -T conditions were checked with anorthite content in plagioclase (An=Ca/(Ca + Na + K)), X Mg (Mg/(Mg + Mn + Fe tot )) in biotite, and Si content in white mica, and are shown on some pseudosections where they assist with further constraining P -T estimates. Fluid content was derived from the loss on ignition (LOI) values unless stated otherwise, and considered as pure H 2 O in all calculations.

Petrography and mineral chemistry
Sample UL248 (Figs. 2, 3, 4b; Table 1) represents a mylonitized migmatitic paragneiss from the upper part of the KNC (Fig. 4a, b). The fine-grained matrix of sample UL248 contains porphyroblasts of garnet (up to 0.9 mm) and fragments (up to 5 mm) of plagioclase, K-feldspar, and quartz-feldspar aggregates (Fig. 4c). Biotite, muscovite, and elongate grains and aggregates of quartz and feldspar define the S 2 mylonitic foliation. Two generations of muscovite are found as (1) rare large (0.2 mm) mica fish parallel to the S 2 foliation (Ms 1 ), and (2) as small grains intergrown with biotite and along garnet and K-feldspar boundaries, within the S 2 foliation (Ms 2 ). Quartz, biotite, and rutile are common as inclusions in garnet (Fig. 4c). Minor rutile is also found in the matrix. Titanite is abundant as elongate, 0.02-0.35 mm long, grains parallel to the S 2 foliation and as inclusions in garnet rims (Figs. 4c,6a). Distinctive garnet zoning displays two generations of garnet. Garnet cores (Grt 1 ) have a relatively flat compositional profile, followed by a transition zone to a ∼ 0.1-0.15 mm thick rim (Grt 2 ; Fig. 6a). The two generations show significantly different compositions. Grt 1 is lower in Grs (Grs 10−12 ), higher in spessartine (Sps 12−15 ), and has a higher X Mg content (0.11-0.13) than Grt 2 . The latter displays the following composition: Grs 36−38 , Sps 3−4 , and X Mg content of 0.08-0.11. Grt 2 (rims) grew together with biotite, titanite, plagioclase (rims), and Ms 2 during top-to-SE S 2associated shearing. X Mg in biotite ranges between 0.36 and 0.45. Plagioclase generally shows zoning with a higher anorthite content in the cores (An 25−28 ) than in rims (An 20−23 ). Early muscovite (Ms 1 ) has a Si content of 3.06-3.08 (a.p.f.u., atoms per formula unit), whereas later muscovite (Ms 2 ) associated with the S 2 foliation has a higher Si content of 3.20-3.27 (a.p.f.u.; Table 4 Grt 1 , plagioclase cores, and large Ms 1 fish indicate that pre-D 2 P -T conditions may be preserved in this sample.

P -T modelling
Garnet cores represent an earlier P -T history, shielded from re-equilibration with the matrix by garnet rims during subsequent metamorphism, and the conditions for garnet core and rim formation were therefore modelled separately. Garnet core formation was modelled using the composition obtained from the bulk rock XRF analysis (Table S1). Since leucosome was observed in less sheared outcrops of the same rock type within 10 m of the sample site, garnet cores in the sample were modelled as part of a migmatite assemblage. The pressure was estimated to be 10 kbar based on the position of garnet isopleths representing measured garnet core compositions in a water-saturated pseudosection. The water content was then estimated from the position of the solidus in a T -X H 2 O pseudosection calculated at 10 kbar. In Fig. 7a modelled X Mg , Grs, and Sps compositions fit measured garnet core compositions within the phase field Grt-Bt-Pl-Kfs-Ms-Rt-Qtz-Melt, constraining pre-D 2 conditions between 705-735 • C and 9.9-10.8 kbar. As diffusion in garnet at upper am-phibolite facies conditions is considered to be fast (e.g. Caddick et al., 2010), care should be taken with interpreting these as absolute P -T conditions as garnet cores could have been modified by diffusion during the overprinting event. However, we interpret the shape of the garnet profile to indicate that diffusion is reflected by the transition zone between Grt 1 cores and Grt 2 rims, and that the flat cores probably represent Grt 1 that escaped diffusion during formation of the S 2 foliation (Fig. 6a). The interpretation is consistent with the presence of rutile and lack of titanite as inclusions in garnet cores (Figs. 4c, 6a). It also agrees well with the anorthite (An 25−29 ) content in zoned plagioclase cores and Si content of large Ms 1 mica fish (Table 4). The P -T estimate for formation of garnet rims and the S 2 foliation (Fig. 7b) was calculated with the XRF bulk rock composition from which garnet cores were subtracted by using the modal Grt 1 proportion estimated from the pseudosection (0.75 modal %). LOI from the bulk rock analysis was used as H 2 O content. Modelled X Mg , Sps, and Grs compositions fit measured garnet rim compositions in the phase field Grt-Bt-Ms-Pl-Kfs-Ttn-Rt-Qtz, and constrain formation of the S 2 foliation between 635 and 690 • C and between 11.5 and 12.3 kbar ( , has a profile that shows some zoning, with mainly Grs and Sps contents displaying a difference between inclusion-rich cores and inclusion-poor rims. rims, X Mg content in matrix biotite, and Si content of Ms 2 agree well with this estimate. Titanite is predicted as part of this assemblage, which is consistent with the large amount of titanite in the matrix and as inclusions in Grt 2 (Figs. 4c, 6a). Based on field structural observations we interpret this P -T estimate for S 2 to represent shearing along the KNC-Vaddas boundary, and titanite growth is directly related to this event.

Petrography and mineral chemistry
Sample AR71 (Figs. 2, 3; Table 1) is a medium-grained garnet-mica-schist with a strong S 2 foliation defined by intergrown biotite and muscovite (Fig. 4f). Muscovite is the dominant mica and occurs as large (1-3 mm long) mica fish. Plagioclase is found as rare porphyroclasts and as single recrystallized grains within quartz layers. Rutile is abundant as elongate grains parallel to the foliation (Fig. 4f). Minor ilmenite occurs as rims on rutile grains. Garnet grains are found as small (0.03-0.1 mm), single, idiomorphic grains and as larger (1-2 mm) fish-shaped clusters of grains paral-lel to the foliation, with a synkinematic geometry displaying a typical Caledonian top-to-SE shear sense (Fig. 4f). Garnet shows an almost flat compositional profile with slight variations between cores, intermediate zones, and rims. In rare cases fish-shaped garnets have a thin rim of a different composition. These rims only occur at the apex of fish-shaped garnets and are never present on the S 2 parallel rims, indicating that they represent growth of garnet during late D 2 shearing (e.g. Fig. 6b, Table 3). In garnets that do not display the thin rim, X Mg varies from 0.17 in the cores to 0.14 in the intermediate zones and garnet edges. Grs is slightly lower in cores (Grs 7−8 ) than intermediate zones and edges (Grs 9−10 ). Sps content varies between 1.4 and 1.6 mol %. Almandine content varies between 68 and 73 mol % (Fig. 6b, Table 3). The thin garnet rims display the following composition: Alm 65−67 Grs 14−16 Sps 1.2 and an X Mg of ∼ 0.11 (e.g. Fig. 6b, Table 3). X Mg in biotite ranges between 0.50 and 0.53, Si content in muscovite is between 3.05 and 3.11 (a.p.f.u.), and plagioclase has an anorthite content of 18-20 mol % (Table 4).
www.solid-earth.net/10/117/2019/ Solid Earth, 10, 117-148, 2019 and Kåfjord (sample AR153). Estimates in the pseudosections are displayed as striped boxes. (a) Garnet core (Grt 1 ) growth conditions were estimated from the X Mg , Grs, and Sps contents of garnet cores. Isopleths for anorthite content are shown, and are in agreement with measured anorthite content in plagioclase cores. The pseudosection was calculated using the bulk composition determined from XRF analysis of the whole rock. (b) S 2 foliation conditions were estimated from X Mg , Grs, and Sps content correlating with measured garnet rims (Grt 2 ). The estimate correlates well with anorthite content in plagioclase rims. The pseudosection was calculated using a bulk composition from which garnet cores were subtracted. (c) Grs, Sps, and X Mg contents in garnet give core and rim at slightly lower and higher pressures, respectively. Si (a.p.f.u) in muscovite was used to further constrain pressures. (d) Grs, Sps, and X Mg contents in zoned garnets from the Kåfjord Nappe give a garnet core estimate at lower pressure conditions, and garnet rim estimate at higher pressure conditions. Chlorite and ilmenite inclusions in garnet cores are consistent with the garnet core estimate.    Modelled X Mg , Grs, and Sps isopleths for cores, intermediate zones, and normal rims corresponding to the observed composition of garnet and biotite intersect in the Grt-Bt-Pl-Ms-Pg-Qtz-Rt-H 2 O phase field (Fig. 7c). This is consistent with the assemblage present in the sample except for paragonite, which was not observed but is predicted to be present only in small amounts (< 2 modal %). The slightly lower Grs and higher X Mg contents of the cores indicate that minor garnet zoning resulted from growth during an increase in pressure (Fig. 7c). Si content in muscovite was used to better constrain the pressure range due to the pressure-insensitive orientation of garnet isopleths in the phase field. Together with garnet compositions, the Si content of muscovite constrains conditions between 630 and 640 • C and between 11.7 and 13 kbar for the formation of the S 2 foliation (Fig. 7c), with cores forming between 11.7 and 12.2 kbar and intermediate zones and normal rims forming at up to ∼ 13 kbar. The higher Grs and lower X Mg contents, and a similar Sps content in the rare thin rims on garnet, together with the growth of ilmenite rims on rutile, indicate retrogression to lower temperature, and probably lower pressure conditions (Fig. 7c).
The lack of an older foliation in the sample or in the outcrop, and the zoning pattern of garnet leads us to interpret that the Vaddas Nappe was only affected by D 2 Caledonian metamorphism and subsequently retrogressed to lower P -T conditions.

Petrography and mineral chemistry
Sample AR153 (Figs. 2, 3; Table 1) is a medium-to finegrained garnet-biotite-zoisite-schist sample taken from the lower Kåfjord Nappe. Biotite, muscovite, zoisite, and elongate quartz-feldspar aggregates define a strong mylonitic S 2 foliation (Fig. 4i). Muscovite generally occurs within the foliation grown together with biotite and as rare grains cross-cutting the S 2 foliation. Garnet is porphyroblastic and idiomorphic and sometimes has inclusion-rich cores and inclusion-poor rims. Cores include quartz, biotite, ilmenite, and chlorite (Fig. 4h). Chlorite is absent in the matrix. Rutile and titanite are both found as elongate grains parallel to the foliation, although rutile is significantly more abundant. Two generations of ilmenite occur as small inclusions in garnet and as thin rims on rutile grains (Fig. 4h, 4i). Garnet shows some compositional change from core to rim ( Fig. 6c; Table 3). Alm content is variable across the garnet (between Alm 53−62 ). Grs content is around 19-21 mol % in the cores and increases towards the rims up to 22-26 mol %. In the profile (Fig. 6c), X Mg does not vary much between cores and rims (0.15-0.18) although spot analyses show that X Mg is often slightly lower in the cores than in the rims (Table 3). The Sps content is higher in the cores (3-5 mol %) than in the garnet rims (1-2 mol %). Plagioclase shows no zoning, and has a composition of An 25 to An 32 . The X Mg in biotite is between 0.56 and 0.59. Si content of the white mica is between 3.16 and 3.2 (a.p.f.u.).

P -T modelling
Zoning in garnet and the difference between garnet inclusion assemblage and matrix assemblage reflect growth at evolving P -T conditions. The presence of chlorite and ilmenite as inclusions in garnet (Fig. 4h, i) indicates that the cores grew at different conditions relative to rims and the matrix, limiting core growth below 620 • C and 8.5 kbar. Garnet cores show a relatively higher Sps content and lower content of Grs than rims while X Mg shows little variation between cores and rims. The sample was modelled under water saturated conditions and the resulting pseudosection is shown in Fig. 7d. The presence of chlorite and ilmenite as inclusions in garnet cores limits core growth conditions below 620 • C and 8.5 kbar. Modelled X Mg , Grs, and Sps isopleths are consistent with measured garnet core composition within the Grt-Bt-Chl-Pl-Ilm-Qtz phase field, and constrains garnet core growth between 590 and 610 • C and between 5.5 and 6.8 kbar (Fig. 7d). The presence of rutile, zoisite, and muscovite in the matrix constrains formation of the S 2 foliation above 9 kbar in the Grt-Bt-Ms-Pl-Zo-Rt phase field, consistent with the matrix assemblage ( Fig. 4i). Modelled X Mg , Grs, and Sps isopleths representing measured garnet rim compositions constrain P -T conditions between 580 and 605 • C and between 9.2 and 10.1 kbar. The morphology of the isopleths is consistent with a lack of significant zoning in pyrope (and X Mg content) in the garnets. A match with model compositions for the rim estimate is also observed with the measured plagioclase composition (Fig. 7d) and measured Si content in muscovite. The X Mg value of biotite in the sample has a large range that is consistent with modelled values for both the garnet core and rim estimates, and which is likely the result of biotite growth and equilibration over evolving P -T conditions. There is no microstructural evidence to suggest that there is more than one generation of biotite. The shape of the garnet compositional profile (Fig. 6c) compared to the morphology of the isopleths on the pseudosection suggests initial garnet growth at low pressures followed by continuous garnet growth with increasing pressure (and slightly increasing temperature). Ilmenite rims on rutile grains are likely a result of later retrogression to lower pressures and temperatures.  as possible so that it contained ∼ 5-10 volume % leucocratic material. The sample has a macroscopic weak foliation defined by biotite and migmatitic banding. The leucosome consists of plagioclase, quartz, and K-feldspar (partly replaced by myrmekite). Kyanite and garnet porphyroblasts (2-3 mm in size) are common along leucosome boundaries and in the restite (Fig. 5f). Kyanite occurs as large (up to 6 mm) porphyroblasts oriented both randomly and parallel to the foliation and almost always in association with biotite, and often with quartz ( Fig. 5f). Kyanite crystals sometimes contain inclusions of quartz and biotite, and are occasionally in contact with garnet. Smaller kyanite crystals are also found as inclusions in quartz in leucosomes. Garnet is usually idiomorphic, and contains inclusions of biotite, quartz, ilmenite, and sillimanite (the latter confirmed by Raman spectra; Figs. 8, 9a). In minor amounts, fine-grained sillimanite is intergrown with biotite or occurs along kyanite and garnet rims. Muscovite occurs as rare small grains overgrowing biotite and is often associated with fine-grained sillimanite. Garnets show a relatively flat compositional profile in the core and gradual changes near the rims, probably indicative of diffusion zoning. X Mg in the cores is 0.23-0.26 and decreases to 0.20-0.23 in the rims. Alm in cores is between 65-69 mol % and 69-71 mol % in the rims. Grs in garnet cores is as low as 4.5 mol % and at the rims 6-7 mol %. Sps content of garnet cores is 3 mol % and 3.5-4 mol % in the rims (Fig. 9a).

P -T modelling
Metamorphic conditions for the migmatitic S 1 foliation in the Nordmannvik Nappe were estimated from sample AR25b (Fig. 10a). To set water content for migmatization, first an approximate pressure for equilibration was estimated from a water-saturated pseudosection using the XRF bulk composition. Isopleths for Grs, Sps, and X Mg contents intersect at ∼ 10 kbar. The water content used in the pseudosection (Fig. 10a) was then determined from the position of the solidus on a T -X H 2 O pseudosection calculated at 10 kbar for the mineral assemblage associated with melting. The pres-ence of kyanite and sillimanite, together with the variation in plagioclase and biotite composition, suggest that the minerals record evolving metamorphic conditions. The flat profile for garnet cores and changes at the rims suggest that the rims may have been modified by diffusion due to disequilibrium (Fig. 9a). The association of kyanite, biotite, and garnet with the leucosome in the S 1 foliation suggests that these minerals formed the stable assemblage during D 1 partial melting. Given the large volume of kyanite associated with the leucosome it can be explained by the following reaction: Ms ± Pl + Qtz = Kfs + Als + Liq. This is consistent with the lack of muscovite in the rock and presence of Kfeldspar. Based on the relationship between garnet rims and the leucosome, garnet rim compositions were taken as representative of the D 1 migmatization conditions. Most measured garnet rim compositions and some intermediate zone compositions correspond with modelled isopleths in the phase field Grt-Bt-Pl-Kfs-Ky-Qtz-Melt on the pseudosection (Fig. 10a). X Mg , Grs, and Sps contents from these analyses constrain D 1 partial melting between 760-790 • C and 9.4-11 kbar. Several measured garnet rims also correspond with modelled isopleths below the solidus within the muscovite-bearing assemblage, consistent with cooling, late muscovite growth, and late minor diffusional resetting of garnet rims. The minor fine-grained sillimanite with biotite and along kyanite and garnet rims is consistent with retrogression to lower pressures. The homogenous garnet core compositions were also plotted on the pseudosection (Fig. 10a) and modelled isopleths for Grs, Sps, and X Mg contents correspond with measured contents in the sillimanitebearing phase field (Grt-Bt-Pl-Kfs-Sill-Qtz-Melt) above the solidus between 790-815 • C and 8.9-9.9 kbar (Fig. 10a). These conditions correspond to the temperatures during D 1 melting, but at slightly lower pressures. The range in measured X Mg content of biotite corresponds with isopleths that plot consistently with both estimates, whereas the measured anorthite content for matrix plagioclase corresponds with isopleths that plot toward the lower pressure conditions (sillimanite-bearing assemblage). It should be noted that the conditions estimated for garnet core formation may be affected by re-equilibration of the cores by later diffusion. Considering this, and given the presence of sillimanite inclusions in garnet, it is possible that garnet core conditions could have been up to ∼ 2-3 kbar lower pressure and ∼ 50 • C hotter than the predicted estimate based on their current composition.  Table 1), which comes from just above the Nordmannvik-Kåfjord boundary, is a finegrained garnet-kyanite micaschist that displays a pervasive mylonitic S 2 foliation defined by biotite, muscovite, and kyanite (Fig. 5b). It has a strong L 2 stretching lineation de-Solid Earth, 10, 117-148, 2019 www.solid-earth.net/10/117/2019/ fined by quartz aggregates. Plagioclase locally forms 0.1-0.3 mm sized sigma clasts. Biotite occurs as two generations: less common large foliation-parallel grains, and as finergrained elongate grains intergrown with muscovite. Muscovite occurs as mica fish parallel to the foliation, and as rare late grains that cross-cut the foliation (Fig. 5b). Garnet is found as large 1-2 mm idiomorphic porphyroblasts (Fig. 9b) that sometimes have inclusion-poor cores (with occasional sillimanite) and more inclusion-rich rims (quartz and biotite inclusions), and as small 0.1-0.4 mm idiomorphic grains that occur singularly (Fig. 9c) or as clusters oriented parallel to the S 2 foliation. The bimodal garnet grain size is also observed in hand specimens and is therefore not an artefact of sectioning. The larger garnets (Fig. 9b) display a slight difference in composition relative to the smaller ones (Fig. 9c). Cores in the larger garnets have slightly lower Sps and Grs contents (Sps 1−1.5 , Grs 4−5 ) than cores in the smaller garnets (Sps 2.5−3.5 , Grs 5−6 ). Almandine content in the larger garnet cores (Alm 71−76 ) is higher than in the smaller garnet cores (Alm 65−72 ; Table 3). The larger garnets display a relatively flat profile with thin rims, while the smaller garnets exhibit more pronounce zoning (Fig. 9b, c). Inclusion-poor cores of the larger garnets appear microstructurally distinct from inclusion-rich rims; however, Grs is the only endmember that reflects this difference (Fig. 9b, grey box). The larger garnets have a thin 10-20 µm wide rim ( Fig. 9b; dark-grey box) with lower X Mg values and Grs contents than cores, whereas smaller garnets have a similar zoning profile, but lack the large flat cores. Based on the similarity in rim compositions, and the lack of a chemical zoning following the microstructure (inclusion free cores in large garnets), it appears that the slight compositional variations in both the large and small garnets probably resulted from diffusion during high-grade metamorphism. The smaller garnets are probably more completely re-equilibrated than the larger ones. The two generations of biotite have different compositions. Older biotite has an X Mg between 0.56 and 0.58, closer to the composition of inclusions in garnet (X Mg ∼ 0.6). Younger biotite displays X Mg values between 0.46 and 0.48. Plagioclase is un-zoned with a composition of 23-24 mol % (Table 4).
www.solid-earth.net/10/117/2019/ Solid Earth, 10, 117-148, 2019 Figure 10. Pseudosections for the Nordmannvik Nappe. (a) Measured garnet core and rim compositions for migmatite (AR25b) agree with modelled compositions that plot above the solidus. Garnet cores give an apparent estimate in the sillimanite stability field. Garnet rims record S 1 foliation-related partial melting in the kyanite stability field. (b) Measured compositions of garnet cores and rims for S 2 mylonitic foliation (AR26) are shown. Garnet cores give an apparent estimate above the solidus in the sillimanite stability field (S 1 foliation -melting). Garnet rims and small garnets give an estimate below the solidus for the S 2 foliation.

P -T modelling
The variations in mineral compositions likely reflect changing P -T conditions with time. The profile shape and compositions of the larger and smaller garnets suggest that both probably started with the same homogenous composition and smaller garnets re-equilibrated more completely than large ones by diffusion during D 2 metamorphism and shearing (e.g. Caddick et al., 2010). Therefore, the cores of larger garnets may reflect the original P -T conditions. Biotite, kyanite, and muscovite are in contact with garnet rims and define a clear S 2 foliation associated with subsolidus shearing. X Mg , Sps, and Grs compositions of the smaller reequilibrated garnets and some larger garnet rims agree with modelled isopleths in the Grt-Bt-Pl-Kfs-Ms-Ky-Qtz phase field, and constrains the formation of the S 2 foliation between 680-730 • C and 9.5-10.9 kbar (Fig. 10b). Measured X Mg values for the younger biotite generation and Si content in muscovite agree well with their respective modelled isopleths for this estimate. The Grs, Sps, and X Mg contents of larger garnet cores (e.g. Fig. 9b, light-grey box) are consistent with modelled isopleths that plot within the Grt-Bt-Pl-Kfs-Ky-Qtz-Melt and Grt-Bt-Pl-Kfs-Sil-Qtz-Melt phase fields between 790-810 • C and 8-10.4 kbar (Fig. 10b). Although garnet core compositions have been affected by later diffusion and the rock has been strongly overprinted by the S 2 foliation, the estimated core conditions are consistent with those recorded by the migmatite sample (AR25b; Fig. 10a). This indicates that garnet cores record the same pre-D 2 history (D 1 melting and/or possibly earlier) reflected in the migmatites. The garnet core estimates agree with the presence of sillimanite inclusions in the cores, the measured X Mg content of biotite inclusions in garnet, and anorthite content of plagioclase sigma clasts in the matrix, suggesting that garnet is not the only mineral that records this earlier P -T evolution.

Geochronology
In order to temporally constrain the metamorphic and deformation history of the RNC, zircon and titanite from several tectonic levels were selected for geochronology. The methods are described in supplement S1 and data are compiled in Supplement Tables 1-3.  Table 1) was taken from an S 1 lens in the Nordmannvik Nappe with clear leucosome segregations associated with kyanite-bearing restite. The sample was cut in an attempt to obtain different zircon populations for the restite and leucosome layers. Zircons from the restite are clear, euhedral, and prismatic crystals 100-350 µm long. All grains have low cathodoluminescence-emission (CL, dark grains) and show weak oscillatory zoning (Fig. 11a). Some grains have brighter rounded cores with discrete boundaries, indicative of a possible xenocrystic origin. Zircon from the leucosome is slightly larger than the grains from the restite. The grains are clear or yellow-orange in colour, and form euhedral prismatic crystals 150-350 µm long. Under CL they are equally dark and show similar weak oscillatory zoning to the zircon from the restite, but inherited cores are absent (Fig. 11d). For the restite, 37 spots were analysed in 26 grains. Five of the analyses are from the bright cores. These cores have Th/U ratios between 0.3-0.53. Two of the core analyses plot on concordia, two plot close by, and one is strongly discordant (Fig. 11b, grey ellipses). The discordant analysis is likely due to mixing of a CL-bright core and darker rim. Of the four analyses that plot on or near concordia, the oldest three give 207 Pb/ 206 Pb dates of 1836 ± 10, 1721±10, and 1599±14 Ma. The fourth core gives a younger 206 Pb/ 238 U date of 585 ± 8 Ma. The remaining 32 analyses are from regions in zircon grains that show low CL-emission and weak oscillatory zoning. They have Th/U ratios between 0.09 and 0.16, overlap on concordia, and give a combined concordia age of 441 ± 2 Ma (Fig. 11c). For the leucosome, 27 spots were analysed in 18 grains from both light and dark CL zones in the zircon crystals. They have Th/U ratios of 0.08 to 0.19, and form a cluster on concordia giving a concordia age of 439±2 Ma (Fig. 11e). The concordia age for the leucosome is slightly younger than, although within error of, the restite concordia age. The oscillatory zoning patterns in the zircons from both the restite and leucosome populations and the shape of the zircons suggest they have crystallized from melt (e.g. Hoskin and Schaltegger, 2003).

Metamorphism in the Nordmannvik Nappe -U-Pb dating of titanite
Sample A01 is a coarse-grained felsic rock from a lowstrain lens in the Nordmannvik Nappe on Arnøya (Figs. 2, 3, 11f, g; Table 1). The sample represents an interaction zone between melt and calc-silicate rocks and is composed of K-feldspar, quartz and plagioclase, and randomly oriented mafic schlieren that include mainly clinopyroxene, amphibole, plagioclase, biotite, quartz, and titanite with minor calcite, ilmenite, magetite, and iron sulphides. Only a weak macroscopic S 2 foliation is present (Fig. 11f). Titanite is abundant and coarse grained (up to 2 mm in size: Fig. 11g). It occurs as large euhedral grains mainly in the boundary zone between the schlieren and the felsic zones and is either associated with magnetite and quartz or as inclusions in the rims of amphibole grains. Sometimes it also occurs as thick rims on ilmenite. These microstructural relationships indicate it formed from interaction between melt and the calc-silicate rocks. Six clear-brown 300-600 µm sub-to euhedral titanite fragments were analysed using the thermal ionization mass spectrometry (TIMS) method. All analyses plot on concordia and four of them overlap (in bold; Fig. 11h) giving a combined concordia age of 432 ± 1 Ma. An older grain with a 206 Pb/ 238 U age of 440 ± 4 Ma and a younger grain with a 206 Pb/ 238 U age of 428 ± 1 Ma do not overlap with these analyses and these were not involved in the calculation of the concordia age. Together, these ages indicate a range between 440 ± 4 and 428 ± 1 Ma, with a peak at 432 ± 1 Ma.

Gabbro intrusion in the Vaddas Nappe -U-Pb SIMS dating of zircon
Sample SK18b was taken from a late-stage pegmatitic gabbro lens within the main medium-grained Kågen gabbro in the Vaddas Nappe (Figs. 2, 12a). It is coarse-grained and consists predominantly of up to 1 cm long amphibole crystals in a plagioclase matrix. Zircons from the sample are euhedral, short-prismatic, and between 100-300 µm long. CL zonation is variable. Twenty spots were analysed in 16 grains. Th/U values are between 0.3 and 0.7, consistent with zircon crystallization from a melt. All analyses are concordant giving a concordia age of 439 ± 1 Ma (Fig. 12c), interpreted as the intrusive age for late-stage gabbroic pegmatites within the Kågen gabbro.

S 2 shearing at the Vaddas-Kalak boundary -U-Pb titanite ages
Microstructural observations and P -T modelling show that titanite crystallization in sample UL248 is associated with the S 2 foliation during garnet rim growth (Figs. 4c, 6a). Titanite from the sample was therefore dated using TIMS. The majority of picked titanite grains have inclusions and the three cleanest fragments were chosen for analysis. They are pale-brown 240-300 µm long inclusion-free grains. All three analyses plot on Concordia with 238 U/ 206 Pb dates between 436 ± 4 and 431 ± 2 Ma (Fig. 12d). The mean 238 U/ 206 Pb age is calculated as 432 ± 6 Ma. We interpret this age to represent the age of D 2 shearing and metamorphism along the KNC-Vaddas boundary, and the likely emplacement age of the Vaddas Nappe over the KNC.

Potential Neoproterozoic metamorphic event in the KNC
Clear evidence for an early metamorphic event is preserved in sample UL248 from the KNC below the Vaddas-KNC boundary. Grt 1 (cores) formed between 705-735 • C and 9.9-10.8 kbar above the solidus (Fig. 7a) at higher temperature and lower pressure than Grt 2 (rims) (Figs. 7b, 13a). Sheared and dismembered leucosome suggests that melt solidified after an early melting event and was subsequently sheared in the solid state during the D 2 overprint. Age constraints for this early melting event are lacking. However, the metamorphic conditions for Grt 1 are similar to those recorded for partial melting in the KNC at Eide nearby (∼ 730-775 • C and ∼ 6.3-9.8 kbar; Gasser et al., 2015;Fig. 13a, b). The Eide melting event is dated at 702 ± 5 Ma, indicating that it is of pre-Caledonian age (Gasser et al., 2015), and Caledonian migmatization has not been recorded in the KNC so far. We therefore suggest that Grt 1 in sample UL248 records the same pre-Caledonian partial melting recognized at Eide (Fig. 13b). Our results indicate that S 1 in the KNC is significantly older than S 1 in the RNC. The S 1 -forming event in the KNC was termed D 1 by Gasser et al. (2015). It is not the same D 1 as we describe here in the RNC. In this paper we consider metamorphism in the KNC as pre-D 1 relative to the RNC.

Migmatization in the Nordmannvik Nappe -D 1
In the Nordmannvik Nappe early high-grade metamorphism is preserved as a S 1 migmatitic foliation in lower strain lenses. Phase equilibrium modelling of this foliation reveals that garnet cores in samples AR25b and AR26 reflect hightemperature mid-to low-pressure partial melting in the sillimanite and kyanite stability fields. The relationship between garnet cores and rims in migmatite sample AR25b and what the conditions estimated using the garnet core compositions might represent is difficult to interpret. The rim compositions are consistent with kyanite-present melting (ca. 760-790 • C and 9-11 kbar), whereas the cores indicate earlier lower pressure, higher temperature conditions (ca. 790-820 • C and 9-10 kbar; Figs. 10a, 13a). The lack of two distinct generations of garnet may suggest that their compositions record a single event. If this were the case, then the estimates constrained from core and rim compositions would record a single migmatization event, beginning at low pressures and continuing to higher pressures. The problem with such a scenario is that the melting reaction is almost completely temperature dependent, and therefore the bulk of melting would have occurred at the lower pressure sillimanite-present conditions with no or little melt produced at higher pressure. This is in contrast to the microstructures that suggest that substanwww.solid-earth.net/10/117/2019/ Solid Earth, 10, 117-148, 2019 tial partial melting produced mainly kyanite. This scenario is only possible if earlier melt formed under sillimanite-grade conditions removed from the system and melting continued under kyanite-grade conditions. However, due to high rates of diffusion in garnet that prevail at the temperatures under consideration during the D 1 partial melting and D 2 shearing (e.g. Caddick et al., 2010), we have to consider other possibilities for the metamorphism related to the sillimanite inclusions: (1) subsolidus prograde metamorphism related to D 1 migmatization, or (2) some unknown pre-D 1 (pre-Caledonian) event either above or below the solidus. Garnet cores in sample AR26 also record conditions above the solidus in the kyanite and sillimanite stability fields (ca. 790-810 • C and 8-10 kbar) prior to the D 2 shearing recorded by garnet rims (Figs. 10b, 13a). Polymetamorphism in the Nordmannvik Nappe has been previously related to Precambrian or Ordovician events (Elvevold et al., 1987;Lindstrøm and Andresen, 1992). Leucosome associated with partial melting and formation of the kyanite-bearing S 1 foliation crystallized between 441 ± 2 and 439 ± 2 Ma (Figs. 11a-e, 13b). This indicates early Silurian melting at ∼ 30-37 km depth with a geothermal gradient of 22-25 • C km −1 . The presence of leucosome veins in the axial planes of Caledonian (F 2 ) folds in the Nordmannvik Nappe (Fig. 5d) indicates that the rocks were partially molten during folding. The rocks are also overprinted by the solid-state S 2 shear foliation (Fig. 5a, c). Considered together, these structures indicate that the Nordmannvik Nappe migmatites underwent initial deformation between 441±2 and 439±2 Ma while still partially molten, and continued shearing as the rocks solidified until the final D 2 overprinting took place in their solid state. This progression represents an anticlockwise P -T path (mainly with temperature decreasing) from D 1 to D 2 conditions (Fig. 13a). There is a large range in titanite dates (440 ± 4 to 428 ± 1 Ma) for the Nordmannvik migmatite sample A01, covering the age range of leucosome formation derived from zircon (between 441±2 and 439±2 Ma), but also extending to much younger dates. This range could be the result of either (1) diffusional Pb loss during cooling above the closure temperature of titanite (e.g. Tucker et al., 2004), or (2) protracted titanite growth (Kohn, 2017). The closure temperature of titanite is debated (e.g. Cherniak, 1993;Frost et al., 2001), although recent well constrained work in the Western Gneiss Region (WGR) suggests Pb diffusion in titanite is slow (Spencer et al., 2013;Kohn et al., 2015) and probably ineffective below 800 • C Solid Earth, 10, 117-148, 2019 www.solid-earth.net/10/117/2019/ (Kohn, 2017). The large grain size (> 300 µm) of the titanites in sample A01, and estimates for both D 1 and D 2 P -T conditions < 800 • C suggest it is unlikely that diffusional Pb loss caused the large age range, and from this we conclude that titanite probably grew over a protracted time period recording long-lived metamorphic processes. The association of titanite and amphibole or titanite, magnetite, and quartz suggest that it formed by rehydration (1: Cpx + Ilm + Qtz + H 2 O = Amph + Ttn) and/or oxidation (2: Cpx + Ilm + O 2 = Ttn + Mag + Qtz) reactions (e.g. Kohn, 2017) during crystallization of melt within calc-silicate layers. Both reactions are consistent with the observation of titanite rims on ilmenite grains. If both reactions were active, then prolonged titanite growth from ∼ 444 to 427 Ma in the Nordmannvik rocks probably resulted from switching between one reaction to the other during evolution of the system, controlled by H 2 O and O 2 availability from injected melt and metamorphic fluids. Due to the chaotic migmatite structures the rocks do not demonstrate a clear S 2 shear foliation, but they record long-lived titanite growth probably driven by melt and/or fluid infiltration during prolonged early and peak Caledonian metamorphism. Both the zircon and titanite ages are consistent with the 439 ± 1 Ma metamorphic age (from zircon overgrowths) reported from the Nordmannvik Nappe at Heia (Augland et al., 2014).

Gabbro intrusion in the Vaddas Nappe
Coeval with the D 1 migmatization in the Nordmannvik Nappe, late-stage gabbroic pegmatites intruded the Vaddas Nappe on Kågen at 439 ± 1 Ma (Fig. 12c) at pressures between 7-9 kbar (Getsinger et al., 2013), corresponding to a depth of 26-34 km (Fig. 13). Getsinger et al. (2013) established the P -T range based on the presence of kyanite in syn-tectonic pegmatites related to gabbro intrusion, and on the composition of igneous biotite and metamorphic hornblende and plagioclase. The temperature range for intrusion is rather large (650-900±50 • C), but the pressures are rather well constrained by the presence of zoisite at lower temperature shearing and kyanite in pegmatites (7-9 kbar). The data indicates a depth of intrusion of ∼ 26-34 km. The gabbroic composition clearly indicates a mantle source of the melt, and the tholeiitic composition, established from other mafic rocks and gabbros in the nappe (Lindahl et al., 2005), suggests melting in an extensional setting. The similar small gabbro bodies in the Kåfjord Nappe are probably also related to this event. The intrusion of the Kågen gabbro at 439±1 Ma is slightly earlier than intrusion of the Heia gabbro in the Nordmannvik Nappe (435 ± 1 Ma; Augland et al., 2014).

Early garnet growth in the Kåfjord Nappe
Garnet cores in the sample from the lower Kåfjord Nappe (AR153) preserve lower amphibolite facies conditions (ca. 600 • C and 6-7 kbar; Figs. 7d, 13a). The P -T estimate is consistent with a slightly elevated geothermal gradient of 24-29 • C km −1 . Without direct age dating it is difficult to give the exact timing of this event. The lack of pre-D 2 structures in the lower Kåfjord metasediments and likely Late Ordovician depositional age constrains the time window for this metamorphism to either immediately prior to D 2 shearing or during early D 2 shearing (Fig. 13b). The P -T ratio indicates a similarity to the geothermal gradients estimated for D 1 migmatization and gabbro intrusion, suggesting that initial garnet growth in the Kåfjord Nappe may record the early Silurian heating event as in the Vaddas and Nordmannvik nappes, but at a shallower depth. Continuous garnet growth records the increase in pressure during D 2 shearing and eventually peak D 2 metamorphism.

Pervasive D 2 shearing
Pervasive amphibolite-facies D 2 shearing with top-to-SE kinematics is recorded throughout the entire RNC. KNC-Vaddas boundary rocks record D 2 shearing at ca. 630-690 • C and 11-12 kbar at similar conditions as D 2 shearing in the Vaddas Nappe (ca. [630][631][632][633][634][635][636][637][638][639][640]Figs. 7c,13a), at depths around ∼ 43-46 km with the P -T ratio giving an estimated geothermal gradient of 14-15 • C km −1 . Garnets in the Kåfjord Nappe show an increase in pressure from 6 to 10 kbar at temperatures between 580 and 610 • C, reflecting a change in geothermal gradient from early high temperature and low pressure conditions to high pressure and low temperature conditions during D 2 shearing (Figs. 7d, 13a). In the Nordmannvik Nappe D 2 shearing is recorded in the mylonitic gneisses near the Kåfjord-Nordmannvik boundary at ca. 680-730 • C and 9-11 kbar (Figs. 10b, 13a). Although it is not possible to constrain uncertainties for individual P -T estimates, considering maximum possible errors, based on geological uncertainties, of ±50 • C and ±1 kbar (e.g. Palin et al., 2016), the P -T estimates for D 2 in the three nappes are similar, but with a trend towards the highest pressure at the base of the RNC, and the coldest temperature in the Kåfjord Nappe in the middle of the nappe stack (Fig. 13a). The age of D 2 shearing along the KNC-Vaddas boundary is constrained by syn-S 2 titanite to 432 ± 6 Ma (Figs. 12d, 13b), recording stacking of the Vaddas Nappe over the KNC. The age of D 2 metamorphism is consistent with ages for Caledonian metamorphism in the equivalent Magerøy Nappe (Andersen et al., 1982;Corfu et al., 2006Corfu et al., , 2011Kirkland et al., 2005Kirkland et al., , 2016, Narvik Nappe Complex (Augland et al., 2014) and upper KNC (Kirkland et al., 2007a;Gasser et al., 2015).

An anticlockwise P -T path for Caledonian metamorphism in the RNC
The Vaddas, Kåfjord, and Nordmannvik nappes all display an increase in pressure and/or decrease in temperature between pre-D 2 or early-D 2 high-temperature and low-pressure conditions and peak D 2 conditions of higher pressure and www.solid-earth.net/10/117/2019/ Solid Earth, 10, 117-148, 2019 lower temperature (Fig. 13a), recording a clear anticlockwise P -T path. Gabbro intrusion is syn-D 1 during high temperature and low pressure conditions at ca. 439 Ma. Although the Vaddas and Kåfjord metasediments do not exhibit partial melting (with the exception of locally around gabbro bodies), they also record this heating event. Anticlockwise metamorphism is recorded as the Nordmannvik Nappe migmatites solidified, between ca. 440-430 Ma, with an increase in pressure and decrease in temperature during Caledonian shearing (Figs. 13,14). Although D 2 pressures in the Kåfjord Nappe are not as high as those in the Nordmannvik and Vaddas nappes, the shift of the Kåfjord Nappe from an environment with a higher geotherm towards one with a lower geotherm (Fig. 13a) is also consistent with an anticlockwise P -T path, recording crustal thickening and subduction of pre-heated rocks (but beginning in a shallower crustal position).
7 Discussion: the tectonic evolution of the RNC

Depositional ages
In order to understand the age-metamorphic-deformation relationships between the different nappes of the RNC and to reconstruct their tectonic evolution, the depositional ages of the metasedimentary rocks within the different nappes have to be known. As some of the metamorphic temperatures are so high that potential basement cover relationships cannot be reconstructed, the distinction between different nappes becomes difficult and can only be made on the basis of age dating or lithological associations. The Vaddas and Kåfjord samples are metasedimentary rocks that so far have not shown any evidence for a pre-Caledonian metamorphic or deformation history. The Magerøy Nappe, considered equivalent to the Vaddas Nappe, was sedimented in the earliest Silurian based on fossil evidence and an age of 438±4 Ma for volcaniclastic rocks, indicating deposition almost synchronously with emplacement of mafic-ultramafic and felsic magmas, and implying rapid burial of sediments (Andersen, 1981;Kirkland et al., 2005Kirkland et al., , 2016. Although no depositional ages for the Vaddas rocks associated with the dated gabbro have yet been obtained, the presence of Late Ordovicianearly Silurian fossils at the base of the upper Vaddas metasediments south of the field area (e.g. Binns and Gayer, 1980) suggests that the upper Vaddas and possibly lower Kåfjord nappes in our study area are Late Ordovician or younger in age (Fig. 13b). The presence of sheared conglomerates at the base of the upper Vaddas Nappe (Figs. 3, 4d;Lindahl et al., 2005) suggests it was deposited on a continental basement; however, it is unclear what this basement was. The similarity of the metasediments in the lower Kåfjord Nappe, upper Vaddas Nappe and early Silurian volcaniclastic rocks of the Magerøy Nappe suggests that they might represent different stratigraphic and/or distal parts of the same Late Ordovicianearly Silurian basin.
The upper Kåfjord rocks, displaying evidence for an S 1 foliation and an age of ca. 440 Ma for anatexis (Dangla et al., 1978) are different from the lower Kåfjord metasediments, and probably represent strongly sheared Nordmannvikderived rocks. We therefore consider that the Kåfjord Nappe is probably comprised of sheared Nordmannvik-derived rocks in its upper part and sheared Vaddas-derived rocks in its lower part; however, the pervasive nature of S 2 foliation makes this difficult to define (Fig. 14b). The tectonostratigraphic position of the Nordmannvik paragneisses and their older zircon cores suggest they could also have been part of a more outboard basement to the Vaddas and Kåfjord metasediments (Fig. 15a-c). The protolith age of the Nordmannvik samples, however, is still unresolved. The presence of marble and calc-silicates suggests they were deposited as shallow water sediments, probably on a continental shelf. The oldest inherited zircon cores (ca. 1800-1600 Ma) may represent detrital ages from the original source rock, and are typical intrusive ages in northern Baltica (e.g. Larson and Berglund, 1992). The youngest inherited core age of 585 ± 8 Ma represents a maximum depositional age for these rocks (Figs. 11a,13b). Farther to the south-west, the Narvik Nappe Complex, at the same tectonostratigraphic level, contains amphibolitefacies marbles with a chemostratigraphic depositional age of 610-590 Ma (Melezhik et al., 2014), indicating that the Nordmannvik rocks could indeed represent Late Neoproterozoic deposits.

Early Silurian (∼ 440 Ma) metamorphic and magmatic event
Prior to the development of top-to-SE S 2 foliation and shear zones, the KNC was the most inboard terrane relative to Baltica, with the lower Vaddas, upper Vaddas, Kåfjord, and Nordmannvik nappes increasingly more outboard (Figs. 14,15). Based on pressure estimates for conditions prior to D 2 shearing, each nappe was also at different depths during this time, with the Kåfjord and Vaddas metasediments at midto lower-crustal levels, respectively, and the Nordmannvik rocks at high temperatures and pressures (lower-crustal level) during early Silurian heating (Fig. 14a)

Similar early Silurian metamorphic and magmatic events in the Scandinavian Caledonides
Early Silurian migmatization in the Nordmannvik Nappe is similar to that observed in the Seve Nappe Complex (SNC), considered part of the outermost Baltica margin (continentocean transition zone). In the SNC, Ordovician high-pressure metamorphism (ca. 460 to 450 Ma) was followed by early Silurian (ca. 445-435 Ma) granulite facies metamorphism and magmatism (Gromet et al., 1996;Majka et al., 2012;Klonowska et al., 2013Klonowska et al., , 2017 at similar conditions to those recorded in the Nordmannvik Nappe. However, no evidence for Ordovician high-pressure metamorphism has been observed in the RNC (sillimanite inclusions in garnet suggest pre-Silurian high temperatures instead). The current lack of evidence for Ordovician high-pressure metamorphism in the RNC indicates that the RNC and SNC may have been in different tectonic positions relative to Baltica at that time. Intrusive rocks of early Silurian age are widely recognized along the length of the Caledonides (e.g. Tucker et al., 1990;Gromet et al., 1996;Vaasjoki and Sipilä, 2001;Andréasson et al., 2003;Kirkland et al., 2005;Corfu et al., 2011;Majka et al., 2012;Klonowska et al., 2013Klonowska et al., , 2017. In north-ern Norway, they are found in the Narvik Nappe Complex (mafic Rånå intrusion; 437 ± 0.5 Ma; Tucker et al., 1990), the Vaddas or Kalak Nappe (mafic Halti Igneous Complex; 434 ± 5 and 438 ± 5 Ma; Vaasjoki and Sipilä, 2001;Andréasson et al., 2003), and granitic and gabbroic rocks on Magerøy (intruded between 440-435 Ma; Kirkland et al., 2005;Corfu et al., 2011). The ages and tectonostratigraphic relationships established in this work support a correlation of the Narvik Nappe Complex with the Nordmannvik Nappe (e.g. Augland et al., 2014). The Magerøy Nappe has been most closely associated with the Vaddas Nappe (e.g. Andresen, 1988;Corfu et al., 2007). Both have similar lithologies and two-phase histories (D 1 and D 2 in the Magerøy Nappe; Andersen, 1981). Gabbro of the Honningsvåg igneous complex, which intrudes syn-D 1 , has a similar age as the Kågen gabbro in the Vaddas Nappe (Corfu et al., 2006). The Skarsvåg Nappe overlies the Magerøy Nappe, is comprised of migmatitic mica-schist and quarzites, and is intruded by granites synchronously with or after D 1 at 436 ± 1 and 435 ± 2 Ma (Andersen, 1984;Corfu et al., 2006). The similarity in lithology, ages, structures, and tectonostratigraphy supports a correlation of the Magerøy Nappe with the Vaddas Nappe and overlying Skarsvåg Nappe with the Nordmannvik Nappe.
Early Silurian mafic magmatism is also recognized south at Sulitjelma and extensively within the Trondheim Nappe Complex (TNC; Fig. 1a; Pedersen et al., 1992;Nilsen et al., 2007;Slagstad and Kirkland, 2018). The age window of this mafic magmatism is between 438-434 Ma (Slagstad and Kirkland, 2018), fitting with the ages determined for the Kågen gabbro and migmatization in the Nordmannvik Nappe. Along the Caledonian orogen most occurrences are attributed to marginal basin settings, and typically considered as part of the Iapetus-derived rocks. In cases such as the TNC, oceanic rocks are recognized tectonostratigraphically below them (e.g. Andersen et al., 2012) in contrast to the RNC.

Large-scale tectonic models for the RNC
Several different tectonic scenarios have previously been proposed to explain the Ordovician high-pressure metamorphism and early Silurian mafic magmatism in the Scandinavian Caledonides (e.g. Pedersen et al., 1992;Northrup, 1997;Andréasson et al., 2003;Kirkland et al., 2005;Corfu et al., 2006;Slagstad and Kirkland, 2018). However, establishing tectonic models that fit rocks in northern Norway is challenging due to the lack of typical arc rocks, and oceanic rocks marking a suture. The only oceanic rocks present belong to the Lyngsfjellet Nappe, which overlies the RNC. The marked discontinuity in metamorphic grade between the Nordmannvik (amphibolite-granulite facies) and Lyngsfjellet (greenschist-low amphibolite facies) nappes suggests that the Lyngsfjellet Nappe was not thrust in sequence over the Nordmannvik Nappe, and that it did not reach the same depths as the RNC or KNC. The nature of the boundary between www.solid-earth.net/10/117/2019/ Solid Earth, 10, 117-148, 2019 the Nordmannvik Nappe and overlying Lyngsfjellet Nappe is still unclear, and it could also be a later detachment. Regardless, its position above the Nordmannvik Nappe means that it must have been somewhere outboard of the Nordmannvik Nappe rocks, with the overlying (possibly Laurentia derived) Nakkedal and Tromsø nappes representing more exotic elements during the early Silurian ( Fig. 15; Corfu et al., 2003;Janák et al., 2012;Augland et al., 2014). The tectonostratigraphic position of the RNC suggests it formed a part of the outer Baltica margin (e.g. Fig. 15a). Any tectonic model of the Caledonian evolution of the RNC in northern Norway needs to account for the tectonostratigraphy; Late Ordovician-early Silurian sedimentation and volcanism (Vaddas and Kåfjord); rapid burial and intrusion of tholeiitic gabbros (Vaddas) concurrent with migmatization in older, deeper rocks further outboard (Nordmannvik), followed by the onset of D 2 shearing; pressure increase and cooling (anticlockwise P -T path); and the development of a lower-crustal nappe stack. Anticlockwise P -T paths are typically observed where an initial heating event is followed by burial or cooling (Wakabayashi, 2004). Nappe stacking it-self is not considered to be capable of providing the amount of heat required to explain the estimated temperatures for the initial heating event (e.g. Johnson and Strachan, 2006). Ridge subduction (e.g. Northrup, 1997;Corfu et al., 2006) has been mentioned as a possible explanation for the apparent coeval deposition of volcaniclastic sediments and tholeiitic mafic intrusions and the relatively short time window of mafic magmatism. However, the effect of ridge subduction on magmatism is usually very localized (e.g. Lomize and Luchitskaya, 2012) and does not explain the widespread nature of the early Silurian mafic magmatism several hundreds of kilometres along strike of the Caledonian orogeny from the Magerøy Nappe in the north to Trondheim in the south. The widespread regional heating event is best explained by magmatic underplating of mantle melts. The gabbros with ∼ 440 Ma intrusion ages (Kågen, Heia, Magerøy) are evidence for mantle melts at this time. Deposition of volcaniclastic sediments and tholeiitic mafic magmatism is often indicative of an extensional setting, such as a back-arc basin (e.g. Pedersen et al., 1992;Kirkland et al., 2005;Slagstad and Kirkland, 2018), and we consider this setting in several Solid Earth, 10, 117-148, 2019 www.solid-earth.net/10/117/2019/ different tectonic scenarios. Considering the RNC in a backarc position on the Baltica margin would require eastwarddipping subduction (e.g. Fig. 15a; Andréasson et al., 2003). However, the KNC lacks the expected magmatism for this model, and there is no evidence for a switch in subduction polarity immediately prior to collision. For a back-arc environment we therefore have to consider that there is a missing suture below or within the RNC or KNC. Suture zones can be cryptic and difficult to identify in highly deformed rocks, with oceanic rocks typically not present along the entire length of a suture (Dewey, 1977). Similar nappes further south in the Caledonides display evidence for underlying oceanic rocks (e.g. Andersen et al., 2012). Models that consider the RNC as formed in an extensional back-arc setting on the upper plate during westward subduction (e.g. Fig. 15b, c; cf. Slagstad and Kirkland, 2018) provide a good explanation for Ordovician-early Silurian deposition of volcanic rocks and sediments (Vaddas and Kåfjord) and intrusion of mafic rocks (e.g. Kågen gabbro; large-scale magma underplating), and provide an early Silurian heat source for an anticlockwise P -T path. The early collision and slab roll-back model of Slagstad and Kirkland (2018) supports short-lived early Silurian deposition, mafic magmatism and anticlockwise P -T path with rapid switching from extension and heating to compressional D 2 shearing (meltfilled D 2 axial planes in the Nordmannvik Nappe). However, palaeomagnetic data from the Honningsvåg igneous complex on Magerøy suggests it intruded the Magerøy Nappe rocks when they were ∼ 1350 km north of Baltica, supporting a possible early Silurian separation between Baltica and the RNC (e.g. Fig. 15b, c; Torsvik et al., 1992;Corfu et al., 2006). Figure 15b and c offer two different scenarios based on the palaeogeographic position of the KNC. Both consider the Vaddas and Kåfjord having deposited in a back-arc basin on a continental basement. Option B places the KNC as an outboard microcontinent, while option C considers the KNC as part of the down-going plate. In this scenario the edge of the Baltica continent or some continent-ocean transition zone would have undergone subduction, consistent with Ordovician high pressure metamorphism in the SNC (Fig. 15a, b; e.g. Dallmeyer and Gee, 1986;Majka et al., 2014). Both scenarios require the presence of arc rocks either between the RNC and KNC or below the KNC. Traces of these arc rocks have not been found yet. It is possible that subduction was relatively short-lived, providing too little time for a proper arc to develop.
There are several problems that arise when considering the RNC rocks as having formed as part of a back-arc basin. The fairly high pressures of D 1 migmatization in the Nordmannvik Nappe (kyanite bearing; Fig. 13a) and deep-seated Kågen gabbro intrusion are not typical for an extensional environment, and intrusion of gabbro at 26 to 34 km depth into sediments in a back-arc basin requires a very deep basin, especially if it was short lived (e.g. Slagstad and Kirkland, 2018). The P -T path for D 2 shearing in the Vaddas and the near-isothermal increase in pressure from early heating to D 2 in the lower Kåfjord Nappe (Figs. 13a, 14) indicate subduction of pre-heated Vaddas rocks (deeper) and Kåfjord (shallower) as part of the lower plate. Some models discussed here fit better with our observations from northern Norway than others. Eastward-dipping subduction (e.g. Fig. 15a; Andréasson et al., 2003) seems unlikely or, at least, evidence for it is lacking so far. It is likely that the KNC (or at least its upper part) formed a microcontinent on the upper plate (e.g. Fig. 15b), placing the suture below it. The RNC as a microcontinent on the upper plate (e.g. Fig. 15c) is also a promising model, assuming that a yet unrecognized suture is present between the KNC and RNC. Both of these models (Fig. 15b, c) are possible in a scenario where early Silurian heating and gabbro intrusion resulted from extensive magma underplating in a back-arc setting with a rapid switch to shortening and collision over ∼ 10 million years. Discrepancies in the models could be due to uncertainties in observations (e.g. palaeomagnetic datasets or P -T estimates), or due to oblique collision, introducing strike-slip aspects to the story, and a possible overprinted pre-collisional accretion history.
The temperature and pressure evolution of the Vaddas and Kåfjord nappes from ca. 440 to 430 Ma can be used to discuss large-scale nappe stacking. Both indicate a clear pressure increase while temperatures remained similar (Fig. 14). An internal shear zone between the Kåfjord and Vaddas units was active to produce the thrust stacking of these two units (the pressure difference between them becomes more pronounced; Figs. 13a, 14). Both units combined show a pressure increase with respect to the Nordmannvik Nappe, which maintained its depth position (Figs. 13a, 14). This relative pressure development corresponds to thrusting of the Nordmannvik Nappe over the Vaddas and Kåfjord nappes or to initial subduction of these units underneath the Nordmannvik. The Vaddas Nappe was at the same depth as the KNC during the onset of this thrusting, together with the shear zone between them. This scenario suggests an initiation of a subduction zone towards the west with Kåfjord or Vaddas units at its leading edge or thrusting of the Nordmannvik Nappe over Kåfjord and Vaddas during continental collision.

Discussion: nappe stacking in continental collision zones
Even though the early Silurian tectonic history in northern Norway could have resulted from several different scenarios, the P -T evolution of the rocks indicates they formed a relatively hot and weak part of the middle to lower crust (especially the Nordmannvik Nappe) prior to Scandian nappe stacking (D 2 ; Fig. 13a). Additionally, D 2 deformation is pervasive and internal nappe strain can be as high as the strain along nappe boundaries. This is in contrast to typical Alpinestyle nappes (e.g. Escher et al., 1993;Escher and Beau-www.solid-earth.net/10/117/2019/ Solid Earth, 10, 117-148, 2019 mont, 1997), where basement and cover relationships can be established on the basis of metamorphic grade or nappe dividers. It is also in contrast with anorthositic and granulitic Baltica crust in the WGR and Lofoten further south in Norway, which displays only very minor partial melting and little internal deformation, despite being subducted to depths > 100 km (Engvik et al., 2000;Hacker et al., 2010;Froitzheim et al., 2016). The anorthositic and granulitic rocks are not fertile or require large amounts of H 2 O for partial melting. Conversely, the Nordmannvik rocks are fertile and show extensive partial melting. Caledonian deformation started under melt-present conditions (ca. 440 Ma) in a relatively high-temperature environment, then pressures increased and temperatures decreased from D 1 to D 2 , constituting an anticlockwise P -T path. Subsequently, pervasive deformation continued into solid-state conditions (ca. 430 Ma). This suggests that initial high temperatures and associated partial melts played a key role in controlling the style of deformation at greater depth in continental collision zones: rock strength is significantly decreased even with melt fractions as low as ∼ 7%, and melt-bearing systems may facilitate strain localization (e.g. Rosenberg et al., 2005;Cavalcante et al., 2016). In northern Norway, the fertile Nordmannvik rocks are significantly different in composition compared to the dry, strong Baltica basement of the WGR and Lofoten. Partial melting in the fertile subducted Nordmannvik rocks led to strain localization and decoupling of nappe units from the crustal part of the lithosphere, and facilitated pervasive deformation (even in the internal parts of the nappes). Continental crust units are underthrusted and stacked as nappes (shown by increasing pressures in the Vaddas and KNC; Fig. 14b) instead of being subducted. This process destroys the integrity of the lithospheric slab and produces thin crustal nappe slivers. The initial high-temperature metamorphism and partial melting are critical to weakening processes in such a case. The response of the lithology to the metamorphic conditions exerts a key control on the strength of the rocks in a continental collision zone and can decide whether rocks become subducted or not and how they subsequently deform. In this way the structure, composition, and rheology of the lower crust during continental collision determines the geometry and kinematics during subduction of continental crust.

Conclusions
Our reported mineral assemblages, structures, and U-Pb dates from the RNC in northern Norway define a continuous sequence of events: (1) Late Ordovician-early Silurian sedimentation and mafic volcanism in the Vaddas and Kåfjord nappes; (2) early high-temperature partial melting in the Nordmannvik Nappe associated with a migmatitic fabric S 1 and mafic magmatism in the Vaddas, Kåfjord, and Nordmannvik nappes at 441 ± 2 to 439 ± 2 Ma, and 439 ± 1 Ma (D 1 ), respectively; followed by (3) pervasive amphibolite-facies shearing, formation of a subhorizontal S 2 foliation and L 2 lineation, and nappe stacking recorded at ∼ 432 ± 6 Ma (D 2 ). The transition from high-temperature medium-pressure D 1 metamorphism to higher-pressure lower-temperature D 2 metamorphism results in an anticlockwise P -T path. Deformation between these two events is continuous, and S 2 formation began while the rocks were still hot. The tectonic model that best accounts for our observations considers early Silurian heating resulting from extensive magma underplating in a back-arc system, followed rapidly (within 10 million years) by collision, crustal thickening, and nappe stacking. Early Silurian high temperature metamorphism in the RNC (particularly the Nordmannvik Nappe) led to weakening in the rocks, promoting extensive deformation, dismembering of the continental crust, and pervasive D 2 shearing during nappe stacking in early collision.
Data availability. Underlying data for this paper can be found in the supplement.
Author contributions. CF wrote most of the manuscript (with input from all co-authors), conducted ∼ 90 % of the field work, and gathered most of the data (XRF, EMP, SEM, SIMS, and TIMS) with assistance and input from the co-authors. CF, HS, DG, and PJ were involved in project design. Field data from Uløya, and XRF analysis, EMP analysis, and earlier phase equilibrium modelling for the sample from there (UL248) was conducted as part of the MSc project of KK. HS, DG, PJ, JK, and KK participated in the field work and discussions. JK and EKR gave advice on the phase equilibrium modelling, and EKR contributed previous knowledge of the rocks on Arnøya. FC assisted with sample preparation and TIMS analyses and contributed significantly to the discussion and the design of Fig. 15.
and Planets for supporting this work. The publication charges for this article have been funded by a grant from the publication fund of UiT The Arctic University of Norway. Very constructive reviews by Torgeir B. Andersen and Thorsten Nagel have improved the manuscript.
Edited by: Federico Rossetti Reviewed by: Torgeir B. Andersen and Thorsten Nagel