Zircon and monazite reveal late Cambrian/early Ordovician partial melting of the Central Seve Nappe Complex, Scandinavian Caledonides

The Seve Nappe Complex (SNC) comprises continental rocks of Baltica that were subducted and exhumed during the Caledonian orogeny prior to collision with Laurentia. The tectonic history of the central SNC is investigated by applying in-situ zircon and monazite (Th-)U–Pb geochronology and trace element analysis to (ultra-)high pressure (UHP) paragneisses in the Avardo and Marsfjället gneisses. Zircons in the Avardo Gneiss exposed at Sippmikk creek exhibit xenocrystic cores with metamorphic rims. Cores show typical igneous REE profiles and were affected by partial Pb-loss. The rims have flat HREE profiles and are interpreted to have crystallized at 482.5 ± 3.7 Ma during biotite-dehydration melting and peritectic garnet growth. Monazites in the paragneiss are chemically homogeneous and record metamorphism at 420.6 ± 2.0 Ma. In the Marsfjället Gneiss exposed near Kittelfjäll, monazites exhibit complex zoning with cores enveloped by mantles and rims. The cores are interpreted to have crystallized at 481.6 ± 2.1 Ma, possibly during garnet resorption. The mantles and rims provide a dispersion of dates and are interpreted to have formed by melt-driven dissolution-reprecipitation of pre-existing monazites until 463.1 ± 1.8 Ma. Depletion of Y, HREE, and U in the mantles and rims compared to the cores record peritectic garnet and zircon growth. Altogether, the Avardo and Marsfjället gneisses show evidence of late Cambrian/early Ordovician partial melting (possibly in (U)HP conditions), Middle Ordovician (U)HP metamorphism, and late Silurian tectonism. These results indicate that the SNC underwent south-to-north oblique subduction in late Cambrian time, followed by progressive north-to-south exhumation to crustal levels prior to late Silurian continental collision.


Introduction
The Scandinavian Caledonides is an excellent natural laboratory for studying Wilson-cycle tectonics. It comprises several allochthons that altogether record Neoproterozoic opening of the Iapetus Ocean, followed by closure of the ocean starting in late Cambrian time, progressing to Silurian to Devonian collision of Baltica and Laurentia (Gee et al. , 2008(Gee et al. , 2020Stephens and Gee 1985;Stephens 1988;Roberts 2003;Corfu et al. 2014;Stephens et al. 2020). An important element for understanding the evolution of the Scandinavian Caledonides is the Seve Nappe Complex (SNC), which is traditionally regarded to be remnants of the Baltican passive margin that were subducted to mantle depths during Iapetus Ocean closure (Andréasson 1994;Gee et al. 2020 1 3 eclogites, peridotites, and pyroxenites recording (ultra-) high pressure metamorphism (e.g., Janák et al. 2013;Majka et al. 2014a;Gilio et al. 2015;Klonowska et al. 2016Klonowska et al. , 2017Bukała et al. 2018Bukała et al. , 2020aPetrík et al. 2019). It is divided into several (U)HP terranes > 1000 km along strike of the Scandinavian Caledonides, and are often grouped as the northern, central and west-central SNC terranes. The record of (U)HP metamorphism has often been discussed to be older for the northern terranes (late Cambrian/early Ordovician metamorphism; Root and Corfu 2012; Barnes et al. 2021a;Fassmer et al. 2021) versus the central and west-central SNC (Middle Ordovician metamorphism; Brueckner and Van Roermund 2007;Fassmer et al. 2017). The apparently different temporal (U)HP records have led to tectonic models involving diachronous subduction of the Baltican margin, either progressing from north to the south through time, or reflecting a promontory in the north that was subducted prior to the southern portions of the margin Bukała et al. 2018;Fassmer et al. 2021). Further understanding and evaluation of these models is critical for unraveling the evolution of the Scandinavian Caledonides and Iapetus Ocean closure.
Zircon and monazite are two important refractory minerals in metasedimentary rocks that can record a wide range of sub and suprasolidus metamorphic reactions (e.g., Engi 2017; Rubatto 2017). The metamorphic records of the central SNC have predominantly been extracted from (U)HP mafic and ultramafic lithologies hosted within metasedimentary rocks. However, recent studies demonstrate that zircon and monazite often record older metamorphic histories than those preserved in the mafic and ultramafic lithologies (Barnes et al. 2019(Barnes et al. , 2021bPetrík et al. 2019;Walczak et al. 2021). In particular, zircon geochronology of the Tväråklumparna Gneiss in the west-central SNC (Walczak et al. 2021) and monazite geochronology of the Marsfjället Gneiss in the central SNC (Petrík et al. 2019) indicate late Cambrian/early Ordovician (U)HP metamorphism and partial melting of the SNC terranes that were previously thought to have only undergone Middle Ordovician (U)HP metamorphism. Following these discoveries, in-situ zircon and monazite geochronology and trace element analysis are applied to paragneisses of the central SNC that host Middle Ordovician mafic and ultramafic (U)HP lithologies in order to expand the tectonic history of these terranes. The results are compared within the broader history of the SNC to test previous models for subduction of the Baltican margin during Iapetus Ocean closure.
The central SNC is further divided into the Western, Central, and Eastern belts (Fig. 1A;Trouw 1973;Williams and Zwart 1977;Zachrisson and Sjöstrand 1990). The Western Belt is referred to as the Svartsjöbäcken Schists, which are dominated by garnet white mica schists, locally kyanitebearing, which host garnet amphibolites (Trouw 1973). Pressure-temperature (P-T) conditions of the schists were calculated to > 1.5 GPa and 670 ± 30 °C, with garnet Sm-Nd geochronology establishing the timing of subduction to 462.0 ± 3.5 Ma (Grimmer et al. 2015).
The Marsfjället Gneiss constitutes the Central Belt that underlies the Svartsjöbäcken Schists (Fig. 1A). The contact between the units has been described as metamorphic facies transition corresponding to the appearance of ubiquitous kyanite within the rocks (Trouw 1973). The Marsfjället Gneiss consists of garnet-and kyanite-bearing pelitic and quartzo-feldspathic metasedimentary rocks that host garnet amphibolites and ultramafic rocks. Paragneisses of the Marsfjället Gneiss exposed near Saxnäs contain metamorphic microdiamonds included in garnet, providing direct evidence of UHP metamorphism. The timing of metamorphism was interpreted to be 472 ± 3 Ma, provided by monazite Th-U-total Pb geochronology (Fig. 1B;Petrík et al. 2019). In contrast, the paragneisses and amphibolites near Kittelfjäll, ~ 30 km northeast of Saxnäs, demonstrate that the Marsfjället Gneiss underwent partial melting at 474.2 ± 9.6 Ma, yielded by zircon U-Pb geochronology (Bukała et al. 2020b). Orogenic spinel peridotites are exposed close to the Kittelfjäll paragneisses and preserve metamorphic conditions of 650-830 °C and 1-2 GPa during subduction, followed by exhumation of the peridotite to 550-600 °C and 0.45-0.60 GPa (Clos et al. 2014). White mica Rb-Sr geochronology resolved exhumation of the Svartsjöbäcken Schists and Marsfjället Gneiss to 433.2 ± 4.8 Ma (Grimmer et al. 2015).
The Central Belt further consists of the Lillfjället Gneiss, the Avardo Gneiss, and the Lejaren Unit ( Fig. 1A; Zachrisson and Sjöstrand, 1990). The Lillfjället Gneiss and the Avardo Gneiss are both characterized by kyanite-sillimanitegarnet metapelites with quartzites and quartzo-feldspathic rocks and show evidence of partial melting (Van Roermund and Bakker 1983). The two gneisses are distinguished as the latter hosts eclogites and ultramafic rocks, whereas the former contains amphibolites. The Lejaren Unit is predominantly composed of quartzo-feldspathic rocks that host amphibolites and ultramafic rocks. The quartzo-feldspathic rocks generally show evidence of partial melting but apparently do not bear abundant metamorphic index minerals (Van Roermund and Bakker 1983). Of the three, the Avardo Gneiss has been the focus of investigation due to several key localities of eclogite as well as orogenic peridotites and pyroxenites, some of which are garnet-bearing (Van Roermund 1989). Mineral assemblages of the metasedimentary rocks indicate two distinct HP events as two different kyanite generations have been identified (Van Roermund and Bakker 1983). First evidence of UHP metamorphism in the Avardo Gneiss (and in the SNC as a whole) was obtained from an eclogitized portion of a garnet pyroxenite dyke within a garnet peridotite exposure in the Avardo Gneiss northeast of Lake Friningen that yielded conditions of ~ 3 GPa and ~ 800 °C . These conditions were reproduced by the garnet peridotite (Gilio et al. 2015). Garnet Sm-Nd geochronology of garnet pyroxenite yielded 452.8 ± 7.5 Ma, 452 ± 20 Ma and 448 ± 17 Ma, interpreted as the timing of peak metamorphism (Fig. 1B;Brueckner and Van Roermund 2007). Eclogite exposures at Sippmikk creek have yielded minimum P-T conditions of ~ 1.8 GPa and ~ 750 °C (Van Roermund 1985). Garnet Sm-Nd geochronology from eclogite indicates its formation at 454.8 ± 9.7 Ma (Brueckner and Van Roermund 2007), similar to the Friningen garnet pyroxenite. White mica Rb-Sr geochronology of the Avardo Gneiss resolved exhumation to 433.2 ± 4.8 Ma (Grimmer et al. 2015) and zircon U-Pb geochronology applied to leucosomes of partially melted metasedimentary rocks yielded 423 ± 5 Ma ( Fig. 1B; Williams and Claesson 1987).
The Central Belt is separated from the underlying Eastern Belt by a distinct mylonitic zone that represents a tectonic contact (Fig. 1A;Trouw 1973;Van Roermund and Bakker 1983). In the northern exposures, the Eastern Belt is also referred to as the Eastern Schist and Amphibolite Belt. Schists of variable composition dominate with subordinate amphibolite layers and lenses throughout (Trouw, 1973). To the southwest, the Eastern Belt also comprises the Blåsjöälvan Unit, the Gakkafjället Unit, and the Sjouten Unit (Zachrisson and Sjöstrand 1990), which collectively overly the Eastern Schist and Amphibolite Belt (Fig. 1A). The Blåsjöälvan Unit is dominated by amphibolites, often  Williams and Claesson (1987); 7 Dallmeyer and Gee (1988); 8 Fassmer et al. (2017); 9 Root and Corfu (2012) garnet-bearing, in addition to less frequent dolerites, paragneisses and mica schists. Mylonitic amphibolites in the upper parts of the unit in contact with the Avardo Gneiss, exposed near Lake Blåsjön, yielded hornblende 40 Ar/ 39 Ar dates in the range of c. 436-432 Ma (Dallmeyer and Gee 1988). The Gakkafjället Unit and underlying Sjouten Unit are both dominated by quartzo-feldspathic gneisses and garnet mica schists, however, they are distinguished as the former hosts amphibolites and metagabbros and typically preserves amphibolite-facies mineral assemblages, whereas the latter contains numerous eclogites and ultramafic lithologies, including garnet peridotite and pyroxenite, indicating (U)HP metamorphic conditions (Van Roermund and Bakker 1983;Van Roermund 1985, 1989Majka et al. 2014b;Klonowska et al. 2016;Fassmer et al. 2017).
Two key localities have been studied in the Sjouten Unit, at Tjeliken Mountain and on the northern shores of Stor Jougdan. At Tjeliken, eclogites and gneisses preserve P-T conditions of 2.5-2.7 GPa and 650-750 °C (Majka et al. 2014b). Garnet Sm-Nd geochronology of the eclogites yielded 463.7 ± 8.9 Ma (Fig. 1B;Brueckner and Van Roermund 2007), similar to the age of (U)HP metamorphism recorded by the Avardo Gneiss pyroxenite and eclogite. Root and Corfu (2012) questioned the age as their results of zircon U-Pb geochronology of the Tjeliken eclogite produced 445.6 ± 1.2 Ma, which was also reproduced by two individual rutile U-Pb dates of 445.0 ± 2.4 Ma and 446.3 ± 3.7 Ma. Subsequent garnet Lu-Hf geochronology of the eclogite and zircon U-Pb geochronology of the host paragneiss provided 458.1 ± 1.0 Ma and 458.9 ± 2.5 Ma, supporting the previous garnet Sm-Nd results for the timing of UHP metamorphism and suggesting the eclogite zircon and rutile may record exhumation (Fassmer et al. 2017). Conditions for UHP metamorphism are also preserved by eclogite and garnet pyroxenite at Stor Jougdan, which collectively constrain the metamorphic conditions to 2.3-4.0 GPa and 750-960 °C (Klonowska et al. 2016). Garnet Sm-Nd geochronology conducted on garnet peridotite from Stor Jougdan, produced 459.6 ± 4.2 Ma (Fig. 1B;Brueckner and Van Roermund, 2007).

Microscopy and mineral chemistry
Three paragneiss samples were investigated for this study (Fig. 1A). Two were obtained from the Sippmikk creek exposure of the Avardo Gneiss (samples SM19-02, and SM19-03) and one was obtained from a roadcut exposure of the Marsfjället Gneiss northwest of the town of Kittelfjäll (sample MJ18-05A), which was previously studied by Bukała et al. (2020b). All paragneisses were investigated using transmitted light microscopy to characterize mineral assemblages, textural relationships, and microstructures. Information gathered from transmitted light microscopy was used as the basis for electron microprobe (EMP) analysis using a JEOL JXA8230 electron microprobe located at AGH University of Science and Technology, Kraków, Poland. Back-scattered electron (BSE) images were acquired with both high-brightness/low-contrast and low-brightness/ high-contrast settings, the latter utilized specifically to image garnet, zircon, and monazite to observe chemical zoning of the minerals. Semi-quantitative wavelength-dispersive spectroscopy (WDS) major element chemical maps of Ca, Fe, Mg, Mn, and Y were acquired for one representative garnet grain in paragneiss SM19-03. The map was acquired with a beam diameter of 4 μm, an accelerating voltage of 15 keV and beam current of 100 nA with a dwell time of 100 ms. The map was coupled with a WDS profile of garnet, acquiring quantitative data for Na 2 O, SiO 2 , Al 2 O 3 , MgO, K 2 O, CaO, TiO 2 , FeO, MnO, and Cr 2 O 3 . The quantitative WDS analyses were acquired using a 1 μm beam diameter, an accelerating voltage of 15 keV, and a beam current of 20 nA. The following natural minerals and synthetic standards were used for calibration: albite (Si, Al, Na), diopside (Ca, Mg), sanidine (K), rutile (Ti), fayalite (Fe), rhodonite (Mn), and Cr 2 O 3 (Cr). Reflective light (RL) images were obtained of all zircons and monazites that were selected for geochronology and trace element analysis.

Zircon and monazite (Th-)U-Pb geochronology and trace element analysis
Zircon and monazite grains were analyzed in situ within standard 30 µm polished thin sections. The BSE and RL images of the grains were used to select the locations for (Th-)U-Pb geochronology and trace element analysis. Locations were selected to analyze different chemical domains within the zircon or monazite grains and also to avoid fractures, inclusions, domain boundaries, and grain edges. The geochronology and trace element analyses were conducted during four sessions (two for zircon and two for monazite) at the at the Vegacenter, Swedish Museum of Natural History (Stockholm, Sweden). Geochronology was conducted using an ESI NWR193 ArF excimer laser ablation system and a Nu Plasma II multi-collector inductively coupled plasma mass spectrometer (MC-ICP-MS). The m/z (mass-to-charge ratios) corresponding to masses 202, 204, 206, 207 and 208 were measured on ion counters and those corresponding to 232, 235 and 238 were measured on Faraday collectors. For zircon, a beam diameter of 15 µm was used for ablation at a frequency of 7 Hz and laser fluence of 2.0 J/cm 2 . For monazite, a beam diameter of 6 µm was used for ablation at a frequency of 5 Hz and laser fluence of 2.7 J/cm 2 . Helium was used as a sample carrier gas (flow rate of 0.45 L/min for zircon and 0.35 L/min for monazite), which was then mixed with Argon (mixed gas flow rate of 0.80 L/min for zircon and 0.62 L/min for monazite). The data collection procedures included a 20 s ablation time, a 20 s washout time, and a 0.5 s integration time. The data for both zircon and monazite were processed using the Iolite add-on 'VizualAge' data reduction scheme (v. 2.5; Petrus and Kamber, 2012).
The zircon trace element (REE) and monazite trace element (Sr, Y, REE) data were collected during subsequent sessions using an Attom high-resolution inductively-coupled plasma mass spectrometer (HR-ICP-MS) and the same ESI NWR193 ArF excimer laser ablation system that was used for geochronology. For zircon, the laser was operated with a beam diameter of 15 µm at a frequency of 8 Hz and laser fluence of 4.2 J/cm 2 . For monazite, the laser was operated with a beam diameter of 6 µm at a frequency of 5 Hz and laser fluence of 2.7 J/cm 2 . Helium was used as a sample carrier gas (flow rate of 0.55 L/min for zircon and 0.35 L/min for monazite), which was then mixed with Argon (mixed gas flow rate of 0.84 L/min for zircon and 0.83 L/min for monazite). The data for both zircon and monazite were processed using an in-house Excel spreadsheet.
After geochronological and trace element analyses, a second set of RL images of all grains were obtained to show the exact locations of the ablation pits for each zircon and monazite grain. Due to the small sizes of some grains/domains, the new RL images were overlaid on the BSE images to identify ablations that overlapped with grain edges, fractures, inclusions, or that potentially mixed different chemical zones. The resulting datasets were also scrutinized for outliers and were compared with the locations of ablation pits, the combination of which was used as the basis to accept or reject each individual analysis. No blanket criteria were used to filter entire datasets (e.g., discordance cutoffs) to prevent discarding individual data points that are outliers due to geological processes (e.g., Pb-loss). Concordia plots of the accepted analyses were constructed using IsoplotR (v. 4.1) with the default decay constants and 238 U/ 235 U value (Vermeesch, 2018). The reported error correlations (rho) for 207 Pb/ 235 U and 206 Pb/ 238 U were used for the Wetherill concordia (for zircon) and Th-U-Pb concordia (for monazite; Vermeesch, 2020) diagrams. Error correlations between other ratios were not used to create the Th-U-Pb concordia diagrams of the monazite geochronological data. All uncertainties are reported at the 2σ level. Rare earth element plots were created according to McDonough and Sun (1995).

Petrography of the Sippmikk paragneiss
The two samples of the Sippmikk paragneiss (SM19-02 and SM19-03) are predominantly composed of quartz, plagioclase, biotite, white mica, garnet, and kyanite with lesser volumes of K-feldspar, zircon, monazite, allanite, apatite, clinozoisite, rutile, magnetite, pyrite, and carbonate (Figs. 2,3,4,5). The rock foliations are defined both by alignment of minerals (i.e., white mica, biotite, and quartz ribbons) as well as sub-mm compositional layering consisting of quartzofeldspathic and micaceous domains. The quartzofeldspathic domains are dominated by fine-grain plagioclase and quartz grains (< 200 µm diameter). Quartz grains show undulous extinction and evidence of subgrain development. Thin layers of fine grain biotite and white mica (< 150 µm length) are found throughout the plagioclase and quartz matrix. The micaceous domains are dominated by white mica porphyroblasts (0.5-4.5 mm length) with finer biotite (0.1-1 mm length) and minor plagioclase and quartz grains ( Fig. 2A). Pockets of coarse-grain biotite (0.2-1.5 mm length of individual grains) are found locally within the paragneiss (Fig. 2B). Porphyroblasts of plagioclase (1-5 mm diameter), subhedral to euhedral garnet (1-2.5 mm diameter), and kyanite (100-250 µm length) are found within both domains. Plagioclase often contain lobate and ovoid quartz intergrowths (Fig. 2C). Two morphologies of kyanite can be observed, the first with quartz intergrowths similar to plagioclase, and the second without intergrowths or abundant inclusions. The grain boundaries of the second type terminate against garnet or matrix quartz (Fig. 2D).
Many of the garnet grains in the Sippmikk paragneiss exhibited three zones in BSE intensity: (1) a dark core, (2) a bright mantle, and (3) a dark rim (Fig. 3). The boundaries between the cores and mantles are indistinct and curvilinear, whereas the boundary between the mantles and rims are much more definitive, as best observed  (Table S1). The chemical maps of the garnet illustrate the chemical variations between the three zones that correspond to the BSE intensities (compare Fig. 3A with Fig. 3B, C). Inclusions in garnets in all paragneiss samples were characterized with BSE imaging. For the Sippmikk paragneiss garnet, inclusions of biotite, white mica, albite, quartz, zircon, and magnetite are found in all three zones (Figs. 3,4,5). Most notably, polymineralic inclusions containing these minerals with negative garnet crystal shapes are abundant in the BSE bright zones (Fig. 3). Biotite can be found throughout garnet zones and are often subrounded in shape (Figs. 2C,4E). Kyanite inclusions are also found in the BSE bright garnet. Rutile inclusions are only found within the BSE dark rims (Fig. 3). No Tibearing phases (i.e., ilmenite nor titanite) are found within the BSE dark cores or bright mantles. The garnet porphyroblasts are locally embayed by coarse-grain polycrystalline quartz.
The garnet porphyroblasts typically show subhedral crystal habits with a low-density of inclusions. However, two instances of anhedral poikiloblastic garnet with abundant, coarse inclusions of quartz and quartz surrounded by plagioclase (< 200 µm diameter) are also present within the quartzofeldspathic domain (Fig. 2E). These garnet porphyroblasts also have embayments of coarse-grain polycrystalline quartz. The subhedral garnets contain inclusions of quartz, rutile, biotite, white mica, zircon, and monazite (Fig. 6E, F). Similar to the Sippmikk paragneisses, some garnet grains also contain polymineralic inclusions with negative garnet crystal shapes (Fig. 3G) and subrounded biotite inclusions (Fig. 2F). High-contrast BSE images showed no zoning of any garnet porphyroblasts, in accordance with the results of Bukała et al. (2020b). The lack of identifiable garnet zoning A White mica porphyroblasts with undulous extinction surrounded by smaller, recrystallized white mica grains. The matrix consists of recrystallized quartz, plagioclase porphyroblasts, fine-grain biotite, and clinozoisite. B A pocket of coarse-grain biotite and recrystallized quartz adjacent to a kyanite porphyroblast within the fine-grain matrix. C Plagioclase porphyroblast with lobate and ovoid quartz intergrowths adjacent to a garnet porphyroblast containing a rounded biotite inclusion. The porphyroblasts are surrounded by fine-grain white mica, biotite, and recrystallized quartz constituting the matrix. D Kyanite porphyroblasts showing different morphological characteristics. The kyanite in the center contains lobate and ovoid quartz intergrowths whereas the other kyanite porphyroblasts do not show intergrowths and their grain abruptly boundaries terminate against garnet (top left) or again recrystallized quartz (lower right). E Poikiloblastic garnet within the quartzofeldspathic domain of the paragneiss. The garnet contains lobate and ovoid quartz inclusions as well as inclusions of quartz surrounded by plagioclase (top left and right). The lower half of the microphotograph shows the boundary of the garnet porphyroblast embayed by coarse-grain quartz adjacent to fine-grain recrystallized quartz with biotite that define the foliated rock matrix. F Irregularly shaped garnet porphyroblasts surrounded and embayed by coarsegrain biotite within the biotite-kyanite domain of the paragneiss. The coarse-grain biotite is associated with plagioclase and fine-grain quartz. G Kyanite porphyroblast within the biotite-kyanite domain of the paragneiss. The kyanite is partially pseudomorphed by sericite and intergrown with coarse-grain biotite, collectively surrounded by a fine-grain matrix of biotite, white mica, and quartz. H Garnet porphyroblast with a round biotite inclusion within the biotite-kyanite domain of the paragneiss. Locally, the garnet boundary has a thin film of quartz in contact with plagioclase, exhibiting myrmekite texture. The matrix is characterized by fine-grain biotite, plagioclase, quartz, white mica, and sericite. A plagioclase porphyroblast is also adjacent to the kyanite. Minerals are abbreviated according to Whitney and Evans (2010)

Zircon U-Pb geochronology
Zircon grains in the Sippmikk paragneiss of the Avardo Gneiss (sample SM19-03) are found amongst the matrix minerals and as inclusions in both garnet and kyanite porphyroblasts. Many large zircon grains exhibit a core-rim texture with BSE bright cores that are chemically heterogeneous, commonly with oscillatory zoning patterns, and BSE dark rims that are relatively chemically homogeneous. Zircon grains with these characteristics are only found within garnet mantles, rims, in kyanite, and predominantly in micaceous domains of the paragneiss matrix (Figs. 3, 4). A thinner, brighter outer rim is observed on some of the grains, but they were not analyzed as they are thinner than the 15 µm beam diameter used for analysis. A total of thirty-six trace element analyses and fortyseven geochronological analyses were accepted from forty-four zircon grains. The analyses were accepted following the protocol outlined in the methodology. The results are summarized in Table 1 and presented in full in Tables S2 and S3. Fourteen trace element analyses of the zircon cores yielded steep chondrite normalized HREE patterns, whereas twenty-two analyses of rims produced significantly shallower HREE patterns ( Fig. 7A; Table 1). LREE were typically below the detection limit for the zircon grains but nine analyses of the cores show more pronounced negative Eu anomalies than four analyses of the rims (Fig. 7A; Table 1). Thirty-two geochronological analyses of zircon cores provide a range of 207 Pb/ 235 U dates from 661.7 ± 13.6 Ma to 1736.4 ± 35.3 Ma, and 206 Pb/ 238 U dates from 639.9 ± 15.3 Ma to 1712.9 ± 33.9 Ma with discordance values ranging from -26.4 to 0.6% (with negative values representing normal discordance). In Wetherill concordia space, the analyses plot as broad discordia line (Fig. 7B). Model-I linear regression of analyses with > 4% normal discordance (n: 21) provides a lower-intercept of 487.5 ± 90.7 Ma (MSWD: 14). Fifteen U-Pb analyses of rims provide 207 Pb/ 235 U dates from 469.5 ± 13.1 Ma to 501.1 ± 10.8 Ma, and 206 Pb/ 238 U dates from 470.9 ± 11.5 Ma to 494.2 ± 12.1 Ma with discordance ranging between -2.5 to 2.0% (Fig. 7B). A Wetherill concordia age for these analyses is calculated to 482.5 ± 3.7 Ma (MSWD: 1.9). The zircon cores are generally higher in U content compared to the rims (Table 1).

Monazite Th-U-Pb geochronology
Monazites were analyzed from two samples of the Sippmikk paragneiss of the Avardo Gneiss (samples SM19-02 and SM19-03). The results are summarized in Table 1 and presented in full in Tables S4 and S5. Monazites in paragneiss SM19-02 are found in the BSE dark garnet rims and in the matrix, whereas monazites in paragneiss SM19-03 are only found within the matrix (Fig. 5). The matrix monazites in both rocks are associated with mica-rich regions. All monazite grains are surrounded by mineral agglomerations consisting of apatite, allanite, clinozoisite and other REEbearing phases. BSE images of monazites demonstrate no obvious chemical zoning (Fig. 5).
For paragneiss SM19-02, a total of twenty geochronological analyses and seventeen trace element analyses were accepted from ten monazite grains. For paragneiss SM19-03, twelve geochronological analyses and ten trace element analyses were accepted from four grains. The analyses were accepted according to the protocol outlined in the methodology. The REE, Y, and Sr content of monazites in both paragneisses are similar (Table 1) and no systematic trends are recognized for either sample. Monazite in paragneiss SM19-02 show higher Th/U values than monazite in paragneiss SM19-03, which overall decrease with younger 208 Pb/ 232 Th dates ( Fig. 8A; Table 1). Monazite in paragneiss SM19-02 yielded a spread of concordant dates with 208 Pb/ 232 Th dates ranging from c. 412.5 to 456.9 Ma, whereas the results from paragneiss SM19-03 cluster on concordia, producing an age of 420.6 ± 2.0 Ma (MSWD: 1.2; Fig. 8B).
Monazites were analyzed from one sample of the Kittelfjäll paragneiss of the Marsfjället Gneiss (sample MJ18-05A; Tables S4 and S5), which was obtained from the same rock sample that Bukała et al. (2020b) used to conduct zircon geochronology. The results are summarized in Table 1 and presented in full in Tables S4 and S5. In the paragneiss, monazite grains are predominantly found in the biotite-kyanite domain, but are also found to be associated with smaller biotite pockets in the quartzofeldspathic domains and also included in garnet and kyanite. Regardless of the textural position, the BSE images of the monazite reveal three domains (Fig. 6): (1) BSE-dark, homogeneous to weakly zoned cores, (2) BSE-medium, heterogeneously zones mantles, and (3) BSE-bright, homogeneous rims. The boundaries between all three domains are irregular, often curvilinear or sharp and angular. In some grains, the mantles appear to be transected by the rims, and the cores transected by both mantles and rims.
A total of sixty trace element analyses and eighty-one geochronological analyses were accepted from thirty-three monazite grains in the Kittelfjäll paragneiss. The analyses were accepted based on the protocol outlined in the methodology. The grains are located in the matrix, in garnet, or in kyanite, but no systematic patterns were recognized according to the textural positions. All three monazite domains were analyzed, and where possible, Th-U-Pb analyses were conducted in the same domains as the trace element analyses. The Th/U ratios of the monazite show a weak correlation with Th content and a strong correlation with U content (Fig. 9A). Fifteen trace element analyses of the monazite cores show distinctly lower Th/U with significantly higher U, Y, HREE, weakly elevated Sr content, and subdued negative Eu anomalies (Eu*) compared to thirty analyses of mantles and fifteen analyses of rims. The mantles and rims are chemically indistinguishable (Fig. 9A, B; Table 1). Seventeen geochronological analyses of the cores cluster in 208 Pb/ 232 Th vs. 206 Pb/ 238 U concordia space, providing an age of 481.6 ± 2.1 Ma (MSWD: 2.0; Fig. 9C). Forty geochronological analyses of the mantles and twenty-four analyses of rims produced indistinguishable results in concordia space and range from c. 493 to 455 Ma with an average concordia age of 476.6 ± 1.5 Ma (MSWD: 4.3; Fig. 9D). To delimit the endmember ages of the collective mantle and rim results, the ten oldest and ten youngest analyses were used to calculate ages, yielding 488.9 ± 1.9 Ma (MSWD: 0.8) and 463.1 ± 1.8 Ma (MSWD: 1.1), respectively.

Metamorphic Evolution of the Avardo Gneiss at Sippmikk Creek
The Avardo Gneiss has previously been described to have undergone partial melting that was responsible for development of compositional layering of the paragneisses (Van Roermund and Bakker 1983). The Sippmikk paragneiss (belonging to the Avardo Gneiss) exhibits layering akin to partial melting. The matrix phases of the paragneiss show evidence of recrystallization (Fig. 2), thereby preventing unequivocal identification of melt-related microstructures throughout the rock volumes. However, the polyphase inclusions with decrepitation cracks that are found within the Sippmikk paragneiss garnet mantles are interpreted as melt inclusions, indicating that garnet is peritectic ( Fig. 3; Cesare et al. 2009;Ferrero et al. 2012Ferrero et al. , 2015Bartoli et al. 2013). The subrounded biotite inclusions in garnet (Fig. 2C, H) and the coarse-grain quartz embayments are further evidence for peritectic garnet growth, attributed to biotite dehydration melting (Braun et al. 1996;Fitzsimons 1996;Waters 2001;Barbey 2007).
There is evidence for white mica-dehydration melting within the kyanite stability field prior to biotite-dehydration melting. The peritectic garnet mantles in the Sippmikk paragneiss envelope garnet cores with bell-shaped Y zoning. This zoning indicates the cores formed during prograde metamorphism with Y preferentially partitioned to garnet according to Rayleigh fractionation (Symmes and Ferry 1992;Lanzirotti 1995;Pyle and Spear 2003). The irregular and curvilinear boundaries between the cores and mantles, together with the high Mn content of the latter (Fig. 3A-F), suggests that the garnet cores were partially resorbed (Perchuk et al. 2005(Perchuk et al. , 2008Liu et al. 2014;Xia et al. 2016;Xia and Zhou 2017). Garnet has been documented to be a reactant during white mica-dehydration melting, producing kyanite and biotite as products (Indares and Dunning 2001;Lang and Gilotti 2007). Thus, the kyanite and biotite inclusions within the garnet mantles may be products of the melting prior to peritectic garnet growth. Furthermore, the kyanite porphyroblasts intergrown with lobate quartz can be of peritectic origin.
Zircon with core-rim structure occur within the garnet mantles, rims and matrix. The zircon structures are interpreted as xenocrystic cores that are overgrown by metamorphic rims (Fig. 4; e.g., Whitehouse and Platt 2003). The xenocrystic cores provide a dispersion of dates and typical igneous REE patterns that is representative of a detrital zircon provenance. The discordance of the core U-Pb dates, producing a lower intercept of 487.5 ± 90.7 Ma (Fig. 7B), demonstrates Pb-mobilization broadly coeval with Caledonian tectonism, requiring high temperature or dissolution-reprecipitation in order to mobilize Pb (Cherniak and Watson 2001;Harley et al. 2007;Rubatto 2017). The shallow HREE slopes of the zircon rims, compared to the zircon cores (Fig. 7A), is a typical pattern for metamorphic zircon that formed either after or coeval with garnet crystallization in a closed system environment (Rubatto 2002;Hermann and Rubatto Fig. 4 BSE images showing the general mineral assemblages of Sippmikk paragneiss sample SM19-03 and highlighting the textural positions of zircons that were analyzed. High-contrast BSE image insets exhibit the zircon core-rim zonation. A-E: Zircon grains located in the paragneiss matrix close to garnet porphyroblasts. D-G: Zircon grains included in garnet porphyroblasts. H: A zircon grain included in a kyanite porphyroblast. The low-contrast BSE images F and G correspond to the high-contrast BSE images in Fig. 2F and E, respectively. Minerals are abbreviated according to Whitney and Evans (2010) ◂ 2003; Whitehouse and Platt 2003;Hokada and Harley 2004). Considering that resorption of garnet cores would release HREE into the reactive bulk rock, the zircon rims are best explained to have formed coeval with peritectic garnet during biotite-dehydration melting. The subdued Eu anomalies of the rims suggest melting may have occurred in (U)HP conditions. This evolution would explain the inclusions of rounded biotite, kyanite, and zircon in peritectic garnet. This would also infer that the zircon rims formed after kyanite, yet they are found within kyanite porphyroblasts. However, the kyanite morphologies within the matrix (Fig. 2D) suggest two generations, which was described for the Avardo Gneiss by Van Roermund and Bakker (1983). The first generation of kyanite was  Whitney and Evans (2010) reported to occur with K-feldspar, in line with white micadehydration melting. The second generation was reported to overgrow sillimanite (not observed in the studied paragneisses) and occur in relation with staurolite, perhaps reflecting Middle Ordovician metamorphism (Brueckner and Van Roermund 2007). Monazites in paragneiss SM19-03 provide a more concise cluster of dates and Th/U values, providing definitive evidence of the late Silurian event at 420.6 ± 2.0 Ma (Fig. 8B), reproducing previous geochronology suggesting partial melting of the Avardo Gneiss at 423 ± 5 Ma (Williams and Claesson 1987). However, there is no obvious textural evidence in the Sippmikk paragneiss to suggest partial melting associated with monazite formation. The monazite in paragneiss SM19-02 are found both within garnet rims and in the matrix. The spread of  Middle Ordovician to late Silurian dates that correlates with decreasing Th/U values (Fig. 8A) indicates that the monazites underwent dissolution-reprecipitation (Seydoux-Guillaume et al. 2003;Weinberg et al. 2020), which partially reset a Middle Ordovician (or older) history during a late Silurian metamorphic event (Fig. 10). As a result, it is possible that garnet rims formed during Middle Ordovician metamorphism and entrapped monazite (associated with staurolite and possibly second generation of kyanite; Van Roermund and Bakker 1983). The monazite then underwent dissolution-reprecipitation in late Silurian time, facilitated by fluid ingress along garnet fractures.

Metamorphic evolution of the Marsfjället Gneiss near Kittelfjäll
The compositional layering of the Kittelfjäll paragneiss is attributed to partial melting, resolved to late Cambrian/early Ordovician time (Bukała et al. 2020b). Similar to the Sippmikk paragneiss, the matrix of the Kittelfjäll paragneiss is recrystallized and partially retrogressed. However, the polymineralic inclusions (Fig. 3G), the embayed and poikiloblastic garnet (Fig. 2E), and subrounded biotite inclusions (Fig. 2H), all indicate biotite-dehydration melting of the paragneiss and peritectic garnet growth (Braun et al. 1996;Fitzsimons 1996;Waters 2001;Barbey 2007). Furthermore, many garnet grains are observed to be partially replaced by biotite (Fig. 2F) and are locally enveloped by quartz films in contact with plagioclase with myrmekitic texture (Fig. 2H). These features indicate garnet porphyroblasts were partially resorbed (Waters 2001), as demonstrated by chemical maps of the garnet (Bukała et al. 2020b). Furthermore, the occurrence of kyanite intergrown with biotite in association with plagioclase (Fig. 2G), suggests that the biotite-dehydration reaction boundary was crossed during cooling to produce biotite, kyanite, and plagioclase at the expense of garnet and K-feldspar, which has been reported for (U)HP rocks (Lang and Gilotti 2007;Kendrick and Indares 2018). The complexly zoned monazite grains are found as inclusions within both garnet and kyanite porphyroblasts and are concentrated within the matrix of biotitekyanite domain. Their textural positions indicate that the  Pb/ 232 Th concordia space. The analyses are presented together as they cannot be distinguished according to U, Th, Th/U, REE, Y, Eu*, nor Sr content. The center white ellipse represents the concordia age for all analyses, whereas the upper and lower ellipses represent the concordia ages for the ten oldest and ten youngest analyses, respectively zoning patterns may be related to interaction with melt (Harlov et al. 2011;Weinberg et al. 2020;Ding et al. 2021).
The age of monazite core formation is established at 481.6 ± 2.1 Ma (Fig. 9C). The high Y and HREE content of the cores with respect to mantles and rims suggest that it formed in the absence of garnet, yet garnet is a major constituent of the paragneiss. Metamorphic zircon from the Kittelfjäll paragneiss show variable HREE patterns, with some suggesting dissolution of garnet during their formation, while other zircon grains indicate HREE depletion due to garnet growth (Bukała et al. 2020b). Although the garnet chemistry does reveal distinct growth zones (Bukała et al. 2020b), it is plausible that garnet was present during the formation of the monazite cores but the relatively high Y and HREE content reflects formation during dissolution of garnet ( Fig. 10; Kelly et al. 2006;Ding et al. 2021). Subsequent growth of peritectic garnet would have then provided a Y sink during mantle and rim formation (Foster et al. 2000;Pyle et al. 2001;Kohn et al. 2005;Rubatto et al. 2013). It may also be possible that monazite formed prior to any  Fig. 10 Schematic summaries depicting the important metamorphic processes and the related response of zircon and monazite for the Sippmikk paragneisses (top) and the Kittelfjäll paragneiss (bottom). The coloring of the minerals is according to the color representation in previous figures, highlighting the results of this study. The coloring of garnet zones for the Sippmikk paragneiss reflects the Mn zoning pattern presented in Fig. 3 with hotter colors equaling higher Mn content. The evolution of the Sippmikk paragneiss summarized as follows: A Possible white mica-dehydration melting at or before 482.5 ± 3.7 Ma. Peritectic kyanite crystallizes with lobate quartz intergrowths. Garnet is partially consumed as a reactant of the melting reaction. B Biotite-dehydration melting at 482.5 ± 3.7 Ma. Zircon rims crystallize on the xenocrystic detrital zircon, which may have been partially resorbed during initial melting. Peritectic garnet is crystallized, trapping melt inclusions (MI), zircon, and remnants of kyanite that were formed from white mica-dehydration. HREE is sequestered by peritectic garnet. C Crystallization of garnet rims and a second generation of kyanite porphyroblasts, possibly occurring at 454.8 ± 9.7 Ma (i.e., the timing of (U)HP metamorphism recorded by eclogite at Sippmikk Creek; Brueckner and Van Roermund, 2007). Monazite is captured by garnet rims and is present in the matrix. Both positions of the monazite undergo dissolution-reprecipitation at 420.6 ± 2.0 Ma. The evolution of the Kittelfjäll paragneiss is summarized as follows: D Possible onset of partial melting at 481.6 ± 2.1 Ma. Partial resorption of pre-existing prograde garnet porphyroblasts provides conditions to crystallize high-Y monazite, and metamorphic zircon with HREE patterns indicating absence of garnet (Bukała et al. 2020b). E Partial melting of the paragneiss, possibly with multiple events from 481.6 ± 2.1 Ma to 463.1 ± 1.8 Ma. Monazites undergo dissolution-reprecipitation, mediated by the melt. Biotite-dehydration melting results in the growth of peritectic garnet, which captures melt inclusions (MI) and sequesters Y and HREE from monazite and zircon. Zircon crystallization (Bukała et al. 2020b) is responsible for sequestration of U from monazite. Pb is partially lost from both mantle (yellow) and rims (blue). F Final crystallization of melt after 463.1 ± 1.8 Ma. Kyanite porphyroblasts are produced and entrap monazite and zircon that recorded partial melting. Garnet is partially resorbed garnet growth, however, this would not explain the HREE patterns of zircon from the paragneiss (Bukała et al. 2020b).
The chemistry of the mantles and rims are distinct from the cores but are themselves indistinguishable (Fig. 9A, B), indicating that the mantles and rims formed in similar conditions (Fig. 10). The lower U content in the mantles and rims, relative to the cores, may be a function of zircon growth in the paragneiss (Bukała et al. 2020b), whereas Y and HREE depletion indicates crystallization of the subhedral peritectic garnet ( Fig. 10; Foster et al. 2000;Pyle et al. 2001;Kohn et al. 2005;Rubatto et al. 2013). Furthermore, the mantles and rims show lower Sr content and stronger Eu anomalies than the cores (Fig. 9B), which suggests crystallization of plagioclase that was not present during core formation (Holder et al. 2015), consistent with cooling during decompression and melt crystallization (Lang and Gilotti 2007;Kendrick and Indares, 2018). However, other factors may control the Eu anomalies (Holder et al. 2020).
The textures of the mantles and rims resemble formation via dissolution-reprecipitation of pre-existing monazite in the presence of melt (Harlov et al. 2011;Weinberg et al. 2020;Ding et al. 2021). The age of 463.1 ± 1.8 Ma reflects the most complete Pb loss from the monazite and likely record the youngest episode of partial melting. Based on the association and textures of kyanite, biotite, and plagioclase (Fig. 2G), it is likely that the Kittelfjäll paragneiss cooled through the biotite-dehydration reaction boundary, consuming garnet as a result (Fig. 2F;Waters 2001;Lang and Gilotti 2007;Kendrick and Indares 2018). However, the monazite grains do not record the melt crystallization as garnet resorption should release Y and HREE to be incorporated in the monazite. Therefore, the monazite mantles and rims provide a maximum age for cooling and final melt crystallization of the paragneiss. The monazite inclusions within peritectic garnet should record an older partial melting event than those encased in kyanite. However, the oldest ages of the monazite mantles and rims are older than the monazite cores, indicating Pb re-distribution in the reprecipitated monazite volume, which obfuscates detection of older events in the inclusions. Although multiple partial melting events is suggested by re-melting of possible peritectic garnet, possibly starting at 481.6 ± 2.1 Ma, the timing of the events are  Barnes et al. (2021a); 19 Dallmeyer and Gee (1986); 20 Barnes et al. (2021b); 21 Barnes et al. (2020a); 22 Dallmeyer and Stephens (1991) Barnes et al. (2020b); 28 Dallmeyer et al. (1990) equivocal due to Pb behavior within the monazite during dissolution-reprecipitation (Seydoux-Guillaume et al. 2003;Weinberg et al. 2020;Varga et al. 2020).

Implications for the evolution of the Seve Nappe complex
The Sippmikk paragneiss (belonging to the Avardo Gneiss) and the Kittelfjäll paragneiss (belonging to the Marsfjället Gneiss) provide the first evidence for metamorphism of the central SNC in late Cambrian/early Ordovician time (Fig. 11). The P-T conditions for the partial melting events are not quantified for the Avardo and Marsfjället gneisses. However, metamorphic microdiamonds were found within the outermost regions of prograde garnet cores enveloped by peritectic garnet rims in a paragneiss near Saxnäs (belonging to the Marsfjället Gneiss), indicating UHP metamorphism prior to partial melting (Petrík et al. 2019). The timing of UHP metamorphism was interpreted to have occurred at 472 ± 3 Ma based on monazite Th-U-total Pb geochronology, yet, the monazites are characterized by a dispersion of dates defining isochron endmember ages of 479 ± 9 Ma and 460 ± 16 Ma, reflecting monazite dates of the Kittelfjäll paragneiss. Provided that partial melting in the Sippmikk paragneiss occurred at 482.5 ± 3.7 Ma, and possibly started at 481.6 ± 2.1 Ma in the Kittelfjäll paragneiss, it seems likely that the Saxnäs paragneiss underwent a similar history, inferring late Cambrian UHP metamorphism. The Sippmikk paragneiss may also have been in (U)HP conditions during or prior to partial melting, evinced by kyanite within peritectic garnet, together with two distinct generations of kyanite that are present elsewhere in the Avardo Gneiss (Van Roermund and Bakker 1983). These records reflect the evolution of the microdiamond-bearing Tväråklumparna Gneiss in the westcentral SNC at 482.6 ± 3.8 Ma, either during or immediately proceeding (U)HP metamorphism (Walczak et al. 2021).
Prior to the investigation of the Tväråklumparna Gneiss, late Cambrian/early Ordovician metamorphism was only identified in the northern SNC terranes as well as the Vestgötabreen Complex on Svalbard ( Fig. 11; Root and Corfu 2012; Andréasson et al. 2018;Barnes et al. 2021a, b;Fassmer et al. 2021;Baird et al. 2022). Therefore, the records of the Avardo, Marsfjället, and Tväråklumparna gneisses indicate that a larger extent of the SNC was subducted in late Cambrian/ early Ordovician time and that the west-central and central SNC may have been deeper than the northern SNC. Eclogites, peridotites, and pyroxenites have established (U)HP metamorphism of the Avardo Gneiss and the Sjouten Unit in Middle Ordovician time (Fig. 11) with general conditions of 2.3-3.0 GPa and 650-800 °C Majka et al. 2014b;Gilio et al. 2015;Klonowska et al. 2016). Spinel peridotite in the Kittelfjäll paragneiss provided P-T conditions of 1-2 GPa and 650-830 °C, representing extraction from the mantle wedge during exhumation of the gneiss (Clos et al. 2014). The youngest monazite in the Kittelfjäll paragneiss indicate metamorphism proceeded until at least 463.1 ± 1.8 Ma. Considering the regional record of Middle Ordovician metamorphism, it is likely that the spinel peridotite approximates the Middle Ordovician P-T conditions of the Kittelfjäll paragneiss. The record of Middle Ordovician (U)HP metamorphism is also documented in the west-central SNC terranes (Majka et al. 2012), however, it is apparently absent from the northern SNC (Fig. 11).
No zircon and monazite investigated in this study (Figs. 7,8,9; Table 1) nor previous geochronological studies of the central SNC Brueckner and Van Roermund 2007;Fassmer et al. 2017;Petrík et al. 2019) indicate Late Ordovician/early Silurian metamorphism, with the exception of zircon and rutile in eclogite belonging to the Sjouten Unit (Root and Corfu 2012; Fig. 11), which may record exhumation of the unit (Fassmer et al. 2017). The sparse Late Ordovician/early Silurian geochronological record of the central SNC contrasts both the west-central and northern SNC terranes. In the former, the rocks bear a strong record of partial melting associated with exhumation of the rocks from (U)HP conditions, recorded by zircon and monazite (Majka et al. 2012;Ladenberger et al. 2014;Klonowska et al. 2017;Walczak et al. 2021). In the latter, the rocks are characterized by deformation during exhumation in crustal conditions, recorded by white mica (Dallmeyer and Gee 1986;Dallmeyer and Stephens 1991;Page 1992;Barnes et al. 2020aBarnes et al. , 2021b, which is also recorded in the Vestgötabreen Complex ( Fig. 11; Dallmeyer et al. 1990; Barnes et al. 2020b). All of the SNC terranes were subsequently involved in continental collision between Baltica and Laurentia in Silurian to Devonian time (Fig. 11).
In summary, the overall P-T conditions of the SNC are the highest in late Cambrian/early Ordovician time. The P-T conditions decrease progressively southwards through time with Middle Ordovician (U)HP metamorphism restricted to central/west-central SNC terranes. Subsequent Late Ordovician/early Silurian partial melting is only evident in the west-central SNC, coeval with possible subsolidus exhumation of the central SNC, and deformation in crustal levels for the northern SNC. Altogether, this pattern reflects south-to-north oblique subduction of the SNC starting in late Cambrian/early Ordovician time, followed by progressive north-to-south exhumation of the (U)HP terranes to crustal levels from early Ordovician to Late Ordovician/early Silurian time. Oblique collision from south-to-north is supported by counter-clockwise rotation of Baltica (Bottrill et al. 2014) that is reflected in its paleomagnetic record in late Cambrian time (Torsvik and Rehnström 2001;Cocks and Torsvik 2002). This evolution contrasts previous models for subduction of the SNC (e.g., Bukała et al. 2018;Fassmer et al. 2021) providing a new basis for evaluating the closure of the Iapetus Ocean leading to continental collision of Baltica and Laurentia.

Conclusions
Zircon geochronology of the Sippmikk paragneiss (belonging to the Avardo Gneiss) resolves the timing of partial melting to 482.5 ± 3.7 Ma, likely related to biotite-dehydration. Monazite geochronology of the Sippmikk paragneiss reveals another metamorphic event at 420.6 ± 2.0 Ma, possibly overprinting a Middle Ordovician history of monazite. In the Kittelfjäll paragneiss (belonging to the Marsfjället Gneiss), relatively high-Y and low Th/U monazite cores provide a crystallization age of 481.6 ± 2.1 Ma and may be a response to resorption of garnet during initial partial melting. The monazite subsequently underwent dissolution-reprecipitation in the presence of melt, which re-equilibrated the trace element content as Y and U were sequestered by garnet and zircon, respectively. The Pb loss from monazite during dissolution-reprecipitation was less efficient, producing a spread of Th-U-Pb isotopic dates from late Cambrian to Middle Ordovician time. The youngest partial melting event that caused dissolution-reprecipitation of monazite occurred at 463.1 ± 1.8 Ma. Subsequently, the rock cooled and final melt crystallization produced biotite, kyanite, and plagioclase. The textural evidence of the paragneiss suggests that several partial melting events occurred.
The Sippmikk and Kittelfjäll paragneisses provide first evidence for late Cambrian/early Ordovician metamorphism of the central SNC. The tectonic evolutions of the paragneisses, together with geochronological and petrological studies, indicates the SNC underwent south-to-north oblique subduction in late Cambrian time, followed by progressive north-to-south exhumation to crustal levels prior to late Silurian continental collision. These results demonstrate the effectiveness for refractory minerals in metasedimentary rocks to record older histories compared to mafic and ultramafic lithologies to establish a comprehensive tectonic evolution.