Rocks from the Blair River inlier of Northern Cape Breton Island (Nova Scotia, Canada) have been correlated with either the Grenville basement of eastern Laurentia or the accreted Avalon terrane. Additional zircon U-Pb dates of spatially associated anorthositic dykes ( Ma) and a metagabbro ( Ma) from the Fox Back Ridge intrusion of the Blair River inlier reveal Late Silurian emplacement ages. Their contemporaneity suggests that they may be members of a larger intrusive complex. The anorthositic rocks have high Eu/Eu values (>2.5), and bulk compositions are similar to the mineral compositions of labradorite (An50-70) and andesine (An30-50). The metagabbro is compositionally similar to alkali basalt and does not seem to have been affected by crustal contamination (; ) although it was metamorphosed. The high Tb/YbN (1.8-1.9) ratios suggest that the parental magma of the metagabbro was derived from a garnet-bearing peridotite. Fractional crystallization and mass balance calculations indicate that the anorthositic rocks can be derived by mineral accumulation from a mafic parental magma similar in composition to the metagabbro of this study. The Late Silurian ages suggest that the rocks were emplaced into the Avalon terrane after the closure of the Iapetus Ocean but before Early Devonian (415-410 Ma) sinistral transpression.
The eastern margin of North America records a complex middle Proterozoic to middle Paleozoic geological history. The most complete record is preserved in a belt of basement inliers with overlying cover sedimentary sequences that extend discontinuously from Alabama to Newfoundland along the western flank of the Appalachian orogen [1–3]. The inliers are composed of middle to late Proterozoic rocks that were either formed or deformed during the Grenvillian orogeny and dated at 1200-980 Ma . These inliers include rocks of easternmost Laurentia that are elsewhere typically included in the Appalachian orogen . They are distinct from the outboard Appalachian terranes (i.e., Avalonia, Meguma, Carolinia) composed of Neoproterozoic and Paleozoic units that accreted to Laurentia during Paleozoic orogenesis [6, 7]. Laurentian basement inliers in the Appalachian orogen include distinctive lithologies such as anorthosites, charnockites, and granulites. These rock types are common in the Grenville Province but are mostly missing from the Appalachian outboard terranes.
Although earlier studies proposed a lithologic correlation between rocks in northern Cape Breton Island (defined as the Blair River inlier) and basement rocks of western Newfoundland (e.g., [8, 9]), subsequent tectonic division models considered Cape Breton Island to be part of the exotic Avalon terrane [10–12]. On the other hand, Barr and Raeside  revived the correlation of the Blair River inlier with the Humber Zone (Laurentian margin) as part of their redefinition of tectonic terranes in the northern Appalachian orogen. At present, the basement rocks of the Blair River inlier in Northern Cape Breton Island have different interpretations. Brown , Barr and Raeside , and Miller and Barr  correlate the inlier with Grenvillian (Laurentian) basement of the Humber Zone, whereas Neale and Kennedy  correlated them with the Gander terrane (peri-Gondwana) and Keppie  and Keppie et al.  interpreted them as the basement rocks of the Avalon terrane (peri-Gondwana), thereby providing a source for the common 1-1.3 Ga detrital zircons in Avalonia (Figure 1). A critical factor in the correlation with Grenvillian basement is the presence of anorthosite bodies, which are restricted to the area of the Blair River inlier. In order to provide a better database for the comparison with anorthosite of the Grenville Province, we have investigated these Cape Breton rocks.
The Blair River inlier is located at the northernmost tip of Cape Breton Island, Nova Scotia, and is comprised of Mesoproterozoic to Carboniferous rocks. It is bounded by the Gulf of St. Lawrence to the west and north and separated from the adjacent Aspy terrane by the Wilkie Brook fault to the east and the Red River fault to the south. The Mesoproterozoic (~1.2 Ga to ~1.0 Ga) rocks consist of composite orthogneiss complexes (i.e., Sailor Brook gneiss, Otter Brook gneiss, and Polletts Cove River gneiss) that are intruded by comparatively less deformed plutonic rocks. The plutonic rocks include gabbro, syenite, anorthosite, granite, and granodiorite [15, 17]. Single-grain analyses of abraded zircon from the anorthosite of the Blair River inlier (Red River anorthosite) indicate that the minimum crystallization U-Pb age was Ma although there is evidence suggesting that the zircons experienced multiple stages of lead loss . The uniqueness of the Blair River inlier relative to the peri-Gondwanan terranes (i.e., Avalon, Gander, and Meguma) in the northern Appalachians suggests that it may represent a fragment of Laurentian basement and that the Red River anorthosite is interpreted to be correlative with Proterozoic massif-type anorthosite found in the Grenville Province [15, 17, 18]. However, recent zircon U-Pb geochronology on the Red River anorthosite indicates that it was emplaced during the Late Silurian ( Ma) and contains Meso- to Neoproterozoic ( Ma to 1044 ± 20 Ma)-inherited zircons . An evaluation of the whole-rock Pb isotopes of the Blair River inlier rocks, including the Red River anorthosite, indicates that it has greater affinity to peri-Gondwana terranes than to Laurentia . Moreover, the absence of ophiolites between the Blair River inlier and the neighbouring Aspy terrane, which are present in southern Newfoundland, argues against a Laurentia-Avalon suture zone [16, 19].
In this paper, we present new zircon U-Pb geochronology and whole-rock geochemistry of metagabbro from the Fox Back Ridge diorite/gabbro intrusion and an anorthositic dyke from the Blair River inlier. Our results provide new limits on the emplacement age of anorthositic rocks within the Blair River inlier, which, in turn, have direct implications for the location of the boundary between the cratonic margin of Laurentia and the accreted Avalon terrane. Moreover, we present a petrogenetic model that attempts to constrain the formation of the anorthositic dykes and the Red River anorthosite.
2. Geological Background
The Blair River inlier includes several anorthosite bodies, the largest of which is the Red River anorthosite suite (Figure 2). This intrusion is about 11 km long and up to 1 km wide. It occurs within a NNE-trending belt of quartz-feldspathic gneiss, which is about 30 km long and 2-8 km wide . The belt is bounded by mylonite zones, and thus, any direct correlation with other units in the area is tenuous. The Red River anorthosite is representative of the Cape Breton anorthosites . It is composed of anorthosites passing in places into interlayered anorthositic gabbro and pyroxene-rich rocks [15, 21]. The anorthosite grades from massive and equigranular to foliated. Its texture ranges from protoclastic and porphyroclastic to recrystallized granoblastic. The dominant mineral is plagioclase (An35-50) with grain size ranging from 1 to 5 mm . The anorthositic rocks typically contain minor amphibole with rare pyroxene relicts in the core. Both the anorthosite and the surrounding gneisses of the belt have experienced polyphase deformation, which has transposed the contacts into parallel features. The rocks within and bounding the Blair River inlier were metamorphosed to amphibolite facies along shear zones. Deformation was accompanied by partial melting and mylonitization in some cases. In between these shear zones, anorthosites, monzodiorite-gabbro, and granitic and syenitic plugs are undeformed.
A suite of eight samples was collected from the Blair River inlier. Two samples from the main body of the Red River anorthosite, three samples from an anorthositic dyke that intrudes the Fox Back Ridge diorite/gabbro intrusion, and three host metagabbroic rocks of the Fox Back Ridge intrusion  were collected. The metagabbro unit was originally mapped together as orthogneiss and paragneiss with minor amphibolite- and kyanite-bearing horizons by Smith and MacDonald  and assigned to the Helikian (850 Ma to 1600 Ma) Pollets Cover River Group by Barr et al. . However, Miller  reexamined the area and identified the region west of the Red River anorthosite as a distinct intrusive unit and named it the Fox Back Ridge intrusion of probable Paleozoic age . The dated samples were collected along the Cabot Trail highway immediately west of the Red River anorthosite (Figure 3). Sample CBCH-1 is an amphibole-rich metagabbro, and CBGR-1 is an anorthositic dyke that cross-cuts the Fox Back Ridge intrusion.
Sample CBGR-1 (46°4818.17N, 60°4157.08W) was collected from an undeformed anorthositic dyke that intruded the host Fox Back Ridge intrusion. It is comprised of 90-95 vol.% plagioclase with minor (~5 vol.%) pyroxene and accessory amounts (<<1 vol.%) of Fe-Ti oxide minerals and apatite. The plagioclase crystals tend to be blocky rather than tabular and display polysynthetic twinning under cross-polars. Generally, the crystals have similar grain size but some, based on the twinning angle, appear to be bent (Figures 4(a) and 4(b)). The pyroxene crystals are small, anhedral, and interstitial to the plagioclase. Although both orthopyroxene and clinopyroxene are identified in the anorthosite, most pyroxene crystals are altered to uralite suggesting that it was likely clinopyroxene . Secondary muscovite is replacing some plagioclases, and a narrow carbonate vein cross-cuts a portion of the thin section.
The metagabbro (CBCH-1; 46°4831.57N, 60°4144.53W) from the Fox Back Ridge intrusion is coarse to medium grained and is comprised primarily of replacement minerals of plagioclase (45-50 vol.%) and clinopyroxene (40-45 vol.%), with abundant apatite (~5 vol.%), and minor amounts (≤5 vol.%) of secondary oxide minerals (Figures 4(c) and 4(d)). The plagioclase is mostly altered to saussurite, but there is a minority of unaltered crystals that retain the euhedral to subhedral texture and polysynthetic twinning. Similar to plagioclase, there is a minority of small, subround, unaltered clinopyroxene crystals but most crystals are altered to uralite. The original rock texture is difficult to discern as there is a weak foliation of the uralite and saussurite, but it appears that the rock may have been ophitic. Apatite is relatively abundant and occurs as either small euhedral crystals or larger tabular crystals. The oxide minerals, including hematite, appear to be exclusively secondary as they are commonly found within or spatially associated with uralite. Short and narrow quartz veins are present but comprised <2 vol.%.
4.1. Zircon U-Pb Geochronology
Zircon separation from samples CBCH-1 and CBGR-1 was processed using magnetic separation and heavy-liquid techniques at the Yu-Neng Rock and Mineral Separation Company (Lanfang, Hebei, China). Cathodoluminescence (CL) images were taken at the Institute of Earth Sciences, Academia Sinica, Taipei, for selecting suitable positions for U-Pb analyses and to examine individual crystal internal structures. Zircon U-Pb isotopic analyses were performed by laser ablation-inductively coupled plasma mass spectrometry (LA-ICP-MS) at the Department of Geosciences, National Taiwan University, Taipei, using an Agilent 7500s Q-ICP-MS and a Photon Machines Analyte G2 193 nm laser ablation system following the method of Chiu et al. . A spot size of 35 μm with a laser repetition rate of 5 Hz was used, and the laser energy density was 3.83 to 5.33 J/cm2. Calibration was performed by using the zircon standard GJ-1 ( Ma; ) and 91500 (1065 Ma; ). The Plešovice ( Ma; ) zircon standard was also used for data quality control. Measured U-Th-Pb isotope ratios were calculated using the GLITTER 4.4.4 software , and the relative standard deviations of reference values for GJ-1 were set at 2%. The common lead was corrected using the common lead correction function of Andersen [30, 31], and the weighted mean U-Pb ages and concordia plots were created using Isoplot v. 4.1 . Only results within 10% discordance were used for the data interpretations.
4.2. Major and Trace Elemental Geochemistry
Whole-rock major and trace elements were analyzed at the Activation Laboratories Ltd. in Ancaster, Ontario, Canada, using lithium metaborate-tetraborate fusion. Major elements were determined by an inductively coupled plasma-optical emission spectrometer, whereas trace element analyses were performed by an inductively coupled plasma-mass spectrometer. Replicate analyses of the reference standard rocks indicate that 1 sigma error is between 2 and 8% of the values cited. The detailed information on the major and trace element analyses is available at the Activation Laboratories web site (https://www.actlabs.com).
The CL images show that some crystals have bright cores with darker rims, whereas others are entirely bright or entirely dark. The crystals have a variety of habits that include euhedral prismatic, subhedral round to subround, and anhedral fragmented (Fig. DR1). Individual zircon grains of the anorthosite sample (CBGR-1) range in size from ~25 μm to ~100 μm in length/width and typically show oscillatory zonation (Fig. DR4). Ninety-five of 97 analyses of zircons are within 10% discordance and yielded a total range of 206Pb/238U ages from Ma (1σ) to Ma (1σ). Four populations are identified based on the clustering of the data (Table DR1). The oldest three zircons yielded ages of Ma (1σ), Ma (1σ), and 794 ± 18 Ma (1σ). The second oldest group consists of 13 middle to late Ordovician zircons. This group appears to cluster into two subgroups: (1) Ma (1σ) to Ma (1σ) and (2) Ma (1σ) to Ma (1σ). The largest group (76) of zircons range in age from Ma (1σ) to Ma (1σ), whereas the youngest group is comprised of two Early Devonian zircons ( Ma (1σ) and Ma (1σ)). The oldest zircons and the Middle to Late Ordovician zircons are interpreted to represent inheritance, whereas the youngest zircons may be related to the effects of regional deformation during the Early Devonian. Consequently, we interpret the crystallization age to be based on the Late Silurian zircons which yielded a weighted-mean 206Pb/238U age of Ma (2σ, Figure 5(a)).
Individual zircon grains of the metagabbro (CBCH-1) range in size from ~50 μm to ~225 μm (Fig. DR2), typically show oscillatory zonation, and have fragmented, anhedral to subhedral textures. The zircons from this sample tend to have uniform CL brightness compared to those in the anorthosite (CBGR-1). Seventy-seven of 82 analyses are within 10% discordance and yielded a total range of 206Pb/238U ages from Ma (1σ) to Ma (1σ). Three distinct zircon age populations are identified based on age gaps (Table DR2). The largest (65 zircons) group ranges from Ma (1σ) to Ma (1σ). The oldest group (6 zircons) ranges from Ma (1σ) to Ma (1), and the youngest (6 zircons) group ranges from Ma (1σ) to Ma (1σ). The oldest age group consists of two subgroups where four of the zircons are between Ma (1σ) and Ma (1σ), and the remaining two are Ma (1σ) and Ma (1σ). This group is considered to be inherited due to their relatively minor abundance but also because middle to late Ordovician rocks are known in Cape Breton and throughout Avalonia [33, 34]. We interpret the crystallization age to be indicative of the weighted-mean 206Pb/238U age ( Ma, 2σ) of the zircons that comprise the middle group (Figure 5(b)). We note that the error ellipses of the zircons that comprise the weighted-mean age are slightly discordant which may be related to Pb loss during the development of the foliation within the metagabbro.
5.2. Whole-Rock Geochemistry
Five anorthosites and three metagabbros from the Blair River inlier were analyzed for whole-rock and trace elemental analyses (Table 1). The anorthosite samples are alkalic and can be divided into two groups (Figure 6(a)). One group (BR-1-1, BR-1-2) from the Red River anorthosite is broadly similar to the mineral chemistry of labradorite (SiO2 ≈53.7 wt%, Al2O3 ≈ 26.2 wt%, CaO ≈9 wt%, Na2O ≈4.8 wt%, and K2O ≈0.6 wt%), although Fe2O3t, MgO, and TiO2 are higher and reflect the presence of oxide and mafic silicate minerals in the mode (Figure 6(b)). In contrast, the second group (BR-3-1, BR-3-2; CBGR-1) from the anorthositic dykes are somewhat similar to andesine (SiO2 ≈56 wt%, Al2O3 ≈ 24 wt%, CaO ≈5 wt%, Na2O ≈6.2 wt%, and K2O ≈3 wt%) but K2O, TiO2, Fe2O3t, and MgO are higher and reflect the presence of oxide minerals, mafic silicates (biotite), K-feldspar, and mica (Figure 6(b)). Including data from Miller et al. (2000), the anorthositic rocks have a relatively restricted SiO2 (wt%) content but the Fe/(Fe + Mg) ratio is variable (~0.3 to ~0.8) (Figure 6(c)). Both types of anorthosite are peraluminous (Al/Ca+Na+K (mol.) ) and have relatively high volatile content as loss on ignition (LOI) is between 1.4 and 2.3 wt%. The quantities of transition elements for the anorthosites are low (<30 ppm) across all elements with Zn having the highest concentrations (50 to 100 ppm). The primitive mantle-normalized incompatible element patterns show enrichment of the large ion lithophile elements (Cs, Rb, Ba, and Sr), and Th and U relative to the rare earth elements and the remaining high-field-strength elements (Figure 7(a)). The chondrite-normalized rare earth element patterns show enrichment of the light rare earth elements with relatively flat heavy rare earth elements. The prominent Eu-anomaly () is typical of anorthositic rocks and is a direct consequence of their feldspar-rich nature (Figure 7(b)).
The gabbroic rocks of this study (BR-2-1, BR-2-4, and CBCH-1) are compositionally similar to alkali basalt as they have low SiO2 (<47 wt%), are alkalic ( wt%), and have moderate Mg# values (~45) and TiO2 (~2.0 wt%) contents (Figure 6(d)). The P2O5 is unusually high (~1.5 wt%), and the loss on ignition values are relatively high ( wt%). The transition metals (Sc ≈ 20 ppm, ppm, ppm, ppm, and ppm) are relatively low although V is high (240 ppm to 310 ppm). The primitive mantle-normalized incompatible element pattern of the gabbro is similar to ocean-island basalt (OIB), but there is enrichment in Cs, Ba, and a minor depletion of Zr (Figure 7(a)). The chondrite-normalized rare earth element patterns of the gabbro are also similar to ocean-island basalt as they have very high La/YbN (22 to 25) values without a significant Eu anomaly (Figure 7(b)).
6.1. Regional Temporal Correlation
The Late Silurian emplacement ages of the metagabbro ( Ma) and anorthosite dyke ( Ma) confirm that the two units are contemporaneous. Consequently, it is possible that the intrusions may be members of the same igneous complex or magmatic system and either directly or indirectly related. There are significant local and regional lithotectonic implications given that the anorthositic rocks are Late Silurian rather Mesoproterozoic as previously interpreted.
Locally, the zircon U-Pb ages are similar to titanite, rutile, and hornblende ages reported for the dioritic rocks of the Fox Back Ridge intrusion ( Ma, U-Pb titanite; Ma, 40Ar/39Ar hornblende) and Red Ravine syenite ( Ma, U-Pb titanite) suggesting that they may be members of a larger igneous complex (Table 2; ). Moreover, in western Cape Breton, there are a number of Silurian plutonic rocks ( Ma; Ma; Ma, Ma; Ma; Ma; Ma; Ma; Ma; and Ma) that have 206Pb/238U ages from Ma to Ma [17, 19, 33, 35]. In addition, middle to late Silurian ( Ma to Ma) postmetamorphic cooling ages are reported for a variety of rocks in the Blair River inlier .
Regionally, middle to late Silurian magmatic rocks are common in the Gander terrane of Newfoundland and New Brunswick ( Ma to Ma), and Silurian ( Ma to Ma) magmatic rocks (Notre Dame arc, Attean pluton, Quimby Formation) were also emplaced into the Laurentian margin crust [36–42]. Silurian magmatism is attributed to the Salinic orogeny (445 Ma to 422 Ma) which corresponds to the subduction and accretion of the leading (western) edge of Gander terrane to Laurentia. It is thought that the Salinic orogeny produced older subduction-related plutons and younger postcollisional (back-arc) plutons via slab detachment after the cessation of subduction [40–43]. The Salinic arc magmatic rocks are located in northern and central New Brunswick and western and central Newfoundland proximal to and within the Laurentian margin, but there are a number of postcollisional plutons, batholiths, and intrusions that were emplaced closer to and along the boundary region of the Gander terrane and the Avalon terrane [38, 42, 44–46].
Silurian 206Pb/238U ages of intrusive rocks (mafic to felsic) are reported from the Avalon terrane of Newfoundland ( 4 Ma), Gander terrane of southern New Brunswick ( Ma; Ma, Ma, Ma, and Ma; and Ma; Ma), and coastal Maine ( Ma; Ma; Ma; Ma and Ma; Ma; and Ma). Early Ludlovian (427.4 Ma to 425.6 Ma) K-bentonite beds were identified at Arisaig (Avalon terrane, Nova Scotia) and are thought to be regionally correlative with ash beds from Gotland (Sweden) and those of the British Isles . Similar aged rocks are found within the Meguma terrane of Western Nova Scotia ( Ma; Ma; and Ma) as well [44–46, 48–58].
Plutonism in the region continued from the Late Silurian to Early Devonian and may be related to either the early stages of the Acadian orogeny (420 Ma to 410 Ma) or postcollisional crustal relaxation. The plutons are found across the trailing (eastern) edge regions of the Gander terrane and the boundary region between the Gander and Avalon terranes [38, 45, 46, 48, 53, 59–63]. Subsequent to the Avalon-Laurentia collision (435 Ma to 430 Ma), the Meguma terrane was thrust over Avalonia (~400 Ma) and both terranes were intruded by middle to late Devonian mafic and silicic plutonic rocks [33, 45, 60, 64–67]. Consequently, the precise terrane correlation of the Blair River inlier rocks is difficult to discern based on just their ages as they are more likely to be related to either the Avalon terrane or the Gander terrane than Laurentia. On the basis of whole-rock Nd isotope model ages and Pb isotopes of rocks from the Blair River inlier, including ~580 Ma mafic dykes, Keppie et al.  suggests that they have more in common with Avalonian crust, specifically Oaxacian crust (a proxy for Amazonia), than Laurentia [9, 10, 15, 16, 33, 46, 48, 68]. Nevertheless, the new age dates indicate that they are contemporaneous with the period of postcollisional magmatism that developed along the boundary between the Gander terrane and Avalon terrane.
An examination of the inherited zircons from the anorthositic rocks may provide additional support for the specific terrane correlation as they indicate the age of the country rock through which the parental magma interacted with during emplacement. The majority of inherited zircons have ages of ~440 Ma and~ 465 Ma with individual ages of ~500 Ma and~ 660 Ma and a slightly discordant age of ~790 Ma (Tables DR1 and DR2). Cambrian-Ordovician ages are common in Avalonian sedimentary and metasedimentary rocks but are also consistent with plutonic rocks in the Avalon and Gander terranes [34, 67, 69, 70]. The oldest concordant inherited zircon from the anorthosite dyke has a 206Pb/238U age of Ma (1σ) which is more indicative of the Avalon terrane than the Gander terrane . Zircons with ~650 Ma ages are identified in the Arisaig Group of the Antigonish Highlands in the Avalon terrane, metasedimentary granulites from the structural basement of Meguma, sedimentary rocks from eastern Newfoundland (Avalonia), and metavolcanic schists from Creignish Hills (Avalonia) of southwestern Cape Breton [67, 69, 71, 72]. Sedimentary rocks of southwest New Brunswick and coastal Maine have a few detrital zircons with ages of ~650 Ma (12 zircons) but four of the six units are south of the Turtle Head Fault . The Turtle Head Fault was considered to be the boundary between Gander and Avalon terrane although this interpretation seems to have fallen out of favour [6, 57]. Presently, the zircon inheritance suggests that the rocks of the Blair River inlier are more likely correlative to the Avalon terrane.
6.2. Postcollisional Origin of the Fox Back Ridge Intrusion?
The metagabbro is chemically similar to within-plate alkali basalt with LREE-enriched chondrite-normalized patterns but dissimilar to gabbroic rocks of the Red River anorthosite suite (Figures 6(d), 7(b), and 8(a)–8(c)). It appears that the rocks were either unaffected or minimally affected by crustal contamination as they do not have negative primitive mantle-normalized Nb-Ta anomalies and their Nb/U (> 24) and Th/NbPM (≤ 1.1) ratios do not trend toward average crustal values (Figure 8(d)). Furthermore, it does not appear that fractional crystallization of plagioclase or Fe-Ti oxide minerals played a role in their formation as the Eu/Eu values are ~1 and they do not have negative anomalies of Sr, Ba, or Ti (Figure 7(a)). This suggests, if early fractional crystallization occurred and produced the metagabbro, that only mafic silicate minerals (olivine, pyroxene) fractionated. Therefore, it is likely that the chondrite-normalized rare earth element patterns are primarily the result of partial melting. The high Tb/YbN (1.8-1.9) and Sm/YbPM (4.2-4.5) ratios of the metagabbro indicate that garnet was likely a residual phase in the mantle source suggesting that the minimum depth of melting was equivalent to the garnet-spinel transition ( km) zone . Farther to the west, two diorite rocks from the Fox Back Ridge intrusion reported by Miller  are interpreted to be within-plate basalt but also indicate that the intrusion likely differentiated from mafic to intermediate.
The chemistry of the metagabbro is typical of mafic rocks from within-plate (Figure 8(b)) settings but the regional tectonic setting just prior (>430 Ma) to emplacement was likely one of compression (i.e., Salinic orogeny), and thus the gabbros should be similar to volcanic-arc basalt or more likely demonstrate subduction-related geochemical signatures [43, 62]. The tectonomagmatic dichotomy is difficult to reconcile. However, compositionally similar rocks to the metagabbro are known to occur in postcollisional settings. The Late Ediacaran (~580 Ma) Kekem (Cameroon) and Guéra Massif (Chad) gabbros have low SiO2 (45.7 to 51.0 wt%) and high TiO2 and P2O5 ( to 3.6 wt%; to 2.0 wt%) contents (Figures 8(b) and 8(c)). They were emplaced 15-25 million years after the terminal collision (600-590 Ma) between the Saharan Metacraton and the Congo-Saõ Francisco Craton. Both the Kekam and Guéra gabbros are considered to be derived from a subduction-modified lithospheric mantle source from the spinel-garnet transition zone [75, 76].
The Late Ediacaran tectonomagmatic setting of Central Africa is broadly similar to that of the Blair River inlier before and during emplacement of the metagabbro. The collision between the Avalon terrane and Laurentia was likely preceded by westward subduction toward Laurentian and the formation of Late Silurian arc and later back-arc (post-collisional) magmatic rocks throughout Newfoundland, New Brunswick, and Nova Scotia [42, 43]. Collision and crustal thickening occurred before 428 Ma which could account for subsequent mantle melting within the garnet stability field to produce the metagabbro. Keppie et al.  suggests that the fault and shear zone pattern of NW Cape Breton is indicative of a positive flower structure and that the Blair River inlier was rapidly exhumed during the Early Devonian (415-410 Ma) probably due to sinistral reverse motion along the Cape Ray suture that is related to sinistral transpression along the Dover fault in Newfoundland. Furthermore, the metagabbro is contemporaneous with the emplacement of the ‘A-type’ granites of the Saint George Batholith of southern New Brunswick. The Saint George Batholith (Welsford, Utopia, and Jake Lee Mountain plutons and Bocabec gabbro) intruded rocks of the southern Gander terrane and is very close to the boundary with Avalonia . In other words, it is possible that the gabbro was emplaced during a period of postcollisional crustal relaxation after the Avalon-Gander collision but prior to Early Devonian sinistral fault motion in Newfoundland.
6.3. Petrogenetic Model of the Anorthositic Rocks
Anorthosite is an uncommon plutonic rock that is almost entirely (≥ 90%) comprised of plagioclase feldspar [77, 78]. There are at least six different varieties of anorthosite identified by Ashwal [77, 79] that include: Archean anorthosite, Proterozoic massif anorthosite, layered intrusion anorthosite, oceanic anorthosite, anorthosite inclusions, and extraterrestrial anorthosite. The most abundant type of terrestrial anorthosite is the Proterozoic massif-type anorthosite [80, 81]. Proterozoic massif-type anorthosite is associated with large (up to ~18 000 km2) composite batholiths that range in age from ~2.6 Ga to ~0.5 Ga. The batholiths also include silicic rocks (mangerite, charnockite, and granitoids) that are coeval but not petrogenetically related to the anorthosite [81, 82]. The plagioclase compositions are commonly between An70 and An30 although some anorthosite massifs may have a restricted range (An45±3) whereas others have a wide range . The formation of massif-type anorthosite is not entirely understood but is ultimately related to efficient segregation and accumulation of plagioclase from a mafic parental magma [77, 80–83].
In order to test the possible petrological relationship between the metagabbro and anorthositic rocks, we apply fractional crystallization modeling using MELTS [84, 85]. For the model, we use a parental magma composition equal to the metagabbro. The precise magmatic conditions of the system at the time of emplacement are unknown; however, given the likely collisional/post-collisional tectonic setting , we based the modeling conditions (H2O content and ƒO2) on such settings (Table 3). We use a relatively moderate water content (~2.0 wt%) for the parental magma and a relative oxidation state equal to the Ni-NiO buffer (ΔFMQ +0.7). The moderate water content is selected for two reasons: (1) the rock has high P2O5 content that suggests it was derived from a volatile-rich system and (2) the rock composition indicates it already experienced fractionation of mafic silicate minerals (olivine, pyroxene) that likely increased the bulk water content of the residual magma. The choice of the redox conditions is based on the fact that mafic magmas associated with collisional settings tend to be +1 to +4 log units (ΔFMQ +1 to +4) above that of mid-ocean ridge basalt although our choice of the NNO (ΔFMQ +0.7) buffer is not anomalously high [86, 87]. The pressure (0.2 GPa) of the model is equal to a depth associated with the lower upper crust or upper middle crust.
The results of the fractionation modeling are summarized in Table DR3 and shows that the system is dominated by feldspar (55.6%), pyroxene (10.9%), and spinel (11.9%) fractionation with ~10% residual liquid remaining. The residual liquid will not match the bulk composition of the anorthositic rocks given that they are ‘cumulate’ rocks, but we selected the mineral compositions generated by the model that are within the compositional range of the minerals of the Red River anorthosites presented by Miller  in order to ‘reconstruct’ the rocks by mass balance. Tables 4 and 5 summarize the results of mass balance and show that the modeled mineral compositions derived from the metagabbro can yield compositions that are very similar to the anorthosite dykes and main body anorthosite. Miller  observed that the anorthositic rocks are 90% to 97% plagioclase with a compositional range from An43-50 to An50-56 in the least altered samples. The amount of plagioclase required in both mass balance calculations are 94.2% and 90.2%, respectively. The compositions of the feldspars used in the calculations are An49 for calculation 1 (Table 4) and An29 for calculation 2 (Table 5). The lower An content reflects the higher SiO2 and Na2O contents and lower Al2O3 content of the ‘sodic’ anorthosite but also the lowest measured An component (An29) observed by Miller . The mafic minerals observed in the anorthosite include biotite, chlorite, epidote, Fe-Ti oxide minerals with recrystallized orthopyroxene, and altered clinopyroxene . In addition to 94.2% plagioclase (An49) for calculation 1, 1.7% titanomagnetite, 4.0% orthopyroxene (Mg# 70), and 0.1% apatite can reproduce the main body anorthosite composition. In contrast, calculation 2 requires 0.4% titanomagnetite, 0.17% apatite, and 9.0% biotite to reproduce the bulk composition of the anorthositic dykes.
The modeled mineral compositions we used in our calculation did not necessarily crystallize at the same time. For example, the modeled orthopyroxene composition used in the first mass balance calculation (Table 4) crystallized at a temperature of 950°C, whereas the modeled plagioclase compositions crystallized at 920°C and the Ti-magnetite at 820°C. Conceptually, this suggests that some minerals (e.g., orthopyroxene) within the anorthosite were entrapped within a plagioclase mush whereas others (e.g., Ti-magnetite) may have crystallized from an interstitial liquid within the plagioclase mush.
6.4. Emplacement Model
Undoubtedly, there are a variety of processes that contribute to the formation of anorthosite massifs, but the principal feature is that plagioclase must separate, probably by density stratification, from a mafic parental magma and accumulate [77, 80, 81]. Our modeling results suggest that a parental magma similar to the metagabbro of this study can produce the requisite minerals and their compositions that comprise the anorthositic dykes and the Red River anorthosite. Although we cannot quantify the processes that led to the concentration of the plagioclase, it was likely due to a combination of density stratification and residual melt segregation (c.f., ). Regarding the likely tectonic setting of this system, as mentioned previously, the Avalon terrane collided with the margin of Laurentia during the Late Silurian [43, 61, 88]. The emplacement of older (>428 Ma) silicic plutons in Cape Breton Island and the southern margin of the Gander terrane prior to collision between the accreted margin of Laurentia and Avalonia are similar to I-type (cordilleran) granites, whereas younger (≤ 428 Ma) plutons are similar to A-type (postcollision) granites and within-plate gabbros [33, 35, 46, 55, 89–91]. Given the emplacement age and the timing of A-type granite plutonism in southern New Brunswick (Saint George Batholith), it is likely that the anorthositic rocks of the Blair River inlier were emplaced at the same postcollisional setting, which may be an important environment for the creation of some Proterozoic and Phanerozoic massif-type anorthosites [81, 92, 93].
The anorthositic and mafic rocks collected from the southern Blair River inlier of Cape Breton are contemporaneous and were emplaced during a period of Late Silurian crustal relaxation. The anorthositic rocks of this study are different whereas one is compositionally similar to labradorite, the other is more similar to andesine although it has an appreciable amount of K2O. The cumulate nature of the rocks is confirmed by their petrography, low concentrations of incompatible trace elements, and their high positive chondrite-normalized Eu anomalies (). The Fox Back Ridge metagabbro is compositionally similar to alkali basalt and was likely generated by partial melting of garnet-bearing peridotite. Fractional crystallization modeling coupled with mass balance calculations indicates that the anorthositic rocks may be derived via density and melt segregation from a mafic parental magma similar in composition to the metagabbro of this study. It is likely that the emplacement of the Red River anorthosite occurred at a postcollision setting within the Avalon terrane and that the Laurentian boundary is located farther to the west.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
We thank Sarah Roeske for editorial handling and comments J.S. Daly and an anonymous reviewer that helped to improve this manuscript. We are grateful to Carol Chuang and Terri Tang for their laboratory assistance. This research was supported by a funding provided by Ministry of Science and Technology (MOST) (Taiwan) grant 107-2628-M-003-003 to JGS and NSERC Discovery grant to JD.