Multi-method U-Pb baddeleyite dating: insights from the Spread Eagle Intrusive Complex and Cape St. Mary’s sills, Newfoundland, Canada.

. Baddeleyite (ZrO 2 ) is widely used in U-Pb geochronology, but analysis and age interpretation are often difficult, especially for samples which have experienced post-intrusive alteration and/or metamorphism. Here, we combine high spatial resolution (secondary ionization mass spectrometry, SIMS) and high precision (isotope dilution thermal ionization mass spectrometry, ID-TIMS) analyses of baddeleyite from the Spread Eagle Intrusive Complex (SEIC) and Cape St. Mary’s sills (CSMS) from Newfoundland. Literature data and our own detailed microtextural analysis suggest that at least seven different 20 types of baddeleyite–zircon intergrowths can be distinguished in nature. These include secondary baddeleyite inclusions in altered zircon, previously unreported from low-grade rocks, and likely the first discovery of xenocrystic zircon inclusions mantled by baddeleyite


Introduction
Baddeleyite, a monoclinic ZrO2 polymorph, is one of the most commonly used minerals in U-Pb geochronology, especially for mafic rocks, which are traditionally difficult to date (e.g., Schaltegger and Davies, 2017).It forms most readily during the late stage of igneous crystallization from a silica-undersaturated residual melt, and can co-exist with zircon at conditions near silica saturation (Heaman and LeCheminant, 1993;Schaltegger and Davies, 2017).Where both minerals co-exist, geochronologists have tended to prefer baddeleyite dates because it is (1) a primary igneous mineral, facilitating age interpretation, (2) rarely inherited from country rock, and (3) more resistant to Pb loss than zircon (e.g., Heaman and LeCheminant, 1993), as it remains crystalline even at high radiation doses (Lumpkin 1999).U-Pb baddeleyite dating has proved powerful in solving numerous problems in earth and planetary sciences (e.g., Olsson et al., 2011;Moser et al., 2013;Wall and Scoates, 2016;Davies et al., 2017;White et al., 2020).
However, many of these studies were able to work with rather large, texturally simple, unaltered and concordant baddeleyite crystals.Such favorable conditions are the exception rather than the rule for large parts of the geologic record.In fact, crystals are often too small for mineral separation (Söderlund and Johansson, 2002), prohibiting high-precision ID-TIMS (isotope dilution-thermal ionization mass spectrometry) analysis.Instead, they can be analyzed in situ by secondary ionization mass spectrometry (SIMS; e.g., Schmitt et al., 2010;Chamberlain et al., 2010) or laser ablation ICP-MS (e.g., Renna et al., 2011).
In situ methods requires for relative sensitivity corrections, which may be complicated by crystal orientation effects (Wingate and Compston, 2000), although these effects for SIMS can be reduced by oxygen flooding (Schmitt et al., 2010;Li et al., 2010).Furthermore, baddeleyite can be intergrown with other Zr-bearing minerals, especially in rocks with a metamorphic or hydrothermal overprint, where fluids with high SiO2 activity often cause partial reaction of baddeleyite to polycrystalline zircon (Heaman and LeCheminant, 1993;Söderlund et al., 2013).In situ dating after careful micropetrographic investigation can resolve such intergrowths, and for ID-TIMS, baddeleyite can be selectively dissolved in hydrochloric acid, leaving zircon rims essentially undissolved (Rioux et al., 2010).However, even baddeleyite that lacks zircon intergrowths is often discordant, and chemical abrasion, which is now routinely used to produce more concordant zircon data, is ineffective for reducing baddeleyite discordance (Rioux et al., 2010).Mechanisms such as alpha recoil (Davis and Sutcliffe, 1985;Davis and Davis, 2018), Pb loss due to fast pathway diffusion (Rioux et al., 2010;Söderlund et al., 2013;Schaltegger and Davies, 2017), and isotopic disequilibrium due to 231 Pa excess (Amelin and Zaitsev, 2002;Crowley and Schmitz, 2009) or 222 Rn loss (Heaman and LeCheminant, 2000) have been proposed to explain baddeleyite discordance, but none of these is universally accepted as the dominant baddeleyite discordance mechanism.
It is necessary to investigate to which degree reliable baddeleyite dating is possible when several of the above-mentioned challenges come together.Here, we present a case study of dating the Spread Eagle Intrusive Complex (SEIC) and Cape St.
Mary's sills (CSMS) of the Avalon Zone of Newfoundland.These early Paleozoic mafic dikes and sills were affected by lowgrade metamorphism and contain abundant texturally complex baddeleyite, as commonly found in the geologic record.
Moreover, the approximate intrusion age of the SEIC is constrained by its stratigraphic context.We applied U-Pb geochronology by SIMS and ID-TIMS on the same baddeleyite crystals, combined with detailed micro-petrographic characterization of baddeleyite by scanning electron microscopy (SEM) before and after SIMS analysis.Our comparison of SIMS and ID-TIMS dating also includes essentially unaltered baddeleyite from the Duluth gabbro (sample FC-4b, Schmitt et al., 2010) and Freetown Layered Complex (sample SL18, Callegaro et al., 2017), leading to a critical evaluation of possibilities and limitations in dating small, texturally complex and/or altered baddeleyite crystals.Based on our micro-petrographic and geochronologic data, we discuss various types of baddeleyite-zircon intergrowths, possible mechanisms of baddeleyite discordance and the reliability of different types of baddeleyite dates (e.g., 206 Pb/ 238 U, 207 Pb/ 206 Pb, concordia upper intercept).
Many of these implications are also significant for unaltered baddeleyite.

Regional geology
The Avalon Peninsula consists of rocks that were formed as part of the microcontinent Avalonia during the Neoproterozoic and early Paleozoic (e.g., Williams, 1979;Murphy et al., 1999; Figure 1).The Cambrian Adeyton and Harcourt groups (Hutchinson, 1962;King, 1988; Figure 2), which unconformably overlie Precambrian rocks, consist of well-preserved marine sediments with intercalated pillow basalts and mafic tuffs.The feeder pipes or dike-like conduits of these volcanic rocks make up a mafic intrusive complex, called the "Spread Eagle Gabbro" (McCartney, 1967) or "Spread Eagle Gabbro and equivalents" (King, 1988).To avoid confusion, we define this unit, consisting of at least eleven dikes, as the Spread Eagle Intrusive Complex (SEIC).Lower Ordovician sedimentary rocks are exposed only on Conception Bay islands northeast of the study area (King, 1988).They are roughly coeval with Avalonia's presumably early Ordovician separation from Gondwana (e.g., Cocks et al., 1997;Murphy et al., 2006;Linnemann et al., 2008Linnemann et al., , 2012;;Pollock et al., 2009Pollock et al., , 2012)).The Cambrian pillow basalts and their feeder dikes are confined to the western Avalon Peninsula, which experienced deformation and pervasive low-grade metamorphism during the Acadian Orogeny (McCartney, 1967), lasting from ca. 420-360 Ma (van Staal and Barr, 2012;Willner et al., 2014).

Spread Eagle Intrusive Complex (SEIC) and Cambrian volcanic rocks
Cambrian shales of the Chamberlain's Brook Formation and Manuels River Formation contain pillow basalt flows and/or basaltic metatuffs in five different localities (McCartney, 1967;Greenough, 1984;Greenough and Papezik, 1985b;Fletcher, 2006).Additionally, our field observations suggest their presence also within the overlying Elliot Cove formation (sensu stricto King, 1988).Several lines of indirect field evidence substantiate that these volcanic rocks were fed by the SEIC dikes or pipes which crop out in their vicinity (Greenough, 1984;Greenough and Papezik, 1985b).The SEIC intrusions are arranged in a N-S trending array (Figure 1), and are usually subcircular pipes with diameters of several hundreds of meters (McCartney, 1967).
Petrologically, they are internally differentiated, ranging from melanogabbros to leucogabbros (Greenough, 1984;Greenough and Papezik, 1985b) and monzonites (this study).These subvolcanic rocks are considerably altered to a greenschist facies assemblage, but better preserved than the volcanic rocks, which almost completely lack primary minerals (Greenough and Papezik, 1985a,b).Geochemical features of the SEIC rocks are similar to rift basalts, suggesting that Avalonia experienced extensional tectonics in the Cambrian (Greenough and Papezik, 1985b;Greenough and Papezik, 1986a; new whole-rock major and trace element data are given in the supplement S7).Early attempts at radiometric dating of the SEIC failed, as alteration hampered Rb-Sr whole-rock methods and zircons were not found (Greenough, 1984).Samples of the current study were collected from several SEIC feeder pipes (Table 1).

Cape St. Mary's sills (CSMS)
In addition to the occurrence of the Cambrian igneous rocks described above, the Cambrian succession is intruded by the CSMS in the southwestern Avalon Peninsula, especially in the upper Cambrian Gull Cove Formation (Fletcher, 2006).The sills are up to 60 m thick and consist mostly of gabbro, but the thickest sills also include up to 3 m thick granophyric dikes and pockets (Greenough and Papezik, 1986b;Fletcher, 2006).The gabbros and granophyres are both unusually rich in hydrous minerals, notably amphibole and biotite, and volatile complexing was probably an important differentiation process of the subalkaline parental magma (Greenough and Papezik, 1986b).Trace element patterns indicate that the magma was generated from a mantle source that had been metasomatically enriched by subduction zone-derived fluids (Greenough et al., 1993;this study, supplement S7).A CSMS granophyre was the first terrestrial rock for which baddeleyite dating was performed (Greenough, 1984 and pers. comm.).These ID-TIMS analyses of multiple-crystal aliquots yielded an early Silurian weighted mean 207 Pb/ 206 Pb date of 441 ± 2 Ma, but all analyses are discordant to 2.0-3.5% (Greenough et al., 1993).Combining this date with a paleolatitude of 32° S ± 8°, the CSMS were interpreted as the result of an igneous event after Avalonia's separation from Gondwana, but before complete closure of the Iapetus (Hodych and Buchan, 1998).Samples of the current study were collected from a 60 m thick sill at the southwestern coast of Lance Cove (Table 1) from the same sill as in Greenough et al. (1993).

U-Pb geochronology methods
Care was taken during field work to collect samples as coarse-grained and unaltered as possible.Such samples are most likely to yield large, high-quality baddeleyite crystals, as baddeleyite typically forms in interstitial melt pockets and can react to zircon in response to Si-bearing metamorphic fluids (Heaman and LeCheminant, 1993).Polished thin sections of selected samples were investigated by SEM (backscattered electron (BSE) and cathodoluminescence (CL) imaging and energy dispersive X-ray spectrometry (EDS)) to localize baddeleyite and study it textural properties.When baddeleyite and/or zircon crystals were too small for mineral separation (<50 µm), they were analyzed in situ in thin section by SIMS for U-Pb dates (Table 1).For samples with larger crystals, polished epoxy grain mounts were prepared from handpicked separates, followed by SEM imaging and SIMS analysis.Every SIMS session was followed by re-imaging of the analysis craters by SEM to identify analyses with contributions from adjacent phases (see supplement S1).Selected crystals were removed from grain mounts for single crystal ID-TIMS analyses.SIMS analyses were performed using a CAMECA ims1280-HR ion probe at Heidelberg University.Oxygen flooding of the sample chamber was employed to mitigate crystal orientation effects and improve secondary ion yields (Schmitt et al., 2010;Chamberlain et al., 2010;Li et al., 2010), using oxygen pressures of 2.0-3.0×10 - Pa.The primary ion beam was focused to a diameter of about 10-15 µm.In cases where an even higher spatial resolution was needed, the field aperture of the secondary beam was closed to a square of 5-8 µm.Analytical procedures for baddeleyite and zircon were adapted from Schmitt et al. (2010).The U/Pb relative sensitivity calibration (RSC) against the UO2/U ratio accounts for differences in Pb ionization caused by spot-to-spot differences in sputtering behavior.For this, Phalaborwa baddeleyite (Heaman, 2009) was always used as primary reference material.For grain mount sessions, FC-4b baddeleyite (Schmitt et al., 2010) was included as secondary reference material.Zircon analyses were calibrated using the reference materials AS3 (Schmitz et al., 2003) for U-Pb ages and 91500 (Wiedenbeck et al., 2004) for U concentrations.Zirconolite from sample FP7G was analyzed to investigate the influence of its matrix on U-Pb baddeleyite dates when the primary ion beam overlaps onto both minerals.
For ID-TIMS analyses of samples FP6D and S2E, baddeleyite dissolution and chemistry were adapted from Rioux et al. (2010).
Baddeleyite crystals were plucked from the SIMS grain mount, spiked with a mixed 205 Pb/ 233 U/ 235 U tracer (ET535) and dissolved in HCl acid.Solutions were pipetted into beakers to separate them from undissolved zircon domains.Pb and UO2 from baddeleyite were loaded onto single rhenium filaments with silica gel without ion exchange cleanup.Isotopic compositions were measured on a Micromass Sector 54 TIMS at the University of Wyoming in single Daly-photomultiplier mode.The model of Stacey & Kramers (1975) at 400 Ma was used for common Pb corrections of SIMS and ID-TIMS analyses.
The decay constants and 238 U/ 235 U ratio are from Steiger and Jäger (1977).Concordia coordinates and uncertainties were calculated using IsoplotR for SIMS (Vermeesch, 2018) and PBMacDAT and IsoplotEX for ID-TIMS (Ludwig, 1988(Ludwig, , 2003)).Callegaro et al. (2017) discussed different age interpretations of SL18 extensively but were unfortunately unable to determine a robust U-Pb crystallization age.We present additional single crystal ID-TIMS data of SL18, obtained with exactly the same methodology as for Callegaro et al. (2017) at the University of Geneva (details in supplement S2), but with baddeleyite crystals from the same separate that are 10-30 times smaller.

Spread Eagle Intrusive Complex (SEIC)
The SEIC rocks have well-preserved igneous textures, but major and accessory minerals are frequently replaced by parageneses indicative of alteration and low-grade metamorphism.Grain size and colour index can vary considerably, ranging from finegrained melanogabbros to coarser-grained, rarely pegmatoidal monzonites and monzosyenites.Plagioclase is always completely altered to albite.Many samples also contain large amounts of K-feldspar.In the less potassic samples, K-feldspar is concentrated in interstitial areas or leucocratic pockets, often together with accessory minerals.Minor quartz is often present in baddeleyite-bearing and baddeleyite-free rocks, but usually in secondary pockets and veins.Clinopyroxene is replaced by chlorite to a variable extent, often forming pseudomorphs.Some samples have essentially unaltered clinopyroxene, but also contain chlorite pseudomorphs.Ilmenite is largely replaced by titanite ± rutile ± magnetite.Other common secondary minerals are calcite, epidote, prehnite and/or pumpellyite, and accessory sulphides (pyrite, chalcopyrite, galena, sphalerite).Apatite is ubiquitous.
All samples contain Zr-bearing accessory phases.Baddeleyite occurs in many samples, usually as <20 µm long, lath-shaped euhedral crystals.Crystals 20-80 µm in length occasionally occur in some samples, and FP6D is the only sample with large (50-200 µm) baddeleyite crystals (Figure 3b-d).SEIC baddeleyite is commonly intergrown with zircon, forming a variety of textures.The most common case is that baddeleyite crystals contain zircon domains mostly at their rim, but commonly penetrating into the core.This largely pseudomorphic replacement texture is accompanied by feather-like zircon coronas around the crystal (e.g., Figure 4d).Baddeleyite preservation tends to be better if the rock is less altered, but also if crystals are large, as the presence of rather modest zircon overgrowths in the strongly altered monzonite FP6D indicates.Baddeleyite inclusions in K-feldspar usually lack zircon intergrowths.In samples FP1F and FP1I, clusters of baddeleyite needles are enclosed by zircon crystals (Figure 3a).The enclosing zircon is sometimes almost euhedral, but often with feather-like zircon overgrowths.A special feature we identified in FP12A is that some baddeleyite crystals have zircon inclusions up to 3 µm wide and at most 12 µm long (Figure 4f-k).Secondary zircon overgrowth on the enclosing baddeleyite was rarely observed.
Besides that, FP12A contains prismatic euhedral baddeleyite crystals essentially free of zircon, but with a notch on one prism plane that penetrates into the crystal's core (Figure 4e, l).Zircon crystals without baddeleyite intergrowth are sometimes present also in baddeleyite-bearing rocks.Some of these crystals have an amoeboid surface and a scarred interior.Baddeleyite-free rocks are often either melanocratic or rich in quartz.They contain zircon as the only Zr-bearing mineral, forming euhedral (≤ 20 µm) or anhedral (≤ 50 µm) crystals.Only one sample (FP7G) contains zirconolite CaZrTi2O7 (all mineral formulae given as stochiometric end-members from Anthony et al., 2001) and rare pyrochlore (Ca,Na)2Nb2O6(OH,F) in addition to baddeleyite and zircon.Zirconolite occurs mainly in the vicinity of baddeleyite crystals of similar size and form (Figure 3f).It sometimes has an altered appearance and/or significant Si contents.In rare cases it is partly replaced by titanite + zircon (Figure 3g).
In contrast to the gabbros, granophyres are strongly leucocratic (> 70 vol.-%albite).Albite crystals are typically 1 × 0.5 cm large in the interior of granophyre pockets, whereas the outermost ca. 5 cm of the pockets are somewhat less leucocratic with smaller albite crystals.Nonetheless, a strong contrast of grain size and mineralogy characterizes the sharp contact between gabbros and granophyre pockets.Granophyre samples S2C (center of a pocket) and S2E (transition of a pocket center to the gabbro contact) are largely identical in major phases (albite, clinopyroxene, ilmenite, chlorite, biotite ± Ti-rich hornblende).
However, the accessory mineral assemblages are highly different.In S2C, zircon is the only Zr-bearing phase with contact to groundmass minerals.Crystals are up to 200 μm large, isometrical and growth-zoned.Shapes vary from euhedral sector-zoned with radial cracks in the outer zones (Figure 5c) over grain clusters (Figure 5f) to fan-shaped forms (Figure 5i).BSE and CL intensities are usually inversely correlated (Figure 5f, g).Some crystals have partly microporous textures and contain anhedral baddeleyite inclusions (mostly < 2 μm, max. 4 × 10 μm in size), usually in the outer zones.Sometimes the inclusions seem to retrace cracks or crystallographical orientations (Figure 5e).In S2E, the pocket interior has large (up to 200 μm) zircon crystals, but without baddeleyite inclusions, whereas baddeleyite sometimes forms isolated crystals in the groundmass.Towards the gabbro contact, baddeleyite becomes more and more predominant until the presence of zircon is mostly confined to baddeleyite replacement textures.Near the gabbro contact, baddeleyite occurs either as euhedral, lath-shaped crystals with dimensions up to 20 × 300 μm (Figure 5j), or as short prismatic crystals up to 30 μm in diameter (Figure 5k).Baddeleyite within albite usually lacks zircon intergrowths.Contrastingly, most baddeleyite within chlorite pseudomorphs is pseudomorphically replaced by zircon, sometimes containing baddeleyite relicts, and surrounded by a feather-like zircon corona (Figure 5m).Near the gabbro contact, zircon with amoeboidal grain boundaries occurs, as well as zirconolite and gittinsite CaZrSi2O7.Gittinsite has not been reported from similar rock types before.Textures include gittinsite-zirconolite intergrowths (Figure 5n) and gittinsitetitanite intergrowths branching along the fissure plains of surrounding chlorite (Figure 5o).
Common Pb contents tend to be lower than in FP7G and FP12A, but are often still significant.Zircon rims on baddeleyite and baddeleyite-free zircon from all SEIC samples yielded a wide range of 206 Pb/ 238 U dates from 142-517 Ma (Figure 6d, f).At least for sample FP6D, 206 Pb/ 238 U dates become younger with increasing U contents.Zircon analyses from SEIC samples have mostly high common Pb contents (<90% radiogenic 206 Pb).Zircon inclusions in baddeleyite (FP12A) yielded 206 Pb/ 238 U date ranges of 470-733 Ma and 297-607 Ma for baddeleyite-and zircon-based RSC, respectively.
For CSMS, baddeleyite of sample S2E (Figure 7; Table 2 9; Table 2; Table S6).The 206 Pb/ 238 U dates from these sessions are in good agreement with the CAMP age of ~201.5 Ma based on worldwide samples using zircon (Blackburn et al., 2013;Davies et al., 2017).

ID-TIMS data
ID-TIMS analyses of baddeleyite from SEIC (sample FP6D) and CSMS (sample S2E) yielded normally discordant data that form linear arrays (Figure 8; Table 3).The upper intercept dates are 498.7 ± 4.5 Ma (FP6D) and 437.0 ± 7.9 Ma (S2E).The weighted mean 207 Pb/ 206 Pb date of S2E is 444.1 ± 4.4 Ma (95% confidence; MSWD = 0.82), within error of the upper intercept date.In contrast, 207 Pb/ 206 Pb dates of FP6D show scatter beyond uncertainty (Figure 8c) with a weighted mean that is statistically meaningless (MSWD = 5.5).For both samples, there is a direct correlation between the 207 Pb/ 206 Pb dates and the percentage of discordance, leading to negative lower intercepts for linear regressions on concordia plots (Figure 8).Like the corresponding SIMS analyses, baddeleyite analyses from FP6D and S2E contained significant common Pb (up to 6 pg).We tested different common Pb isotopic compositions and blank/non-blank proportions however, and no geologically realistic choice produces consistent 207 Pb/ 206 Pb dates for sample FP6D (see supplement S5 for details), thus the increase in 207 Pb/ 206 Pb dates with increasing discordance must be attributed to the radiogenic Pb in FP6D baddeleyite.Furthermore, the negative lower intercepts are controlled by the most discordant analyses, as a regression of the 4 most discordant analyses from FP6D, the ones with the most total Pb loss, leads to a lower intercept of -462 ± 430 Ma (supplement S5), negative outside of error and similar to the regression using all 6 analyses.Additional ID-TIMS data for very small baddeleyite crystals of sample SL18 (Figure 9c; Table 3) yielded 206 Pb/ 238 U dates that are younger than those of the larger crystals of SL18 published in Callegaro et al. (2017).All analyses overlap with Concordia within uncertainty and have low common Pb contents (~0.5 pg) which are attributed to laboratory blank rather than initial common Pb or common Pb that has been added during a secondary alteration process.

Occurrence, textures and interrelations of accessory minerals
Zirconium-bearing accessory minerals in mafic magmas form during late stages of crystallization in more differentiated interstitial melt (Heaman and LeCheminant, 1993;Schaltegger and Davies, 2017).In our study, abundance and crystal size of accessory minerals lack a strong correlation with whole rock Zr content (supplement S7), but the more coarse-grained samples tend to contain larger baddeleyite and zircon crystals.Regardless of crystal sizes, baddeleyite and zircon in SEIC and CSMS rocks form various types of intergrowths.Baddeleyite in mafic rocks is typically of igneous origin, but metamorphic processes can cause it to react to polycrystalline zircon (Heaman and LeCheminant, 1993).Metamorphic zircon is therefore expected to be the most common type of zircon intergrown with baddeleyite in SEIC and CSMS rocks, which all have experienced lowgrade metamorphism.Although probably less common, magmatic zircon overgrowths on pre-existing baddeleyite can also form during late-stage igneous crystallization due to an increased SiO2 activity in the melt (e.g., Renna et al., 2011).Such igneous zircon overgrowths on baddeleyite have rather euhedral crystal faces and straight interfaces with baddeleyite (Renna et al., 2011).By contrast, metamorphic zircon shares more irregular crystal interfaces with baddeleyite and has an anhedral exterior, described as "raspberry texture" (Heaman and LeCheminant, 1993) or "frosty appearance" (Söderlund et al., 2013).
For SEIC and CSMS, the typical textural features of igneous zircon overgrowth on baddeleyite are rarely displayed (Figure 5l), whereas features of metamorphic replacement zircon are frequently observed.Zircon seems to pseudomorphically replace baddeleyite, accompanied by feather-like zircon coronas (e.g., Figure 4d, 5l), probably due to volume enlargement by the addition of silica during metamorphism.The presence of such coronas can therefore help to distinguish zircon with a baddeleyite precursor from primary zircon in altered igneous rocks.
The extent of baddeleyite replacement by zircon in this study often depends on the host minerals.Baddeleyite surrounded by chlorite shows replacement by zircon more commonly than baddeleyite in albite or epidote group minerals, and baddeleyite inclusions in K-feldspar mostly lack zircon.We attribute this to local variations in the SiO2 activity during metamorphism: the chloritization of pyroxenes liberates large amounts of Si, whereas alteration of alkali feldspars has a neutral Si balance.Si release sometimes also causes replacement of zirconolite by titanite + zircon (Figure 3g), or titanite + gittinsite (Figure 5o).
Alteration by fluids with high SiO2 activity causes baddeleyite replacement by zircon, but fluids poor in Si and rich in Ca can induce the opposite effect even in siliceous rocks (Lewerentz et al., 2019).In sample S2C, multiple µm-sized baddeleyite inclusions are hosted in the outer zones of zircon, which shows porous domains (Figure 5a), cracks (Figure 5c), and high contents of common Pb, which are typical alteration features (e.g., Corfu et al., 2003;Rayner et al., 2005;Aranovich et al., 2017).Whereas other fluid-mediated processes may also be capable of forming secondary baddeleyite inclusions in altered zircon (Lewerentz et al., 2019), the former presence of fluids with high Ca/Si ratio in S2C is likely due to widespread albitization of plagioclase.Previously reported occurrences of secondary baddeleyite inclusions in zircon are from rocks that experienced high temperature (mostly amphibolite facies) alteration (Barth et al., 2002;Aranovich et al., 2013Aranovich et al., , 2017;;Lewerentz et al., 2019), and experiments reproducing this texture were conducted at 600°C and 900°C (Lewerentz et al., 2019).
However, Cape St. Mary's sills have experienced only subgreenschist facies (Greenough and Papezik, 1986b), or at most, lower greenschist facies conditions.Hence, we present the first evidence that secondary baddeleyite inclusions in zircon can also form at low temperatures, and low-temperature reactions of zircon to baddeleyite and vice versa can occur within the same dike.
A peculiar texture in sample FP12A is baddeleyite with zircon inclusions (Figure 4e-k).BSE and CL images provide only two-dimensional petrographic information, so it can be argued whether these are only apparent inclusions as part of metamorphic zircon rims that locally penetrate the crystal.However, most of these baddeleyite crystals lack any visible zircon overgrowth, and the baddeleyite mantle is coherent even when as thin as 1 µm.Furthermore, we observed this texture repeatedly in sample FP12A, but not in other samples.Thus this texture most likely represents actual zircon inclusions of xenocrystic origin.An attempt to test this hypothesis with SIMS dating of these zircon inclusions yielded an older 206 Pb/ 238 U date only for one crystal (Figure 4e; Table S3, analysis 1_3), but the primary beam overlapped both the baddeleyite and zircon.
Furthermore, difficulties with the RSC (see Sect. 6.2.1) and/or Pb loss from possibly metamict zircon may have biased the results.Lastly, the xenocrysts may be derived from an assimilated Cambrian or Ediacaran sedimentary country rocks, whose detrital zircon is often not much older than ca.500 Ma (Pollock et al. 2009), which is roughly the age of FP12A.
In a hot, low SiO2 activity magma such as the primary magma of FP12A, xenocrystic zircon is undersaturated and dissolves (see e.g., Boehnke et al., 2013).Slow Zr diffusion in the melt limits the zircon dissolution rate, so that the melt adjacent to the zircon will develop an exponentially decreasing Zr concentration gradient (e.g., Harrison and Watson, 1983).Hence, partially dissolved xenocrystic zircon will be surrounded by a halo of elevated Zr concentration in the zircon-undersaturated magma.Such a halo of elevated Zr concentrations is a preferential location for baddeleyite nucleation, even if the bulk of the magma remains undersaturated with regard to baddeleyite.Once a nucleus is formed, baddeleyite will grow preferentially where Zr concentration is highest (at the dissolution interface).If a coherent baddeleyite mantle is formed, the zircon xenocryst becomes shielded from further dissolution.
Although such xenocrystic inclusions can be overlooked easily, the lack of previous reports suggests that this is a rare texture.This may be a consequence of numerous factors: 1) A melt with low SiO2 activity is needed, being zircon-undersaturated, but close to baddeleyite saturation.
2) The country rock must have zircon, but should not be too siliceous, because otherwise zircon would be stabilized and baddeleyite destabilized.If the country rock is only weakly consolidated, zircon liberation is facilitated.
3) High temperatures and low crystal fraction of the magma are necessary to assimilate country rock effectively.However, this also favors rapid dissolution of zircon or xenocrysts or their entrapment in major phases before baddeleyite saturation is achieved.
4) In specific cases, baddeleyite may fail to enclose zircon before the latter dissolves completely.This may leave a notch on the baddeleyite crystal, such as in Figures 4e and 4k.
Our detailed micropetrography revealed that baddeleyite and zircon can be intergrown in various ways, especially for samples with complex metamorphic histories.Including our newly proposed xenocrystic zircon inclusions in baddeleyite, at least seven different types of baddeleyite-zircon intergrowths have to be considered (compiled in Table 4).Three of these types are not observed in our samples: (1) granular baddeleyite droplets can rim zircon that decomposed to baddeleyite + SiO2 as a result of impact shockwave heating (e.g., El Goresy, 1965;Wittmann et al., 2006).(2) The inversion of this reaction was found in a shergottite sample, where primary baddeleyite is often partially rimmed by polycrystalline zircon in the vicinity of impact melt (Moser et al., 2013;Darling et al., 2016).This is the only baddeleyite-zircon intergrowth type where baddeleyite is affected by shockwave metamorphism (Moser et al., 2013;Darling et al., 2016).(3) Feather-like polycrystalline baddeleyite reaction rims on mantle-derived zircons in kimberlites were found to be the result of desilicification reactions (Kresten, 1973).These seven intergrowth types can occur in combination, complicating textural interpretation.Dating by SIMS or LA-ICP-MS provides the high spatial resolution that is required to unravel the age relationships of baddeleyite-zircon intergrowths.

Challenges in baddeleyite geochronology by SIMS
Despite many examples of good agreement between SIMS and ID-TIMS data (e.g., the SIMS 207 Pb/ 206 Pb date agrees with the ID-TIMS 207 Pb/ 206 Pb and upper intercept dates of FP6D, the SIMS in situ dates agree with both ID-TIMS dates of S2E, and all dates of SL18 and FC-4b agree), some SIMS sessions yielded significantly deviating dates, notably in case of grain mount sessions and/or 206 Pb/ 238 U dates.Although SIMS dates of baddeleyite can be more accurate than ID-TIMS dates in certain cases (see Sect. 6.2.3 and 6.2.4), baddeleyite poses analytical challenges that are specific for in situ methods: beam overlap on adjacent phases, and possible bias in the U/Pb relative sensitivity calibration (RSC).
Small crystal sizes of baddeleyite (< 10 µm) often result in primary beam overlap on adjacent phases during SIMS sessions.This does not necessarily affect the accuracy of baddeleyite dates if the adjacent minerals are U-and Pb-free (e.g., chlorite), but otherwise, especially in case of intergrowths with zircon and zirconolite, accuracy of baddeleyite dates can be severely affected.But even so, estimations of the baddeleyite crystallization ages are possible if the extents of potential resulting inaccuracies can be assessed.This requires knowledge about (1) the approximate U-Pb crystallization ages of baddeleyite and the contaminating phase -we suppose that the zircon rims of SEIC samples formed at ca. 400 Ma during the Acadian Orogeny, being ca.20% younger than SEIC baddeleyite, and that some of the zircon rims experienced Pb loss; (2) possible differences in U content of the involved minerals -baddeleyite tends to have the same (this study) or higher (Heaman and LeCheminant, 1993) U contents than zircon, and lower U contents than zirconolite (Rasmussen and Fletcher, 2004); and (3) matrix effects that lead to different U/Pb relative sensitivities.As observed from the RSC, baddeleyite 206 Pb/ 238 U dates shift to younger ages when computed against a zircon RSC (-54 to +3% total, -23% average), and zircon dates become older vice versa (+1 to +121%, +43% average).Similar bias results from baddeleyite-zirconolite beam overlap, although presently this bias cannot be quantified due to lack of well-characterized zirconolite reference materials.Notably, when SIMS baddeleyite analyses contain metamorphic zircon contributions, their younger age is partly compensated by the RSC bias, so their influence on data quality is less severe than in cases where zircon is coeval or older than baddeleyite.
Precision and accuracy of SIMS 206 Pb/ 238 U dates strongly depend on the quality of the RSC.Variable degrees of reverse discordance in different SIMS grain mount sessions of the same sample most likely reflect difficulties in quantification of matrix effects.RSC accuracy depends on numerous factors, but crystal orientation effects have been long recognized as particularly important for SIMS U-Pb dating of baddeleyite (Wingate and Compston, 2000).Although oxygen flooding of the sample chamber proved to be effective for reducing crystal orientation effects, residual bias remains (Schmitt et al., 2010;Li et al., 2010).In case of grain mounts, tabular crystals will be preferentially oriented with their c axis parallel to the sample surface.This may be a cause of matrix mismatch in grain mount analyses, whereas in situ mounts, which lack reverse discordance in this study, are expected to have more random distributions of crystal orientations.Analysis of a sufficiently large number of crystals (> ca.25) and randomizing crystal orientations of grain mounts may lead to more accurate SIMS 206 Pb/ 238 U baddeleyite dates.However, crystal orientation effects, if present, would also be expected for the secondary reference baddeleyite FC-4b, which nevertheless yielded 206 Pb/ 238 U dates in agreement with published ID-TIMS data (Schmitt et al., 2010) although it was analyzed during the same grain mount sessions as FP6D and S2E baddeleyite.Hence, other factors such as alteration may be more significant for the RSC than crystal orientation.

Common Pb in Zr-bearing minerals
Ideal U-Pb mineral geochronometers exclude common Pb during crystallization, a behavior that baddeleyite, zircon, and zirconolite approximately fulfill (e.g., Heaman and LeCheminant, 1993;Rasmussen and Fletcher, 2004).However, in our samples, all three minerals often contain abundant common Pb, as evident from SIMS and ID-TIMS analyses of SEIC and CSMS (Tables 3 and S1-S6).Similar abundances of common Pb were also detected in baddeleyite from other mafic dikes (e.g., Olsson et al., 2011).Adjacent phases, surface contamination (SIMS), laboratory blank (ID-TIMS) and mineral inclusions (SIMS and ID-TIMS) can be external sources of common Pb, but steady 204 Pb counting rates even in some SIMS analyses with the spot entirely on baddeleyite demonstrate that some amount of common Pb is intrinsic to these crystals.In case of zircon, it is known that common Pb can be incorporated during alteration (e.g., Watson et al., 1997;Rayner et al., 2005;Geisler et al., 2007).By analogy, the interplay of radiation damage and interaction with fluids may trigger a similar process in baddeleyite and zirconolite.Alteration could also modify the chemical composition of baddeleyite and may therefore bias the RSC, but a correlation of 206 Pb/ 238 U date and common Pb content is detectable only for sample FP7G, which also shows a correlation of common Pb content and the (RSC-independent) 207 Pb/ 206 Pb date.Moreover, common Pb contents of baddeleyite in FP7G and FP12A increase with decreasing U content, making it difficult to explain common Pb incorporation with radiation damage, unless U was mobilized as well.Uranium mobilization can cause normal or reverse discordance and may therefore be detectable by further ID-TIMS analyses.Furthermore, baddeleyite and zircon can be expected to alter differently, as baddeleyite is more resistant to radiation damage (Schaltegger and Davies, 2017).Oxygen isotope analysis by SIMS has been suggested as a tool to detect baddeleyite alteration (Davies et al., 2018), although its routine use is still in its infancy.

Baddeleyite discordance
Although volume diffusion of Pb in baddeleyite is thought to be extremely slow (Heaman and LeCheminant, 2000), baddeleyite analyses often show a few percent of discordance.Naturally observed discordance patterns range from linear and complex arrays to uniform clusters beneath concordia (e.g., Söderlund et al., 2013;Schaltegger and Davies, 2017;Heaman & LeCheminant, 2000;Ibañez-Mejia and Tissot, 2019).This underscores that discordance is not always created by the same mechanism, and it is crucial for age interpretation to identify the correct discordance mechanism(s) in the specific sample by means of the discordance pattern (Figure 10).If contributions from metamorphic zircon rims have caused discordance, a threecomponent mixing model can be applied, with end-members defined by (1) the igneous crystallization age of baddeleyite, (2) metamorphic formation of zircon, and (3) recent Pb loss of zircon and/or baddeleyite (Söderlund et al., 2013).The oldest 207 Pb/ 206 Pb date would then yield a minimum estimate of (1).In contrast, 231 Pa excess would influence only the 207 Pb/ 235 U decay system, making 206 Pb/ 238 U ages most accurate (Amelin and Zaitsev, 2002;Crowley and Schmitz, 2009), and vice versa in case of 222 Rn loss (Heaman and LeCheminant, 2000).Pb loss due to alpha recoil or fast pathway diffusion (e.g., along twinplanes) affects both systems similarly, making 207 Pb/ 206 Pb ages most accurate (Davis and Davis, 2018).
Metamorphic zircon overgrowth was absent in baddeleyite used for ID-TIMS analysis of S2E and SL18.Even in FP6D, where it is petrographically evident, the ID-TIMS data appear to be free of a significant isotopic component of metamorphic zircon.This confirms that baddeleyite and zircon can be separated successfully with the method of Rioux et al. (2010), using only hydrochloric acid for dissolution.Consequently, discordance in our samples should be attributed rather to baddeleyite itself than to zircon intergrowths.Our ID-TIMS analyses show linear arrays that are typical for varying degrees of recent Pb loss.
The portion of radiogenic Pb lost by alpha recoil can be predicted based on crystal shapes (Davis and Davis, 2018) and is generally <0.3% for the crystals of this study.However, Pb loss indicated for our samples exceeds this extent often by more than one order of magnitude (Table 3).Hence alpha recoil can explain only a small part of the discordance observed for FP6D, S2E and SL18, depending on the U zonation of the crystals.Fast pathway and/or volume diffusion is therefore likely to dominate baddeleyite discordance here.The radiation dose does not seem to be crucial: except for the in situ session of S2E, negative correlations between U content and 206 Pb/ 238 U dates appear to be absent in SIMS and ID-TIMS data of our samples (cf.Söderlund et al., 2013).
Intriguingly, the discordia trends for samples FP6D and S2E have negative lower intercepts and show a positive correlation of the 207 Pb/ 206 Pb date with the percentage of discordance (Figure 8).Similar discordance patterns have been observed from baddeleyite ID-TIMS analyses of the Duluth gabbro (Hoaglund, 2010;Ibañez-Mejia and Tissot, 2019) which plot essentially within the concordia error swath, but yielded older 207 Pb/ 206 Pb dates, younger 206 Pb/ 238 U dates and undistinguishable 207 Pb/ 235 U dates compared to coexisting zircon.These studies interpreted 231 Pa excess to be responsible for this pattern, which cannot be excluded as evidence for modest amounts of excess 231 Pa was recently obtained from young baddeleyite (Sun et al., 2020).
However, the negative lower intercepts for our Newfoundland samples would require that the magnitude of 231 Pa excess increases with discordance, which is difficult to rationalize.For the Duluth gabbro analyses, the overlap of baddeleyite and zircon 207 Pb/ 235 U dates is impossible to produce with 231 Pa excess alone.We therefore propose preferential loss of 206 Pb, possibly due to 222 Rn mobility, as a more likely explanation of these peculiar patterns (see, e.g., Heaman and LeCheminant, 2000), as such excess 206 Pb loss could simply increase with overall Pb loss (see supplement S6 for more details).Negative lower intercepts have rarely been reported in previous baddeleyite dating studies, suggesting that they are either caused by a rare process or are often masked by uncertainties and/or additional discordance mechanisms.
Our data also imply that baddeleyite rims are more strongly affected by Pb loss than cores.The ID-TIMS data of SL18 (Figure 9c) show younger 206 Pb/ 238 U dates for the smallest crystals, which have larger surface to volume ratios than the larger crystals.
Furthermore, the ID-TIMS analyses of large SL18 baddeleyite yielded younger 206 Pb/ 238 U dates than the SIMS analyses (Table 3 vs.Table S6).The SIMS spots were typically placed in the centers of the grains, avoiding the apparently less discordant rims, whereas dissolution of the plucked grains for ID-TIMS analysis included them, resulting in more discordant ID-TIMS analyses.Admittedly, it cannot be excluded that SIMS 206 Pb/ 238 U data of the relatively pristine SL18 baddeleyite may be also affected by RSC bias, similar to the altered Newfoundland samples.However, this concern is irrelevant for the 207 Pb/ 206 Pb data, which agree within error between SIMS and ID-TIMS for SL18, but 207 Pb/ 206 Pb dates from the more discordant Newfoundland samples tend to be older for ID-TIMS than for SIMS.This may imply that baddeleyite rims in FP6D have undergone preferential 206 Pb loss, which would drive the more discordant ID-TIMS analyses towards older 207 Pb/ 206 Pb dates.
In any case, fast pathway and/or volume diffusion can both explain intensified Pb loss from crystal rims, but to differentiate between these processes, both the U zonation within the crystals and the post emplacement thermal history of the samples are not sufficiently known.

Implications for baddeleyite dating approaches
Targeting the least altered and most coarse-grained samples in the field is highly advantageous for recovery of pristine baddeleyite which often yields concordant ages.However, as the examples of SEIC and CSMS rocks show, alteration can be ubiquitous in certain geologic units, so that usage of texturally complex, discordant baddeleyite is inevitable.But even baddeleyite from unaltered samples, such as SL18, often shows evidence for minor Pb loss (see also Söderlund et al., 2015).
With Pb loss as a dominant discordance mechanism, 206 Pb/ 238 U dates of baddeleyite often underestimate intrusion ages whereas 207 Pb/ 206 Pb dates are more accurate.In case of SIMS, another advantage of 207 Pb/ 206 Pb dates is their independence from the RSC.This favors the use of 207 Pb/ 206 Pb rather than 206 Pb/ 238 U ages even for Phanerozoic baddeleyite, in spite of its typically lower precision.For Mesozoic samples such as SL18, however, comparatively low precisions of SIMS data make it difficult to use 207 Pb/ 206 Pb dates meaningfully, although improved 207 Pb/ 206 Pb precisions have been achieved with a multi-collection SIMS method (Li et al., 2009).Another limitation for 207 Pb/ 206 Pb dates is possible bias by preferential 206 Pb loss, as suggested by the correlation of 207 Pb/ 206 Pb dates with the percentage of discordance (Figure 8).If this bias is important, the 207 Pb/ 206 Pb dates of the least discordant analyses are most accurate, and more discordant analyses should be excluded from weighted mean calculations.Alternatively, if the excess 206 Pb loss is a linear function of total Pb loss, upper intercept dates are the most accurate.This can be shown for sample FP6D, where the upper intercept date from ID-TIMS (499 ± 5 Ma) and the 207 Pb/ 206 Pb date from SIMS (497 ± 7 Ma) are indistinguishable, but most of the ID-TIMS 207 Pb/ 206 Pb dates are considerably older (504 to 530 Ma, Table 3, Figure 8).If baddeleyite cores are less discordant than rims, the cores should be preferentially targeted for analysis, which is easier to achieve with SIMS.In the case of ID-TIMS, mechanical abrasion is potentially helpful in this respect, and there are documented cases where discordance was reduced by mechanical abrasion (e.g., Corfu and Lightfoot, 1996;cf. e.g., Greenough et al., 1993).Alternatively, baddeleyite cores can be cut out with a focused ion beam before ID-TIMS analysis (White et al., 2020).

Intrusion ages of the SEIC, CSMS and FLC
Stratigraphic constraints require that the SEIC is coeval with the Cambrian volcano-sedimentary succession on the western Avalon Peninsula (Sect.2.1).Trilobite biostratigraphy suggests that the age span of deposition of the Manuels River Formation roughly equals that of the Drumian stage (Hildenbrand, 2016;Austermann, 2016), which is currently calibrated from 504.5 to 500.5 Ma with uncertainties on these bounds of ca. 2 Ma (Peng et al., 2012).The occurrence of pillow lavas in the Chamberlain's Brook Formation and Elliot Cove formation expands the age range of Cambrian volcanism on the Avalon Peninsula to both pre-and post-Drumian (Figure 2).Our SIMS and ID-TIMS data are thus in good agreement with these stratigraphic constraints.We regard the ID-TIMS concordia upper intercept date of FP6D (498.7 ± 4.5 Ma) as the best available estimate of the intrusion age of the corresponding feeder pipe, as it is more likely unbiased by the sample-specific baddeleyite discordance mechanisms than the 206 Pb/ 238 U and 207 Pb/ 206 Pb dates.As indicated by the stratigraphic distribution of volcanic rocks, the other feeder pipes may be significantly older and/or younger, but this first radiometric age for the SEIC indicates that SEIC magmatism clearly predates the opening of the Rheic Ocean (see Sect. 2).
For the CSMS, zircon from granophyre sample S2C is considerably altered and the degree of secondary Pb loss is often very large (supplement S4), limiting its use for geochronology.Baddeleyite from sample S2E, although potentially also altered, better preserves the igneous age.Our ID-TIMS 207 Pb/ 206 Pb baddeleyite date (444.1 ± 4.4 Ma) agrees within error with that of Greenough et al. (1993), and the samples of both studies are derived from the same granophyric dike (Greenough, pers. comm.).Greenough et al. (1993) analyzed bulk separates of baddeleyite, and the large number of crystals per aliquot in their analyses enhanced precision due to stronger Pb signals, but may have obscured the discordance patterns that we observed in our single crystal analyses.Combining the data of both studies, we regard a weighted mean ID-TIMS 207 Pb/ 206 Pb age of 439.4 ± 0.8 Ma (95% conf.; MSWD = 0.94; Figure 11) as the best available estimate of the intrusion age of the Lance Cove sill, using only the analyses with <3% discordance to minimize bias due to possible preferential 206 Pb loss.If this bias is significant even for the least discordant analyses, this 207 Pb/ 206 Pb date may be an overestimate, but the ID-TIMS upper intercept date of 437 ± 8 Ma, which would be more accurate in this case, does not lead to a better age estimate due to inferior precision.The date that we report does not rule out the possibility that other sills of Cape St. Mary's are significantly older or younger.
For the FLC, our new ID-TIMS and SIMS data suggest that the intrusive age based on a weighted mean 206 Pb/ 238 U date reported in Callegaro et al. (2017) from large baddeleyite crystals is likely too young and reflects some degree of Pb loss.The small baddeleyite crystals from SL18 yielded younger ID-TIMS 206 Pb/ 238 U dates due to more intense Pb loss, likely due to fast pathway or volume diffusion.Therefore, we regard the weighted mean207 Pb/ 206 Pb ID-TIMS date (201.07 ± 0.64 Ma) as the best estimate for the intrusion age of SL18.This age is in agreement with the SIMS dates and all other age constraints from both the FLC and the CAMP (Davies et al., 2017;Callegaro et al., 2017).This new age does not change the overall interpretations of Callegaro et al. (2017).The good agreement of the SIMS and ID-TIMS dates from sample SL18 validates that baddeleyite dating by SIMS, which has been applied so far mostly for Precambrian and early Paleozoic samples (e.g., Schmitt et al., 2010), can also yield accurate ages for Mesozoic samples.Table 4. Compilation of all known baddeleyite (Bad)-zircon (Zrn) intergrowth types, and guideline for their interpretation.

(
207 Pb/ 206 Pb).Many baddeleyite analyses show surprisingly high contents of common Pb, those with <90% radiogenic 206 Pb were generally excluded from weighted mean calculations.During grain mount sessions, FC-4b baddeleyite was analyzed as a secondary reference in addition to Phalaborwa baddeleyite.Weighted mean 206 Pb/ 238 U dates of FC-4b calculated with Phalaborwa reference are 1118 ± 39 Ma (MSWD = 0.63; n = 28), 1101 ± 44 Ma (MSWD = 0.41; n = 29), 1124 ± 56 Ma (MSWD = 0.91; n = 9) and 1117 ± 23 Ma (MSWD = 2.42; n = 18; session with sample S2E).Therefore, in all grain mount on new and published microtextural observations, at least seven different types of baddeleyite-zircon intergrowths can be identified.Secondary baddeleyite inclusions in zircon can also form during low-temperature alteration.Furthermore, we interpret a previously undocumented texture as xenocrystic zircon cores mantled by igneous baddeleyite overgrowths.2) SIMS 206 Pb/ 238 U data, at least from grain mounts, show bias in the U/Pb relative sensitivity calibration, possibly due to crystal orientation effects and/or alteration.3)Unusually high common Pb contents in baddeleyite are evident from SIMS and ID-TIMS data of the same crystals and probably a consequence of alteration.4) For our samples, we identified secondary Pb loss, due to fast pathway and/or volume diffusion, as the dominant discordance mechanism.Correlations between 207 Pb/ 206 Pb dates and degree of discordance imply that there may be a component of preferential206 Pb loss, possibly due to 222 Rn mobility.Baddeleyite rims appear to be more affected by Pb loss than cores.5) Any kind of Pb loss makes 207 Pb/ 206 Pb ages more reliable than 206 Pb/ 238 U ages, even for Phanerozoic baddeleyite.In case of preferential 206 Pb loss, the most accurate age is either the concordia upper intercept date or the weighted mean 207 Pb/ 206 Pb date of the least discordant analyses.Improved concordance of baddeleyite dates may be achieved when preferentially targeting the less discordant crystal cores, either by SIMS spot analysis, or possibly by applying mechanical abrasion prior to ID-TIMS analysis.6) We present new and refined intrusion ages for the Spread Eagle Intrusive Complex, Newfoundland (498.7 ± 4.5 Ma, ID-TIMS concordia upper intercept), which formed during rifting before the opening of the Rheic Ocean, as well as for Cape St. Mary's sills, Newfoundland (444.1 ± 4.4 Ma, ID-TIMS weighted mean 207 Pb/ 206 Pb date of the least discordant analyses) and the Freetown Layered Complex, Sierra Leone (201.07 ± 0.64 Ma, ID-TIMS weighted mean

Figure 6 .
Figure 6.SIMS U-Pb baddeleyite and zircon results for SEIC samples.Ellipses in red represent analyses of baddeleyite with zircon inclusions (e) or zircon analyses with <90% radiogenic 206 Pb (d, f).Error ellipses of individual analyses are 1σ whereas the weighted mean ellipse (blue) is enlarged to 2σ.

Figure 7 .
Figure 7. SIMS U-Pb baddeleyite (a, b) and zircon (c) results for CSMS samples.Error ellipses of individual analyses (only those used for weighted mean date calculations are shown) are 1σ whereas the weighted mean ellipse (blue) is enlarged to 2σ.

Figure 8 .
Figure 8. ID-TIMS U-Pb results for samples FP6D and S2E. 207Pb/ 206 Pb dates are less scattered than 206 Pb/ 238 U dates, but show a linear correlation with the percentage of discordance.All error bars and ellipses represent 2σ uncertainties.

Figure 9 .
Figure 9. U-Pb baddeleyite data of sample SL18.For the SIMS analyses (a, b), error ellipses of individual analyses are 1σ whereas the weighted mean ellipse (blue) is enlarged to 2σ.For the ID-TIMS data (c) of large crystals (blue; Data from Callegaro et al., 2017) and small crystals (red; this study), all ellipses represent 2σ uncertainties.

Figure 10 .
Figure 10.Schematic model of the consequences of different discordance mechanism on the isotopic composition of baddeleyite.