<<

J OURNAL OF Journal of Petrology, 2018, 1–29 doi: 10.1093/petrology/egy094 P ETROLOGY Advance Access Publication Date: 15 October 2018 Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 Original Article

Mechanical Mixing of and in the Orogenic Peridotite Lenses of the Tvaerdal Complex, Liverpool Land, Greenland Caledonides Hannes K. Brueckner1,2*, L. Gordon Medaris, Jr3, William L. Griffin4, Scott M. Johnston5, Ebbe H. Hartz6, Norman Pearson4, Yue Cai1 and Arild Andresen7

1Lamont-Doherty Earth Observatory of Columbia University, 61 US-9W Palisades, NY 10964, USA; 2The CUNY Graduate Center and School of Earth and Environmental Sciences, Queens College, 65-21 Main St. Flushing, NY 11367, USA; 3Department of Geoscience, University of Wisconsin-Madison, Madison, WI 53706, USA; 4Department of Earth and Planetary Sciences, Australian Research Council Centre of Excellence for Core to Crust Fluid Systems (CCF)/(GEMOC), Macquarie University, NSW 2109, Australia; 5Physics Department, California Polytechnic State University, San Luis Obispo, CA 93407, USA; 6Centre for Earth Evolution and Dynamics, Department of Geosciences, Oslo University, P.O. Box 1028 Blindern, 0316 Oslo, Norway; 7Department of Geosciences, University of Oslo, P.O. Box 1047 Blindern, 0316 Oslo, Norway

*Corresponding author. P.O. Box 1000, 61 Route 9W, Palisades, NY 10964-1000 USA, Telephone: (845) 359-2000. Fax: (845) 365-8150. E-mail: [email protected]

Received April 9, 2018; Accepted September 28, 2018

ABSTRACT The Tvaerdal Complex is an -bearing metamorphic terrane in Liverpool Land at the southern tip of the Greenland Caledonides. It is a Baltic terrane that was transferred to Laurentia during the Scandian . It exposes a few small garnet and lenses, some containing par- allel layers of garnet pyroxenite and peridotite (including ). Sm–Nd ages from the indicate recrystallization occurred at the same time (405 Ma) as eclogite recrystallization in the enclosing gneiss. Geothermobarometry indicates these and pyroxenites shared a simi- lar pressure-temperature history. This congruent evolution suggests pyroxenite-bearing peridotite lenses were introduced from a mantle wedge into subducted Baltic and subsequent- ly shared a common history with this crust and its eclogites during the Scandian orogeny. Some gar-

net peridotite samples contain two garnet populations: one Cr-rich (35–62wt% Cr2O3)andtheother Cr-poor (02–14 wt %). Sm–Nd analyses of two such garnet define two sets of apparent ages: one older (>800 Ma) for Cr-rich and the other younger (<650 Ma) for Cr-poor garnets. We propose that the younger Cr-poor garnets were derived from fractured and disaggregated garnet

pyroxenite layers (i.e. are M2) and were mixed mechanically with older (i.e. M1) garnets of the host peridotite during intense Scandian shearing. Mechanical mixing may be an important mantle process.

Key words: orogenic peridotite; mechanical mixing; pyroxenite; peridotite, Greenland Caledonides

INTRODUCTION rocks in the south (Fig. 1). The metamorphic rocks in- Liverpool Land (LL) at the southern end of the clude numerous eclogites (Augland et al., 2010; Greenland Caledonides is a crystalline terrane com- Johnston et al., 2010; Corfu & Hartz, 2011) and a few, prised of plutonic rocks in the north and metamorphic small (decametre) peridotite lenses that locally are

VC The Author(s) 2018. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: [email protected] 1 2 Journal of Petrology, 2018, Vol. 0, No. 0

21°30’W Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018

4

2 ? 5 1 3 TJ

6 ? 70°30’N

0 1:1000001:100000 5 kmkm Scoresbysund

Glaciers and Glacial Deposits Permian Sediments Hurry Inlet Intrusive Suite Hurry Inlet Eclogitized Undifferentiated Gneiss Body LL Metamorphic Complexes Caledonides Eclogite Hurry Inlet Detachment Ultramafic Body } Jaettedal Sequence Marble Ittoqqortoormiit Shear Zone LL Tvaerdal Sequence 22°W

Fig. 1. Map of the Liverpool Land (LL) metamorphic terrane showing the Tvaerdal and Jaettedal metamorphic complexes, major tectonic features (from Johnston et al., 2010) and locations of peridotite lenses (circles) and eclogites (). Sample numbers refer to peridotite localities in Table 1. Geographic coordinates of most peridotite lenses are given in Table 1. T and J are the valleys called Tvaerdal and Jaettedal. garnetiferous (Augland et al., 2011). Augland et al. Orogeny is generally modeled by the westward (pre- (2010) proposed that the entire metamorphic terrane be sent coordinates) of Baltica (Scandinavia) called the Liverpool Land Eclogite Terrane (LLET). They beneath Laurentia (Greenland) (Krogh, 1977; Andersen proposed further that the LLET was a Baltic, rather than et al., 1991; see also Brueckner & Van Roermund, 2004; Laurentian, terrane (see also Smith & Cheeney, 1981) Hacker et al., 2010; Gee et al., 2012 for recent reviews). that became stranded on the Greenland side of the If the Tvaerdal Complex were a Laurentian terrane, the Caledonides during the opening of the Atlantic. peridotite bodies within it would have had to move up- Subsequent work has confirmed the Scandian (411–398 ward into the overlying plate. But if the Tvaerdal Ma) evolution of the eclogites and the Baltic origin of Complex was part of Baltica when it subducted beneath the gneisses that host them (Johnston, et al., 2010; Laurentia, it could have picked up the peridotite frag- Augland et al., 2011; Corfu & Hartz, 2011; Brueckner ments from the overlying Laurentian mantle wedge et al., 2016). However, mapping by Johnston et al. (Brueckner & Medaris, 2000). A contrary view to the (2010) showed that the metamorphic rocks are not a sin- mantle subduction model is, however, held by one of gle high pressure/ultrahigh pressure (HP/UHP) complex, us (see Hartz et al., 2005, 2007). Models in which the but rather two complexes, with the eclogites, perido- peridotites were derived from the Baltic rather than tites and their host rocks comprising a Baltican HP/UHP Laurentian mantle have been presented by Terry et al. terrane, the Tvaerdal Complex and the rest of the meta- (1999) and Majka et al. (2014). morphic rocks forming a Laurentian granulite facies ter- To test the possibility that the Tvaerdal Complex is a rane, the Jaettedal Complex. Baltic terrane, and the peridotites within it have a This division into a Baltic Complex (Tvaerdal) and a Laurentian affinity, seven of the known peridotite Laurentian Complex (Jaettedal) resolves the problem of bodies in Liverpool Land, as well as peridotite and pyr- how the peridotite lenses came to be embedded in the oxenite blocks in glacial till, were investigated. The southernmost Greenland Caledonides. The Scandian major, minor and trace element compositions of Journal of Petrology, 2018, Vol. 0, No. 0 3 peridotite and pyroxenite whole- samples, as well As noted above, Johnston et al. (2010) concluded as the compositions of their constituent were that the eclogites and peridotites within the meta- Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 determined. Electron microprobe scans were recorded morphic complex are restricted to a tectonically lower across key minerals. Pressure-temperature (P–T) equili- orthogneiss sequence (the Tvaerdal Complex) sepa- bration conditions and Sm–Nd mineral ages were deter- rated by another shear zone (the Ittoqqortoormiit shear mined from selected samples and compared with zone) from a tectonically higher eclogite-free unit con- published pressure-temperature-time (P–T–t) informa- taining predominant orthogneiss and subordinate tion from the associated Tvaerdal eclogites and host granulite, pelitic , calc-silicate rocks and marble gneisses. Re–Os isotopic data from sulfides are also (the Jaettedal Complex). Zircon cores from the presented for selected samples. Jaettedal Complex give and Mesoproterozoic The study became complicated when we observed U–Pb ages, while zircon cores from the Tvaerdal that some garnet peridotites contain garnets with vari- Complex give Mesoproterozoic ages only. The able Cr concentrations, which, in at least one case could Jaettedal Complex underwent granulite facies meta- be divided into two distinct populations: Cr-rich and Cr- morphism and melting between 460 and 410 Ma poor. Sm–Nd results from two other peridotites indicate (Johnston et al., 2010, 2015), whereas the peak meta- that their garnets formed or re-equilibrated during two morphic history of the Tvaerdal Complex occurred be- separate metamorphic events. Thus some peridotites tween 412–395 Ma (see below). Thus the two contain two garnet populations and we present evi- complexes had different Precambrian histories and dence that they might have been mixed together mech- were metamorphosed at different grades and at differ- anically. Mechanical mixing may be an under- ent times, providing compelling evidence that the two appreciated process that could play an important role in complexes originated from Laurentia and Baltica, re- mantle evolution. spectively. This subdivision requires further evaluation, but in the meantime we suggest that the two meta- morphic complexes should be called collectively the LL REGIONAL SETTING Metamorphic Terrane, with eclogite-facies HP/UHP The Greenland Caledonides restricted to the Tvaerdal Complex. The southern half of the Greenland Caledonides is com- Together the Hurry Inlet Plutonic Terrane and posed of west vergent nappe complexes in the west Jaettedal Complex represent deep levels of a and crystalline complexes, including LL, to the east Laurentian continental arc complex that was active dur- (Gee et al., 2008, 2012; Higgins et al., 2008; Gasser, ing the closure of Iapetus, while the Tvaerdal Complex 2014). The allochthons are separated from LL by 100 km was a fragment of the approaching Baltic passive mar- of post-orogenic rift deposits, but LL is generally consid- gin. Igneous activity within this arc complex ended, and ered to be tectonically positioned below these alloch- metamorphism of the Tvaerdal Complex began, when thons and to represent the deepest tectonic level of the Baltica and Laurentia collided and the Tvaerdal Greenland Caledonides, similar in tectonic position to Complex was subducted into the mantle beneath the Western Gneiss Complex (WGC) of the Norwegian Laurentia. There, the Tvaerdal Complex underwent HP/ Caledonides. UHP metamorphism and collected fragments of perido- tite from the overlying mantle wedge. The Tvaerdal and The Liverpool Land Metamorphic Terrane Jaettedal Complexes were subsequently juxtaposed Liverpool Land is the southernmost exposed crystalline along the Ittoqqortoormiit shear zone as the Tvaerdal terrane of the Greenland Caledonides (Fig. 1), covered Complex was exhumed from the . This to the south by the fjord Scorsbysund and by Mesozoic juxtaposition occurred in the lower to middle crust sedimentary sequences. Cheeney (1985) divided the (Johnston et al., 2015) at approximately 395 Ma. Finally, crystalline rocks into a northern, largely igneous com- as a composite metamorphic terrane (i.e. the LL plex, the Hurry Inlet Plutonic Terrane and a southern Metamorphic Terrane), they shared a common retro- metamorphic complex (Fig. 1). The shallowly north- grade and melting history while they were exhumed to- dipping Gubbedalen Shear Zone (aka Hurry Inlet wards the surface along the Grubbedalen shear zone Detachment) with top to the north kinematic indicators between 395–360 Ma (Augland et al., 2010). separates the two complexes (Augland et al., 2010; Johnston, et al., 2010). The Hurry Inlet plutonic Complex is composed of calc-alkaline intrusives with Eclogites of the Tvaerdal complex screens and enclaves of Proterozoic supracrustal rocks Eclogites in Liverpool Land were first reported by (Coe, 1975; Johnston et al., 2010; Corfu & Hartz, 2011; Sahlstein (1935). They occur as metre-scale boudins Augland et al., 2012). The plutons have been dated be- and boudin trains and as portions of decametre to tween 475 to 415 Ma (Corfu & Hartz, 2011; Augland kilometre-scale mafic lenses and layers. Most eclogites et al., 2012; Brueckner et al., 2016) and probably repre- are massive and medium grained and contain ompha- sent deep intrusive levels of a continental arc complex cite þ garnet 6 quartz while some are layered, coarse that developed on the eastern margin of Laurentia dur- grained and commonly contain orthopyroxene 6 biotite ing the closure of Iapetus. in addition to garnet and . Further details can 4 Journal of Petrology, 2018, Vol. 0, No. 0 be found in Augland et al. (2010, 2011), Johnston et al. this fabric is believed to be the result of intense shear (2010) and Corfu & Hartz (2011). Brueckner et al. (2016) and/or shortening (Fig. 2 c, d and f). -bearing Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 calculated peak metamorphic conditions of 800–900C kelyphite reaction rims around garnets are locally and 35–40 kbar. Peak metamorphism was followed by stretched out into this planar fabric (Fig. 2f) indicating nearly isothermal decompression to 16 kbar and that some deformation occurred within the crust during 800C. Dating by a variety of techniques indicate that retrogression under facies conditions. UHP metamorphism occurred between 412–395 Ma fol- Folds were not found in any of the peridotite bodies. lowed by retrograde metamorphism and anatectic melt- Penetrative, parallel, thin veins crosscut ing between 395–375 Ma (Augland et al. 2010; Corfu & the earlier compositional banding in most lenses. The Hartz, 2011; Brueckner et al., 2016). Taken together, the serpentinite veins consist of serpentine (largely chryso- very consistent ages confirm that HP/UHP metamorph- tile) and trains of opaque minerals, presumably Fe- ism of the Tvaerdal Complex was essentially simultan- oxides, which initially segment grains into eous with peak HP/UHP metamorphism in the WGC of aligned micro-ellipses and rectangular blocks before the Norwegian Caledonides during the Scandian they are completely replaced by serpentine. Some orogeny. lenses contain additional cross-cutting planar fabrics defined by serpentine þ oxides. TVAERDAL PERIDOTITES AND PYROXENITES Field relationships Petrography Peridotite and pyroxenite were first reported in Initial mineral proportions were calculated for most sam- Liverpool Land as blocks in glacial deposits by Smith & ples (Table 1) based on their bulk composition and the Cheeney (1981) and Cheeney (1985). Smith & Cheeney assumption that the assemblages consisted of variable (1981) noted the unusually high Cr contents of garnet proportions of olivine, orthopyroxene, clinopyroxene and clinopyroxene in one of them. There were no and garnet. These abundances were used to divide the reports of peridotite in outcrop until field parties led by samples into different peridotite and pyroxenite classes. Ebbe Hartz in 2005 and 2006 and Hannes Brueckner and The observed abundances of retrograde minerals, such Scott Johnston in 2009 located seven small (m–dm) as amphibole and opaque minerals are noted in paren- peridotite bodies (open circles in Fig. 1). Most occur- theses in Table 1. The calculated modes reveal a spec- rences are moderately to pervasively serpentinized and trum of rock types including dunite, harzburgite, only three (Fig. 1 and Table 1: locations 1, 2 and 5) pro- lherzolite, , olivine websterite, websterite and cli- vided samples with sufficient unaltered clinopyroxene nopyroxenite (Table 1). We suspect that some olivine- and garnet to warrant Sm–Nd analyses. Augland et al. rich assemblages, such as the , resulted from (2011) described another peridotite lens (map location mechanical mixing rather than through some chemical 6) with unaltered clinopyroxene and garnet and sam- or magmatic process. We present evidence for this sus- ples from that lens (LEA 08–04 and AA 08–47) were picion throughout the manuscript. donated for this study. We also collected additional and display a 05–20mmgrano- relatively unaltered samples from talus (HKB-6P, LM- blastic assemblage of olivine6 orthopyroxene6 clino- 42A LM-42B and LM-43) and from till (HKB Camp 1–1, 6amphibole6 and scattered, larger HKB Camp 2–1; Table 1). Some important data are from (2–30 mm) grains of garnet (Fig. 2 b, g and h). Lattice pre- these samples. ferred orientation (LPO) fabrics for LM-42B (Luc Mehl, All peridotite bodies exposed in outcrop are relative- personal communication) give [001](100) (C-type) pat- ly small concordant lenses a metre to a few metres terns for olivine and a [001](100) pattern for orthopyrox- wide and metres to tens of metres long (Fig. 2a). They ene. These patterns indicate olivine and orthopyroxene are (or were prior to serpentinization) dominantly fine- underwent dynamic recrystallization while undergoing to medium-grained dunite or harzburgite (olivine 6 large strains (Wang et al.,2013), while mechanically orthopyroxene) commonly containing dispersed, easily stronger garnet resisted recrystallization. visible (2–30 mm), partially to completely kelyphitized Garnet lherzolite, olivine websterite and wehrlite garnets (Fig. 2b). A few bodies are strikingly layered, form 05 to 10 cm thick layers within dunite and harz- defined by parallel garnet-rich and garnet-poor dunite burgite (Fig. 2d, e and h) and are texturally similar, with (Fig. 2c), pyroxenite and lherzolite layers (Fig. 2d and e). a relatively fine grained olivine matrix (05–10 mm) Some garnet-rich layers are thin, locally one or two enclosing larger grains of orthopyroxene, clinopyrox- grains wide, with the garnets separated from each other ene (10–30 mm) and garnet (2.0–30 mm). The size by recrystallized olivine (Fig. 2d, left side; see also range for the garnet is particularly striking (Fig. 2 g and Augland et al., 2011, Fig. 3). Three lenses contain 05–10 h). Especially notable are garnets in some samples (e.g. cm thick garnet and clinopyroxene-rich layers of garnet LM-42B; Fig. 2g) that show considerable colour varia- pyroxenite, garnet wehrlite and garnet lherzolite (Fig. 2c tions in plane light, ranging from a deep red or violet to and d). These layers are parallel to the overall foliation, orange to light pink. The differently coloured garnets which, in turn is concordant or nearly concordant to the are randomly distributed in some samples (LM-42B) enclosing gneissic foliation. The striking parallelism of while in others the lighter garnets are concentrated in Journal of Petrology, 2018, Vol. 0, No. 0 5

Table 1: Peridotite and pyroxenite localities and calculated modes, Tvaerdal Complex, Liverpool Land Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 Map Sample Grt bearing Coordinates ol opx cpx grt amp op other comment # # rock type

1. LM-38 Clinopyroxenite N70.56, W22.18 <1 0 40 60 (40) (5) tr green sp, bi fragmented grt 2. HKB 2A Dunite Tvaerdal 92 5 0 3 (10) 5% bi grt completely kely HKB 2B Websterite N70.56, W22.18 0 37 22 41 (2) (5) dark & pale grts 7 cm thick HKB 2C Ol-Websterite Tvaerdal 25 16 29 30 cpxite-dunite contact 3. HKB 3 Lherzolite Nu´ kaitsoq 46 14 19 20 (7) 37% serp kely largely amph N70.54, W22.25 4. HKB 4 Dunite Gubbedal 87 8 1 4 (8) (2) 3% kely N70.60, W22.22 5. EH-21A Lherzolite Tvaerdal 50 11 13 26 (15) 5% kely 30% serp EH-21B Harzburgite N70.55, W22.21 69 27 0 4 1.5 cm grt pclast EH-22 Wehrlite Tvaerdal 48 2 31 19 (10) dark & pale grts 30% serp 6. LEA 08–04 Peridotite N70.503, W22.34 nd nd nd nd nd see Augland et al. 2011 AA 08–47 Pyroxenite N70.503, W22.34 nd nd nd nd nd see Augland et al. 2011 — Camp 1–1 Harzburgite block in till 68 31 0 1 (8) sp 90% serp — Camp 2-1A Lherzolite block in till 50 15 29 7 tr tr sp 4 cm thick in dunite — Camp 2-1B Lhrz-Dun contact – – – – partially disaggretated — Camp 2–2 Dunite block in till 83 9 0 8 — LM-42A Harzburgite talus below 1 78 15 2 6 (tr) (2) 2% kely 40% serp LM-42B* Lherzolite (N70.56, W22.17) 57 10 12 10 7 2 2% sp, 8% kely Cr rich & Cr poor grt — LM-43 Websterite Talus below 1. 0 27 19 55 (10) 5cm thick vein in dunite — HKB-6P Harzburgite talus below 5 84 12 0 3 grt in single grain layers

The first column lists the sample locations shown in Fig. 1. Most modes were calculated from the anhydrous rock compositions recalculated to 100% (Table 2). Numbers in parenthesis are hydrous minerals and oxides present in thin section. *Whole-rock not analysed; abundance estimated from thin section. ol, olivine; opx, orthopyroxene; cpx, clinopyroxene; grt, garnet; amp, amphibole; op, opaques; sp, spinel; kely, kelyphyte; serp, serpentinite.

thin (mm scale) zones (Fig. 2h). These colour differen- Pyroxenites include garnet websterite and clinopyr- ces are also obvious in mineral separates viewed be- oxenite and are composed of medium grained (20–50 neath a binocular microscope. The colour variations are mm) granoblastic garnet þ clinopyroxene 6 orthopyr- primarily the result of different Cr contents, which range oxene. Most pyroxenites contain little or no olivine (0– from 02to63 wt % in sample LM-42B. Some garnet 10%; Fig. 2j and l) except for olivine websterite sample grains are blocky, angular or irregular in outline rather HKB-2C which has 25%. Amphibole occurs in most than equant granoblastic (Fig. 2g and h), as if the gar- samples, ranging from 2 to 40% and two samples con- nets had been mechanically fragmented prior to kely- tain biotite, one with 5%. Textures of the hydrous phitization. Many grains contain internal planar cracks phases indicate that they are retrograde. that are locally offset by other fractures. Garnets in other thin sections of both pyroxenite (Fig. 2i, j and k) and peridotite (Fig. 2m) are fractured and locally Whole-rock chemistry detached and separated, again suggesting mechanical Major elements breakdown. Whole-rock major element compositions of Tvaerdal All garnet grains in the peridotites are partly to com- peridotite and pyroxenite samples were determined by pletely replaced by kelyphite consisting of radially XRF. Analytical techniques are described in arranged fibres of orthopyroxene, spinel and amphibole Supplementary Data Table S1 (Supplementary Data are (Fig. 2n and o). Some kelyphites display a thin, discon- available for downloading at http://www.petrology. tinuous outer rim of orthopyroxene against the kely- oxfordjournals.org). Averaged compositions are listed phite, which in turn is enclosed by a zone of coarse in Tables 2 (major elements) and 3 (trace elements; grained amphibole (similar to textures described by complete data in Supplementary Data Table S2). Obata, 2011). The amphibole locally has a vermicular Serpentinization and the resultant production of structure where it borders the kelyphite and in some oxides has perturbed the chemistry of most samples; cases is optically continuous with amphibole threads thus the analyses are recalculated to 100% on an anhyd- within the kelyphite. Amphibole also occurs as isolated rous basis, with all Fe reported as FeO. Based on the grains within the rock matrix. One garnet in LM-42B proportions of olivine, orthopyroxene, and clinopyrox- (Fig. 2o) displays two reactions rims around a Cr-rich ene, the analysed Tvaerdal ultramafic suite includes grain: an inner one composed of Cr-poor garnet and three samples of dunite, four harzburgites, three lherzo- spinel (bottom of grain) and an outer one of orthopyrox- lites, one wehrlite, one olivine websterite, two webster- ene, spinel and amphibole. It is likely the strong devel- ites and one clinopyroxenite (Fig. 3). opment of the outer rim destroyed evidence for the A plot for selected major element abundances vs inner rim in most garnets. MgO concentration is presented in Fig. 4 in which 6 Journal of Petrology, 2018, Vol. 0, No. 0 Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018

Fig. 2. (a) Typical peridodite lens within gneisses of the Tvaerdal Complex. (b) Loose block of garnet dunite showing scattered M1 garnets. (c) Peridotite lens EH-22 with strikingly parallel garnet dunite, garnet pyroxenite and garnet-train layers of varying thick- ness. (d) Close up of the peridotite shown in (c). Garnet wehrlite (EH-22) was sampled from the thick layer on the right. The thin par- allel dark layers are garnet-rich trains. The dark colour is the result of kelyphite formation and preferential serpentinization. Note the colour variation in garnets. (e) Loose boulder showing pyroxenite layer in garnet dunite with scattered M1 garnets. (f) Pyroxenite layer within dunite. Sheared out kelyphite reaction rims around M1 garnets define a parallel planar fabric. (g) Thick sec- tion of lherzolite LM-42B showing two garnet generations; one Cr-rich (dark), the other Cr-poor (light). The yellow lines show scans GM 1 and GM 2 in Fig. 9. (h) Polished slab of wehrlite EH-22 with partially disaggregated garnet-clinopyroxene rich layer. Again, note colour variations in the garnets. (i) Thin section of garnet pyroxenite LM-38 showing evidence of fracturing and disaggrega- tion. (j) Thin section of garnet websterite HKB-2B showing evidence of fracturing and disaggregation. (k) Backscattered electron image of a garnet in HKB-2B showing asymmetrical Cr zoning. Two small grains appear to be in the process of breaking away from the larger grain. (l) Thin section of a portion of garnet websterite HKB-2B showing a relatively coherent granoblastic fabric. (m) Garnet harzburgite EH-21B showing fragmented and partially detached garnets in an olivine-rich matrix. (n) Kelyphite reaction rim composed of amphibole and spinel around garnet. The kelyphite in turn is partially enclosed by coarser amphibole grains. (o) Electron microprobe image of garnet 1 in LM-42B (Fig. 2g) showing an inner kelyphite composed of spinel and low-Cr garnet (bot- tom) and an outer kelyphite composed of spinel, orthopyroxene and amphibole. Journal of Petrology, 2018, Vol. 0, No. 0 7 Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018

Fig. 2. Continued compositions of the Tvaerdal peridotites and pyroxen- MgO, three samples of lherzolite and one of wehrlite ites are compared to those from the WGC of the that cluster around the composition of Primitive Mantle Norwegian Caledonides (Supplementary Data Table (McDonough & Sun, 1995) and four samples of pyrox- S2a). The LL peridotite and pyroxenite suite shows a enite that contain <30% MgO. pattern of chemical variation with respect to wt % MgO Although the lherzolite and wehrlite samples gener- that is common in such suites worldwide (Bodinier & ally plot in the vicinity of Primitive Mantle, their compo-

Godard, 2004); namely, an increase in MgO is accompa- sitions are highly variable, particularly for TiO2,Al2O3, nied by a pronounced decrease in TiO2,Al2O3, and CaO, Cr2O3 and CaO, which have relative standard deviations with a lesser decrease in SiO2 and FeO, and an increase of 72, 29, 75, and 59%, respectively. Interestingly, the in NiO and Cr2O3. In detail, the LL suite consists of three data from the WGC shows, with two exceptions, a gap chemical groups of rocks: seven samples of strongly from 31 to 39% MgO between the pyroxenites and peri- depleted dunite and harzburgite that contain >45% dotites. The three LL lherzolites (EH-21A, HKB-3 and 8 Journal of Petrology, 2018, Vol. 0, No. 0 Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018

Fig. 2. Continued

Camp 2–1) and one wehrlite (EH-22), fall into this gap Whole-rock Rare Earth Element (REE) patterns for (Fig. 4). peridotites and pyroxenites are given in Table 3 and plotted in Fig. 6a (normalized to primitive mantle; Trace elements McDonough & Sun, 1995). Dunites and harzburgites are The compatible element Ni decreases, and the incom- strongly depleted in HREE (01–02 x PM) and one sam- patible elements V, Y and Sr increase systematically ple of dunite displays LREE enrichment, as do two sam- with decreasing MgO (Fig. 5); such patterns are consist- ples of harzburgite. In contrast, lherzolites and wehrlite ent with enrichment of highly depleted mantle through have HREE contents ranging from 04to14 x PM, and the introduction of a pyroxenitic component. Zr and Ba garnet pyroxenites have more enriched HREE contents, also show a tendency to increase with a decrease in generally between 27 and 50 x PM. The three pyroxen- MgO, although with much greater scatter than in the ite samples exhibit a relative depletion in LREE, where- other incompatible elements. as clinopyroxenite LM-38 is LREE enriched. The Journal of Petrology, 2018, Vol. 0, No. 0 9

Ol garnet are 878 and 720, respectively. In contrast to the HKB-2A limited variation in Mg-numbers among the minerals, Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 HKB-04 HKB-6P Dunite Cr-numbers in , and garnet vary wide- Camp 2-2 ly among rock types, being higher in peridotite than in LM-42A pyroxenite. Cr-numbers for enstatite in peridotite and Camp 1-1 EH-21A Wehrlite pyroxenite are 78–169 and 47, respectively; for diop- EH-21B side are 191–330 and 26–29; and for garnet are 65– Harzburgite HKB-03 EH-22 140 and 05–15. FeO and Na2O contents in clinopyrox- Camp 2-1 ene also vary widely, as illustrated in Fig. 7b. The con- Lherzolite centrations of Cr2O3 in clinopyroxene are especially noteworthy, in view of the mixed assemblage model HKB-2C discussed subsequently. Clinopyroxenes in pyroxenites (HKB-2B and LM-38) have very low Cr concentrations, Olivine Websterite whereas clinopyroxenes in dunites and harzburgites

Clinopyroxenite (HKB-04, EH-21B and EH-22) have high Cr concentra- tions (Fig. 7b). Lherzolite and wehrlite clinopyroxenes Websterite have intermediate values and appear to define a single Opx HKB-2B LM-43 LM-38 Cpx population, which contrasts with the two populations shown by garnet (see below). Whole Rock Compositions Representative EMP scans across garnet, orthopyr- Fig. 3. Ternary diagram showing the relative proportions of ini- oxene, and clinopyroxene in pyroxenite HKB-2B and tial olivine (ol), orthopyroxene (opx) and clinopyroxene (cpx) garnet in harzburgite EH-21B are illustrated in Fig. 8a (garnet is understood to be present). The position on the dia- gram dictates the rock type identification listed in Table 1 and and b. Most minerals are unzoned to weakly zoned for used throughout this study. MgO and FeO, except for some garnets, which display a decrease in MgO and increase in FeO very close to the boundaries with kelyphite rims. Cr O and CaO, how- pyroxenites lack Eu anomalies, suggesting they - 2 3 ever, are strongly zoned in garnets especially near the lized initially as garnet- rather than - rims, as discussed further below. Orthopyroxenes in all pyroxenites. Extended trace element diagrams, normal- scans, including lherzolite LM-42B (Fig. 8 a, c and d) ized to primitive mantle, are plotted in Fig. 6b. They show a characteristic increase in Al O from core to rim show positive anomalies for Ta and Pb and negative 2 3 (05% to 20%), reflecting re-equilibration during anomalies for Nb, Zr, Hf and Ti. Sr exhibits negative exhumation. anomalies in some samples, but positive anomalies in When garnet grains were separated from crushed others. Overall the pyroxenite patterns suggest a sub- rocks, it was discovered that different colours of garnet duction signature. exist in four samples: dunite HKB-04 (violet and or- ange), lherzolite LM-42B (violet and pink), wehrlite EH- Mineral chemistry 22 (violet, red and pink) and websterite HKB-02B (red Major elements and orange). These colour differences in some samples

The major element concentrations of major phases reflect core to rim variations in Cr2O3 content (e.g. HKB- within the Tvaerdal peridotites and pyroxenites were 2B); however, in other samples (e.g. LM-42B) they rep- measured by EMP at Macquarie University, The resent two populations of garnet containing different

American Museum of Natural History and The concentrations of Cr2O3, with higher amounts of Cr2O3 University of Wisconsin-Madison. Table 4 presents occurring in dark red or violet garnets and lower average mineral concentrations based on point counts amounts of Cr2O3 in the lighter coloured garnets. of adjacent minerals in polished thin sections for some Lherzolite, sample LM–42B is particularly informative samples and epoxy mounted grain separates for others, with respect to compositional variation in garnet and so and averages of the flat portions of scans for selected its mineral chemistry was studied in more detail than in samples. Complete mineral data are presented the other samples. It contains one set of violet, high–Cr Supplementary Data Table S3a, including scans and garnet and a second set of pink, low–Cr garnet (Fig. 2g). averages for cores and rims. Additional mineral analy- Figure 9a illustrates the striking difference in Cr2O3 con- ses for peridotite #6 are given in Augland et al. (2011). centrations between the two sets. Also, in contrast to Primary minerals in the Tvaerdal ultramafic suite, like the garnets measured in most samples, the garnets in minerals in orogenic peridotite occurrences elsewhere, LM-42B are compositionally zoned (Fig. 9), although are magnesian, consisting of , enstatite, diop- EMP traverses across grains do not reveal the entire, side and -rich garnet (Fig. 7a). Mg-numbers for original compositional profiles because the grain mar- forsterite (915–926), enstatite (893–930), diopside gins were destroyed by subsequent growth of kelyphite (907– 941), and garnet (769–820) largely overlap rims. The Cr2O3 content in a high-Cr garnet (GM 1) among the different rock types, except in clinopyroxen- decreases from 6% in the core to 5% in the inner ite, LM-38, where the Mg-numbers for diopside and rim, and abruptly to 15% at the contact with kelyphite 10

Table 2: Peridotite and pyroxenite whole-rock major element chemistry, Tvaerdal complex, Liverpool Land

Sample HKB-4 Camp 2-2 HKB-2A HKB-6P EH21A LM42 Camp 1-1 EH21B HKB-3 EH22 Camp2-1 HKB-2C HKB-2B LM 43 LM-38

Lith. Dunite Hrz Dunite Hrz Lhrz Pd Hrz Hrz Lhrz Lhrz Whrl Wbs Cpxite Wbstr Cpxite

Oxides Arranged in order of decreasing Mg#

n¼2n¼2n¼2n¼2n¼2n¼2n¼2n¼2n¼2

Wt %

SiO2 430427418428478434455453457448466471481466452 TiO2 0012 0023 0020 0020 0098 0024 00175 0038 0161 0046 0046 0200 0205 0305 0619 Al2O3 0709 209 0768 119 497 180 0811 131 607 610 312 891 102124132 Cr2O3 0335 0256 0347 0184 0419 0243 0203 0300 1035 0542 0207 0354 0366 0102 0129 FeOt 701 708 745 736 604 736 735 732 609 730 783 689 688 873 106 NiO 0338 0322 0371 0348 0205 0322 0377 0358 0286 0236 0216 0115 0082 0083 0039 MnO 0096 0112 0071 0109 0174 0111 0055 0091 0171 0169 0151 0205 0227 0239 0255 MgO 481469490475385458455452358338360290257231173 CaO 0396 0548 0093 0470 144 0975 0137 0080 419 672 563 700 816 760 118 Na2O 0001 0016 0000 0031 0246 0004 0000 0000 0505 0282 0135 0096 0084 0818 0511 K2O 0000 0000 0026 0000 0040 0000 0000 0000 0029 0000 0000 0096 0005 0034 0134 P2O5 0002 0000 0052 0000 0000 0000 0063 0002 0000 0004 0000 0024 0023 0000 0182 Total 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 Mg# 924922921920919917917917913892891883869825744 Cr# 240758 232943 536 830 144133103563 425 260 236 0546 0651 Si 201199195200224203213212214210218220225218211 Al 0375 111 0407 0627 263 0954 0429 0694 321 323 165 471 538 658 697 Mg 290283296287232276274272216204217175155139104 Ca 0283 0392 0067 0336 103 0697 0098 0057 299 481 402 500 583 544 842 Mg/Si 144 142 151 143 104 136 129 129 101 0972 100 0796 0688 0639 0494

Ca/Al 0753 0353 0164 0536 0391 0730 0229 0082 0932 149 243 106 108 0826 121 Petrology of Journal Al/Si 0019 0056 0021 0031 0118 0047 0020 0033 0150 0154 0076 0214 0239 0302 0330 n ¼ 2: average of two analyses (see Supplementary Data Appendix 2). All samples contain garnet or garnet relicts. Hrz, Harzburgite; Lhrz, Lherzolite; Whrl, Wherlite; Wbs, Websterite; Cpxite, Clinopyroxenite.

08 o.0 o 0 No. 0, Vol. 2018, , Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 December 04 on user Wisconsin-Madison of University by https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 from Downloaded Journal of Petrology, 2018, Vol. 0, No. 0 11

Major Element Patterns of Tvaerdal Peridotites and Pyroxenites Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 LL Peridotite LL Pyroxenite WGR pd WGR px Primitive Mantle 52 0.8

0.6 48 2 2 0.4 SiO TiO 44 0.2

40 0.0 20 30 40 50 20 30 40 50 1.2 16

12 0.8 3 3 O O 2 2 8 Cr Al 0.4 4

0 0.0 20 30 40 50 20 30 40 50 0.4 14 0.3

10

FeO 0.2 NiO

6 0.1

2 0.0 20 30 40 50 20 30 40 50 16 1.2 12 O 2 0.8 8 Na CaO

4 0.4

0 0.0 20 30 40 50 20 30 40 50 MgO MgO

Fig. 4. Plots of major element concentrations of peridotites (large circles) and pyroxenites (large triangles) from the Tvaerdal Complex, Liverpool Land. The small symbols are analyses of peridotites (þ) and pyroxenites (x) from the Western Gneiss Complex, Norway. WGC data are from Brueckner (HKB unpublished); Medaris (unpublished); Eskola, 1921; Mercy & O’Hara, 1965; Carswell, 1966 unpublished PhD dissertation, 1968a, 1968b, 1973, 1983; Carswell et al., 2006; Lappin, 1974; Moore & Kvale, 1977; Al- Samman, 1985 unpublished PhD dissertation; Osland, 1997 and Beyer et al., 2006. pd, peridotite; px, pyroxenite.

(Fig. 9a). In contrast, the Cr2O3 content in a low-Cr gar- the opposite direction. LM-42B may be different from net (GM 2, Fig. 9a) increases from 06% in the core to most other studied samples in that some garnets are 10% at the rim. Most other garnet grains in this sam- systematically zoned with respect to major elements ple show either a flat pattern or a core to rim increase in (Fig. 9b). Most zoning is subtle, for example Cr-rich GM Cr2O3 and in several cases (e.g. HKB 1, GM 3 and HKB 1 shows a slight decrease from core to rim in pyrope 3) this increase is asymmetric with an increase in Cr2O3 and grossularite components and an increase in alman- in one direction but no or only a very weak increase in dine component. The zoning of major elements in some 12 Journal of Petrology, 2018, Vol. 0, No. 0

Trace Element Patterns of Tvaerdal Peridotites and Pyroxenites Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 LL Peridotite LL Pyroxenite Primitive Mantle

2500 300 Ni V 200 1500 ppm 100 500

0 0

30 200 Y Sr 20

ppm 100 10

0 0

80 25

60 Ba Zr 15 40 ppm

20 5

0 0 10 20 30 40 50 10 20 30 40 50 MgO (wgt %) MgO (wgt %)

Fig 5. Plots of selected trace element concentrations from whole-rock Tvaerdal Complex peridotites, Liverpool Land. low-Cr garnets differs from that in high-Cr garnet, with Trace elements (REE) pyrope and grossularite slightly increasing from core to Selected peridotite minerals were analysed for trace rim, and almandine decreasing. Other garnets show flat elements by LA-ICP-MS and their REE abundances rela- patterns. These patterns are enhanced by plotting the tive to chondritic mantle are plotted in Fig. 11a (clino- variation in Mg-number across the garnet (Fig. 9b), pyroxene) and 11 b (garnet). REE patterns for each where it can be seen that only the high-Cr garnet GM 1 separate sample are plotted in Supplementary Data shows a marked decrease in Mg-number from core to Table S3B. Clinopyroxenes from most samples show rim. The rest of the garnets show flat patterns or only the classic ‘hooked’ pattern, with positive slopes for the minor variation. Nevertheless, the figure does illustrate Light Rare Earth Elements (LREE) followed by negative that there is significant variation in Mg-number be- slopes for the Middle Rare Earth Elements (MREE) and tween grains. Heavy Rare Earth Elements (HREE). In contrast, the Amphibole is localized in and around the kelyphite clinopyroxene REE pattern for wehrlite EH-22 displays a rims and consists of pargasite in harzburgite, lherzolite, negative slope throughout. and wehrlite, and tschermakite in websterite. The spinel REE patterns for garnet display LREE depletion, associated with garnet shows a large range in compos- which is a common feature in garnet from orogenic ition between (Fig. 10a) and within samples (Fig. 10b). peridotites. Garnets from four rock samples have com- in LM-42B within garnet and orthopyroxene, parable MREE to HREE patterns (i.e. Sm to Lu), which and the olivine-rich matrix are relatively Cr-rich, where- are either flat or gently sloping. However, slopes for the as those within the kelyphite reaction rims are more LREE vary among the four samples, depending on rock enriched in Al (Fig. 10b), although there is a large over- type and the colour of garnet. The LREE slope is steep lap in compositions between the domains. for a single garnet separate from lherzolite LM–42B and ora fPetrology of Journal

Table 3: Peridotite and pyroxenite whole-rock trace element chemistry (in ppm), Tvaerdal Complex, Liverpool Land

Sample HKB-4 Camp 2-2 HKB-2A HKB-6P EH-21A LM42 Camp 1-1 EH-21B HKB-3 EH-22 Camp 2-1 HKB2C HKB2B LM 43 LM38 Lithology Dunite Hrz Dun Hrz Lhrz Pd Hrz Hrz Lhrz Hhrz Whrl Wbstr Cpxite Wbstr Cpxite 08 o.0 o 0 No. 0, Vol. 2018, , Li 218 128 250 0961 626 0980 467 438 52251 694 471 329 771 Be 0007 0047 0047 0026 0145 0016 0308 0355 0727 Sc 597 864 375 727 130916 665 598 30 2012233013 349364614 Ti 987 172 103 188 528 224 886 232 901 308 329 1145 1307 1846 3564 V3023952597 300579400248326 133 804590978 119 174 248 Cr 1954 1501 1468 1130 2387 1501 1192 1361 5914 2663 1266 1787 1789 689 531 Mn 736 701 479 761 1221 745 378 670 1200 1268 993 1528 1810 1817 2008 Co 109 915 112 102 683 107 107 954818799556334409 952444 Ni 2478 2172 2520 2454 1402 2285 2542 2274 1909 1618 1520 864 557 649 289 Cu 591 128245 711 174127363 787 977416586 39352 117113 Zn 376357576415313408486363348258351249171198177 Ga 0548 123 0712 089 353 129 074 119 338 362 207 560 610 864 789 Rb 0073 0058 0686 0089 179 0098 0240 011 141 011 0049 267 053 0713 767 Sr 333 306 435 400 122640 561 110 6063981805396017 508 227 Y0460 0741 0674 0687 410 0973 0750 092 60251 164 120148154319 Zr 104 140 0420 136 828 200 0632 298 109297 122 993 981 215128 Nb 0057 0092 0169 0035 0062 0089 0210 0042 036 007 0100 072 062 0379 100 Mo 0107 0355 0067 0246 0150 0132 0076 0186 016 017 0147 033 016 022 030 Cd 0008 0013 0013 0012 0026 0015 0008 0011 007 003 0027 005 006 0076 009 Cs 0030 0026 0063 0011 170017 0202 0044 056 001 0054 060 038 014 051 Ba 448 112279661 300703 338411 710200145 670105118375 La 0079 0242 0450 0092 0208 0139 0521 0190 137 066 0359 061 056 068 897 Ce 0170 0634 136 0239 0710 0377 142 0416 466 106 0669 243 244 221 262 Pr 0028 0093 0192 0042 0147 0064 0187 0077 093 013 0083 046 049 041 362 Nd 0133 0408 0820 0206 0949 0327 0779 0451 496 054 0352 256 278 237 165 Sm 0056 0100 0152 0070 0380 0116 0178 0156 080 015 0126 109 126 087 466 Eu 0023 0037 0041 0026 0128 0043 0062 0045 023 006 0047 047 053 034 136 Gd 0071 0119 0140 0095 0499 0147 0196 0167 081 023 0201 147 177 152 576 Tb 0012 0020 0020 0016 0087 0024 0030 0025 013 004 0038 025 031 0298 101 Dy 0072 0127 0110 0105 0581 0156 0154 0146 088 035 0269 168 210 218 556 Ho 0015 0028 0022 0024 0131 0035 0026 0030 021 009 0062 039 049 052 102 Er 0042 0082 0058 0070 0384 0099 0060 0083 063 027 0181 115 146 154 260 Yb 00447 0087 0053 0076 0387 0093 0046 0078 061 030 0185 116 148 153 223 Lu 00073 0014 00083 0012 0061 0015 0007 0013 010 005 0027 018 023 024 034 Hf 00223 0040 00052 0036 0212 0054 0011 0086 039 010 0049 022 021 0664 030 Ta 0278 0030 00051 0025 0194 0107 0026 0014 025 022 0178 118 004 0817 008 Pb 241 139 208 142 322 266 Th 00136 0039 0020 0007 0008 0012 0035 0016 003 004 0043 005 005 0060 061 U00049 0009 00166 0003 0002 0009 0006 0307 001 001 0006 001 001 0028 016

13 Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 December 04 on user Wisconsin-Madison of University by https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 from Downloaded 14 Journal of Petrology, 2018, Vol. 0, No. 0

LM-38 Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 10 (a) Clinopyroxenite HKB-2A

Pyroxenite HKB-2A & 1 Camp 1-1 Lherzolite EH-21A DMM Camp 2-1 EH-22 Wehrlite

0.1 HKB-4 Dunite & Harzburgite Normalized to Primitive Mantle

0.01 La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu 100 (b)

10 Clinopyroxenite

Pyroxenite 1 Lherzolite

Wehrlite

0.1 Dunite & Harzburgite Normalized to Primitive Mantle

0.01 Cs Rb Th U Nb Ta La Ce Pb Pr Sr NdZr Sm Eu Ti Gd Tb Dy Ho Y Yb Lu

Fig. 6. (a) Rare Earth Element and (b) Extended trace element patterns of Tvaerdal ultramafic rocks, normalized to Primitive Mantle (McDonough & Sun, 1995). comparably steep for violet, red and orange garnets exchange geothermometers for olivine and garnet (Wu from dunite HKB–04 (notwithstanding the ‘spike’ in Ce & Zhao, 2007) and orthopyroxene and garnet (Harley, for one of the violet garnet separates). In contrast, in 1984), and the Al-in-opx geobarometer (Nickel & Green, wehrlite EH–22, the LREE slope for violet garnet is 1985). Application of these methods to the core compo- steep, whereas the LREE slope for red and orange gar- sitions of adjacent minerals (Fig. 12) yields P–T esti- nets is flat. Interestingly, this wehrlite sample is the only mates of 369 kbar/804C for dunite (HKB–04), 318 kbar/ one among the four that shows whole-rock LREE en- 848C for harzburgite (EH21B) and 244 kbar/818C for richment. The LREE slope for garnet from olivine web- lherzolite (HKB–03). Websterite (HKB-2B) yields 329 sterite HKB–2C is intermediate between those of the kbar/864C, which is comparable in temperature to the other samples. Nevertheless, the flat or even slightly peridotites, but lower in pressure. Wehrlite EH-22 and negative slopes for the LREE shown by some garnets clinopyroxenite LM-38 did not give reliable pressure suggest LREE enrichment. estimates, but an assumed pressure of 30 kbar yielded temperatures (900C and 830C, respectively) con- P–T calculations sistent with the other measurements. These values A combination of geothermobarometers (Supplemen- yield an overall estimate of 305 kbar/801C, which lies tary Data Table S3C) was utilized to estimate P–T condi- just below the very high maximum average P–T (351 tions for the LL ultramafic suite, including Fe–Mg kbar/877C) calculated for an external orthopyroxene- Journal of Petrology, 2018, Vol. 0, No. 0 15

Table 4: Mineral major element chemistry, Tvaerdal Complex, Liverpool Land Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 SiO2 TiO2 Al2O3 Cr2O3 FeO MnO MgO CaO Na2O K2O NiO Total Mg# Cr# Ca# LM-42B Lherz ol scan (20) 4050010 0001 0011 826 0119 5010008 0008 0454 994915869 opx scan (90) 5660044 0962 0140 564 0117 3530220 0024 0089 991918889 cpx scan (25) 5320108 223 0845 294 0065 1612141 146 0051 983907203465 hi-Cr grt scan (109) 4090070 200412 995 0501 191475 0029 0011 994773122122 lo-Cr grt scan (200) 4170033 2320355 901 0432 209373 0019 0014 994805102 937 LM-38 Cpxite cpx disc (2) 5460056 149 0067 403 0057 1632240 128 0011 0068 1004878293 465 amph disc(1) 4350055 1710212 1320421 2680579 0017 0009 0036 1019783083 120 grt disc (2) 4220076 2330183 1290463 186458 0021 0009 0017 818720053 113 HKB-2B Websterite ol no ol opx scan (6) 5650071 155 0115 609 0122 3430267 0010 0009 0093 9909 909474 cpx scan (7) 5350115 141 0055 230 0063 1732441 0282 0003 0042 9941 930254 486 amph spot (1) 4650395 1480324 327 0060 1771295 127 0029 0033 9745 906144 322 grt scan (32) 4140069 233043 1030513 192499 0005 0006 0009 1002769123 013 red grt spot (1) 4250070 2370577 971 0420 198511 0002 0006 0037 1020784161 127 orng grt spot (1) 4260077 2420086 992 0427 203432 0010 0009 0000 1019785024 107 HKB-3 Lherz ol spot (3) 4180008 0005 0009 794 0169 5050017 0004 0010 0409 1008919564 opx disc (2) 5750061 157 0263 540 0178 3560269 001 0014 0061 1009922101 cpx spot (2) 5460134 276 148 184 0068 1632251 134 0009 0022 1011941264482 grt disc (1) 4230133 227238 885 0447 204435 0003 0006 0003 1016805655 110 amph disc (2) 4530436 127180 320 0020 1861253 255 0177 0099 974912872 307 HKB-4 Dunite ol spot (3) 4160006 0006 0010 727 0116 5100014 0019 0011 0396 1005926515 opx disc (2) 5860027 0552 0167 487 0110 3620204 0026 0015 0104 1009930168 cpx spot (3) 5470049 174 128 202 0054 1642246 149 0000 0060 1003936330479 violet grt disc (2) 4160086 206495 908 0463 194533 0015 0000 0000 1015792139135 orng grt disc (2) 4260142 231159 934 0434 207420 0029 0000 0012 1021798441 104 EH-21B Harzburg ol no analyses opx scan (14) 5780053 0956 0121 502 0108 3560187 0009 0003 0087 9994 927782 cpx spot (3) 5430193 351 154 189 0037 1532097 236 0004 0026 1002935227480 amph spot (2) 4400746 148118 354 0027 1791183 323 0468 0133 9782 900507 299 grt scan (32) 4260060 237105 844 0405 216350 0010 0006 0009 1013820289 871 grt spot (1) 4240179 221305 897 048 206403 0025 0005 0004 1018803847 101 grt spot (1) 4240097 212438 874 047 212339 0001 0000 0000 1020812122852 EH-22 Wehrlite ol spot (3) 4100002 0004 0010 989 018 4890011 0002 0001 041 1004898602 opx spot 1 5690006 263 0092 724 0206 3380297 0009 0002 012 1013893228 cpx scan (10) 5420067 169 0669 232 005 172229088 0006 004 1001930210472 amph spot (7) 4540302 133143 427 006 183122281 0079 010 982884672 297 grt scans (21) 4220062 229157 998 048 196487 0010 0003 001 1016778441 122 violet grt disc (1) 4250072 219354 853 0480 202496 0008 0006 0000 1022808977 125 red Grt disc (1) 4270102 233159 938 0503 202441 0015 0003 0000 1022793439 111 pink grt disc (1) 4230082 234139 932 0457 202448 0000 0000 0018 1017795382 112

Spinel compositions vary too much for averaging. See Appendix 3. Mineral abbreviations as in Table 1. Parenthesis, number of analyses; scan, grain traverse; spot, single analysis of grain thin section. disc, single grains mounted and polished in epoxy and analysed. bearing eclogite (EH-15B) in the nearby gneisses Application of the same geothermobarometers to (Brueckner et al., 2016). rim compositions of minerals yields pressures of 15–25 Two measurements, from two different laboratories, kbar and temperatures of 750–800C(Fig. 12). The ori- of high-Cr garnets from LM-42B yielded similar values ginal rims of the garnet were replaced by kelyphite, so (297 kbar/784C and 321 kbar/789C). These pressures these P–T values do not represent the values of the ori- are similar to the pressures defined by the ultramafic ginal rims, which were probably lower than the core samples discussed above, but the temperatures are values. Regardless, the present P–T rim values indicate lower. A low-Cr garnet, when combined with the same that the initial stage of decompression was largely iso- forsterite and enstatite compositions as the high-Cr gar- thermal, similar to the pattern shown by the external net, yields the same P–T (386 kbar/897C) as the ortho- eclogite. pyroxene eclogite. However, if this sample contains The core P–T values of most samples fall more or mixed peridotite/pyroxenite assemblages, as we dis- less along the retrograde path defined by the external cuss further below, all calculations from LM-42B should eclogite (Fig. 12). We suggest that these samples re- be treated with caution. equilibrated to varying degrees during post-Scandian 16 Journal of Petrology, 2018, Vol. 0, No. 0

(a) C which is seemingly confirmed by Cr-profiles across 50 Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 Cpx grains (Figs 8 and 9). Because the REE are, like Cr, triva- clinopyroxenite lent, it is likely that REE concentrations vary from core 60 websterite wehrlite to rim as well, precluding the possibility that a given lherzolite sample would yield a single recrystallization age. 70 harzburgite dunite Nevertheless, garnet grains were separated by colour 80 under a binocular microscope with the hope of obtain- ing at least some time constraints on core vs rim forma- 90 Grt tion. Due to the nature and amount of separated garnet, Opx it was only possible to divide garnet separates into two M F fractions and it proved difficult to obtain absolutely 90 80 70 60 50 EH-22 pure end-member separates. core Parts of websterite sample HKB-2B display clean (b) Cr2O3 granoblastic fabrics with little evidence of fracturing (Fig. 2l). Nevertheless, repeated analyses of garnets Dunite HKB-04 from this websterite plot as scattered points and fail to Lherzolite LM-42B GM Lherzolite HKB-03 Lherzolite LM-42B-HKB Websterite HKB-02B show a consistent relationship between age and garnet Harzburgite EH-21B Clinopyroxenite LM-38 colour (Fig. 13a). Despite this scatter, the data give inter- Wehrlite EH-22 pretable results in that the garnet that generates the steepest slope gives a clinopyroxene–whole-rock–gar- net isochron age of 432 6 5Ma (2 r, MSWD ¼ 0003) whereas the garnet that generates the shallowest slope gives a near-isochron age of 409 6 6 Ma (MSWD ¼ 13). Thus, despite the scatter, the ‘ages’ from this webster- ite, indicate that garnet formation in the pyroxenites FeO Clinopyroxene Na2O probably began during the earliest stages of the Compositions Scandian Orogeny (430 Ma Sm–Nd isochron ages were obtained from garnet pyroxenites of the WGC by Fig. 7. (a) Representative compositions of orthopyroxene (opx), clinopyroxene (cpx) and garnet (grt) in the system Ca– Spengler et al., 2009). The 409 Ma age of the ‘youngest’ Mg–Fe (C–M–F), identified by rock type. (b) Representative garnet is within error or very close to the Scandian ages compositions of clinopyroxenes in the Cr2O3–FeO–NaO sys- obtained from Tvaerdal eclogites (Brueckner et al., tem. Note the positions of LM-42B and EH-22 clinopyroxenes 2016), suggesting a shared history. This conclusion is with Cr2O3 contents intermediate between clinopyroxenes from dunites and harzburgites and clinopyroxenes from pyrox- supported by the results from pyroxenite LM-43 enites. These intermediate compositions may be the result of (Fig. 13b) from which several garnet fractions were ana- re-equilibration between the two end-members (see discussion lysed along with clinopyroxene, amphibole and whole- for further details). rock. The different fractions are slightly scattered, but cluster around an age of 399 6 18 Ma (MSWD ¼ 18). retrogression. The garnets from HKB-04 and high-Cr Some of this scatter is attributable to the different col- LM-42B fall off this trend and plot at lower tempera- ours of the garnet fractions (the exact colour of grt2 was tures. Both garnets are high-Cr garnets, which are not noted). The most precise age of 4036 6 75 (MSWD believed to have formed during a pre-Scandian event ¼ 056) is defined by grt2 and the lighter coloured (see discussion below). Their P–T values may reflect (‘pale’) garnet fraction, whereas the darker red garnet conditions during this older metamorphism. defines a younger age (3878 6 75 Ma; MSWD ¼ 084). Alternatively the standard geothermobarometers may Despite the lack of a true isochron age, the results from not apply accurately to garnets with high Cr-contents or both pyroxenite samples are consistent with pyroxenite none of the garnets were in equilibrium with the adja- HP/UHP equilibration during the Scandian Orogeny. cent minerals as a result of mechanical mixing. The data from the garnets in wehrlite EH-22 (Fig. 13c) and peridotite LEA 08–04 (Fig. 13d) complicate this rela- tively straightforward picture. As in websterite HKB-2B, Geochronology the darker garnet fractions in both samples collectively Mineral separates from websterites HKB-2B and LM-43, define steeper best-fit lines than those defined by pale wehrlite EH-22 (Table 1) and pyroxenite LEA 08–04 were garnets. The clinopyroxene–whole-rock–garnet best-fit analysed for their Sm, Nd and Sr concentrations and Sr lines with the steepest slopes generate apparent ages and Nd compositions (Table 5, Fig. 13). of 808 6 48 Ma in EH-22 and 729 6 18 Ma in LEA 08–04. Obtaining reliable ages was complicated by the fact that These ages do not match pre-Scandian orogenic events all four samples contain garnet separates with mixtures recorded in the associated gneisses (i.e. the 16Ga of different colours, including violet, dark red, pink, or- Labradoran orogeny and the 10 Ga Grenville orogeny) ange and pale red. These colour variations were initially and should be considered minimum estimates that interpreted as reflecting core to rim chemical variations, were partially re-set during Scandian crystallization. Journal of Petrology, 2018, Vol. 0, No. 0 17

MgO FeO CaO Al O Cr2O3 2 3 Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 Olivine Orthopyroxene Clinopyroxene Garnet 40 40

30 (a) HKB-2B Grt Websterite 30 (b) EH-21B Grt Harzburgite 20 20

10 10

6 Grt Grt Grt

5 5 Weight % Weight

4 4 Opx Cpx Grt 3 3

2 2

1 1

0 10 50 100 140 0 3500 0 8000 Steps (Total Scan 8000 µ) Distance (microns) 50 50

30 30 (c) LM-42B Grt Lherzolite 10 6 10 (d) LM-42B Grt Lherzolite 4 Cpx Opx Weight % Weight 6 2 Opx Ol Opx Ol Opx 4 1 2

0 0 0 1000 2000 3000 0 400 800 1200 Distance (microns) Distance (microns)

Fig. 8. Rim-to-rim EMP scans for major oxide concentrations for selected minerals from Tvaerdal peridotites and pyroxenites. (a) Typical scan across orthopyroxene, clinopyroxene and garnet in pyroxenite HKB-2B, showing flat patterns for most elements ex- cept for Al in orthopyroxene and Cr and Ca in garnet (Cr concentrations deleted for clinopyroxene and orthopyroxene for clarity). Also note different Cr contents for different garnets. (b) Scans across garnets in EH-21B. Again, note the flat patterns for most ele- ments except Cr and Ca. (c, d). Scans across olivine, orthopyroxene and clinopyroxene in lherzolite LM-42B. Major elements show flat patterns except for Al2O3 in orthopyroxene in contact with olivine (c) and clinopyroxene (d). Some measurements have been deleted for clarity where the data are too crowded. Data for these scans are given in Supplementary Data Table S3a.

The pale garnets with the shallowest slopes give appar- indicate that the pale garnets in this sample formed dur- ent ages of 647 6 13 Ma in EH-22 and 384 6 32 Ma in ing the Scandian orogeny, consistent with the results LEA 08–04. The pre-Caledonian 647 Ma age in EH-22 is from the garnet pyroxenites. considered a maximum age of the younger of two gar- net generations, which proved impossible to complete- ly separate from the older darker garnets. The two pale Sr–Nd isotope variations garnets in LEA 08–04, however, appear to have been 87Sr/86Sr and 143Nd/144Nd values of clinopyroxenes and successfully separated from older garnets. They give from nine Tvaerdal garnet peridotites and Scandian ages of 384 6 32 Ma and 425 6 15 Ma (not pyroxenites (Table 5) and two Tvaerdal eclogites plotted) when regressed separately with two clinopyr- (Brueckner et al., 2016) are plotted in Fig. 14, where they oxenes and the whole-rock. The apparent age calcu- are compared with the isotopic compositions of clino- lated using both pale garnets (414 6 44 Ma) is not an from WGC peridotites, pyroxenites isochron (MSWD ¼ 14), but together the apparent ages (Brueckner et al., 2010) and eclogites (Brueckner, 18 Journal of Petrology, 2018, Vol. 0, No. 0

LM-42B Grt Profiles Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 (a) 6 6 GM 1 5 5 HKB 1 3 O

2 4 4 HKB 2 HKB 4 3 3 Wt % Cr 0.9 GM 2 GM 3 2 Scale GMfor 1 only

0.5 1 HKB 3 0.1 HKB 5 2000 4000 6000 8000 Distance (microns) 84 (b) Cr-rich 7.7 mm (GM-1) Cr-poor 4.3 mm (GM-2) 82 Cr-poor 6.9 mm Intermediate Cr 3.6 mm 80 Intermediate Cr 2.4 mm

78

Mg# GM 2 76

74

72 GM 1 70 0.0 0.5 1.0 Distance, normalized

Fig. 9. (a) Rim to rim Cr2O3 traverses across garnets within lherzolite LM-42B. Note large gap in Cr2O3 concentrations between Cr- rich and Cr-poor garnets. Grt 1 is the only garnet scanned in LM-42B that shows a decrease for core to rim. The remaining garnets show increasing Cr concentrations towards the rims, many of which are asymmetric, suggesting that they are fractured grains. (b) Rim-to-rim scans showing changes in Mg-number across garnet grains. The use of Mg-number and a restricted scale (70–84) enhances the relatively slight variations across the grains, which nevertheless are flat for many grains (GM 2) except for the Cr-rich garnet GM 2. The significant variations in Mg content between grains also suggests that the garnets are fragments of larger grains that had significant gradients in Mg concentration. HKB and GM refer to the analysts.

unpublished data). Most WGC ultramafic clinopyrox- fluids derived either from a subducting slab or from the enes define a steep array with variable 143Nd/144Nd and surrounding crustal rocks. restricted, low 87Sr/86Sr. Such an array implies that most WGC clinopyroxenes formed or equilibrated with- in the mantle (Brueckner et al., 2010). In contrast, clino- Re–Os pyroxenes from the Tvaerdal peridotites and The results of a LA-ICP-MS Re–Os study of sulfides pyroxenites plot along the trend defined by WGC eclo- from the Tvaerdal peridotite are listed in Table 6 (the gites, though not to the very high 87Sr/86Sr ratios in complete data set is included in Supplementary Data some enriched WGC eclogites. This trend suggests that File 4). Only four samples, EH-21 and 22, HKB-3 and ultramafic Tvaerdal clinopyroxenes with initial mantle Camp 2–1 contained sulfides with measurable Os and values (high 143Nd/144Nd, low 87Sr/86Sr) equilibrated to relatively low Re concentrations. Figure 15a is a Re–Os varying degrees with high 87Sr/86Sr, low 143Nd/144Nd isochron diagram that plots all measured sulfides with Journal of Petrology, 2018, Vol. 0, No. 0 19

incl. in grt incl. in opx blebs in grt LM-42B HKB-04 HKB-03 Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 Grt Lherz Grt Dunite Grt Lherz X outside kely + inside kely 1.0 1.0 (a) (b) 0.9 0.9

0.8 0.8 LM-42B 0.7 0.7

0.6 0.6

0.5 0.5 X +X Cr/(Cr+Al) Cr/(Cr+Al) 0.4 0.4 +X +X 0.3 0.3 + X X++ X+ 0.2 0.2 ++ + 0.1 0.1

0.0 0.0 1.0 0.8 0.6 0.4 0.2 0.0 1.0 0.8 0.6 0.4 0.2 0.0 (Mg)/(Mg+Fe2+) (Mg)/(Mg+Fe2+)

Fig. 10. Mg# versus Cr# in spinels from (a) peridotites LM-42B, HKB-3 and HKB-4 and (b) within peridotite LM-42B. See text for discussion.

187Re/188Os ratios less than 5 (the less 187Re, the more enriched at that time (Beyer et al., 2004, 2006, 2012) and likely the measured 187Os/188Os ratio is reliable). The the 16 Ga scatterchron defined by some sulfides data are widely scattered (i.e. 13Ga6 0.4, MSWD ¼ (Fig. 15a). But the evidence for a 16 Ga event in the 325) both for individual samples and when all data are Tvaerdal peridotites is not strong and a younger enrich- considered. A significant number of analyses (35 of ment would have resulted in the same reduction in the 53) fall on or near a best-fit line that defines an age of scatter of the model ages. A few pyroxenite sulfides do 16 Ga, which is essentially the same time a mantle en- generate reasonable model ages without assuming Re- richment event re-fertilized WGC peridotites (Beyer addition, but only those from EH-21 define straightfor- et al., 2006). A few sulfide analyses (10) plot along a ward model ages (CHUR) around 1.6 Ga. The data from shallower trend, probably related to crustal fluid intro- Camp 2–1 sulfides are particularly interesting (Table 7) duction during Scandian recrystallization. The remain- in that the Re/Os ratios are very low and the Os signals ing (8) analyses have high 187Os/188Os ratios and lines very strong. One analysis of sulfide in Camp 2–1-A has fitted through each point yield meaningless ages older a particularly low 187Os/188Os ratio (009973, fraction- than the Earth. Several of these analyses, particularly ation corrected) and generates very old Archean ages from garnet lherzolite Camp 2–1, are robust with strong regardless of assumptions about Re-addition or deple- Os signals and resultant small errors (shown by small tion. The remaining analyses from this sample also gen- error ellipses) and very low 187Re corrections. They de- erate Archean ages. fine reasonable ages if it is assumed that the samples underwent a Re enrichment event (see below). Model ages (Supplementary Data Table S4) are DISCUSSION widely scattered, from ages older than the Earth to fu- Origin of the Liverpool Land peridotites and ture ages. All future ages are eliminated when analyses pyroxenites with 187Re/188Os >05 and Os concentrations <10 ppb The Tvaerdal peridotite lenses have all the characteris- are excluded, but ages older than 455 Ga remain. tics of mantle wedge peridotites (Brueckner & Medaris, These ages are modified if it is assumed that Re was 2000) and belong to the Mg–Cr type garnet peridotite introduced into Re-free Archean sulfides during an en- class initially proposed by Carswell et al. (1983) for the richment event (Table 6). Figure 15b plots ages that re- WGC of Norway and orogenic peridotites worldwide sult if enrichment occurred at 16 Ga. The ages range (Bodinier & Godard, 2004). Most lenses are highly between 21–41 Ga, with possible clusters between 24– depleted garnet harzburgite and garnet dunite, the resi- 28 and 35–41 Ga. Assigning enrichment at 16Ga is dues of extensive melt extraction. We suggest that based on evidence that peridotites in the WGC were melting removed most of the clinopyroxene, but not the 20 Journal of Petrology, 2018, Vol. 0, No. 0

100 Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 (a) Clinopyroxene

10 pyroxenite dunite lherzolite

1

0.1 HKB-2B 1

REE/Chondrite HKB-04 1 HKB-2B 2 HKB-04 2 HKB-2C 1 HKB-04 3 0.01 HKB-2C 2 EH-22 1 HKB-03 1 EH-22 2 HKB-03 2 EH-22 3 LM-38 1 0.001 La Ce Pr Nd Sm EuGd Tb Dy Ho Er Tm Yb Lu 100 (b) Garnet

pyroxenite 10 dunite lherzolite

1

HKB-2B red

REE/Chondrite 0.1 HKB-04 orng HKB-2B orng HKB-04 red HKB-2C HKB-04 violet HKB-3 HKB-04 violet EH-21 1 0.01 EH-21 2 EH-22-pink LM-38 1 EH-22 red LM-38 2 EH-22 violet LM-42B 0.001 La Ce Pr Nd Sm EuGd Tb Dy Ho Er Tm Yb Lu

Fig. 11. REE patterns for (a) clinopyroxene and (b) garnet separated from Liverpool Land peridotites and pyroxenites, normalized to chondritic mantle (McDonough & Sun, 1995). Most clinopyroxenes show similar hooked patterns except for the straight pattern dis- played by EH-22, a garnet wehrlite, which also shows a weak positive Eu anomaly for one grain. The garnet patterns are variable with some showing evidence of LREE enrichment and others none. Light coloured garnets from wehrlite EH-22 (filled triangles) show evidence of LREE enrichment but the dark one does not.

Cr-rich garnet that is scattered throughout the peridotite Ga errorchron generated by some sulphides on a Re– host (Fig. 2b and f) or occurs as thin trains separated by Os isochron diagram (Fig. 15a). More recent Re intro- fine grained, recrystallized olivine (Fig. 2c and d). Re–Os duction could also generate reasonable Re–Os model results indicate that melt extraction episodes occurred ages and the LREE enrichment patterns of whole-rock throughout the Archean and Early Proterozoic (Fig. 15a and garnet and ‘crustal’ Sr–Nd isotope ratios of some and b) and it is possible that the garnets date back to clinopyroxenes (Figs 11 and 14) are consistent with a these ancient events. Subsequent refertilization involv- Scandian refertilization event that could have occurred ing Re addition is required to avoid Re–Os model ages in the mantle wedge from subduction zone fluids or older than 455 Ga. One such event could have been from crustal fluids after the peridotite lenses were through the introduction of pyroxenite melts, as has inserted into Tvaerdal Complex. been proposed for the peridotites in the WGC (Beyer et al., 2006). WGC garnet pyroxenites give Proterozoic recrystallization ages (Brueckner et al., 2010 and refer- Pyroxenite evolution ences therein), but evidence for Proterozoic refertiliza- The evolution of most Tvaerdal pyroxenites is relatively tion of the Tvaerdal peridotites is restricted to the 16 straightforward, albeit with some caveats. Some Journal of Petrology, 2018, Vol. 0, No. 0 21

50 Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018

45 150

40 HKB-04

Eclogite 120 core P-T 35 HKB-2B

LM-42B EH-21B 30 Hi-Cr EH-22 LM-38 = Dunite 90 25 = Harzburgite HKB-03 = Wehrlite P (kbar) = Lherzolite near

20 Depth (km) rim Eclogite = Websterite P-T rim P-T 60

15 = Clinopyroxenite = Eclogite 10 30

5

0 400 600 800 1000 1200 T°C

Fig. 12. P–T estimates for selected Tvaerdal peridotites and pyroxenites based on major element measurements from mineral cores and rims. Pyroxenite HKB-2B and harzburgite EH-21B plot close to the P–T estimate for an orthopyroxene-bearing eclogite in the host gneisses (Brueckner et al., 2016). The pressures of pyroxenite LM-38 and wherlite EH-22 were set at 30 kbar in order to obtain temperatures. Peridotites, some of which are believed to contain two garnet generations, plot as more scattered points. See Supplementary Data Table S3C for complete data. garnets in websterites HKB-2B and LM-43 give the (Figs 2k and 8a) as well as in most garnets from LM-42B same Scandian Sm–Nd ages (Fig. 13) as the eclogites in (Fig. 9a). This increase is consistent with the progres- the enclosing gneisses (Brueckner et al., 2016), indicat- sive destruction of Cr-rich spinel to form Cr-garnet as P– ing that the garnet pyroxenites recrystallized at the T conditions increased during subduction through the same time as the eclogites. This age correspondence reaction (O’Neill, 1981) suggests that the pyroxenites were inserted into the 1. orthopyroxene þ clinopyroxene þ Cr-spinel ! oliv- Tvaerdal crust as passengers within their peridotite ine þ Cr-rich garnet hosts and continued recrystallizing at progressively higher pressures and temperatures during the subduc- Trivalent Cr diffuses at a significantly slower rate in tion of the Tvaerdal Complex into the mantle (the ‘pro- garnet than bivalent elements (diffusion coefficients are grade type’ of Brueckner & Medaris, 2000). However, 05–15 log10 units smaller, Carlson, 2012). Therefore the only pyroxenite (websterite HKB-2B) analysed for the re-equilibration that erased the prograde patterns of P-T (Fig. 12) plots at a slightly lower pressure and tem- bivalent Mg, Fe, and Mn in most garnets was not in- perature than that defined by the eclogites (Fig. 12), tense enough to homogenize the profile defined by Cr. suggesting that it may have re-equilibrated as it was Caþþ would normally be expected to re-equilibrate to exhumed from mantle depths. Several peridotite sam- the same flat profiles defined by the other bivalent cati- ples plot close to the exhumation path defined by the ons, but instead it replicates the Cr pattern (Fig. 8a and Tvaerdal eclogites, suggesting that they also equili- b). Reaction (1) forms garnet that is progressively richer brated during exhumation. in Ca at the expense of clinopyroxene. The rimward in- Re-equilibration at peak P–T and/or during subse- crease in Ca in garnet is probably the result of this quent retrogression appears to have erased most evi- transfer. Also, Ca diffuses more slowly in garnet than dence for the prograde evolution of the pyroxenites, the other bivalent cations (diffusion coefficients for Ca þþ consistent with the flat concentration patterns across are 05 log10 units smaller than for Mg, Fe , and Mn; garnet grains for most major elements (Figs 8 and 9). Cr Carlson, 2012; Fig. 10). If this reasoning is correct, the Cr and Ca however are exceptions. Their concentrations and Ca profiles are the only elements that reflect the increase from garnet core to rim in websterite HKB-2B prograde evolution of the Tvaerdal garnet pyroxenites. 22 Journal of Petrology, 2018, Vol. 0, No. 0

Table 5: Isotope data from garnet peridotite and pyroxenite Tvaerdal Complex, Liverpool Land Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 Sr ppm 87Sr/86Sr 2 sigma Sm ppm Nm ppm 147Sm/144Nd 143Nd/144Nd 2 sigma

Camp 2-1 cpx LDEO 7270704206 0000010 113 0513308 0000018 HKB-2C cpx LDEO 734 0709850 0000010 242 827 0177 0512253 0000009 HKB-3 cpx LDEO 289 0705026 0000007 313 2040093 0511414 0000007 HKB-4 cpx LDEO 130 0705963 0000010 211 698 0183 0513038 0000011 LM-42B cpx SB nm 0704272 0000002 1031 365 0171 0513094 0000008 LM-42B amph SB nm 0704444 0000013 nm nm nm 0513201 0000016 AA 08-47 cpx SB nm 07072713 0000008 212 1113 0115 0512105 0000008 HKB 2B cpx 1 LDEO 197 0706722 0000010 256 946 0164 0512234 0000011 cpx2 LDEO 184 0706708 0000008 238 789 0182 0512267 0000007 wr1 LDEO nm nm nm 128 354 0218 0512388 0000009 orange grt1 LDEO nm nm nm 113 0607 113 0514904 0000008 orange grt 2 LDEO nm nm nm 1127 0589 116 0514956 0000007 red grt 2 LDEO nm nm nm 117 0690 102 0514536 0000021 red grt1 LDEO nm nm nm 120 0685 106 0514764 0000009 LM-43 amph SB nm 0705033 0000002 0592 1486 0241 0513478 0000006 cpx SB nm 0703971 0000010 2033 6822 0180 0513257 0000006 wr SB nm nm nm 0756 2066 0221 0513397 0000006 grt 2 SB nm nm nm 0380 0204 113 0515804 0000060 grt 3 light SB nm nm nm 0425 0226 114 051817 0000009 grt 3 dark SB nm nm nm 0356 0192 112 0515685 0000010 EH-22 cpx1 LDEO 181 0703961 0000008 0430 199 0130 0512976 0000017 cpx2 LDEO 170 0703902 0000008 0427 195 0132 0513012 0000010 dark grt LDEO nm nm nm 0169 0120 0855 0516804 0000019 dark grt2 LDEO nm nm nm nm 0102 0516717 0000017 pale grt LDEO nm nm nm 0146 0100 0888 0516192 0000015 pale grt 2 LDEO nm nm nm 0211 0149 0856 0516301 0000027 wr LDEO nm nm nm 0149 0579 0155 0513031 0000106 wr2 LDEO nm nm nm 0218 0846 0156 0513064 0000015 LEA 08-04 cpx LDEO nm nm nm 1762 7193 0148 0512632 0000008 cpx2 LDEO 171 070751 0000010 1649 6754 0148 0512601 0000007 wr LDEO nm nm nm 0623 220 0171 0512693 0000008 pink grt2 LDEO nm nm nm 0499 0364 0830 051434 0000050 pink grt LDEO nm nm nm 0545 0292 1130 051535 0000022 orange grt LDEO nm nm nm 0448 0269 1010 0516759 0000017 orange grt2 SB nm nm nm 0559 0333 1018 0516726 0000017

Analyses at Lamont Doherty Earth Observatory (LDEO) and Stony Brook (SB). nm, not measured.

Peridotite evolution not be dated directly by Sm–Nd because of the general In contrast to the apparent solely Scandian history of paucity of clinopyroxene in dunite and harzburgite. the LL garnet pyroxenites, the mixed garnet assemb- However, the Re–Os evidence presented above sug- lages within some garnet peridotites appear to require gests they originated during the Archean or two widely separated metamorphic events and a mech- Proterozoic. Refractory garnets in dunites and harzbur- anism for mixing these assemblages. There are two gites give mixed, but generally high Cr contents (11to garnet generations in all three peridotite samples which 50wt% Cr2O3, Table 4). Thus we propose the high-Cr we studied in detail, LM-42B, EH-22 and LEA 08–04. garnets in the mixed assemblages are also refractory

They all contain mixtures of Cr-rich and Cr-poor garnet garnets and we label them M1. We further propose that and the two dated samples define two age clusters with the M1 garnets in all dunites, harzburgites and lherzo- Cr-rich dark garnets much older than Cr-poor garnets lites were originally larger porphyroclasts that had (LM-42B lacked sufficient clinopyroxene to obtain a steep core to rim gradients in Cr content. Fracturing Sm–Nd mineral age). A Cr-poor garnet aliquot of one these porphyroclasts into smaller dispersed fragments peridotite (LEA 08–04) gave an age (384 6 32 Ma; Fig. would explain the large variation in Cr content between 13d) that is within error of the Scandian ages obtained individual grains (note for example variations between from garnet pyroxenites (Fig. 13a and b) and eclogites different coloured garnets in Table 4). It might also ex- (Brueckner et al., 2016). Therefore we provisionally plain the variations in Mg-number between garnet label these assemblages M2. The high-Cr garnets are grains in sample LM-42B (Fig. 9b). This model is dis- suggested to be the refractory garnets that occur within cussed further below. dunite and harzburgite (discussed above) as suggested We suggest that fragmentation and dispersal can by their older (minimum) ages in the mixed assemb- also provide a mechanism for the mixing of Cr-poor lages (Fig. 13c and d). These refractory garnets could and Cr-rich garnets observed in some samples. We Journal of Petrology, 2018, Vol. 0, No. 0 23

orange grts (a) (b) Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 0.516 red grts LM-43 grt2 Garnet Pyroxenite pale grt 0.514 dark grt

0.515 Nd

144 0.513 432 ± 5 409 Ma ±(MSWD=0.0030) 6 Ma (MSWD=1.3) Nd/ wr 0.514 143 HKB-2B 403.6 ± 7.5 Ma (MSWD = 0.56) 0.512 Garnet Websterite amph all 399 ± 18 Ma (MSWD = 4.4) 2 cpxs wr 0.513 cpx

0.511 0.0 0.4 0.8 1.2 0.0 0.4 0.8 1.2 0.517 0.518 dark grt (c) (d)

0.516 pale grts 2 orange grts 0.516 Nd pale grt1

144 0.515

Nd/ pale grt2 143 0.514 808 ± 48 Ma (MSWD=2.4) 0.514 729 ± 18 MA, MSWD = 8.4 647 ± 13 Ma (MSWD=1.1)EH-22 cpx Garnet Wehrlite 0.513 wr 384 ± 32 MA, MSWDLEA = 2.608-04 2 wrs 2 cpxs Garnet Lherzolite 0.512 0.512 0.0 0.4 0.8 0.0 0.4 0.8 1.2 147Sm/144Nd 147Sm/144Nd

Fig. 13. Garnet–clinopyroxene–whole-rock Sm–Nd isochron diagrams. (a) Pyroxenite HKB-2B and (b) LM-43 give a relatively narrow age span for garnets of different colours, suggesting that a single garnet generation grew or re-equilibrated over 20 million years during the Scandian orogeny. The much larger age difference between different coloured garnets in (c) wehrlite EH-22 and (d) LEA 08–04 indicates two separate garnet-forming events. The youngest age for a pale garnet from LEA 08–04 (d) represents a valid Scandian date, but with a large error. propose the Cr-poor garnets were mixed in from dis- could also be the result of extreme shearing; (8) thicker membered pyroxenites that also fragmented as a result garnet pyroxenite layers that have remained coherent of extreme shearing during the Scandian orogeny. (i.e. Fig. 2l) nevertheless show locally fractured garnets Support for the mixing model is as follows: (1) the two with pieces that have moved apart (Fig. 2i, j and k). garnet generations have different Cr concentrations A possible problem with this mixing model is that M2 (Fig. 9a); (2) they show a bimodal distribution in appar- clinopyroxene and orthopyroxene from dismembered ent ages (Fig. 13c and d), indicating two garnet-forming pyroxenites should also have been mixed in along with events; the lighter garnets should give younger ages garnets. We measured the compositions of five clino- than the darker garnets if they were derived from young pyroxenes from sample LM-42B (Fig. 7b) and none of (Scandian) pyroxenites; (3) many garnet grains are them have the Cr-poor compositions expected of clino- blocky, angular and irregular and show a range of sizes, pyroxenes from pyroxenites. However, neither do they suggesting they are fragments of earlier larger gains have the Cr-rich compositions of clinopyroxenes from (see further discussion below); (4) large garnets with dunite, harzburgite and another lherzolite (Fig. 7b). symmetric Cr profiles would break into smaller grains Instead they occupy an , with different Cr contents and the observed asymmetric which suggests that they might have re-equilibrated profiles (Figs 8 and 9); (5) rotation and flattening during during Scandian recrystallization. shearing and/or shortening could explain the striking Mixed assemblages with M2 clinopyroxene and Cr- parallelism of garnet pyroxenite and peridotite layers poor garnet and M1 Cr-rich garnet and olivine (no or and M1 garnet trains (Fig. 2c–f) in all peridotite lenses; few M1 clinopyroxenes) would lead to spurious mineral (6) the C-type fabric determined from LM-42B (Luc isochron ages. If the clinopyroxene is solely from pyrox-

Mehl, personal communication) is consistent with such enite (i.e. M2) the mixing of it with garnets of two differ- large strains; (7) some garnet-rich trains are remarkably ent ages (M1 and M2) would most likely result in a thin, one or two minerals wide (Fig. 2c, d and h), which minimum age for the older event. However, depending 24 Journal of Petrology, 2018, Vol. 0, No. 0

0.516 Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018

Liverpool Land eclogite 0.515 pyroxenite dunite/harzburgite lherzolite

WGC Western Gneiss Cmplx

0.514 Peridotite Trend peridotite & pyroxenite external eclogite

Nd internal eclogite 144 Nd/

143 0.513

WGC Eclogite Trend 0.512

0.511 0.700 0.705 0.710 0.715 0.720 0.725 0.730 87Sr/86Sr

Fig. 14. 143Nd/144Nd-87Sr/86Sr covariance diagram for clinopyroxenes and two amphiboles (smaller symbols) from Liverpool Land eclogites (diamonds), pyroxenites (stars), dunites (circles) and lherzolites (triangles). Values for eclogites (Brueckner, unpublished) and garnet peridotites and pyroxenites from the WGC (Brueckner et al., 2010) are plotted for comparison. The LL clinopyroxenes show evidence of crustal contamination (high 87Sr/86Sr, low 143Nd/144Nd) relative to most peridotites in the WGC. The few clinopyr- 87 86 143 144 oxenes from WGC peridotites with high Sr/ Sr and low Nd/ Nd are from M3 (Scandian) assemblages.

Table 6: Re–Os data for peridotite and pyroxenites, Tvaerdal Complex, Liverpool Land

description 187Os/188Os 6 187Re/188Os 6 CHUR RD at RE at Os 185Re Sample 16Ga 16Ga ppb ppb

CHUR 01270 0423 0 000 169 LL06 EH-21 s1A 01179 00006 0014 0003 140 135 140 280008 “ s1B 01218 00004 0149 0012 121 077 135 116013 “ s1C 01216 00008 0212 0002 168 080 164 738 012 “ S4A 01162 00021 0245 0008 400 159 257 537 010 LL06 EH-22 3A 01143 00003 0148 0007 293 187 247 390 004 “1B01092 00024 0275 0013 788 260 371 240 004 “1C01192 00033 0225 0013 259 115 205 195 003 Camp 2-1-A fractionated 01086 00004 0294 0009 944 269 387 199021 Camp 2-1-A corrected 00997 00012 0036 0012 431 394 414 746015 Camp 2-1-A “ 01083 00003 0298 0008 994 274 393 195022 Camp 2-1-H “ 01083 00002 0172 0003 469 273 343 186012 Camp 2-1-H “ 01080 00003 0161 0003 457 278 344 219013 Camp 2-1-H “ 01081 00005 0185 0003 502 276 352 148010 Camp 2-1-I “ 01079 00001 0028 0001 299 279 293 137 015 Camp 2-1-I “ 01082 00001 0038 0002 302 274 292 998013 Camp 2-1-I “ 01082 00003 0020 0004 289 274 285 511 010 Camp 2-1-I “ 01096 00003 0067 0002 304 254 283 301008 Camp 2-1-I “ 01083 00004 0030 0003 295 273 288 132 008 Camp 2-1-5 “ 01182 00005 0245 0016 327 130 228 158? Camp 2-1-15 “ 01076 00005 0314 0012 1203 283 409 210? Camp2-bl1-1 “ 01084 00020 0851 0036 254 272 603 580 ? Camp2-bl1-8 “ 01076 00009 0223 0009 618 283 374 116? Journal of Petrology, 2018, Vol. 0, No. 0 25

(a) Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018

0.22

0.18 Os

188 1.6 Ga ± 0.2

/ 0.14 Os

187 Caledonian? 0.10 HKB-03 Dunite EH-21 & 22 Mixed Lherzolite & Wehrlite Camp 2-1 Lherzolite 0.06 012345 187Re/188Os 0.14 (b) Liverpool Land Mantle Evolution EH-21 EH-22 Camp 2-1A 0.13 Camp 2-1H Camp 2-1I CHUR Camp 2-1-5 Camp 2-bl1-4 Camp 2-1-15 Camp 2-bl1-8 Os 0.12 188

Os/ 5 0.11 4 187 I I H H 8 A 0.10 15 A

0.09 Re addition at 1.6 Ga

0.08 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 Age (Ga)

Fig. 15. (a) Re–Os isochron diagram for LA-ICP-MS analysed sulfides from Liverpool Land peridotites and pyroxenites. Symbols with dashed borders were excluded from calculating the 16 Ga errorchron. (b) 187Os/188Os evolution diagram for sulfides with low 187Re/188Os ratios and present day 187Os/188Os ratios below CHUR, assuming Re enrichment at 16 Ga. The letters and numbers at the intercepts of Camp 2–1 sulfides with the CHUR line refer to separate sulfides within the dunite, which otherwise would give ages older than the Earth. Repeated letters refer to maximum and minimum ages from individual sulfides (see Supplementary Data Table S4).

on the relative proportions of M1 and M2 garnet, the ap- orthopyroxene with a M1 garnet or vice versa (this prob- parent ages could be older, younger and even, by ser- lem was largely avoided by analysing cores and rims of endipity, the same age as the earlier garnet-forming adjacent grains in several domains in the rock sample). event. Therefore the two sets of Sm–Nd ages defined The mechanical mixing of two garnet-bearing by minerals from EH-22 and LEA 08–04 (Fig. 13c and d) assemblages of different age resolves many of the yield no precise information on their time of formation, problems that arise if the garnets had formed during a except for the Scandian age of the pale garnet in LEA single metamorphic cycle. However, certain conditions

08–04, which appears to be an aliquot of M2 garnets al- would have had to be in place for which there is little most completely separated from M1 garnets. Similar direct evidence. For example, postulating the existence arguments suggest that the P–T estimates from LM-42B of large parental M1 grains would explain the large could be invalid. For example, the Al-in-orthopyroxene inter-grain variation in the Cr content (and Mg-number) geobarometer could potentially pair a M2 of garnets in sample LM-42B (Fig. 9). It would also 26 Journal of Petrology, 2018, Vol. 0, No. 0 explain why the Cr-rich (darker) garnets in EH-22 and subsequent re-equilibration and, therefore, tend to pre- LEA 08–04 give older ages than the Cr-poor (i.e. lighter) serve older ages. The more completely equilibrated Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 garnets. Cr content increases from core to rim in most areas near the rim would give younger ages. But the other measured samples (Fig. 8) and so paler cores garnets from the parental cores would be at least par- should give older ages than darker rims, but the oppos- tially re-equilibrated, resulting in minimum M1 ages, ite is true. The solution is two metamorphic events, with and the garnets from the rims would not equilibrate both events producing core to rim gradients in Cr con- completely, resulting in maximum M2 ages. Dating gar- centrations. The first event (M1) generated negative net fragments produced by this mechanism would give core to rim Cr gradients in garnet as a result of decreas- older ages for the cores (dark garnets) and younger ing P–T conditions, perhaps as a result of decompres- ages for the rims (pale garnets), which is what is sion within a mantle plume. The core to rim decrease in observed. pyrope and increase in almandine content of the Cr-rich The mixing of M1 and M2 assemblages appears to garnet in LM-42B (GM 1 in Fig. 9b) is consistent with have generated some peridotites with whole-rock com- this hypothesis. The result would be a large difference positions within or near the lherzolite field (Table 3, in Cr concentration between the core and rim, but with- Fig. 3). Lherzolites are classically interpreted as samples in large garnets, resulting in Cr gradients with relatively of primitive mantle (Bodinier & Godard, 2004). More re- shallow slopes (the larger the grain, the shallower the cently some lherzolites have been shown to be the re- slope). The scan of the Cr-rich garnet in LM-42B (GM 1 sult of (Su et al., 2016), even in the in Fig. 9a) might be from the core of such a garnet as it classic type locality (i.e. Lherz; La Roux et al., 2007). The is the only measured garnet that displays a convincing mechanical mixing model offers yet another hypothesis symmetrical core to rim decrease in Cr content. to explain the compositions of LL lherzolites, olivine Fragmenting these garnets would produce smaller websterites, and . grains, each displaying a fraction of the original profile (Fig. 9b). Thus different grains would have different Cr concentrations (and colours) depending on where they A COMPARISON OF TVAERDAL PERIDOTITES originally occurred in the parent grain (i.e. dark frag- WITH WGC PERIDOTITES ments from the core, light fragments from the rim). If the Tvaerdal Complex is a stranded Baltic terrain, These fragments would also inherit the shallow Cr- identical in age, tectonic position and history to the profiles of the parental grain and, therefore, not show WGC of Norway, it follows that the peridotites within much internal colour change, while still showing visible them might have had a common origin and history. colour differences between grains. There are similarities, but also significant differences.

The second required recrystallization event (M2), pre- Re–Os data suggest both ultramafic assemblages origi- sumably Scandian metamorphism, imposed new zon- nated from a mixed Archean/Proterozoic mantle (Beyer ing patterns for Cr at the same time the parent garnets et al., 2004), consistent with a Laurentian mantle wedge broke into smaller fragments as a result of shearing. above subducted Baltica (Brueckner et al., 2010, Fig. 7; Each fragment would have preserved the overall Cr Beyer et al., 2012). The WGC garnet assemblages that content of its position in the original porphyroclast, but are most similar to those in LL may be the spectacular with a new imposed M2 gradient. Most grains show a ‘relict-type’ garnet peridotites and garnet pyroxenites core to rim increase in Cr concentration, so the new gra- (Brueckner & Medaris, 2000) of the northwestern WGC dient probably formed during prograde metamorphism (Spengler et al., 2006; 2009; Van Roermund, 2009a,b; (i.e. Equation 1). Further fragmentation would break Brueckner et al., 2010; Scambelluri et al., 2010). These these grains into still smaller grains, many with incom- peridotites and their host gneisses were subducted to plete (asymmetric) profiles. Stated succinctly, grain- the deepest level of the mantle during the terminal colli- to-grain variations in Cr content reflect gradients sion of Baltica and Laurentia where they generated generated in the precursor porphyroclast during M1 Scandian-age garnets (correlated here with the M2 gar- metamorphism, but intra-grain variations reflect gra- nets in the Tvaerdal pyroxenites). Additionally, the dients in fragments of these porphyroclasts during M2 WGC ultramafics also contain older garnet assemb- (Scandian) metamorphism. The large parental porphyr- lages; M1 and M2 garnets generated in the mantle dur- oclasts hypothesized in this model (>10 cm) have not ing the Archean and Proterozoic, respectively. The WGC been observed in the peridotites of the Tvaerdal M1 garnets formed initially as very large porphyro- Complex. The largest M1 garnet measured is 2 cm wide. blasts, as has been proposed above for the M1 Cr-rich However, such porphyroclasts >10 cm wide have been Tvaerdal garnets. However the maximum errorchron described in the peridotites of the northwestern WGC age generated by M1 garnets in LL peridotites is only (Spengler, et al., 2009; Van Roermund, 2009a,b). 808 Ma (Fig. 13), much younger than the oldest ages

The trivalent REE elements presumably developed (16 Ga) defined by M2 garnets in WGC peridotites and core to rim patterns similar to those of Crþþþ. These pyroxenites (Brueckner et al., 2010). Sm–Nd clinopyrox- gradients could explain the variations in the obtained ene isotope variations and Re–Os sulfide model ages Sm–Nd mineral ages (Fig. 13). REE in the cores of the from the Tvaerdal peridotites (Figs 14 and 15; Tables 5 parental porphyroclasts would be least affected by and 6) suggest they experienced a mantle fertilization Journal of Petrology, 2018, Vol. 0, No. 0 27 event, similar to that which affected the peridotites of shared history. Recrystallization under peak HP/UHP the WGC (Beyer et al., 2004, 2006). However the conditions resulted in the nearly complete equilibration Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018

Tvaerdal peridotites lack convincing evidence that this of M2 garnet-bearing assemblages in the pyroxenites: event occurred in the mantle at 1–6 Ga as was the only Cr and Ca profiles may record a prograde history. case in the WGC. A final conclusion on the relationships The exhumation of the Tvaerdal Complex brought it between these ultramafic occurrences requires consid- against the lower crustal Jaettedal Complex, resulting erably more research, but the overall impression so far in a transfer from Baltica to Laurentia. Subsequently is each association has its own unique petrogenetic his- both complexes moved towards the surface along the tory, suggesting either the (Laurentian?) mantle wedge Grubbedalen Shear Zone. from where they were derived was strikingly heteroge- Mechanical peridotite/pyroxenite mixing may be a neous or there was more than one mantle wedge be- more important process than generally realized. A crit- neath potentially different overriding terranes. ical question is where such mixing occurred: within the mantle or later, after introduction into the subducted crust? If in the mantle, the mechanical mixing model CONCLUSIONS has consequences for understanding the geochemical Six garnet peridotite lenses occur within the gneisses of evolution of the mantle where most variations in mantle the HP/UHP Tvaerdal Complex in Liverpool Land, composition are assumed to be the result of melting Greenland Caledonides. The peridotites are predomin- and, or, chemical processes. antly dunite and harzburgite, but include layers of strik- ingly parallel lherzolite, websterite and wehrlite. Several samples include two garnet populations, one ACKNOWLEDGEMENTS older and Cr-rich (M1) and the other younger (400 Ma) and Cr-poor (M2). We propose the two garnet popula- Hannes Brueckner and Scott Johnston are grateful to tions in these samples were mechanically mixed Ebbe Hartz, who invited us to join field excursions in LL through intense shearing during the Scandian orogeny. in the summers of 2006 and 2007 with equipment pro- Evidence for this intense shearing includes: garnet vided by him and expenses covered by Statoil ASA. We grains that are fractured, separated and locally disag- thank Luc Mehl, Kenjo Agustsson and Richard Bubbico gregated; garnets that have asymmetric edge-to-edge for assistance with fieldwork and Aka Lynge for assist- gradients in Cr and Ca content; garnet grains that are ance with field logistics. Luc Mehl also measured LPO strung out in thin layers; striking variations in Cr content olivine and orthopyroxene fabrics and Richard Bubbico in neighboring garnet grains; and the parallelism of all assisted in sample preparation and analysis. We are pyroxenite, peridotite and garnet-rich layers. The M2 also very grateful to many people who supported this garnets were probably derived from disaggregated pyr- research and helped us acquire and manipulate the oxenite layers. The source of the M1 garnets is less cer- data presented herein; specifically Suzanne O’Reilly tain, but we suggest they are disaggregated fragments and Norman Pearson (MU), Juliane Gross and Kennet of refractory garnet porphyroclasts that were left behind Flores (AMNH), and Terry Plank and Louise Bolge after Archean and/or Proterozoic melting events in the (LDEO). Finally we thank Shuguang Song, Iwona mantle. The intense shearing that fragmented garnet Klonowska and an unnamed reviewer for helpful com- porphyroclasts and thin pyroxenite veins, and locally ments and suggestions. mixed their assemblages, simultaneously thinned the This is contribution 8242 of the Lamont-Doherty Earth layers and rotated them into parallelism. Observatory of Columbia University, contribution 1196 Clinopyroxenes within these mixed rocks, unlike the of the ARC Centre of Excellence for Core to Crust Fluid garnets, do not plot as two populations. Either most M1 Systems (http://www.ccfs.mq.edu.au) and contribution clinopyroxene was removed during Precambrian melt- 1247 of the GEMOC Key Centre (http://www.gemoc.mq. ing or mixed M1 and M2 clinopyroxenes re-equilibrated edu.au) and is related to IGCP-662. with each other to form clinopyroxenes with intermedi- ate compositions. Garnets and other minerals within thick (>5 cm) pyr- oxenite layers appear to have remained coherent and FUNDING contain only M2 assemblages. P–T–t results from these This research was supported by The City University of pyroxenites, combined with earlier studies of LL eclo- New York Faculty Research Award Grants 66287–00 35, gites within the adjacent gneisses, are consistent with 67663–00 36, and 68556–00 37 to HKB and by National the subduction of the Tvaerdal Complex around 400 Ma Science Foundation EAR 08–38530 to HKB and SJ. ago to mantle depths of 110–120 km and temperatures of 850–900C. The similar Scandian history of eclogite and garnet pyroxenite, which had very different origins, indicates that the peridotite lenses were inserted from SUPPLEMENTARY DATA the mantle wedge into the Tvaerdal Complex relatively Supplementary data are available at Journal of early, during subduction, resulting in the subsequent Petrology online. 28 Journal of Petrology, 2018, Vol. 0, No. 0

REFERENCES Carswell, D. A., Harvey, M. A. & Al-Samman, A. (1983). The petrogenesis of contrasting Fe–Ti and Mg–Cr garnet perido- Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 Al-Samman, A. H. (1985). Mineralogy and of tite types in the high-grade gneiss complex of western ultramafic rocks in the west Norway basement gneiss ter- Norway. Bulletinin de Mine´ ralogie 106, 727–750. rain. Ph.D. thesis, University of Sheffield, Sheffield. Carswell, D. A., Van Roermund, H. L. M. & Wiggers de Vries, D. Andersen, T. B., Jamtveit, B., Dewey, J. F. & Swensson, E. F. (2006). Scandian ultrahigh-pressure metamorphism of (1991). Subduction and eduction of continental crust: major Proterozoic basement rocks on Fjørtof and Otrøy, Western mechanisms during continent-continent collision and oro- Gneiss Region, Norway. International Geology Review 48, genic extensional collapse, a model based on the southern 957–977. Norwegian Caledonides. Terra Nova 3, 303–310. Cheeney, R. F. (1985). The plutonic igneous and high-grade Augland, L. E., Andresen, A. & Corfu, F. (2010). Age, structural metamorphic rocks of southern Liverpool Land, central East setting, and exhumation of the Liverpool Land eclogite ter- Greenland, part of a supposed Caledonian and Precambrian rane, East Greenland Caledonides. 2, 267–286. complex. Grønlands Geologiske Undersøgelse, Rapport Augland, L. E., Andresen, A. & Corfu, F. (2011). Terrane transfer 123, 1–39. during the Caledonian orogeny: Baltican affinities of the Coe, K. (1975). The Hurry Inlet granite and related rocks of Liverpool Land Eclogite Terrane in East Greenland. Journal Liverpool Land, East Greenland. Bulletin Grønlands of the Geological Society 168, 15–26. Geologiske Undersøgelse 115, 1–34. Augland, L. E., Andresen, A., Corfu, F. & Daviknes, H. K. (2012). Corfu, F. & Hartz, E. H. (2011). U–Pb geochronology in Liverpool Late Ordovician to Early Silurian ensialic magmatism in Land and Canning Land, East Greenland—the complex re- Liverpool Land, East Greenland: new evidence extending cord of a polyphase Caledonian orogeny. Canadian Journal the northeastern branch of the continental Laurentian mag- of Earth Sciences 48, 473–494. matic arc. Geological Magazine 149, 561–577. Eskola, P. (1921). On the eclogites of Norway. Videnskapsselskapets Beyer, E. E., Brueckner, H. K., Griffin, W. L. & O’Reilly, S. Y. Skrifter. I. Mathematisk–naturvidenskabelig Klasse 8, 1–121. (2012). Laurentian provenance of Archean mantle fragments Gasser, D. (2014). The Caledonides of Greenland, Svalbard and in Proterozoic Baltic crust of the Norwegian Caledonides. other arctic areas: status of research and open questions. In: Journal of Petrology 53, 1357–1383. Corfu, F., Gasser, D. & Chew, D. M. (eds) New Perspectives Beyer, E. E., Brueckner, H. K., Griffin, W. L., O’Reilly, S. Y. & on the Caledonides of Scandinavia and Related Areas, Graham, S. (2004). Archean mantle fragments in Proterozoic Volume 390. London: Geological Society Special Publications, crust, Western Gneiss Region, Norway. Geology 32, 609–612. pp. 93–129. Beyer, E. E., Griffin, W. L. & O’Reilly, S. Y. (2006). Gee, D. G., Fossen, H., Henriksen, N. & Higgins, A. K. (2008). Transformation of Archaean lithospheric mantle by refertili- From the early Paleozoic platforms of Baltica and Laurentia zation: evidence from exposed peridotites in the Western to the Caledonide Orogen of Scandinavia and Greenland. Gneiss Region, Norway. Journal of Petrology 47, 1611–1636. Episodes 31, 44–51. Bodinier, J.-L. & Godard, M. (2004). Orogenic, ophiolitic, and Gee, D. G., Jana´ k, M., Majka, J., Robinson, P. & van Roermund, abyssal peridotites. In: Holland, H. D. and Turekian, K. K. H. (2012). Subduction along and within the Baltoscandian (eds) Treatise on Geochemistry. Amsterdam: Elsevier, pp. margin during closing of the Iapetus Ocean and 103–170. Baltica-Laurentia collision. Lithosphere 5, 169–178. Brueckner, H. K., Carswell, D. A., Griffin, W. L., Medaris, L. G., Hacker, B. R., Andersen, T. B., Johnston, S., Kylander-Clark, A. Jr, Van Roermund, H. L. M. & Cuthbert, S. J. (2010). The R. C., Peterman, E. M., Walsh, E. O. & Young, D. (2010). mantle and crustal evolution of two garnet peridotite suites High-temperature deformation during continental-margin from the Western Gneiss Region, Norwegian Caledonides: subduction and exhumation: the ultrahigh-pressure an isotopic investigation. Lithos 117, 1–19. Western Gneiss Region of Norway. Tectonophysics 480, Brueckner, H. K. & Medaris, L. G., Jr (2000). A General Model 149–171. for the Intrusion and Evolution of “Mantle” Peridotites in Harley, S. L. (1984). An experimental study of the partitioning of High-Pressure and Ultrahigh-Pressure Metamorphic Fe and Mg between garnet and orthopyroxene. Terranes. Journal of Metamorphic Geology 18, 119–130. Contributions to Mineralogy and Petrology 86, 359–373. Brueckner, H. K., Medaris, L. G., Jr, Belousova, E. A., Johnston, Hartz, E. H., Condon, D., Austrheim, H. & Erambert, M. (2005). S. M., Griffin, W. L., Hartz, E. H., Hemming, S., Ghent, E. & Rediscovery of the Liverpool Land Eclogites (Central East Bubbico, R. (2016). An orphaned baltic terrane in the Greenland): a post supra-subduction UHP province. Abstract, Greenland Caledonides: a Sm-Nd and detrital zircon study Mitteilungen der O¨ sterreichischen Mineralogischen of a high-pressure/ultrahigh-pressure complex in Liverpool Gesellschaft 150. Land. The Journal of Geology 124, 541–567. Hartz, E. H., Podladchikov, Y. Y. & Dabrowski, D. (2007). Brueckner, H. K. & Van Roermund, H. L. M. (2004). Dunk Tectonic and reaction overpressures: theoretical models Tectonics: a multiple subduction/eduction model for the evo- and natural examples. Abstract, European Geophysical lution of the Scandianvian Caledonides. Tectonics 23, 1–20. Union General Assembly, EGU2007-A-10430. Carlson, W. D. (2012). Rates and mechanism of Y, REE, and Cr Higgins, A. K., Gilotti, J. A. & Smith, M. P. (eds) (2008). The diffusion in garnet. American Mineralogist 97, 1598–1618. Greenland Caledonides: Evolution of the Northeast Margin Carswell, D. A. (1966). Some ultrabasic bodies and related rocks of Laurentia, Volume 202. Geological Society of America in Sunnmo¨ re, South Norway. Ph.D. thesis, University of Memoir, pp. 1–368. doi:10.1130/2008.1202. Sheffield, Sheffield. Johnston, S. M., Hartz, E. H., Brueckner, H. K. & Gehrels, G. E. CarswelI, D. A. (1968a). Picritic -residual dunite relation- (2010). U–Pb zircon geochronology and tectonostratigraphy ships in garnet peridotite at Kalskaret near Tafjord, southern of southern Liverpool Land, East Greenland: implications for Norway. Contributions to Mineralogy and Petrology 19, 97–124. deformation in the overriding plates of continental colli- Carswell, D. A. (1968b). Possible primary upper mantle perido- sions. Earth and Planetary Science Letters 297, 512–524. tite in Norwegian basal gneiss. Lithos 1, 322–355. Johnston, S. M., Kylander-Clark, A. & Brueckner, H. K. (2015). Carswell, D. A. (1973). Garnet pyroxenite lens within Ugelvik Migmatite gneiss of the Jættedal Complex, Liverpool Land, layered garnet peridotite. Earth and Planetary Science East Greenland: protracted high-T metamorphism in the Letters 20, 347–352. Journal of Petrology, 2018, Vol. 0, No. 0 29

overriding plate of the Caledonian orogen. Journal of Scambelluri, M., Van Roermund, H. L. M. & Pettke, T. (2010).

Metamorphic Geology 1, 22. Mantle wedge peridotites: fossil reservoirs of deep subduc- Downloaded from https://academic.oup.com/petrology/advance-article-abstract/doi/10.1093/petrology/egy094/5131633 by University of Wisconsin-Madison user on 04 December 2018 Krogh, E. J. (1977). Evidence for a Precambrian continent- tion zone processes Inferences from high and continent collision in western Norway. Nature 267, 17–19. ultrahigh-pressure rocks from Bardane (Western Norway) Lappin, M. A. (1974). Eclogites from the Sunndal-Grubse ultra- and Ulten (Italian Alps). Lithos 120, 186–201. mafic mass, Almklovdalen, Norway and the T-P history of Smith, D. C. & Cheeney, R. F. (1981). A new occurrence of the Almklovdalen masses. Journal of Petrology 15, 567–601. garnet-ultrabasite in the Caledonides; a Cr-rich –- La Roux, V., Bodinier, J.-L., Tommasi, A., Alard, O., Dautria, J.- garnet-lherzolite from Tvaerdalen, Liverpool Land, East M., Vauchez, A. & Riches, A. J. V. (2007). The Lherz spinel Greenland. Terra Cognita 1, 74. lherzolite: refertilized rather than pristine mantle. Earth and Spengler, D., Brueckner, H. K., van Roermund, H. L. M., Drury, Planetary Science Letters 259, 599–612. M. R. & Mason, P. R. D. (2009). Long-lived, cold burial of McDonough, W. F. & Sun, S.-S. (1995). The composition of the Baltica towards 200 km depth. Earth and Planetary Science Earth. Chemical Geology 120, 223–253. Letters 281, 27–35. Majka, J., Rose´ n, A., Jana´ k, M., Froitzheim, N., Klonowska, I., Spengler, D., van Roermund, H. L. M., Drury, M. R., Ottolini, L., Manecki, M., Sasinkova´ , V. & Yoshida, K. (2014). Mason, R. D. & Davies, G. R. (2006). Deep origin and hot Microdiamond discovered in the Seve Nappe (Scandinavian melting of an Achaean orogenic peridotite massif in Caledonides) and its exhumation by the “vacuum-cleaner” Norway. Nature 440, 913–917. mechanism. Geology 42, 1107–1110. Su, B., Chen, Y., Guo, S. & Liu, J. (2016). Origins of orogenic Mercy, E. L. P. & O’Hara, M. J. (1965). Chemistry of some dunites: petrology, geochemistry and implications. garnet-bearing rocks from south Norwegian peridotites. Gondwana Research 29, 41–59. Norsk Geologisk Tiddskrift 445, 323–332. Terry, M. P., Robinson, P., Brueckner, H. K. (1999). Sinking intru- Moore, A. C. & Kvale, H. (1977). Three varieties of alpine-type sion model for the emplacement of garnet-bearing perido- ultramafic rocks in the Norwegian Caledonides and Basal tites into orogens: comment and reply. Gneiss Complex. Lithos 10, 149–161. Geology 27, 447–448. Nickel, K. G. & Green, D. H. (1985). Empirical geothermobarom- Van Roermund, H. L. M. (2009a). Recent progress in Scandian etry for garnet peridotites and implications for the nature of UHPM in the northernmost domain of the Western Gneiss the lithosphere, and diamonds. Earth and Complex, SW Norway: continental subduction down to Planetary Science Letters 73, 158–170. 180–200 km. Journal of the Geological Society of London Obata, M. (2011). Kelyphite and symplectite: textural and mineral- 166, 1–13. ogical diversities and universality, and a new dynamic view of Van Roermund, H. L. M. (2009b). Mantle wedge garnet perido- their structural formation. In: Sharkov, E. V. (ed.) New Frontiers tites from the northernmost ultra-high pressure domain of in Tectonic Research, London: IntechOpen, pp. 93–122. the Western Gneiss Region, SW Norway. European Journal O’Neill, H. Ste. (1981). The transition between spinel lherzolite of Mineralogy 21, 1085–1096. and garnet lherzolite, and its use as a geobarometer. Wang, Q., Xia, Q.-K., O’Reilly, S. Y., Griffin, W. L., Beyer, E. E. & Contributions to Mineralogy and Petrology 77, 185–194. Brueckner, H. K. (2013). Pressure- and stress-induced fabric Osland, R. (1997). Modeling of variations in Norwegian olivine transition in olivine from peridotites in the Western Gneiss deposits, causes of variation and estimation of key quality Region (Norway): implications for mantle seismic anisot- factors. Doktor Ingeniør thesis, Norwegian University of ropy. Journal of Metamorphic Geology 31, 93–111. Science and Technology, pp. 189. Ph.D. thesis, Norges Wu, C. M., & Zhao, G. C. (2007). A recalibration of the Tekniske Høgskole, pp. 1–189. garnet-olivine geothermometer and a new geobarometer Sahlstein, T. G. (1935). Petrographie der Eklogiteinschliisse in for garnet peridotites and garnet-olivine-plagioclase-bear- den Gneisen des sudwestlichen Liverpool-Landes in ing granulites. Journal of Metamorphic Geology 25, Ost-Gro¨ nland. Meddelelse Gro¨ nland 95, 1–43. 497–505.