Abstract
Avalonian rocks in northern mainland Nova Scotia are characterized by voluminous 640–600 Ma calc-alkalic to tholeiitic mafic to felsic magmas produced in a volcanic arc. However, after the cessation of arc activity, repeated episodes of felsic magmatism between ca. 580 Ma and 350 Ma are dominated by A-type geochemical characteristics. Sm–Nd isotopic data, combined with zircon saturation temperature estimates, indicate that these magmas were formed by high temperature (800–1050 °C) melting of the same anhydrous crustal source. Regional tectonic considerations indicate that A-type felsic magmatism was produced (1) at 580 Ma in a San Andreas-type strike slip setting, (2) at 495 Ma as Avalonia rifted off Gondwana, (3) at 465 and 455 in an ensialic island arc environment and (4) at 360–350 Ma during post-collisional, intra-continental strike-slip activity as Avalonia was translated dextrally along the Laurentian margin. These results attest to the importance of crustal source, rather than tectonic setting, in the generation of these A-type magmas and are an example of how additional insights are provided by comparing the geochemical and isotopic characteristics of igneous suites of different ages within the same terrane. They also suggest that the shallow crustal rocks in northern mainland Nova Scotia were not significantly detached from their lower crustal source between ca. 620 Ma and 350 Ma, a time interval that includes the separation of Avalonia from Gondwana, its drift and accretion to Laurentia as well as post-accretionary strike-slip displacement.
Similar content being viewed by others
Introduction
A-type granites are classified on the basis of geochemical characteristics such as higher SiO2, FeO/FeO+MgO, Ga/Al, Zr, Nb, Y and Ta than either I-type or S-type granites (Loiselle and Wones 1979; Collins et al. 1982; Whalen et al. 1987; Frost et al. 2001). Their tectonic setting and petrogenesis are controversial. Although their high Nb, Y and Ta means that they plot as within-plate granites on tectonic discrimination diagrams (e.g. Pearce et al. 1984), it is evident that they occur in a wide variety of tectonic settings (Whalen et al. 1987; Bonin 2007).
On the one hand, A-type granites are thought to be the result of anhydrous partial melting of a halogen-enriched dry, granulitic lower crustal residue formed after extraction of a felsic magma (Collins et al. 1982; Clemens et al. 1986; Whalen et al. 1987; Landenberger and Collins 1996). The heat required to induce partial melting may arise from rising mantle-derived mafic magma during lithospheric extension. Similarly, Creaser et al. (1991) proposed that some A-type magmas may be derived from melting of tonalitic to granodioritic source rocks. However, Bonin (2007) maintained that leucosomes of A-type composition have not been documented in migmatitic terranes and that no experiments have produced A-type magmas by partial melting of the crust. Instead, Bonin (2007) proposed that they are derived by fractionation either of alkaline mafic magmas or from hybrids produced by mixing between coeval mafic and felsic magmas. Martin (2006) claimed that A-type granites were formed in a source region fertilized by hydrothermal activity above a rising mantle plume and that such events may also be connected to the formation nepheline syenite and carbonatite.
Eby (1990, 1992) showed that granitic magmas classified as A-type by their high Ga/Al also have wide ranges in Nb/Y, a ratio which is often used as a proxy for the degree of alkalinity (Winchester and Floyd 1977; Pearce 1996). On this basis, Eby (1990, 1992) defined two sub-types: (1) A1 which has high Nb/Y (~0.8), may be derived from fractionation of alkalic mafic magma derived from enriched mantle sources, and is emplaced in continental rifts or intra-plate settings; and (2) A2 which has lower Nb/Y, represents magmas derived from continental crust or underplated crust that has been through a previous cycle of magmatism and commonly occurs in a late to post-tectonic environment at the end of a long period of magmatism.
Most petrological studies of A-type magmas are confined to individual suites with a limited range in age. In this paper, we propose that important insights into the genesis of these magmas and their sources can be gained by focusing on terranes where such magmatism occurred in repeated episodes over a long period of time. Four episodes of A-type magmatism have been recognized in Avalonia of northern mainland Nova Scotia, between the late Neoproterozoic and the Late Devonian and the geochemistry of each of these suites has been described individually (Murphy et al.1996a, 2012; Piper et al. 1999; Archibald et al. 2013; Escarraga et al. 2012). Here we show that these suites, when considered collectively, imply that A-type magmas can be generated in different tectonic environments and that they display geochemical and isotopic patterns indicating a remarkable degree of inheritance from an older (Neoproterozoic) crustal source. We propose that when the crustal source becomes dehydrated, a terrane will repeatedly produce A-type magmas irrespective of the tectonic environment. More generally, this study illustrates the importance of comparing the geochemical and isotopic characteristics of igneous suites with different ages in the same terrane in order to constrain the extent to which their signatures are inherited from a common source.
Regional tectonic setting
Regional tectono-stratigraphic comparisons, together with faunal and palaeomagnetic data, suggest that Avalonia is one of several peri-Gondwanan terranes, that were located along the northern periphery of Gondwana in the Late Neoproterozoic–Early Cambrian (Fig. 1). Five important phases of magmatism have been identified in Avalonia (Fig. 2). The oldest phase (760–640 Ma) produced “proto-Avalonia” as a series of oceanic arcs, some possibly underlain by thin slivers of Baltica crust, that formed within the peri-Rodinian (Mirovoi) Ocean and accreted to the northern margin of Gondwana between 660 and 640 Ma (Murphy et al. 2000, 2008; Henderson et al. 2016).
After accretion, subduction zones stepped outboard, producing the second and main tectonic phase (640–570 Ma) of arc-related magmatism and basin formation that characterizes the Avalonia and related peri-Gondwanan terranes (Pe-Piper and Piper 1989; Murphy et al. 1990). Arc magmatism was terminated diachronously between 600 and 540 Ma by the propagation of a San Andreas-style transform fault (Murphy and Nance 1989; Nance et al. 1991), producing the third phase of magmatism during an arc to platform transition, i.e. after the cessation of arc magmatism, but before deposition of Cambrian–Ordovician platformal successions. This strike- slip activity was accompanied by localized intra-plate magmatism and coeval basin development, possibly analogous to the Basin and Range setting in the western United States (e.g. Nance et al. 2002). As the termination of arc magmatism was not due to continental collisions, Cambrian platformal successions were located along the northern margin of Gondwana facing the Mirovoi Ocean.
A combination of paleomagnetic, geochronological and geochemical data indicate that Avalonia was one of several peri-Gondwanan terranes that drifted from the Gondwanan margin in the Late Cambrian or Early Ordovician as the Rheic Ocean opened (Cocks and Torsvik 2002; Stampfli and Borel 2002; Landing 2005; van Staal et al. 2009, 2012; Fig. 1). The next phase of magmatism began with the separation of Avalonia from Gondwana which occurred by ca. 485 Ma (Prigmore et al. 1997; Sanchez-Garcia et al. 2003; Keppie et al. 2006; Murphy et al. 2006). By 460 Ma, Avalonia was about 2000 km north of the Gondwanan margin (Hamilton and Murphy 2004; Cocks and Torsvik 2002) and defined both the northern margin of the Rheic Ocean and the southern flank of the Iapetus Ocean. Magmatism in this time interval is attributed to ensialic arc activity, analogous to northern New Zealand (Murphy et al. 2012).
Avalonia collided with Baltica by the early Silurian and with Laurentia by the early Devonian (e.g. Keppie 1985; 1989; Murphy et al. 1996b; van Staal et al. 1998, 2009, 2012). Between the Late Silurian and Late Carboniferous, Avalonia was located along the northern margin of a contracting Rheic Ocean, the closure of which eventually resulted in the amalgamation of Pangea (Fig. 1; Murphy et al. 2006, 2008). By the middle Devonian, the Meguma terrane was located outboard of Avalonia, and Late Devonian–Early Carboniferous tectonic evolution was dominated by intra-continental dextral transcurrent motion along the boundary between these two terranes, producing the youngest phase of magmatism which is predominantly bimodal (Keppie 1982; Pe-Piper and Piper 1998; Murphy et al. 2011).
Magmatism in northern mainland Nova Scotia
Avalonian rocks in northern mainland Nova Scotia are exposed in the Antigonish Highlands and the Cobequid Highlands (Fig. 3). Geological syntheses indicate that the geology of these highland regions was contiguous from the Late Neoproterozoic until the Late Carboniferous when the Stellarton Basin formed between them (e.g. Murphy et al. 1992; Waldron 2004).
Neoproterozoic arc magmas The oldest rocks occur are a suite of ca. 750–730 Ma granitic to tonalitic gneisses in the eastern Cobequid highlands (Doig et al. 1993; MacHattie et al. 2014). These rocks are interpreted as local representatives of proto-Avalonia that formed outboard of Gondwana within the peri-Rodinian ocean (e.g. Murphy et al. 2000; Henderson et al. 2016).
Both highland regions are predominantly underlain by arc-related ca. 630–608 Ma low-grade volcanic and sedimentary rocks of the Georgeville (Antigonish) and Jeffers (Cobequids) groups that are intruded by ca. 607–600 Ma syn- to post-kinematic mafic to felsic plutonic complexes (Donohoe and Wallace 1982; Keppie et al. 1990; Murphy et al. 1990, 1991, 1996a; Pe-Piper and Piper 1989, 2002; Pe-Piper et al. 1996; MacHattie and White 2015). These igneous rocks are typical of Late Neoproterozoic voluminous arc magmatism that characterizes Avalonia. In both localities, felsic volcanic rocks are interbedded calc-alkalic and tholeiitic mafic to intermediate flows as well as turbidites, thought to reflect local extension and basin generation within the arc.
Published geochemical data from the older arc suite (the ca. 734 ± 2 Ma Economy River orthogneiss; Doig et al. 1993) indicate the rocks are intermediate to felsic in composition, with SiO2 contents ranging from 59.0 to 65.0 wt%, and Fe2O3/MgO from 1.5 to 2.6. The Late Neoproterozoic felsic rocks are more silicic with SiO2 ranging from 62.1 to 76.2 wt% and Fe2O3/MgO from 1.8 to 2.5. These are values typical of calc-alkalic intermediate to felsic rocks. The Economy River and Late Neoproterozoic sequences are both characterized by moderate enrichment in large-ion lithophile elements (K, Rb, Ba, and Th) relative to high-field-strength elements such as Nb, Zr and Y, a signature typical of volcanic arc environments (Pearce et al. 1984). In comparison to typical A-type granites, all these sequences contain significantly lower abundances of Nb (<20 ppm), Y (<26 ppm), Ga (<20 ppm) and low ΣREE. Sm–Nd isotopic data for these sequences are typical of Avalonian crust, with TDM model ages ranging from 0.8 to 1.2 Ga (Nance and Murphy 1996; Murphy et al. 2000, 2008) and ε tNd values ranging from +1.29 to +4.09 (t = 734 Ma) for the Economy River Gneiss and from +0.1 to +3.0 (t = 610 Ma) for the Late Neoproterozoic sequences. These Late Neoproterozoic sequences are the youngest phase of calc-alkalic magmatism to have affected Avalonia in northern mainland Nova Scotia.
Ediacaran The oldest rocks in this region with an A-type chemistry belong to the ca. 580 Ma post-tectonic Georgeville Pluton, which intrudes the lower greenschist to sub-greenschist Georgeville Group in the northernmost Antigonish Highlands adjacent to the Hollow Fault (Murphy et al. 1991, 1998). This granite is characterized by high SiO2, Th, Nb, Y, Zr, Ta and Ga/Al; very low CaO, TiO2, MgO, FeO and MnO and by positively sloped rare earth element (REE) profiles generated by extreme light REE depletion (Figs. 4, 5, 6, 7). Nb-Ga-Th and Ce abundances straddle the divide between A1 and A2-type granites.
Compared to both older and younger felsic rocks, ε tNd values of the Georgeville pluton, are significantly lower (ranging from −1.0 to −7.25, Fig. 8), whereas Sm/Nd (0.6–1.5) and 147Sm/144Nd (0.3–0.6) ratios are significantly higher. In addition, calculated depleted mantle (TDM) model ages (−266 to + 231 Ma) are much younger than the emplacement age, and, therefore, have no geological meaning. Anderson et al. (2008) attributed both the anomalous LREE depletion and low Nd isotopic values of the granite to subsolidus re-equilibration and hydrothermal alteration of the granite, including the metamict zircon.
Early Ordovician Voluminous A-type magmas occur at various times in the Ordovician in the Antigonish Highlands. An early Ordovician plutonic complex has recently been documented in the Cobequid Highlands (MacHattie et al. 2014), but published geochemical data are unavailable. In the central and southern Antigonish Highlands, A-type granitoid magmas, with zircon U–Pb age dates ranging from 495 Ma to 465 Ma, are part of the bimodal West Barneys River Plutonic Suite, a shallow crustal (sub-greenschist) composite intrusive body that also includes several plutons, many of which contain coeval alkalic gabbros (Escarraga et al. 2012; Archibald et al. 2013). These studies show that the granitoids are relatively enriched in Na2O, K2O, FeOt, Ga/Al, Zr, Y and REE, low in Al2O3 and CaO (Fig. 4, 5, 6, 7) and are classified as ferroan A-type granitoids as described by Frost et al. (2001) and Frost and Frost (2011). The granite samples have similar ε tNd (t = 470 Ma) values ranging from +2.4 to +4.9 (Fig. 8), with T DM model ages ranging from 0.8 and 1.1 Ga. The geochemistry of the coeval gabbroic rocks is transitional from alkalic to tholeiitic and is typical of a within-plate setting (Archibald et al. 2013) and their ε tNd (t = 470 Ma) values ranging from +0.8 to +4.7 with TDM model ages between 0.7 and 1.1 Ga.
In the northern Antigonish Highlands, a sample from the felsic volcanic rocks of the bimodal Dunn Point Formation (DPF) yields a zircon U/Pb age of 460 ± 3.4 Ma (Hamilton and Murphy 2004). Mafic volcanic rocks in the formation have continental tholeiitic, within-plate affinities (Keppie et al. 1979; Murphy and Dostal 2007; Murphy et al. 2012). A sample from the conformably overlying McGillivray Brook Formation (MBF) yields a 454 ± 0.7 Ma age (Murphy et al. 2012). The MBF is dominated by felsic tuffs and flows. Although their chemistry is quite distinct from one another, the felsic rocks of the DPF and MBF both resemble A-type (alkalic) SiO2-rich magmas (Murphy et al. 2012). The MBF rocks are more highly fractionated than DPF felsic rocks and have higher SiO2 (76.0–81.8 wt%), lower Al2O3, and very low MgO, CaO, TiO2 and P2O5. They are also characterized by significantly higher concentrations of Zr (745–1965 ppm), Y (65–213 ppm), Nb (57–185 ppm), Ta (2.1–6.0 ppm), Yb (5.6–16. 7 ppm) and Ga/Al (Figs. 4, 5, 6, 7).
ε tNd values for Dunn Point rhyolites range from +2.9 to +3.7 (t = 460 Ma) and overlap with McGillivray Brook Formation felsic rocks (t = 455 Ma) which range from +1.5 to + 3.9 (Fig. 8). The DPF analyses, however, have a much lower range in 147Sm/144Nd (0.128 to 0.148) compared to the MBF (0.140 to 0.348). TDM ages vary from 0.9 to 1.2 Ga for DPF and 0.9 to 1.1 Ga for MBF.
Late Devonian–Early Carboniferous A-type silicic rocks in the Cobequid Highlands include voluminous Late Famennian to Early Tournaisean bimodal volcanic (Fountain Lake Group) and coeval shallow crustal plutonic rocks which are predominantly restricted to the area north of the Rockland Brook Fault (Donohoe and Wallace 1982; Pe-Piper and Piper 2002; MacHattie et al. 2014). Zircon U/Pb (TIMS) ages range from ca. 365 to 354 Ma (Dunning et al. 2002). MacHattie (2011) discovered significant high-field-strength/rare earth element (HFSE/REE) prospects in the vicinity of the granites, an association typical of A-type granites. Regional relationships suggest that these rocks formed in an intra-continental tectonic environment, with much of the magmatism spatially associated with coeval strike-slip motion along major faults (Pe-Piper et al. 1998; Koukouvelas et al. 2002).
In the eastern highlands, the Fountain Lake Group can be subdivided into a ca. 4000 m sequence of predominantly rhyolitic volcanic rocks (Byers Brook Formation) conformably overlain by ca. 1500 m of predominantly basalt (Diamond Brook Formation). According to Piper et al. (1999), the rhyolites have high abundances of petrogenetic-discriminating trace elements (e.g., Nb, Y, Zr, Ga and Rb) that are typical of A-type magmas and can be subdivided into a lower member dominated by relatively low-Zr rhyolites and an upper member dominated by relatively high-Zr rhyolites. The geochemistry of the low-Zr rhyolites is similar to that of the coeval plutons and is characterized by significantly lower Zr (<600 ppm), low εNd (<+1.0), Ga/Al and La/Th compared to the high-Zr rhyolites in which εNd is >+ 3.0 and has Zr > 1000 ppm (Figs. 4, 5, 6, 7, 8). The overlying basalts have within-plate, continental tholeiitic affinities, characterized by a high Ni, Cr and Mg number, and relatively high ε tNd (>+ 4), which have been interpreted to reflect high potential mantle temperatures influenced by the presence of a mantle plume (Dessureau et al. 2000).
Late Devonian–Early Carboniferous plutonic rocks are predominantly granite and gabbro in composition, with field evidence of magma mixing and mingling between co-existing felsic and mafic magmas that generated intermediate compositions (Pe-Piper et al. 1996). U–Pb dating of plutonic rocks shows that they were emplaced in the upper crust between ca. 363 Ma and 350 Ma (Doig et al. 1996; Pe-Piper et al. 2004). Structural studies show that major shear zones (including the Rockland Brook Fault) underwent dextral strike-slip motion at this time (Miller et al. 1995). Local transtensional zones within this fault system facilitated emplacement of the magmas into the upper crust (Pe-Piper et al. 1998; Koukouvelas et al. 2002). Most plutons are composite, displaying evidence of repeated injections of felsic and mafic magma with mixing and mingling between phases (Pe-Piper and Piper 2002).
The granite plutons are subalkaline to alkaline and all have a geochemical signature typical of A-type granites (Pe-Piper et al. 1991; Papoutsa and Pe-Piper 2014; Papoutsa et al. 2016). According to Pe-Piper and Piper (1998) and Koukouvelas et al. (2002), the A-type geochemistry combined with Nd and Pb isotopic data indicates that the granites were derived by anhydrous melting of Neoproterozoic lower crust. Most of the granite bodies have chemical compositions similar to that of the low-Zr rhyolite in the Fountain Lake Group. However, the Whirley Brook intrusion has a composition similar to the coeval Fountain Lake high-Zr rhyolite (Piper et al. 1993). The youngest intrusive bodies include a voluminous 354–350 Ma gabbro in the southwestern portion of the Wentworth Pluton, and some granite bodies are attributed to partial melting of older granite during intrusion of the mafic magma (Koukouvelas et al. 2002). All mafic intrusive rocks have a continental tholeiitic, within-plate signature, indistinguishable from the mafic rocks of the Diamond Brook Formation.
Zircon saturation thermometry estimates
The temperature at which zircon saturation occurs in a felsic melt can be used to discriminate between felsic magmas with different petrogenetic as well as thermal histories (King et al. 1997; Miller et al. 2003; Schmitt and Simon 2004; Shellnutt et al. 2011; Dostal et al. 2015; Mahdy et al. 2015; Xia et al. 2016). In general, crustally derived calc-alkalic felsic magmas are typically generated by water-fluxed melting and, therefore, form at lower temperatures (<800 °C) than A-type felsic magmas which are generated by relatively high temperature (>800 °C) anhydrous melting (Weinberg and Hasalová 2015; Collins et al. 2016).
The temperature of zircon saturation (T Zr) can be estimated by relating the concentration of Zr to the bulk composition of the magma (Watson and Harrison 1983; Hanchar and Watson 2003; Boehnke et al. 2013). This relationship is expressed by
where the M value is defined by [Na + K + 2Ca]/[Al*Si], D Zr is the distribution coefficient of Zr in zircon (500 000) and T Zr is the temperature in degrees Kelvin. The experimental results are based on compositions that are metaluminous to peraluminous and have M values between 1.3 and 1.9 (Watson and Harrison 1983).
The Economy River Gneiss is not sufficiently silicic to use this method. Only felsic rock samples that yielded M values between 1.3 and 1.9 were considered meaningful for this study as it was the range for the original experimental results on which Eq. 1 is based (Hanchar and Watson 2003). The data are shown in Table 1 and are plotted on Fig. 9. The large temperature range in some suites could be related to their evolution in the shallow crust (fractionation, host rock assimilation, alteration) for Zr saturation as the magma crystallizes. As such, only the maximum TZr is considered diagnostic of the temperature at which melting commences.
The arc-related ca. 618 Ma Georgeville Group felsic volcanic rocks yield the lowest T Zrn (oC) values, ranging from 658 to 780 °C with one outlier at 526 °C. The 580 Ma Georgeville granites yield a relatively restricted range, but higher T Zr (°C) values (790 to 840 °C) compared to the Georgeville Group. The ca. 495 Ma West Barneys River silicic rocks produced a bimodal distribution with two T Zr (°C) values of 850 and 860 °C and an additional two of 720 and 620 °C. The ca. 460 Ma Dunn Point felsic rocks have T Zr (°C) values ranging from 860 to 875 °C and the ca. 455 Ma McGillivray Brook rhyolite produced only one T Zr (°C) value (1050 °C). The silicic rocks from ca. 360 to 350 Ma Fountain Lake Group yield T Zr (°C) values ranging from 770 to 1020 °C.
With the exception of the calc-alkalic Georgeville Group felsic rocks, most of the T Zr (°C) values for the plutons and volcanic rocks (arranged in chronological order in Fig. 9) are greater than 800 °C, the temperature used by Collins et al. (2016) to distinguish between predominantly calc-alkalic and predominantly A-type granitic magmas. The maximum TZr value for each suite is considered the best indicator of the melting temperature of the source rock and appears to increase from the Late Neoproterozoic to the Devonian-Carboniferous. These results, together with the geochemical and Sm–Nd isotopic evidence, are consistent with the interpretation that the felsic magmas were each derived by melting the lower crust and the temperatures correspond to those expected for the formation of a lower crust hot zone associated with crustal melting (Huppert and Sparks 1988; Murphy et al. 1998; Annen and Sparks 2002; Annen et al. 2006, 2015; Annen 2009). The high T Zr (°C) values of the 580–350 Ma felsic magmas may reflect episodes of dry partial melting of lower crustal rocks.
Summary and discussion
Geochemical and Sm–Nd isotopic data show that Avalonian magmatism in northern mainland Nova Scotia produced repeated episodes of A-type felsic magmatism starting at ca. 580 Ma and continued into the Carboniferous. With the exception of the highly altered Georgeville Granite, all felsic suites plot within the Nd isotopic envelope for crustally derived rocks in West Avalonia, supporting the interpretation that the magmas were derived by anatexis of the crust.
The T Zr (°C) estimates show that the A-type silicic rocks are recording repeated episodes of a relatively hot melting of anhydrous West Avalonian crust between the Ediacaran to Early Carboniferous related to the repeated injection of basaltic magma (Annen et al. 2006). Each magmatic pulse represents a separate tectonomagmatic event but it appears that the same crust is being reworked as the silicic rocks have ε tNd values from ~0 to ~+5 and TDM ages from 0.7 to 1.1 Ga, which are typical of Avalonian basement (Murphy et al. 1998, 2012; Piper et al. 1999; Murphy and Nance 2002; Archibald et al. 2013; Shellnutt and Dostal 2015).
Regional tectonic considerations suggest that A-type felsic magmas were generated in different tectonic environments. The 580 Ma Georgeville Granite was intruded within a continental transform fault system, about 20 Ma after the cessation of subduction. The ca. 490 Ma West Barneys River Plutonic Suite was intruded coevally with respect to the opening of the Rheic Ocean and generation of the Avalonian ribbon continent. The felsic volcanic rocks of the ca. 460 Ma Dunn Point and ca. 455 McGillivray Brook formations were extruded when Avalonia was a microcontinent. Although regional syntheses suggest an ensialic island arc environment, the high T of melting indicates that magmatism was not proximal to regions affected by subduction zone fluids. Regionally extensive hot regions behind arcs are common (Hyndman et al. 2005) and the voluminous felsic magmas with A-type composition produced in ensialic island arcs in the southeastern Pacific, including northern New Zealand (Smith et al. 1977), may provide a modern analogy for ca. 460–455 Ma A-type magmatism in Avalonia. The Fountain Lake Group and coeval plutonic rocks were emplaced in an intra-continental, dextral transtensional environment, up to 60 million years after Avalonia had collided with Laurentia.
The common denominator to the generation of these A-type magma suites is their crustal source. With the exception of the Georgeville Granite, whose Sm–Nd isotopic signature was profoundly affected by alteration of U-rich zircon (Anderson et al. 2008), all suites lie within the envelope for crustally derived Avalonian felsic rocks (Murphy et al. 1996a). In the highly fractionated suites, accessory phase fractionation resulted in LREE depletion, thereby rendering T DM model age calculations geologically meaningless for those suites. However, T DM model age calculations for the West Barneys River Plutonic Suite, Dunn Point Formation and Fountain Lake Group are typical of crustally derived Avalonian felsic rocks. Despite the alteration within the Georgeville granites, if a typical Sm/Nd crustal ratio is assumed (Fig. 9), it is also clear that they were derived from older crust, presumably Avalonian basement.
The inferred high-T metamorphic conditions during magma emplacement of the earliest (~580 Ma) Georgeville plutonic suite (max ~840 °C) are typical for biotite breakdown during anhydrous melting of the lower crust (Collins et al. 1982; Whalen et al. 1987; Creaser et al. 1991). A-type felsic magmatism was preceded by ca. 618 Ma subduction-related, water-fluxed melting of Avalonian crust with maximum temperatures consistently <800 °C (see Collins et al. 2016).
By contrast, the younger magmatic phases typically formed at much higher temperatures (>900 °C) under ultra-high temperature (UHT) metamorphic conditions. The higher melting temperatures suggest that the Avalonian crust had lost much (most) of its hydrous mineral component (especially biotite) during the ~580 Ma event. It was, therefore, probably charnockitic in composition before the onset of Ordovician melting (cf., Landenberger and Collins 1996). The UHT conditions during melting were associated with emplacement of alkalic to tholeiitic mafic magmas, and the high TZr melting temperatures reflect the absence of a hydrous, subduction-related mantle component. The elevated εNd values of the felsic rocks indicate melting of the same Avalonian basement source. The loss of hydrous minerals from the crust resulted in the loss of a buffering mechanism to keep crustal temperatures at ~850 °C or less during melting.
Magmas similar in composition to the Dunn Point and McGillivray Brook A-type rhyolites occur in the modern subduction zone systems of the southwest Pacific (Smith et al. 1977) and may serve as a modern analogue for the tectonic setting in the middle Ordovician. The Late Devonian–Early Carboniferous Fountain Lake Group silicic volcanic rocks formed under UHT conditions, with peak temperatures extending to >1000 °C, consistent with the crustal source having been progressively refined by prior melting events so that all hydrous minerals within that crust have been removed. Like the Ordovician events, emplacement of hot, dry mafic magmas were the ultimate cause of the extreme conditions required to melt the dehydrated crust.
This interpretation extends the conclusion of Pe-Piper and Piper (1998) and Koukouvelas et al. (2002) for Devonian-Carboniferous A-type magmas to include those as old as 580 Ma. It suggests that Avalonian crust in mainland Nova Scotia was removed from the effects of subduction-related, hydrous fluxing by 580 Ma, reflecting the transition from arc to transtensional intraplate regime, possibly similar to the present Basin and Range tectonic setting where “A-type” high-silica rhyolites predominate (e.g., Bachmann and Bergantz 2008). It also shows that, irrespective of changing tectonic settings, the thermal conditions during crustal melting were controlled by a pre-conditioning of the crust associated with progressive removal of hydrous minerals during successive melting events, ultimately resulting in crustal melting under UHT conditions.
References
Anderson AJ, Wirth R, Thomas R (2008) The alteration of metamict zircon and its role in remobilization of high field-strength elements in the Georgeville Granite, Nova Scotia. Can Mineral 46:1–18
Annen C (2009) From plutons to magma chambers: thermal constraints on the accumulation of eruptible silicic magma in the upper crust. Earth Planet Sci Lett 284:409–416
Annen C, Sparks RSJ (2002) Effects of repetitive emplacement of basaltic intrusions on thermal evolution and melt generation in the crust. Earth Planet Sci Lett 203:937–955
Annen C, Blundy JD, Sparks RSJ (2006) The genesis of intermediate and silicic magmas in deep crustal hot zones. J Petrol 47:505–539
Annen C, Blundy JD, Leuthold J, Sparks RSJ (2015) Construction and evolution of igneous bodies: towards an integrated perspective of crustal magmatism. Lithos 230:206–221
Archibald DB, Barr SM, Murphy JB, White CE, MacHattie TG, Escarraga EA, Hamilton MA, MacFarlane CRM (2013) Field relationships, petrography, and tectonic setting of the Ordovician West Barneys River Plutonic Suite, southern Antigonish Highlands, Nova Scotia. Can J Earth Sci 50:727–745
Bachmann O, Bergantz GW (2008) Rhyolites and their source mushes across tectonic settings. J Pet 49:2277–2285
Boehnke P, Watson EB, Trail D, Harrison TM, Schmitt AK (2013) Zircon saturation revisited. Chem Geol 351:324–334
Bonin B (2007) A-type granites and related rocks: evolution of a concept, problems and prospects. Lithos 97:1–29
Clemens JD, Holloway JR, White AJR (1986) Origin of an A-type granite: experimental constraints. Am Mineral 71:317–324
Cocks LRM, Torsvik TH (2002) Earth geography from 500 to 400 million years ago: a faunal and palaeomagnetic review. J Geol Soc London 159:631–644
Collins WJ, Beams SD, White AJR, Chappell BW (1982) Nature and origin of A-type granites with particular reference to southeastern Australia. Contrib Mineral Petr 80:189–200
Collins WJ, Huang HQ, Jiang X (2016) Water-fluxed crustal melting produces Cordilleran batholiths. Geology 44:143–146
Creaser RA, Price RC, Wormald RJ (1991) A-type granites revisited: assessment of a residual-source model. Geology 19:163–166
Dessureau G, Piper DJW, Pe-Piper G (2000) Geochemical evolution of earliest Carboniferous continental tholeiitic basalts along a crustal-scale shear zone, southwestern Maritimes basin, eastern Canada. Lithos 50:27–50
Doig R, Murphy JB, Nance RD (1993) Tectonic significance of the Late Proterozoic Economy River Gneiss, Cobequid Highlands, Avalon composite terrane, Nova Scotia. Can J Earth Sci 30:474–479
Doig R, Murphy JB, Pe-Piper G, Piper DJW (1996) U–Pb geochronology of Late Paleozoic plutons, Cobequid Highlands, Nova Scotia, Canada: evidence for Late Devonian emplacement adjacent to the Meguma-Avalon terrane boundary. Geol J 31:179–188
Donohoe HV, Wallace PI (1982) Geological Map of the Cobequid Highlands, Nova Scotia. Scale 1:50,000. Halifax, Nova Scotia, Canada
Dostal J, Owen JV, Shellnutt JG, Keppie JD, Gerel O, Corney R (2015) Petrogenesis of the Triassic Bayan-Ulan alkaline granitic pluton in the North Gobi rift of central Mongolia: implications for the evolution of the Early Mesozoic granitoid magmatism in the Central Asian Orogenic Belt. J Asian Earth Sci 109:50–62
Dunning G, Barr S, Giles PS, McGregor DC, Pe-Piper G, Piper DJW (2002) Chronology of Devonian to early Carboniferous rifting and igneous activity in southern Magdalen Basin based on U–Pb (zircon) dating. Can J Earth Sci 39:1219–1237
Eby GN (1990) The A-type granitoids: a review of their occurrence and chemical characteristics and speculations on their petrogenesis. Lithos 26(115):134
Eby GN (1992) Chemical subdivision of A-type granitoids: petrogenetic and tectonic implications. Geology 20(641):644
Escarraga E, Barr SM, Murphy JB, Hamilton MA (2012) Ordovician A-type plutons in the Antigonish Highlands, Nova Scotia. Can J Earth Sci 49:329–345
Frost CD, Frost BR (2011) On ferroan (A-type) granitoids: their compositional variability and modes of origin. J Petrol 52:39–53
Frost BR, Barnes CG, Collins WJ, Arculus RJ, Ellis DJ, Frost CD (2001) A geochemical classification for granitic rocks. J Petrol 42:2033–2048
Hamilton MA, Murphy JB (2004) Tectonic significance of a Llanvirn age for the Dunn Point volcanic rocks, Avalon terrane, Nova Scotia, Canada: implications for the evolution of the Iapetus and Rheic oceans. Tectonophysics 379:199–209
Hanchar JM, Watson EB (2003) Zircon saturation thermometry. Rev Mineral Geochem 53:89–112
Henderson BJ, Collins WJ, Murphy JB, Gutierrez-Alonso G (2016) Gondwanan basement terranes of the Variscan-Appalachian orogen: Baltican, Saharan and West African hafnium isotopic fingerprints in Avalonia, Iberia and the Armorican Terranes. Tectonophysics 681:278–304
Huppert HE, Sparks RSJ (1988) The generation of granitic magmas by intrusion of basalt into continental crust. J Petrol 29:599–624
Hyndman R, Currie CA, Mazzotti SP (2005) Subduction zone backarcs, mobile belts, and orogenic heat. GSA Today 15(2):4–10
Keppie JD (1982) Tectonic Map of Nova Scotia. Scale 1:500,000. Nova Scotia Department of Mines and Energy, Halifax, Nova Scotia, Canada
Keppie JD (1985) The Appalachian Collage. In: Gee DG, Sturt B (eds) The Caledonide orogen, Scandinavia, and related areas. Wiley, New York, pp 1217–1226
Keppie JD 1989. Northern Appalachian terranes and their accretionary history. In: Dallmeyer RD (ed) Terranes in the circum-Atlantic Paleozoic orogens. Geological Society of America Special Paper 20, pp 159–192
Keppie JD, Dostal J, Zentilli M (1979) Early Silurian volcanic rocks at Arisaig, Nova Scotia. Can J Earth Sci 16:1635–1639
Keppie JD, Dallmeyer RD, Murphy JB (1990) Tectonic implications of 40Ar/39Ar hornblende ages from Late Proterozoic-Cambrian plutons in the Avalon Composite Terrane in Nova Scotia, Canada. Geol Soc Am Bull 102:516–528
Keppie JD, Nance RD, Fernandez-Suarez J, Storey CD, Jeffries TE, Murphy JB (2006) Detrital zircon data from the eastern Mixteca terrane, southern Mexico: evidence for an Ordovician-Mississippian continental rise and a Permo-Triassic clastic wedge adjacent to Oaxaquia. Internat Geol Rev 48:97–111
King PL, White AJR, Chappell BW, Allen CM (1997) Characterization of origin of aluminous A-type granites from the Lachlan fold belt, Southeastern Australia. J Petrol 38:391
Koukouvelas I, Pe-Piper G, Piper DJW (2002) The role of dextral transpressional faulting in the evolution of an early Carboniferous mafic-felsic plutonic and volcanic complex: Cobequid Highlands, Nova Scotia, Canada. Tectonophysics 348:219–246
Landenberger B, Collins WJ (1996) Derivation of A-type granites from a dehydrated charnockitic lower crust: evidence from the Chaelundi Complex, eastern Australia. J Petrol 37:145–170
Landing E (2005) Early Paleozoic Avalon-Gondwana unity: an obituary-response to ‘‘Palaeontological evidence bearing on global Ordovician-Silurian continental reconstructions’’ by RA Fortey, LRM Cocks: discussion. Earth Sci Rev 69:169–175
Loiselle MC, Wones DR (1979) Characteristics and origin of anorogenic granites. Abstr Progr Geol Soc Am 11:468
MacHattie TG (2011) Nature and setting of late Devonian-early carboniferous rare earth element mineralization in the eastern Cobequid Highlands, Nova Scotia. Mineral Resources Branch, Report of Activities 2010. NS Depart Nat Res Rep ME 2011–1:75–92
MacHattie TG, White CE, Beresford V, Reid M (2014) An Update of Bedrock Mapping in the Eastern Cobequid Highlands, Northern Mainland Nova Scotia. In: Mineral Resources Branch, Report of Activities 2013. Nova Scotia Department of Natural Resources, Report ME 2014-001, pp 145–156
Mahdy NM, El Kalioubi BA, Wohlgemuth-Ueberwasser CC, Shalaby MH, El-Afandy AH (2015) Petrogenesis of U- and Mo-bearing A2-type granite of the Gatter batholith in the Arabian Nubian shield, northeastern Desert, Egypt; evidence for the favorability of host rocks for the origin of associated ore deposits. Ore Geol Rev 71:57–81
Martin RF (2006) A-type granites of crustal origin ultimately result from open-system fenitization-type reactions in an extensional environment. Lithos 91:125–136
Miller BV, Nance RD, Murphy JB (1995) Kinematics of the Rockland Brook Fault, Nova Scotia: implications for the interaction of the Meguma and Avalon terranes. J Geodyn 19:253–270
Miller CF, McDowell SM, Mapes RW (2003) Hot and cold granites? Implications of zircon saturation temperatures and preservation of inheritance. Geology 31:529–532
MacHattie TG, White, CE (2015) An Update of Bedrock Mapping in the Western Cobequid Highlands, Northern Mainland Nova Scotia. In: Geoscience and Mines Branch, Report of Activities 2014; NS Dept of Nat Res Rep ME 2015-001, pp 51–61
Murphy JB, Dostal J (2007) Continental mafic magmatism of different ages in the same terrane: constraints on the evolution of an enriched mantle source. Geology 35:335–338
Murphy JB, Nance RD (1989) A model for the evolution of the Avalonian–Cadomian belt. Geology 17:735–738
Murphy JB, Nance RD (2002) Sm-Nd isotopic systematics as tectonic tracers: an example from West Avalonia in the Canadian Appalachians. Earth-Sci Rev 59:77–100
Murphy JB, Keppie JD, Dostal J, Hynes AJ (1990) Late Precambrian Georgeville Group: A volcanic arc rift succession in the Avalon Terrane of Nova Scotia. Geol Soc Lond Spec Publ 51:383–393
Murphy JB, Keppie JD, Hynes AJ (1991) Geology of the Antigonish Highlands. Geol Sur Can Pap 89–10:114
Murphy JB, Pe-Piper G, Keppie JD, Piper DJW (1992) Correlation of Neoproterozoic III sequences in the Avalon Composite Terrane of mainland Nova Scotia: tectonic implications. Atl Geol 28:147–153
Murphy JB, Keppie JD, Dostal J, Cousens BL (1996) Repeated lower crustal melting beneath the Antigonish Highlands, Avalon Composite Terrane, Nova Scotia: Nd isotopic evidence and tectonic implications. In: Nance RD, Thompson MD (eds) Avalonian and related peri-Gondwanan terranes of the circum North Atlantic, vol 304. Geological Society of America Special Paper, pp 109–120
Murphy JB, Keppie JD, Dostal J, Waldron JWF, Cude M-P (1996b) Geochemical and isotopic constraints on the accretion of the Avalonia in the Appalachian-Caledonide orogen: evidence from Early Silurian clastic sequences in Antigonish Highlands, Nova Scotia, Canada. Can J Earth Sci 33:379–388
Murphy JB, Anderson AJ, Archibald DA (1998) Postorogenic alkali feldspar granite and associated pegmatites in West Avalonia: petrology of the Neoproterozoic Georgeville pluton, Antigonish Highlands, Nova Scotia. Can J Earth Sci 35:110–120
Murphy JB, Strachan RA, Nance RD, Parker KD, Fowler MB (2000) Proto-Avalonia: A 1.2 to 1.0 Ga tectonothermal event and constraints for the evolution of Rodinia. Geology 28:1071–1074
Murphy JB, Gutierrez-Alonso G, Damian Nance R, Fernandez-Suarez J, Keppie JD, Quesada C, Strachan RA, Dostal J (2006) Origin of the Rheic Ocean: Rifting along a Neoproterozoic suture? Geology 34(5):325
Murphy JB, McCausland PJA, O’Brien SJ, Pisarevsky S, Hamilton MA (2008) Age, geochemistry and Sm-Nd isotopic signature of the 0.76 Ga Burin Group: compositional equivalent of Avalonian basement? Precambrian Research 165:37–48
Murphy JB, Waldron JWF, Kontak DJ, Pe-Piper G, Piper DJW (2011) Minas Fault Zone: Late Paleozoic history of an intra-continental orogenic transform fault in the Canadian Appalachians. J Struct Geol 33(3):312–328
Murphy JB, Hamilton MA, LeBlanc B (2012) Tectonic significance of Late Ordovician silicic magmatism Avalon terrane, northern Antigonish Highlands, Nova Scotia. Can J Earth Sci 49:346–358
Nance RD, Murphy JB, Strachan RA, D’Lemos RD, Taylor GK (1991) Tectonostratigraphic evolution of the Avalonian–Cadomian belt and the breakup of a late Precambrian Supercontinent: an interpretive review. Precambr Res 53:41–78
Nance RD, Murphy JB (1994) Contrasting basement isotopic signatures and the palinspastic restoration of peripheral orogens: example from the Neoproterozoic Avalonian-Cadomian belt. Geology 22(7):617
Nance RD, Murphy JB (1996) Basement isotopic signatures and Neoproterozoic paleogeography of Avalonian-Cadomian and related terranes in the circum North Atlantic. In: Nance RD, Thompson MD (eds) Avalonian and related peri-Gondwanan terranes of the circum North Atlantic, vol 304. Geological Society of America Special Paper, pp 333–346
Nance RD, Murphy JB, Keppie JD (2002) Cordilleran model for the evolution of Avalonia. Tectonophys 352:11–31
Papoutsa A, Pe-Piper G (2014) Geochemical variation of amphiboles in A-type granites as an indicator of complex magmatic systems: Wentworth pluton, Nova Scotia, Canada. Chem Geol 384:120–134
Papoutsa A, Pe-Piper G, Piper DJW (2016) Systematic mineralogical diversity in A-type granitic intrusions: control of magma source and geological processes. Geol Soc Am Bull 128:487–501
Pearce JA (1996) A user’s guide to basaltic discrimination diagrams. In: Wyman DA (ed.) Trace element geochemistry of volcanic rocks: applications for massive sulphide exploration: Geological Association of Canada Short Course Notes 12, pp 79–113
Pearce JA, Harris NBW, Tindle AG (1984) Trace element discrimination diagrams for the tectonic interpretation of granitic rocks. J Petrol 25:956–983
Pe-Piper G, Piper DJW (1989) The upper Hadrynian Jeffers Group, Cobequid Highlands, Avalon zone of Nova Scotia: a back-arc volcanic complex. Geol Soc Amer Bull 101:364–376
Pe-Piper G, Piper DJW (1998) Geochemical evolution of Devonian-Carboniferous igneous rocks of the Magdalen basin, Eastern Canada: Pb and Nd isotope evidence for mantle and lower crustal sources. Can J Earth Sci 35:201–221
Pe-Piper G, Piper DJW (2002) A synopsis of the geology of the Cobequid Highlands, Nova Scotia. Atlantic Geol 38:145–160
Pe-Piper G, Piper DJW, Clerk SB (1991) Persistent mafic igneous activity in an A-type granite pluton, Cobequid Highlands, Nova Scotia. Can J Earth Sci 28:1058–1072
Pe-Piper G, Piper DJW, Koukouvelas I (1996) The Precambrian plutons of the Cobequid Highlands, Nova Scotia, Canada. In Nance RD, Thompson MD (eds) Avalonian and Related Peri-Gondwanan Terranes of the Circum-North Atlantic: Boulder, Colorado. vol 304. Geological Society of America Special Paper, pp 121–132
Pe-Piper G, Koukouvelas I, Piper DJW (1998) Synkinematic granite emplacement in a shear zone: the Pleasant Hills pluton, Canadian Appalachians. Geol Soc Am Bull 110:523–536
Pe-Piper G, Reynolds PH, Nearing J, Piper DJW (2004) Early Carboniferous deformation and mineralization in the Cobequid shear zone, Nova Scotia: an Ar/ Ar geochronology study. Can J Earth Sci 41(12):1425–1436
Piper DJW, Pe-Piper G, Loncarevic BD (1993) Devonian-Carboniferous igneous intrusions and their deformation, Cobequid Highlands, Nova Scotia. Atl Geol 29(3):219–232
Piper DJW, Dessureau G, Pe-Pier G (1999) Occurrence of Early Carboniferous high-Zr rhyolites, Cobequid Highlands, Nova Scotia: temperature effect of a contemporaneous mafic magma. Can Mineral 37:619–634
Prigmore JK, Butler AJ, Woodcock NH (1997) Rifting during separation of Eastern Avalonia from Gondwana: evidence from subsidence analysis. Geology 25:203–206
Sanchez-Garcia T, Bellido F, Quesada C (2003) Geodynamic setting and geochemical signatures of Cambrian–Ordovician rift-related igneous rocks (Ossa Morena zone, SW Iberia). Tectonophysics 365:233–255
Schmitt AK, Simon JI (2004) Boron isotopic variation in hydrous rhyolitic melts: a case study from Long Valley, California. Contrib Mineral Petr 146:590–605
Shellnutt JG, Dostal J (2015) Granodiorites of the South Mountain Batholith (Nova Scotia, Canada) derived by partial melting of Avalonia granulite rocks beneath the Meguma terrane: implications for the heat source of the Late Devonian granites of the North Appalachians. Tectonophysics 655:206–212
Shellnutt JG, Jahn BM, Zhou MF (2011) Crustal-derived granites in the Panzhihua region, SW China: implications for felsic magmatism in the Emeishan Large Igneous province. Lithos 123:145–157
Smith IEM, Chappell BW, Ward GK, Freeman RS (1977) Peralkaline rhyolites associated with andesitic arcs of the southwest Pacific. Earth Planet Sci Lett 37:230–236
Stampfli GM, Borel GD (2002) A plate tectonic model for the Paleozoic and Mesozoic constrained by dynamic plate boundaries and restored synthetic oceanic isochrones. Earth Planet Sci Lett 196:17–33
van Staal CR, Dewey JF, Mac Niocaill C, McKerrow WS (1998) The Cambrian-Silurian tectonic evolution of the Northern Appalachians and British Caledonides: history of a complex, west and southwest Pacific type segment of Iapetus. In: Blundell D, Scott AC (eds.) Lyell: the past is the key to the present. Geological Society of London Special Publication 143, pp 199–242
van Staal CR, Whalen JB, Valverde-Vaquero P, Zagorevski A, Rogers N (2009) Pre-Carboniferous, episodic accretion-related, orogenesis along the Laurentian margin of the northern Appalachians. In: Murphy JB et al. (eds.) Ancient orogens and modern analogues. Geological Society of London Special Publication 327, pp 271–316
van Staal CR, Barr SM, Murphy JB (2012) Provenance and tectonic evolution of Ganderia: constraints on the evolution of the Iapetus and Rheic oceans. Geology 40:987–990
Waldron JWF (2004) Anatomy and evolution of a pull-apart basin, Stellarton, Nova Scotia. Geol Soc Amer Bull 116:109–127
Watson EB, Harrison TM (1983) Zircon saturation revisited: temperature and composition effects in a variety of crustal magma types. Earth Planet Sci Lett 64:295–304
Weinberg RF, Hasalová P (2015) Water-fluxed melting of the continental crust: a review. Lithos 212–215:158–188
Whalen JB, Currie KL, Chappell BC (1987) A-type granites: geochemical characteristics, discrimination, and petrogenesis. Contrib Mineral Petr 95:407–419
White CE, Archibald DB, MacHattie TG, Escarraga EA 2011. Preliminary geology of the southern Antigonish Highlands, northern mainland Nova Scotia. In: MacDonald DR (ed) Mineral Resources Branch, Report of Activities 2010. Nova Scotia Dept Natural Res Rept ME 2011-1, pp 145–164
Winchester JA, Floyd PA (1977) Geochemical discrimination of different magma series and their differentiation products using immobile elements. Chem Geol 20:325–343
Xia Y, Xu X, Liu L (2016) Transition from adakitic to bimodal magmatism induced by the paleo-Pacific plate subduction and slab rollback beneath SE China: evidence from petrogenesis and tectonic setting of the dike swarms. Lithos 244:182–204
Acknowledgements
JBM acknowledges the continuing support of N.S.E.R.C. and a Hadyn Williams Fellowship at Curtin University. JGS wishes to thank the Ministry of Science and Technology (Taiwan) through grant #102-2628-M-003-001-MY4. WJC acknowledges support from a Curtin Senior Research Fellowship and ARC DP 120104004. We thank Jarda Dostal and an anonymous reviewer for constructive reviews that significantly improved the manuscript.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Murphy, J.B., Shellnutt, J.G. & Collins, W.J. Late Neoproterozoic to Carboniferous genesis of A-type magmas in Avalonia of northern Nova Scotia: repeated partial melting of anhydrous lower crust in contrasting tectonic environments. Int J Earth Sci (Geol Rundsch) 107, 587–599 (2018). https://doi.org/10.1007/s00531-017-1512-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00531-017-1512-7