Articles | Volume 12, issue 4
Research article
15 Apr 2021
Research article |  | 15 Apr 2021

A tectonic carpet of Variscan flysch at the base of a rootless accretionary prism in northwestern Iberia: U–Pb zircon age constrains from sediments and volcanic olistoliths

Emilio González Clavijo, Ícaro Dias da Silva, José R. Martínez Catalán, Juan Gómez Barreiro, Gabriel Gutiérrez-Alonso, Alejandro Díez Montes, Mandy Hofmann, Andreas Gärtner, and Ulf Linnemann

The allochthonous complexes of Galicia–Trás-os-Montes Zone (NW Iberia) are part of a rootless tectonic stack which preserves part of a Variscan accretionary prism. They are formed by individual tectonic slices marked by specific tectonometamorphic evolutions, which were piled up in a piggy-back mode onto its relative autochthon, the Central Iberian Zone (CIZ). Allochthony decreases from the structurally upper thrust sheets towards the lower ones. The lowermost unit of the stack is known as the Parautochthon or Schistose Domain. It is characterized by a low metamorphic grade in contrast with higher temperatures and/or pressures estimated for the overlying allochthonous units and shares the stratigraphic sequence with the underlying autochthon. The Parautochthon is divided in two structural and stratigraphic sub-units: (i) the Lower Parautochthon (LPa) is made of synorogenic flysch-type sediments with varied turbiditic units and olistostrome bodies, showing Upper Devonian–lower Carboniferous age according to the youngest zircon populations and fossiliferous content; (ii) the Upper Parautochthon (UPa) is composed of highly deformed preorogenic upper Cambrian–Silurian volcano-sedimentary sequence comparable with the nearby autochthon and to some extent, also with the high-P and low-T Lower Allochthon laying structurally above. The UPa was emplaced onto the LPa along the Main-Trás-os-Montes Thrust, and the LPa became detached from the CIZ relative autochthon by a regional-scale structure, the Basal Lower Parautochthon Detachment, which follows a weak horizon of Silurian carbonaceous slates.

A review on the detrital zircon studies on the synorogenic LPa complemented by zircon dating of 17 new samples is presented here. The results support the extension of the LPa underneath the NW Iberian allochthonous complexes, from Cabo Ortegal, to Bragança and Morais massifs. Its current exposure follows the lowermost tectonic boundary between the Galicia–Trás-os-Montes (allochthon) and Central Iberian (autochthon) zones. The youngest zircon age populations point to a maximum sedimentation age for the LPa formations ranging from Famennian to Serpukhovian and supports the piggy-back mode of emplacement of the Galicia–Trás-os-Montes Zone, of which it represents the latest imbricate.

The zircon age populations in the LPa allow the sedimentary provenance areas to be constrained, showing the intervention of nearby sources (mostly the UPa) and/or multiply recycled and long-transport sediments with a typically north-central Gondwana age fingerprint, also found in the Lower Allochthon, UPa and Autochthon. Complementary geochronology of volcanic olistoliths trapped in the LPa sediments and of late Cambrian to Upper Ordovician rhyolites from the UPa is also presented. It shows a direct relationship between the major blocks source area (UPa) and the setting place (LPa). Old zircon age patterns show that the LPa sedimentary rocks were recycled from detrital rocks of the allochthon (advancing wedge) and the nearby autochthon (peripheral bulge).

1 Introduction

Synorogenic marine basins encompass most of the known marine geodynamic settings, from active to passive earthquake-prone margins (Dickinson and Valloni, 1980; Garzanti et al., 2007; DeCelles, 2012). They are found to be associated with Archean to Phanerozoic orogens (e.g., Liang and Li, 2005; Wilmsen et al., 2009; Mulder et al., 2017; Kusky et al., 2020) and hold key evidence to understand the geographic and geodynamic evolution of modern and ancient orogenic belts.

A common sedimentary feature of all synorogenic marine basins is the presence of earthquake-triggered turbiditic flows that promote a variety of sedimentary facies, from cohesive rhythmic flysch sequences to chaotic large-scale mass-wasting bearing heterometric-sized blocks or olistoliths (Franke and Engel, 1986; Coleman and Prior, 1988; Eyles, 1990; Festa et al., 2020), also denominated block-in-matrix formations (hereafter referred as BIMFs) (Festa et al., 2016). Because the most probable source areas of the sediments and olistoliths that fed these basins are in the surrounding orogenically active highs, their stratigraphy can provide important clues on the orogen relief variation along space and time (e.g., Ducassou et al., 2014; Chiocci and Casalbore, 2017).

The synorogenic basins formed during continental convergence are gradually incorporated into the active orogenic edifice as tectonic slices in the accretionary complex, especially at the base, forming a tectonic carpet (Festa et al., 2019; Kusky et al., 2020, and references therein). The progradation of the tectonic front and the basin depocenter favors a systematic intrabasinal sedimentary recycling (i.a. wild flysch) and mixing of synorogenic sediments with other external sources (e.g., Franke and Engel, 1986; Bütler et al., 2011) producing mixed signals of not straightforward paleogeographic interpretation. The key to addressing this problem rests in the regional study of the basin stratigraphy, including the flysch sequences, the mass-wasting deposits and the petrography of the olistoliths (Festa et al., 2019, 2020). Complementary, the detailed geological recognition of the basement and surrounding areas of the synorogenic basins is crucial to identify discriminatory aspects that can help to constrain different variables, such as possible sources, sediment transport distance, regional and local tectonic settings, and paleogeographic limitations (e.g., Alonso et al., 2015; Festa et al., 2016; Krastel et al., 2019).

The sedimentological models can be refined using detrital zircon geochronology. This tool is commonly used to trace source-to-sink relationships, through different statistical approaches that compare the zircon age populations present in detrital rocks (e.g., Meinhold et al., 2011, 2013; Linnemann et al., 2012, 2014). One of the most rigorous procedures for zircon age populations similarity–dissimilarity analysis is multidimensional scaling (MDS) (Vermeech, 2018), which provides a graphic output of the Kolmogorov–Smirnov test when using a large number of samples with individual zircon age populations (Gutiérrez-Alonso et al., 2020; Pereira et al., 2020a).

In the southwestern edge of the European Variscan Belt (Fig. 1), the Iberian Massif preserves some of the best examples of Phanerozoic synorogenic marine basins, which reflect different tectonic settings along the belt during the Upper Devonian–late Carboniferous collision of Laurussia and Gondwana to form Pangea (e.g., Pereira et al., 2012a, b, 2017; Oliveira et al., 2019b). In SW Iberia, the Late Devonian–late Carboniferous flysch basins appear on both sides of the oceanic suture that separates Laurussia from Gondwana (Silva et al., 1990; Braid et al., 2011; Pereira et al., 2012a; Pérez-Cáceres et al., 2017). On the Gondwana-side, the Late Devonian–early Carboniferous marine sedimentation had Gondwana-type sources with massive contribution of intrabasinal volcanism in the Tournaisian–Visean period (Pereira et al., 2012a, 2020a). On the Laurussian side, sediments that filled the synorogenic marine basins resulted from intrabasinal recycling processes and source areas located on both continents. In this case, the sediments were systematically imbricated at the base of the advancing orogenic front, towards inland Laurussia from the Late Devonian (Pulo do Lobo Zone) to the upper Carboniferous (southwestern South Portuguese Zone) (Pereira et al., 2012a, 2020a; Pérez-Cáceres et al., 2017; Braid et al., 2011; Jorge et al., 2013; Rodrigues et al., 2015).

Figure 1Map of the European Variscan belt at the end of the Carboniferous. Modified from Martínez Catalán et al. (2007). Acronyms: CZ – Cantabrian Zone; WALZ – West Asturian–Leonese Zone; GTMZ – Galicia–Trás-os-Montes Zone; CIZ – Central Iberian Zone; OMZ – Ossa-Morena Zone; PLZ – Pulo do Lobo Zone; SPZ – South Portuguese Zone. Blue rectangle area is represented in Fig. 2.

This study presents new data and analysis of synorogenic rocks from marine basins located in the hinterland (internal zones) of the Variscan belt in NW Iberia. While synorogenic basins in the external zones of the orogen (the foreland fold and thrust belt) have been classically classified as foreland basins laying on top of north Gondwana Cambrian to Upper Devonian passive margin sequences (Marcos and Pulgar, 1982; Pastor-Galán et al., 2013; Gutiérrez-Alonso et al., 2015), the basins studied in this work are located in the hinterland and were deposited over rocks showing pervasive strain and that were metamorphosed to different degrees (Martínez Catalán et al., 2004, 2008, 2016; Dias da Silva et al., 2015). The sedimentary sources of the synorogenic flysch and BIMF deposits are related to the development and unrooting of a Variscan accretionary prism (the Galicia–Trás-os-Montes Zone) onto Gondwana, and to the development of a peripheral bulge affecting the extensive passive margin of Gondwana (e.g., González Clavijo and Martínez Catalán, 2002; Keller et al., 2008; Dias da Silva et al., 2015). Some of the Variscan hinterland synorogenic basins have been incorporated into the base of the allochthonous wedge as a parautochthonous unit and then emplaced onto the autochthonous terrain of NW Iberia (Dias da Silva et al., 2015, González Clavijo et al., 2016). The basin depocenter migrated towards inland Gondwana from the Late Devonian to the late Carboniferous following the progression of the orogenic front, from the Galicia–Trás-os-Montes Zone (GTMZ) to the Central Iberian Zone (CIZ), where synorogenic deposits are preserved in the San Clodio Series (Martínez Catalán et al., 2016), possibly following to the West Asturian–Leonese zone (WALZ), although no synorogenic deposits linked to this basin are preserved on it, and finally to the foreland Cantabrian Zone (CZ) (Merino-Tomé et al., 2017; Gutiérrez-Alonso et al., 2020).

In both NW and SW Iberian synorogenic basins, zircon geochronology has been used to constrain sedimentary provenance based on the fingerprint of sources and basin stratigraphic units (Braid et al., 2011; Pereira et al., 2012a; 2014; Jorge et al., 2013; Pastor-Galán et al., 2013; Rodrigues et al., 2015; Martínez Catalán et al., 2016; Pérez-Cáceres et al., 2017). While in the southwest and in the CZ recent works have demonstrated the importance of the MDS statistical approach in identifying the relationships between source and sink (Gutiérrez-Alonso et al., 2020; Pereira et al., 2020b), in NW Iberian Variscan hinterland this approach has not been applied to date.

In this work we present new field observations that revise previous interpretations, supported by new structural data and U–Pb geochronology of igneous and detrital zircon grains, which enables a new vision of the parautochthonous units of the GTMZ here referred to as Upper Parautochthon (UPa; preorogenic) and Lower Parautochthon (LPa; synorogenic) (Dias da Silva et al., 2015). Our findings support the extension of these units to different sectors of the GTMZ, extending beyond the area covered by the UPa and LPa to virtually below all the allochthonous complexes of NW Iberia. Field work has revealed the diverse types of flysch complexes and mélanges present in the LPa, commonly obscured by Variscan polyphasic and pervasive deformation. Tectonic and sedimentary mélanges have been recognized; they combine to produce polygenetic mélanges using the terminology by Festa et al. (2019, 2020). Moreover, the detrital zircon age fingerprinting offers a new general view of the LPa geotectonic setting at Variscan times and sets new constrains on the possible source areas of the synorogenic sediments, olistoliths and blocks. We consider that this review led to a better understanding on the paleogeography and geodynamic setting of the Late Devonian–lower Carboniferous flysch basins in NW Iberia and its incorporation into a general model for the synorogenic basins of the Iberian Variscan Massif.

2 Geological setting

Within the rootless Variscan accretionary wedge which forms the Variscan Massif of NW Iberia – the so-called GTMZ (Figs. 2 and 3) (Ribeiro, 1974; Schermmeron and Kotsch, 1984; Martínez Catalán et al., 2009; Ballèvre et al., 2014; Martínez Catalán et al., 2014; Azor et al., 2019) – a significantly distinct structural unit has been identified and interpreted as a remnant of former oceanic realms (Lower Ordovician Rheic Ocean and an Early Devonian suprasubduction ophiolite). This tectonic unit marks the putative suture zone of the Laurussia–Gondwana continental collision that partially led to the formation of Pangea; this unit is known as a Middle Allochthon (MA) or Ophiolitic unit (Gómez Barreiro et al., 2007; Martínez Catalán et al., 2009; Stampfli et al., 2013; Ballèvre et al., 2014; Arenas and Sánchez-Martínez, 2015; Azor et al., 2019) and currently separates the so-called Upper Allochthon (upper units) and the Lower Allochthon (basal units). The Upper Allochthon (UA) is considered a far-traveled ribbon-shaped terrane that drifted away in the Lower Ordovician from the Gondwanan margin during the opening of the Rheic Ocean and accreted to Laurussia in the Silurian (Gómez-Barreiro et al., 2007). It includes two tectonically stacked units, presenting early-Variscan (ca. 390–380 Ma) HP/HT and Cambro-Ordovician IP/M to HT metamorphism respectively (Martínez Catalán et al., 2019). The Lower Allochthon (LA) is made of a set of nappe folds and tectonic slices (Farias et al., 1987; Díez Fernández et al., 2010; Dias da Silva et al., 2014, 2015) considered to represent the most seaward rim of continental Gondwana (Murphy et al., 2008) which underwent continental subduction (HP/L to MT metamorphism) and obduction (retrogression to amphibolite and greenschist facies) recording the inception of Variscan continental collision at ca. 370–360 Ma (Munhá et al., 1984; Gil Ibarguchi and Dallmeyer, 1991; Arenas et al., 1995, 1997; Gil Ibarguchi, 1995; Rubio Pascual et al., 2002; Rodríguez et al., 2003; López-Carmona et al., 2010, 2014). This latter unit, (LA) despite being interpreted as part of the Gondwanan passive margin ensemble, and being part of the lower plate in the collisional edifice, is classically considered as part of the allochthonous realm in the region but not belonging to the exotic terrains (Ribeiro, 2013; Ribeiro and Sanderson, 1996; Ribeiro et al., 2007). Structurally below the LA, a tectonic unit displaying low metamorphic grade separates the above-mentioned allochthons from their relative autochthon, the Central Iberian Zone (CIZ). This unit, named Schistose Domain (Farias et al., 1987) or Parautochthon (Pa) (Ribeiro et al., 1990; Martínez Catalán et al., 1997), was considered as a thick Silurian sequence, lacking correlation with the condensed Silurian graphite-rich sequences of the underlying autochthon (CIZ). Nevertheless, the paleogeographic affinity of both domains was highlighted by identical north Gondwana Silurian graptolite and conodont faunas (Sarmiento et al., 1998; Piçarra et al., 2003, 2006a, b). However, the stratigraphy and structural features of the lowermost tectonic sheets of the Pa, directly above the CIZ, were interpreted as a synorogenic basin with possible (Middle–Late) Devonian age (Antona and Martínez Catalán, 1990; González Clavijo and Martínez Catalán, 2002; Martínez Catalán et al., 2004; Pereira et al., 2009; Rodrigues et al., 2013).

Figure 2Simplified geological map of the NW Iberian Variscan Massif modified from Martínez Catalán et al. (1997). The limits between Lower Allochthon, Upper and Lower Parautochthon, and the autochthonous unit have been modified according to this work.

Figure 3Representative cross sections of the Variscan belt in the region studied including the inferred geometry of the synorogenic tectonic carpet as proposed in this study. See Fig. 2 for location and legend. Cross sections modified from (1) Marcos et al. (1984), Arenas (1988); (2) Martínez Catalán et al. (2004); (3) Pereira et al. (2006), Ribeiro (1974), Díez Montes (2006); (4) Pereira et al. (2006), Ribeiro (1974), González Clavijo and Martínez Catalán (2002); (5) González Clavijo and Martínez Catalán (2002); and (6) Pereira et al. (2006), Dias da Silva (2014).

The attempt to better understand the tectonostratigraphy of the Pa led to a later division in two tectonically stacked units, Upper and Lower Parautochthon (UPa and LPa), in the sense firstly proposed by Rodrigues et al. (2006, 2013) and updated by Dias da Silva et al. (2014, 2015, 2016). This division limits the UPa to a pre-Variscan upper Cambrian–Silurian sequence comparable with the CIZ and LA that was affected by Variscan recumbent folds and thrusts, and it defines the LPa as an imbricated thrust sequence bearing slices of a foreland synorogenic basin, with the younger slices in the transition to the CIZ (Martínez Catalán et al., 2016). The thrust fault structures bounding the lower tectonic sheets of the GTMZ are as follows (Figs. 2 and 3):

  • i.

    The LA Basal Thrust (LABT, Figs. 3 and 4) or basal thrust of the Centro-Transmontano thrust complex (in the meaning of Ribeiro et al., 1990) represents the roofing thrust of the Parautochthon.

  • ii.

    The UPa–LPa thrust system (Main Trás-os-Montes Thrust, MTMT; Ribeiro, 1974; Ribeiro and Ribeiro, 2004; Meireles et al., 2006; Pereira et al., 2006) is a gently dipping low-grade shear-zone up to 1000 m thick, interpreted to be caused by the thrusting of the upper Cambrian–Silurian (preorogenic UPa) sequence onto the synorogenic LPa, producing significant crustal thickening during the Tournaisian–Visean stage (Dias da Silva et al., 2014, 2015, 2016, 2021; Azor et al., 2019).

  • iii.

    At the base of the LPa another major bedding-parallel fault named the Basal Lower Parautochthon Detachment (BLPD), also gently dipping, is the sole fault separating the synorogenic imbricated slices from the structurally underlying nonimbricated autochthon (Dias da Silva et al., 2014). The BLPD was developed using a slip-favorable stratigraphic unit, the autochthonous Silurian carbonaceous–siliceous slates (SCSSs) (González Clavijo and Martínez Catalán, 2002; Dias da Silva et al., 2014).

The northern CIZ autochthonous domain consists of an Ediacaran to Lower Devonian preorogenic sequence, disturbed by two regional unconformities, and including ca. 490 to 460 Ma felsic to intermediate, locally mafic magmatism; Floian Armorican-type quartzites; a Middle–Upper Ordovician mostly detrital sequence; and the SCSSs (e.g., Sousa, 1984; Valladares et al., 2000; Gutiérrez-Marco et al., 2019; Sánchez García et al., 2019). In the autochthon laying immediately below the BLPD, tectonic overburden mainly occurred due to the action of thin-skinned imbricated thrust-duplexes rooted in the BLPD, developed in the LPa tectonic units (Fig. 2 and sections 4 and 5 in Fig. 3) (Dias da Silva et al., 2021). In the sectors where the LPa is present, thickening in both CIZ and LPa was rapidly attenuated by the succeeding synorogenic extensional processes (Dias da Silva et al., 2021).

Figure 4Regional correlation of the NW Iberian LPa stratigraphic units, and their structural relationship with the tectonically overlying (UPa and allocthonous units) and underlying geologic domains (autochthon, CIZ). See Fig. 5 for more information.


In the area studied, the UPa and CIZ underwent regional Barrovian metamorphism (M1) through the early Variscan compressive events (C1+ C2 on the Martínez-Catalán et al., 2014 proposal; ca. 360–330 Ma) which were followed by a complex extensional (E1 M2; ca. 340–320 Ma) and compressional (C3 M3; ca. 318–300 Ma) tectonothermal history (Dallmeyer et al., 1997; Azor et al., 2019; Dias da Silva et al., 2021, and references therein).

The LPa synorogenic ensemble was also affected by anchizone to low-grade (chlorite zone) metamorphism. Attempts to discriminate if the metamorphic grade was lower in the synorogenic units than in the preorogenic units using illite crystallinity (Antona and Martínez Catalán, 1990) and the Colour Alteration Index in conodonts (Sarmiento and García-López, 1996; Sarmiento et al., 1997) were inconclusive. Petrographic observations made by Matte (1968) in the synorogenic deposits of the San Clodio in the CIZ to the southeast of Monforte (Fig. 2) and in the underlying Ordovician sequence, and by Dias da Silva et al. (2020) in the LPa and CIZ in the eastern rim of Morais Complex, show similar low-grade epizone metamorphism in both pre- and synorogenic sequences. However, the San Clodio flysch rests unconformably above the reverse limb of a large C1 recumbent syncline whose axial planar cleavage is more evolved than that of the flysch above (Martínez Catalán et al., 2016). And it is also older: ca. 360 Ma (Dallmeyer et al., 1997), while detrital zircons are as young as 324 Ma in the San Clodio Series and 340 Ma in the synorogenic deposits of Trás-os-Montes (Martínez Catalán et al., 2004, 2008, 2016), the age of emplacement of the allochthonous complexes during C2. The first foliation in the preorogenic metasediments of the UPa and CIZ is axially planar to recumbent folds of the C1 event and predates the main foliation in the synorogenic deposits. But a second, low-grade penetrative foliation was developed in the UPa, LPa and CIZ during the emplacement of the Allochthon (C2). This second regional foliation is the one showing similar aspect and metamorphic conditions in both UPa and LPa ensembles, but in the latter it represents the first tectonic fabric.

3 Review of the synorogenic marine sequences in NW Iberian Variscan hinterland

The internal zones of the orogenic belts are considered areas with scarcely preserved related synorogenic sequences because of the subsequent denudation caused by the orogenic relief (Martínez Catalán et al., 2008), but they might be preserved in the core of synclines or below post-depositional thrust. In other parts of the Variscan belt, Franke and Engel (1986) described tectonic slices carrying synorogenic sedimentary units from more internal areas. In NW Iberia, while there is a complete record of the synorogenic deposits in the foreland fold and thrust belt (CZ; e.g., Marcos and Pulgar, 1982; Merino-Tomé et al., 2017), in the deeply eroded internal part of the chain, the first identified synorogenic sequence was the San Clodio Series (Matte, 1968), which is preserved at the core of the late-Variscan Sil Syncline in northern Iberia (Figs. 2, 4 and 5, and cross section 2 in Fig. 3). It consists of a rhythmic turbiditic sequence of pelite and greywacke (Riemer, 1966; Pérez-Estaún, 1974) including Upper Devonian or even older fossil plant fragments (Pérez-Estaún, 1974). Previous detrital zircon studies (samples SO-1 and SO-2; location of all samples from previous studies are plotted in the map of Fig. SF1.1 of the Supplement) have reinforced the synorogenic character and supported an Upper Mississippian maximum depositional age (Martínez Catalán et al., 2004) according to the youngest ages found in the detrital zircon population of ca. 324 Ma. Towards the base it displays exotic lithic blocks and pebbles, that is, extrabasinal rocks derived from the basement exposed in the basin surrounding highs. These include carbonaceous chert, quartzite, slate, gneiss and granite. Intrabasinal (“native”) lithified intraclasts and soft pebbles also occur (Riemer, 1966). The San Clodio Series is often separated from the underlying Ordovician formations by a few meters of Silurian black shales (SCSSs) which were mylonitized forming a basal detachment (Barrera Morate et al., 1989). However, the unconformity was preserved from reactivation at a few places (Martínez Catalán et al., 2016). Following Festa et al. (2019, 2020) terminology, the San Clodio Series represents a coherent primary succession above a sedimentary block-in-matrix unit (olistostrome).

Figure 5Synorogenic lithostratigraphic units of the Lower Parautochthon at the different sectors studied in this work. The sketch displays a simplified structure avoiding the minor tectonic slices repeating every stratigraphic unit. The sequences have been prepared considering all the data referred in the text and Fig. 4 for every sector and using the division and names that better fit with our field survey.


The presence of a variably deformed SCSS unit below the synorogenic lithostratigraphic units in the LPa, as described for the Sil Syncline, is a constant feature in all the areas incorporated into this study. The BLPD is a first-order structure that forms a complex arrangement of stacked tectonic slices depicting diverse patterns. Locally, this shear band only deforms the lower part of the SCSSs, preserving the sedimentary unconformity at the base of the synorogenic unit (Alcañices syncline; González Clavijo, 2006). In other places the BLPD involves the entire SCSSs with upper and lower shear bands that limit lower-order shear band structures merging to create a first-order S/C structure (eastern Morais Complex; Dias da Silva, 2014). In many sections, deformation along the BLPD involves the BIMF deposits placed at the base of the LPa, thus resulting in a polygenetic mélange in the Festa et al. (2019) meaning. Also common is the presence of lower-order thrust duplexes within the LPa that are rooted in the BLPD. In these cases (San Vitero Formation in the Alcañices Syncline, described in this section) the LPa is internally repeated by a stack of imbricated tectonic slices each with the SCSSs at the base. The lensoid shape of the SCSSs along the BLPD and associated structures also suggests that it was submitted to a strong tectonic pinching (thinning) and swelling (thickening) during thrusting. However, deformation is not always pervasive, as several size lenses of comparatively undeformed SCSS preserved fossiliferous content at many localities (Romariz, 1962, 1969; Quiroga de la Vega, 1981; González Clavijo et al., 1997; Piçarra et al., 2006b), defining all the Silurian graptolite biozones but lower Rhuddanian and upper Ludfordian stages as well as the Pridoli series (González Clavijo, 2006; Piçarra et al., 2006b).

The most extensive outcrop of the synorogenic rocks studied occurs in the periphery (structurally below) of the Bragança Complex (Figs. 2 and 4). To the east of this allochthonous complex, the core of the Late Variscan Alcañices Synform (Figs. 2, 4 and 5; cross sections 4 and 5 in Fig. 3) is formed by several LPa synorogenic units tectonically piled up in a number of imbricated thrust units (C2) folded by a train of NW–SE-trending upright folds (C3 M3) (González Clavijo and Martínez Catalán, 2002; González-Clavijo et al., 2012). From the top (more internal) structural position to the bottom (more external) the synorogenic units are named: Gimonde, Rábano, San Vitero and Almendra formations. According to the maximum depositional age (MDA) obtained through detrital zircon U–Pb geochronology (Upper Devonian to uppermost Mississippian), these formations include progressively younger rocks from the more internal to the more external ones, thus suggesting migration of the depocenter toward the relative autochthon coeval to the Parautochthon stacking (González-Clavijo et al., 2012; Martínez Catalán et al., 2016).

The structurally highest Gimonde Formation (Pereira et al., 1999; Meireles et al., 1999a, b) is formed by finely bedded phyllites and metagreywackes, and scarce polymictic microconglomerate lenses containing exotic clasts and native intraclasts. Its age has been considered Upper Devonian on base of fossil plant debris (Teixeira and Pais, 1973) and palynomorphs (Pereira et al., 1999). However, detrital zircon ages (samples SO-7, SO-8, SO-9, SO-12; Martínez Catalán et al., 2016) imply an early Carboniferous (Tournaisian–Visean) MDA.

The structurally underlying Rábano Formation (González Clavijo and Martínez Catalán, 2002), also called External Gimonde in Martínez-Catalán et al. (2016), comprises diverse lithologies being the most abundant a BIMF sequence displaying large exotic olistoliths of deformed rhyolites and dacites, felsic metatuffs, epiclastic rocks, white and grey quartzite, Silurian lydite (black radiolarite) and ampelite (carbonaceous shale), greywacke, phyllite, and limestone (Fig. 6a). The ages of these blocks (sometimes hundreds of meters in length) based on fossils and U–Pb zircon ages range from Furongian to Emsian (González Clavijo et al., 2016). At the uppermost Rábano Formation a flyschoid sequence made of phyllite, quartzlitharenite, and local polygenic microconglomerate holds deformed exotic clasts and lithified intraclasts (Fig. 7a, b, and c), including plagioclase and volcanic quartz mineraloclasts and quartz shards indicating a nearby volcanic source (González Clavijo, 2006). Detrital zircon ages performed in this wild flysch (sample SO-6) support a synorogenic nature and point to a Visean MDA (González Clavijo et al., 2016; Martínez Catalán et al., 2016).

Figure 6Field aspects of the synorogenic sediments of the LPa: (a) Lower Ordovician white quartzite olistoliths (arrows) in Rábano Formation; (b) centimeter-sized foliated rhyolite in an Almendra Formation microconglomerate; (c) metric glided block of limestone in an Almendra Formation grey phyllite bed; (d) slump folds in flyschoid facies at the Meirinhos synorogenic LPa; (e) broken beds in a flyschoid facies of the Meirinhos Formation; (f) centimeter-thick flysch facies in the LPa unit exposed in the tectonic window western Mirandela; (g) quartzite, lydite and rhyolitic tuff olistoliths in the LPa northern Bragança Complex; (h) blocks of quartzite (under the hammer) and lydite (above) in the LPa northern Bragança Complex.


Figure 7Microphotographs of the sedimentary textures in the LPa formations, surrounding clasts of different natures. (a) Bioclast in a Rábano Formation quartzlitharenite; (b) plagioclase mineraloclast in a Rábano Formation quartzlitharenite; (c) broken up volcanic quartz crystals in a Rábano Formation quartzlitharenite; (d) lithoclast made of black phyllite in a San Vitero Formation litharenite; (e) rounded clast displaying tectonic foliation almost normal to the surrounding external C1 foliation (San Vitero Formation); (f) partial view of a rounded clast displaying mylonitic banding in an Almendra Formation microconglomerate; (g) randomly oriented clasts bearing previous foliation in an Almendra Formation microconglomerate; (h) rounded foliated clast displaying a microfold in an Almendra Formation microconglomerate.


Figure 8Examples of tectonic mélanges at the base of the synorogenic units in the Alcañices synform. (a) San Vitero duplex slices showing a Lower Ordovician volcanic olistolith above the Silurian fossiliferous sequence. (b) Rábano Formation with rhyolite olistoliths of unknown age and displaying sedimentary block-in-matrix facies lately incorporated to a tectonic mélange.


The San Vitero Formation (Martínez García, 1972) is a flysch made of up to meter-thick phyllite and quartzlitharenite rhythms and local lenses of polygenic microconglomerates holding exotic clasts and lithified intraclasts (Figs. 7d, e and 8a, b) (González Clavijo and Martínez Catalán, 2002). This unit was considered Upper Devonian or younger based on fossil plant debris (Teixeira and Pais, 1973) but detrital zircon studies (samples SO-4, SO-5 and SO-13) support a Tournaisian MDA (Martínez Catalán et al., 2016).

The structurally lower imbricated thrust system hosts the Almendra Formation (Vacas and Martínez Catalán, 1987) which is a calciturbidite made of phyllite and calcarenite rhythms up to several meters thick (González Clavijo, 2006). Local lenses of polygenic conglomerates and microconglomerates holding exotic clasts and pebbles, and lithified intraclasts were firstly described by Aldaya et al. (1976). Their lithologies include phyllite, sandstone, litharenite, quartzite, limestone, orthogneiss, rhyolite and felsic volcanic tuff (Fig. 7f, g and h). Major blocks of lydite with Silurian graptolites and limestones (Fig. 6b and c) containing abundant fossils (bioclasts of corals, sciphocrinoides, bivalves, gastropods and tentaculites) have been identified (González Clavijo and Martínez Catalán, 2002; González Clavijo et al., 2016 and references therein). Conodonts found in calcarenite yielded a Lower Devonian age (Sarmiento et al., 1997). Detrital zircon studies in the Almendra Formation indicated a Visean MDA (sample SO-14), thus supporting the Variscan synorogenic origin of this unit (Martínez Catalán et al., 2016).

The LPa emplaced at the eastern rim of the Morais Complex (Figs. 2, 4 and 5; cross section 6 in Fig. 3) was described by Dias da Silva (2014) as a turbiditic synorogenic sequence comprising two stratigraphic units (Travanca and Vila Chã Formations). They consist of coherent primary units, broken bed units and lenses of block-in-matrix units displaying hectometer-sized olistoliths. Clasts include native intraclasts and soft clasts as well as exotic deformed lydite, ampelite and quartzite. These units lack fossil-based ages. A palynomorph study performed by Gil Machado (in Dias da Silva, 2014) was unfruitful because of poor pollen preservation owing to a Variscan thermal imprint. Detrital zircon studies (samples VC-21ZIR, VC-45ZIR and VC-57ZIR) suggest a Devonian MDA (ca. 390 Ma; Dias da Silva et al., 2015).

In the Marão Range, west of Vila Real (Fig. 2), the westernmost extent of the rocks studied, tectonic slices appertaining to the Pa comprise several turbiditic sequences (Mouquim, Canadelo and Santos Formations) displaying rhythms of phyllite and greywacke with some intercalations of volcanic tuffs towards the top (Pereira, 1987). As they lie concordantly above the SCSSs dated by graptolites (Piçarra et al., 2006b), a Devonian age for the flyschoid units was proposed, without fossiliferous evidence (Pereira et al., 2006). González Clavijo (2006) supported a correlation between these units and the San Vitero flysch, in the Alcañices synform, on the basis of the lithologies and the stratigraphic position, and for that reason here a possible Tournaisian MDA is proposed. According to our proposal of a tectonostratigraphic scheme, all the stacked pile must be considered as belonging to the LPa.

4 LPa synorogenic units: youngest detrital zircon age populations of the LPa units

The known existence of synorogenic LPa units, partially encircling the GTMZ to the east, fostered the present research and aimed to recognize their possible extension to other areas of the GTMZ (Fig. 2). In this work 17 samples from areas surrounding the Bragança, Morais and Cabo Ortegal complexes were used for zircon U–Pb geochronology (location of all samples in the map of Fig. SF1.1; coordinates in Table SF1).

The detailed description of the new zircon geochronology study is presented in the Supplement and in Figs. SI-1 to SI-9. The complete dataset with the new U–Pb isotopic analyses is given in Tables S1–S17 in the Supplement. The reference and complementary U–Pb zircon age datasets used in the multidimensional scaling (MDS) process and other statistical procedures are in Tables S18 and S19.

In this section, the youngest zircon grain and zircon population ages of our new results are discussed for the different parts of the synorogenic carpet and compared with published data.

The Meirinhos area, to south of the Morais Complex (Figs. 2, 4 and 5; cross section 4 in Fig. 3), was divided in two stratigraphic units: Meirinhos and Casal do Rato following Pereira et al. (2009) and Rodrigues et al. (2003). They were considered synorogenic flyschoid deposits bearing olistoliths of quartzite, phyllite, greywacke, felsic and mafic volcanic tuffs, limestone, ampelite and lydite by Pereira et al. (2009). However, different age, stratigraphic features and structure were proposed (Sá et al., 2014) based on the reappraisal of the Lower Ordovician trilobite Cruziana ichnofossils in Armorican-type quartzites originally described by Ribeiro (1974) in this sector. Field data allow us to confirm the earlier proposal and the recognition of new synorogenic features as slump folds, disrupted beds, and olistostromes including hectometer-sized lydite and quartzite olistoliths (Fig. 6d and e). In both stratigraphic units, blocks, cobbles and pebbles of native (lithified intraclasts, soft pebbles) and exotic (quartzites, ampelite and lydite, felsic and mafic volcanic rocks, limestone) sources have been identified. Biostratigraphic ages of the olistoliths range from early Ordovician (trilobite tracks, Ribeiro, 1974; Sá et al., 2014) to Silurian (graptolites, Pereira et al., 2009). Two samples of medium-grained greenish greywacke belonging to the turbiditic sequence (CR-ZR-01 and MEI-ZR-01; see location, coordinates and geochronology of all geochronology samples in the Supplement) were collected in this area with a youngest single zircon (YZ) with Upper Ordovician age (439±8 Ma for CR-ZR-01 and 443±6 Ma for MEI-ZR-01) and MDA (i.e., concordant age given by the youngest zircon population in the 90 % or 95 % concordance interval; it may not include the YZ) of 467 ± 4 Ma (CR-ZR-01) and 486 ± 3 Ma (MEI-ZR-01) (Supplement; Fig. SI-1A and B). These ages do not support a synorogenic character nor a Lower Ordovician age for these siliciclastic rocks. However, the sedimentary features and the cartographic and structural continuity and correlation with the previously described LPa unit east of the Morais Complex (i.e., Travanca Formation in Dias da Silva et al., 2015) make it possible to propose a Mississippian MDA for these synorogenic siliciclastic rocks, with contribution of early Silurian, Middle Ordovician and Tremadocian (possibly magmatic) zircon sources.

West of Mirandela, around the village of Suçães, there is a tectonic window (Rodrigues et al., 2010) partially controlled by late- or post-Variscan NNE–SSW subvertical faults (Figs. 2, 4 and 5; cross sections 3 and 6 in Fig. 3) which displays a turbiditic sequence attributed to the Devonian (Ribeiro, 1974) or to the Silurian with small patches of Lower Devonian siliciclastic rocks in the upper stratigraphic positions (Rodrigues et al., 2010). The MTMT was firstly mapped in this area by Rodrigues et al. (2010), tectonically separating the UPa overturned fold sequence from the LPa imbricated thrust complex. The proposed Silurian–Devonian age was based on graptolite assemblages preserved in the Silurian lydites (Piçarra et al, 2006b). Field and geochronological data allow us to infer that the LPa in this sector is also a synorogenic sequence, presenting the classical flyschoid features (Fig. 6f) (made of centimeter to meter beds of pelite and greywacke; Rodrigues, 2008). Scarce Variscan detrital zircon grains (Famennian) were found in the greywacke layers of the flysch sequence, in the Upper Schists Formation (MIR-41: YZ = 369 ± 7 Ma, MDA = 497 ± 5 Ma, Fig. SI-4A; AD-PO-49: YZ = 468 ± 39 Ma, MDA = 494 ± 27 Ma, Fig. SI-4B; AD-PO-55: YZ = 444 ± 26 Ma, MDA = 488 ± 16 Ma, Fig. SI-5A) and in the culminating slate and greywacke formation (sample AD-PO-57: YZ = 372 ± 6 Ma, MDA=488 ± 16 Ma, Fig. SI-5B). The combination of field and geochronology data suggest that the graptolite-rich lydites belong to exotic olistoliths or to the SCSSs at the base of the synorogenic tectonic slices. Detrital zircon ages and the lithostratigraphy of these stratigraphic units enable us to consider them as coherent primary units with olistoliths, overlying a tectonic mélange developed in the SCSSs.

To the west, in the Vila Pouca de Aguiar area (Fig. 2), the LPa comprise several tectonically stacked flysch units limited by thrust planes (Ribeiro, 1974; Ribeiro et al., 1993; Noronha et al., 1998; Ribeiro, 1998; Rodrigues, 2008). All units are formed by turbidite sequences, with millimeter- to meter-thick interbedded pelites and greywackes (quartzwacke towards the base). The coarse-grained layers are very rich in plagioclase, thus suggesting a near-source (felsic) volcanic input into the basin (Rodrigues, 2008). Differently sized lens-shaped bodies of grey quartzite, black limestone, felsic metavolcanic rocks and highly sheared SCSSs are described within these units (Rodrigues, 2008). Our field research permitted us to define these LPa formations as coherent primary units with block-in-matrix deposits exposing exotic fragments in a synorogenic convolute sediment tectonically repeated by imbricated thrust faults. These units are currently considered as Silurian–Devonian based on Silurian graptolites preserved in the lydites (Piçarra et al., 2006b) and by lithological comparison with similar formations to the east of the Bragança and Morais Complexes (Ribeiro, 1974; Noronha et al., 1998; Ribeiro, 1998; Pereira et al., 2000). The detrital zircon geochronology of one sample of lithoclastic coarse-grained arkose (AD-PO-48B) returned a Famennian–Tournaisian youngest single zircon and maximum depositional age (YZ: 355 ± 34 Ma; MDA: 364 ± 22 Ma; Fig. SI-3B).

At the northern edge of the Bragança Complex (Figs. 2, 4, 5 and 9; cross section 3 in Fig. 3), the Upper Allochthon high P/high T rocks tectonically overlay the LPa. Field work confirmed the absence of the UPa stratigraphic units and the discontinuous cartographic outline of the Middle and Lower Allochthon (Meireles et al., 1999a, b). These units were tectonically thinned by an extensional shear zone or truncated by an out-of-sequence thrust system. In the upper structural part of the LPa, Ribeiro and Ribeiro (1974) noticed the presence of differently sized rock fragments in the stratigraphic sequence. They describe (i) epizonal fragments such as phyllite, quartz-phyllite, quartzite, felsic tuffs and rhyolites, ampelite, and lydite, and (ii) meso-catazonal fragments of paragneiss (albite, chlorite and K feldspar), blastomylonite, and biotite–garnet gneiss. We have verified that these rocks appear as centimeter- to hectometer-sized olistoliths, being dispersed in a wider area than previously estimated (Figs. 6g, h and 10a). The olistoliths usually cluster within large mass-wasting deposits (BIMFs) or occur as isolated bodies in the flysch sequence. The native rock blocks consist of fragments of consolidated fine- to coarse-grained greywacke beds among other siltstone–sandstone intraclasts and pelitic soft pebbles. The exotic blocks include highly deformed lydites and ampelites, rhyolites, felsic tuffs, quartzites, and limestones. Regional work in this area (Meireles, 2013) divides this sector LPa in four formations: Coroto, Rio de Onor, Soutelo and Gimonde, all of them bearing flyschoid characteristics and olistoliths (Figs. 4 and 5). In some of these units (Coroto and Soutelo) the limestone blocks yield upper Silurian to Lower Devonian crinoids (Meireles, 2013). Inherited Cambrian–Middle Ordovician acritarchs and lower–middle Silurian palynomorphs were also found in the Rio de Onor Formation (Pereira et al., 1999). In the Soutelo Formation, we also report the presence of a sandy–quartzitic olistolith with poorly preserved brachiopods of the “Lingula” genus (Sofia Pereira and Jorge Colmenar, personal communication, 2017), which are particularly common in the Lower Ordovician Armorican quartzite of the CIZ (Marão Range; Coke et al., 2001). Three samples collected in the Gimonde Formation (Meireles et al., 1999a, b) were used for geochronology of detrital zircon grains: GIM-ZR-01 (microconglomerate), EC-PO-293 and AD-PO-66 (quartz-lithic sandstones). Although two of these samples have yielded Silurian (EC-PO-293: YZ = 426 ± 44 Ma and MDA = 435 ± 29 Ma; Fig. SI-2B) and Furongian (AD-PO-66: YZ = 431 ± 31 Ma and MDA = 483 ± 18 Ma; Fig. SI-3A) youngest zircon ages, one sample has Variscan detrital zircon ages (GIM-ZR-01: YZ=327 ± 9 Ma and MDA = 354 ± 3 Ma; Fig. SI-2A) confirming the results obtained in previous studies (samples SO-9 and SO-12 in Martínez Catalán et al., 2016), which allowed the characterization of the Gimonde Formation as a Tournaisian–Visean synorogenic stratigraphic unit.

The Picón Beach at Cabo Ortegal Complex is placed east of the Ortigueira locality, in the Galicia northern coastline (Figs. 2, 4 and 5; cross section 1 in Fig. 3). There, the Loiba unit of the Rio Baio thrust sheet (Marcos et al., 2002) is structurally on top of the tectonically sheared (autochthonous) SCSSs that define the BLPD. This tectonic unit is formed by low metamorphic grade flysch sequences and discontinuous lens-shaped BIMFs. Its field aspects and structural relationships with the overlying and underlying tectonostratigraphic units led us to consider this unit as part of the Variscan synorogenic marine basin. The detrital geochronology of a fine-grained quartzite sample (PICON-2) has yielded a Tournaisian maximum depositional age (YZ: 350 ± 7 Ma; MDA: 357 ± 4 Ma; Fig. SI-6) thus supporting the extent of the LPa from northeast Portugal to the northern Spanish coast.

5 Magmatic zircon ages of the LPa olistoliths

In all lithostratigraphic units of the LPa exotic (extrabasinal) and native (intrabasinal) grains, clasts, pebbles and large olistoliths have been identified. Exotic rock fragments show presedimentary mild to high deformation and metamorphic aspects, contrasting with the usually poorly deformed, low-metamorphic flyschoid sequence. The native fragments include soft pebbles and intraclasts, often showing synsedimentary deformation features (slump folds, boudinage, bedding disruption and convolute bedding of turbiditic aspect: Fig. 6d and e). The most consolidated sedimentary fragments suggest that they were recycled within the basin, as a wild flysch (Ribeiro and Ribeiro, 1974; Aldaya et al., 1976; González Clavijo and Martínez Catalán, 2002; Martínez Catalán et., 2016). These fragments were considered a proof of its synorogenic character, and also evidence that the basin was fed from areas of the Variscan belt already deformed and metamorphosed (Antona and Martínez Catalán, 1990; González Clavijo and Martínez Catalán, 2002; Martínez Catalán et al., 2004, 2008). Complementarily, the fossil flora, fauna, and ichnofossil findings in the exotic olistoliths display ages from Lower Ordovician to Middle Devonian (Fig. 5). This wide range of biostratigraphic ages seems to indicate that fossil findings belong to rock blocks within a synorogenic turbiditic unit, which is confirmed by stratigraphic features and detrital zircon ages (González Clavijo et al., 2016). Trying to confirm this hypothesis, a U–Pb zircon geochronology study on volcanic rocks considered to represent olistoliths from the LPa was performed (Fig. SF1.1). We have complemented our study with the integration of already published U–Pb zircon age data from other volcanic olistoliths in the LPa (Farias et al., 2014; González Clavijo et al., 2016). The complete description of the samples, their location and their geochronology is presented in the Supplement.

5.1 New ages from magmatic olistoliths

As a complementary study of the detrital zircon research in the LPa flyschoid sequences, four volcanic rock olistoliths from the LPa in the Alcañices synform and in the northern edge of the Bragança Complex were sampled for U–Pb geochronology.

Sample EC-PO-337 was picked northeast of Bragança in an olistolith made of low-metamorphic-grade foliated green metadacitic pyroclastic tuff from the Rábano Formation. This rock consists of volcanic quartz crystals and plagioclase fragments surrounded by a recrystallized tuffaceous matrix composed of fine-grained quartz, sericite and white micas defining the tectonic foliation in the olistolith. The youngest single zircon age is ca. 435 ± 40 Ma (Telychian), and the magmatic concordia age is 442 ± 22 Ma (three ages, 95 % concordant; Fig. SI-8A) with important age populations defining inherited concordia ages at 471 ± 14 Ma (eight ages), 484 ± 16 Ma (six ages) and 494 ± 13 Ma (nine ages) (Floian to Furongian).

Sample EC-PO-419 was collected in the northern limb of the Verín–Alcañices Synform, to NW of the Bragança Complex in a low-grade foliated medium-grained felsic metatuff representing an olistolith associated with others made of rhyolite, quartzite, lydite, quartzlitharenite and limestone (Fig. 6g) in a mass-wasting slide within the siliciclastic synorogenic sequence (Soutelo Formation in Meireles, 2013). This sample presents a porphyritic texture composed of a highly foliated recrystallized aphanitic matrix made of sericite, fine-grained quartz, white micas, chlorite, biotite, surrounding volcanic quartz and plagioclase phenocrysts and shard fragments. The youngest single zircon age is 439 ± 8 Ma and the magmatic concordia age is 442 ± 12 Ma (five ages, 95 % concordant). Other age populations are Darriwilian (465 ± 11 Ma, six ages) and Floian (474 ± 13 Ma, five ages) (Supplement; Fig. SI-8B).

Sample PET-01 was grabbed at the Spanish–Portuguese border in the Rábano Formation, in an olistolith cluster extending from the northern edge of the Bragança Complex to the northern limb of the Alcañices Syncline. The sampled volcanic body is a grey rhyolitic tuff with disseminated sulfides, which dye the rock with reddish color when weathered. It shows a porphyritic texture made of a roughly foliated recrystallized matrix of sericite, fine-grained quartz, white micas, chlorite retrogressed from biotite and irregular and cubic opaques. The phenocrystals are of sericitized plagioclase and volcanic quartz crystals showing embayments, broken crystals and shards. A concordia age of 494.1 ± 1.1 Ma was attained (lower Furongian) (Fig. SI-9A), supporting an olistolithic nature of this sample, as it was collected in an olistostrome of the synorogenic Rábano Formation (González Clavijo, 2006).

Sample RAB-01, located at the north of the Alcañices village, was hand-picked in a meter-sized block of a weathered intensely foliated rhyolitic tuff representing an olistolith from the lower part of the synorogenic San Vitero Formation. Phenocrysts of plagioclase, sometimes fragmented, and volcanic quartz crystals and shards are surrounded by a fine-grained recrystallized matrix of quartz, sercite, white mica and opaques. Structurally it belongs to a horse with the sheared SCSSs at the base (Fig. 8a). This sample yields a magmatic concordia age of 476.0  ± 1.5 Ma (Floian) (Fig. SI-9B). Its position within an olistostrome stratigraphically higher than the graptolite-rich SCSSs (Fig. 8a) excludes other plausible explanations as an interlayered pyroclastic flow, or a sill (González Clavijo, 2006; Juan Carlos Gutiérrez Marco, Artur Sá, and José Piçarra, personal communication, 2008).

5.2 Published ages from magmatic olistoliths

In the Alcañices synform several olistoliths of felsic volcanic rocks have been identified, all of them of rhyolite to dacite composition, and often forming large clusters elongated NW–SE.

Previous research (González Clavijo et al., 2016) obtained an age of the Nuez olistolith (NUEZ; see location of all samples in the Supplement), one of the major blocks forming a cluster several kilometers long included in an olistostrome inside the synorogenic Rábano Formation, towards the southern limb of the Alcañices synform. This block contains two volcanic facies: dacite lava and dacitic quartz-eyed tuff (Ancochea et al., 1988). The LA-ICP-MS U–Pb isotope analysis of magmatic zircons of a dacitic tuff sample returned a concordant magmatic age of 497 ± 2 Ma (lowermost Furongian).

In the northern limb of the same synform, the Figueruela dacite (COS-8) was dated by SHRIMP-II U–Pb analysis (Farias et al., 2014) yielding a magmatic concordia age of 488.7 ± 3.7 Ma (around the limit Furongian/Tremadocian). This igneous rock was interpreted as a dacitic lava flow interlayered in the Paraño Group of the Schistose Domain or Parautochthon sensu lato. Our field work in the area has revealed that the Figueruela dacite belongs to a major cluster of olistoliths in a large mass wasting deposit, mainly composed of blocks of felsic lavas (dacite and rhyolite) and tuffs, also containing large quartzite and lydite lenses. In our reinterpretation the Figueruela dacite is an olistolith contained in a basal block-in-matrix unit placed below the San Vitero Formation coherent primary unit and above the sheared SCSSs. Thus, here we support that the Figueruela dacite belongs to the synorogenic LPa as previously stated by González Clavijo et al. (2016).

At the northern edge of the Bragança Complex, another significant volcanic body, the Soutelo rhyolite (COS-7), was dated by SHRIMP-II U–Pb analysis (Farias et al., 2014) yielding a 499.8 ± 3.7 Ma (upper Miaolingian) concordia age. The aforementioned authors have included the Soutelo rhyolitic lava in the so-called Paraño Group and considered it as volcanic event in the preorogenic sedimentary sequence of the (Upper) Parautochthon. Our field study disclosed the existence of an important cluster of olistoliths of diverse lithologies such as lydite, grey quartzite, greywacke, limestone, rhyolitic lavas and acidic pyroclastic tuffs, the last two being the most abundant types. Complementarily, this major block-in-matrix unit is placed on top of the mylonitized SCSSs that defines the BLPD. For all these reasons we consider that the Soutelo rhyolite is also an olistolith inside the synorogenic LPa.

5.3 The possible sources of the magmatic olistoliths are in the UPa

The UPa unit, structurally below the Lower Allochthon as defined by Dias da Silva et al. (2014) in the eastern rim of the Morais Complex (Figs. 2 and 3), contains a late Cambrian to Silurian detrital sequence with minor limestones and voluminous volcanism (Pereira et al., 2000, 2006). The main volcanic events are, from bottom to top, the Mora felsic to mafic volcanic rocks (Mora Volcanics; Dias da Silva, 2014; Díez-Montes et al., 2015); the Saldanha gneiss (Ribeiro, 1974; Ribeiro and Ribeiro, 2004; Pereira et al., 2006, 2008); and a large felsic and mafic volcanic–sedimentary complex (Volcano Siliceous Complex of Ribeiro, 1974 and Pereira et al., 2006), latterly renamed to Peso Formation (Dias da Silva et al., 2016; Díez-Montes et al., 2015). The Saldanha gneiss is a rhyolitic dome composed of fine- to coarse-grained porphyritic lavas and tuffs, intercalated in the Cambro-Ordovician terrigenous succession below the Armorican Quartzite (Algoso Formation; Dias da Silva et al., 2014), while the Peso Formation lays above it, in the highest stratigraphic positions of the UPa (Dias da Silva et al., 2016).

Previous zircon U–Pb geochronology of some representative bodies of those volcanic rocks (MOR-18ZIR; SAL-1ZIR; PR-1; PR-2, location given in the Supplement) at the eastern fringe of Morais Complex yielded Furongian (Mora Volcanics, 493.5 ± 2 Ma), Tremadocian (Saldanha Volcanics, 484 ± 2.5 Ma) and Upper Ordovician (Peso Formation, 455–460 Ma) magmatic ages (Dias da Silva et al., 2014; 2016), evidencing episodic voluminous vulcanism along the stratigraphic record of the UPa. These ages are coherent with their stratigraphic position, thus supporting that the polyphasic pervasive deformation undergone by the UPa has not disrupted the original sedimentary architecture as the interlayered volcanic rocks keep the putatively primary chronological order (Díez-Montes et al., 2015; Dias da Silva et al., 2016). Although none of these U–Pb analyses included a study of the inherited zircon ages in these rocks, their magmatic ages were crucial to confirm the presence of rocks with ages similar to those found in the olistoliths, pointing to the UPa as the likely source of the Cambro-Ordovician volcanic olistoliths in the LPa (González Clavijo et al., 2016), such as the Nuez, Soutelo (COS-7) and Figueruela (COS-8) hectometer-sized magmatic olistoliths.

The two new samples from the NW Morais Complex UPa (P-381 and P-385, location in the Supplement) were analyzed for magmatic and inherited zircon ages (Supplement; Fig. SI-7). Sample P-381 is an intensely folded rhyolite with quartz phenocrysts in a foliated sericite and white mica matrix. Its youngest single zircon age is 456 ± 9 Ma and the youngest concordia (magmatic) age is 461 ± 4 Ma (Darriwilian; six ages, 95 % interval; Fig. SI-7A). Although it was collected in a similar structural or stratigraphic position, sample P-385 (foliated dacite with quartz, feldspar and plagioclase phenocrysts in a micaceous highly foliated matrix) seems to be older, with the youngest single zircon at 468 ± 13 Ma and a concordia magmatic age of 475 ± 5 Ma (Floian–Arenigian; five ages, 95 % concordant; Fig. SI-7B). Both samples have inherited Tremadocian (ca. 477–485 Ma) and Furongian (ca. 485–497 Ma) concordia ages, which are the main zircon age populations in these rocks. Ages of both rocks confirm the results obtained in previous studies, which attribute a Middle–Upper Ordovician age for the Peso Formation of the UPa (Dias da Silva et al., 2016). These new results also show that the UPa is likely the major source of the 440–475 Ma magmatic olistoliths and of Cambro-Ordovician–Silurian detrital zircons in the LPa synorogenic basin, also supported by the age distribution plots and the MDS diagrams presented below (see also the Supplement).

6 Discussion

6.1 Structural and stratigraphic meaning of the Lower Parautochthon synorogenic basins

The samples studied contain detrital zircons coeval with the Variscan orogeny time span (400–320 Ma). These results permit us to extend the LPa to areas with flyschoid sequences, some of them with broken beds, slump folds, BIMFs, olistostromes and olistoliths, previously not recognized as Variscan synorogenic deposits. From rocks underlying eastern parts of the Bragança and Morais complexes, where the LPa was described by Dias da Silva et al. (2014, 2015, 2020), it may be continued to the south following the GTMZ boundary through the Meirinhos (MEI-ZR-01 and CR-ZR-01) and western Mirandela (MIR-41, AD-PO-49, AD-PO-5 and AD-PO-57) zones to end in the Vila Pouca de Aguiar imbricated thrust system (AD-PO-48B) as shown in Figs. 2, 4 and 5. Farther south, west of Vila Real, in the Marão Range (Fig. 2) some structurally staked units contain sequences described by Pereira (1987, 1989) as likely to be correlated with the LPa, with a lower SCSS unit (Campanhó Formation) overlain by turbiditic sequences (Santos, Canadelo and Mouquim Formations) which have been lithologically correlated to some of the LPa stratigraphic units in the Alcañices Syncline area (González Clavijo, 2006).

To the west of the Alcañices synform, north of the Bragança Complex, three samples of the flyschoid sequence (EC-PO-293, GIM-ZR-01 and AD-PO-66), two new volcanic olistoliths ages (EC-PO-337 and EC-PO-419) and data from another olistolith (COS-7, Farias et al., 2014) support the interpretation of the LPa rocks in several stacked slices as the Variscan synorogenic sequence (Tournaisian age or younger). These slices contain glided blocks of upper Cambrian to Ordovician–Silurian volcanic rocks (Fig. 10a), and also graptolite-rich Silurian lydites (Meireles et al., 1999a, b; Piçarra et al., 2006a, b) which preclude an in situ interpretation for the volcanic rocks. A similar arrangement was unveiled in the Alcañices Synform from two new samples picked in olistoliths (PET-01 and RAB-01) plus two more from the literature (NUEZ, González Clavijo et al., 2016 and COS-8, Farias et al., 2014) yielding upper Cambrian to Ordovician magmatic ages. These rock blocks appear among other Lower Ordovician–Lower Devonian fossil-bearing (meta)sedimentary olistoliths, which occur in a sedimentary unit of turbiditic nature containing Upper Devonian to Mississippian detrital zircon grains.

Figure 9Geological sketch of one area north of the Bragança Complex showing slices of synorogenic units bounded by tectonics mélanges involving the synorogenic and the black Silurian rocks. The multiplex was openly folded during the Late Variscan event (C3). Coordinates system: UTM-WGS84.


Figure 10Field aspects of the block in matrix formations and olistoliths in the LPa. (a) Rhyolite block in the LPa northern Bragança Complex inside a quartzlitharenite; (b) flyschoid sequence in the Verín synform SW limb; (c) thrust band involving black Silurian ampelite, and rhyolite (whitish) and lydite (black) blocks in the NE limb of the Verín synform; (d) thrust band deforming rhyolite (yellowish and whitish), quartzlitharenite (brown in the right side) and ampelite (black) at the NE limb of the Verín synform; (e) black Silurian condensed facies in a tectonic slice inside a tectonic mélange at the NE limb of the Verín synform; (f) flyschoid sequence with sedimentary and load structures in the LPa at the easternmost part of the Verín synform SW limb; (g) olistolith made of UPa rocks displaying the characteristic arrangement of tectonic foliations in the LPa of the easternmost part of the Verín synform SW limb; (h) tectonic foliation array of the UPa lower part at the southwest of the Morais Complex.


The Gimonde Formation exposed in the Alcañices Synform and in the northern edge of the Bragança Complex (together some other local names: Soutelo, Rio de Onor and Coroto in Meireles, 2013) has cartographic continuity with the Nogueira Group at the Verín Synform and, according to our field observations, also with the lower tectonic unit of the Paraño Group underlying a thrust structure positioned at the base of the “Quartzitic” middle unit of the Paraño Group (Marquínez, 1984; Farias, 1990) as displayed in Fig. 2. Also, a coherent primary unit of decimeter-scale interbedded phyllites and greywackes was newly identified in several zones (Fig. 10b), including the presence of thick beds of greenish fine-grained lithic sandstone close to the town of Verín. Block-in-matrix phacoid bodies enclosing native and exotic blocks were also observed (Fig. 10c and d), some of them large enough to be considered olistoliths (mainly made of quartzite, lydite, black limestone, and felsic lava and tuffs). Towards the base, Silurian lydite and ampelite beds are frequent and strongly sheared, thus forming a tectonic mélange involving synorogenic sediments and the SCSSs (Fig. 10e). As in the other nearby sectors, the Verín Synform area can be also envisaged as a complex imbricated thrust system, where the horses repeat the SCSSs and the Nogueira Group. In the eastern part of the SW limb of the Verín Synform in Portugal (Fig. 2 and cross section 3 in Fig. 3), the Nogueira Group is correlated with the Lower Schists Formation (Pereira et al., 2000), which maintains the sedimentary aspects described above, but with coarser-grained (meta)sandstones very rich in angular quartz and plagioclase grains, thus supporting a more proximal volcanic-rich source area. Also, the olistoliths identified in this area within the synorogenic sequence are made of UPa rocks, as they present polyphasic pervasive deformation characteristic of the tectonically overlaying unit (Fig. 10f, g and h), contrasting with the low deformation observed in the rest of the LPa sequence. According to the terminology adopted in this work (Festa et al., 2019, 2020) the basal detachment zone of this stratigraphic unit (possibly the BLPD) is made of imbricated tectonic slices mixing the BIMF deposits and the mylonitized SCSSs, forming a polygenetic mélange.

The revision of the stratigraphy of the Nogueira Group and the lower unit of the Paraño Group in the Verín Synform allow us to propose a new interpretation for rest of the Paraño Group (Quartzitic and Upper units), which occupies the core of the synform. It may be correlated with the UPa according to the following points:

  • i.

    There is cartographic continuity with the UPa exposures that surround the Bragança Complex (Figs. 2, 3 and 4).

  • ii.

    A thrust fault underlying the Quartzitic Middle unit of the Paraño Group has been identified in this work in coincidence with places where Farias (1990) identified protomylonites and crenulation cleavage, which is here interpreted as the MTMT following Dias da Silva (2014).

  • iii.

    The upper and quartzitic units of the Paraño Group are made of low-grade pervasively deformed detrital rocks like the sequence forming the UPa under the Morais and Bragança complexes (Nuño Ortea et al., 1981; Alonso Alonso et al., 1981; Farias, 1990), and they contain volcanic interbedded bodies with similar geochemistry and U–Pb zircon ages (Nuño Ortea et al., 1981; Alonso Alonso et al., 1981; Farias, 1990; Valverde Vaquero et al., 2007).

  • iv.

    The Paraño quartzitic unit (Farias, 1990), which delineates the hinge zone and limbs of the Verín Synform (Fig. 2), shows spatial continuity and can be lithologically correlated with the Algoso Formation in Portugal, which is considered an Armorican-type early Ordovician quartzite in the UPa (Dias da Silva, 2014; Dias da Silva et al., 2016).

  • v.

    A volcanic body placed above the Lower Ordovician quartzite in the Verín Synform, the Navallo trachyte, has a radiometric age of 439.6 ± 5 Ma (uppermost Ordovician to Llandovery; Valverde Vaquero et al., 2007), younger but coherent with that of felsic metavolcanic rocks in the Peso Formation of the UPa around the eastern and northern rim sections of the Morais Complex.

To the north, in the Cabo Ortegal Complex, the Rio Baio thrust sheet (Marcos et al., 2002) is structurally located under the allochthonous units and has been correlated to the Schistose Domain (Pa) below the Órdenes, Bragança and Morais complexes (Farias et al., 1987; Ribeiro et al., 1990; Martínez Catalán et al., 1997). The internal structure of the Rio Baio thrust sheet is complex, deforming a greenschist facies detrital sequence which includes quartzites and volcanic rocks (Arce Duarte and Fernández Tomás, 1976; Arce Duarte et al., 1977; Fernández Pompa and Piera Rodríguez, 1975; Fernández Pompa and Monteserín López, 1976; Marcos and Farias, 1999). Among the latter, the main bodies are the Loiba dacites, Costa Xuncos rhyolites and Queiroga rhyolites (Arenas, 1984, 1988; Ancochea et al., 1988). The stratigraphic sequence of the Rio Baio thrust sheet was considered Silurian by the fossiliferous content of some beds (Matte, 1968; Romariz, 1969; Iglesias and Robardet, 1980; Piçarra et al., 2006b). Nevertheless, a field reappraisal considered those Silurian levels to be placed at the base of the Rio Baio thrust sheet (Valverde Vaquero et al., 2005), following the basal tectonic contact of the Schistose Domain in the Cabo Ortegal Complex. The base of the Rio Baio thrust sheet was detached from the autochthonous CIZ by a fault developed preferably in the SCSSs (here interpreted as the BLPD). Immediately above the BLPD, a low-grade turbiditic sequence is exposed at the coastline, where the PICON-2 sample was collected (Figs. 2, 3, 4 and 5 and the Supplement) giving support for a Variscan synorogenic origin of this sequence and the ascription of the lower part of the Rio Baio thrust sheet to the LPa. Thus, under the Cabo Ortegal Complex a Parautochthon comparable with that of the Morais and Bragança areas exists, with a small deformed thin LPa unit overlying a tectonic mélange developed in the SCSSs, and the upper part of the Rio Baio thrust sheet above, representing the preorogenic UPa. An isotopic age obtained by Valverde Vaquero et al. (2005) in the Queiroga alkaline rhyolite (U–Pb TIMS, 475 ± 2, Ma – Floian) supports this ascription by comparison with felsic volcanic rocks of the UPa around the Morais Complex (Dias da Silva, 2014; Dias da Silva et al., 2014, 2016; this work).

Based on the geochronology results of the 17 magmatic and detrital rock samples presented here and previously published zircon age data we interpret that the LPa Variscan synorogenic sedimentary and structural unit forms a continuous tectonic carpet underlying the GTMZ separating it from the CIZ. The LPa is not observed in some reaches of the limit between the two zones because if it exists, it has been hidden by late Variscan transcurrent faults or due to the intrusion of Variscan granitoids. Between the Cabo Ortegal and the Bragança complexes and in the northern Porto sector the available data from the literature does not conclusively support the existence of synorogenic sequences which could be endorsed to the LPa, and no detrital zircon studies have been performed until now. Only the San Clodio Series may represent a link between the LPa of Cabo Ortegal and Trás-os-Montes, although it is younger and, although imbricated, it is not fully allochthon.

The strongly deformed SCSSs present at the base of every LPa tectonic slice and frequently separated from the synorogenic sequence by a thrust fault rooted in the BLPD must be considered a tectonic mélange in the meaning proposed by Festa et al. (2019, 2020), as it also incorporates tectonic blocks and olistoliths from the base of the synorogenic sequences. So, when the related shear zone incorporates glided blocks (Figs. 8a, b and 9) it could be considered as a polygenetic mélange (Festa et al., 2019, 2020). We envisage the SCSSs as a mixing unit sharing rocks of the LPa and the Autochthon, where the Silurian components were scraped off from the CIZ local uppermost sequence during the emplacement, a mechanism suggested by Ogata et al. (2019), Smeraglia et al. (2019) and Hajna et al. (2019) for mélange formation. This offscraping mechanism is supported by the presence of Silurian rocks in the Autochthon (CIZ) at the eastern part of the Alcañices Synform (González Clavijo, 2006) and at the Sil Syncline area (Martínez Catalán et al., 2016).

6.2 Provenance of the siliciclastic rocks and olistoliths in the Lower Parautochthon

We have used multidimensional scaling (MDS) (ISOPLOT-R by Vermeech, 2018) to compare the new and the already published data on the zircon age populations of the NW Iberian synorogenic basins, with potential sources within the Iberian Massif terranes (see the Supplement for more details, and complete U–Pb age datasets in Tables S18 and S19). This approach has proven successful in the improvement of paleogeographic reconstruction models in the SW Iberian Variscan belt (Pereira et al., 2020a, b) where the zircon age fingerprinting of possible sedimentary sources feeding the Devono-Carboniferous synorogenic marine basins indicate that sediments were derived from the continental margins of Laurussia and Gondwana, and from a “missing” Variscan volcanic arc (Pereira et al., 2012a). In contrast to the SW Iberian case, the potential source areas of sediments and olistoliths of the NW Iberian synorogenic basin can be easily found in the nearby tectonostratigraphic domains (Autochthon and Allochthon). Because synorogenic basins are usually fed from the neighboring surrounding (orogenically active) highs, there is no need to recur to other more distant Variscan sectors like the Laurussian domains of Meguma, Avalonia or Baltica, to accurately determine sedimentary provenance.

In this study, the zircon age data used for estimating possible source areas of NW Iberian synorogenic marine basins are a selection of published detrital zircon U–Pb ages of pre-Upper Devonian siliciclastic rocks of the NW Iberian Autochthon and Allochthon. The data compiled by Puetz (2018) and Stephan et al. (2018) were combined with other more recent literature (Abati et al., 2007; Albert et al., 2015; Díez Fernández et al., 2010, 2012, 2013; Dinis et al., 2012; Fernández-Suárez et al., 2000, 2002, 2003, 2014; Gutiérrez-Alonso et al., 2003, 2015; Linnemann et al., 2008, 2018; Martínez Catalán et al., 2004; Naidoo et al., 2018; Pastor-Galán et al., 2013; Pereira et al., 2011, 2012a, b; Shaw et al., 2014; Talavera et al., 2012, 2015; Teixeira et al., 2011; Zimmermann et al., 2015) (reference samples in Table S18). We have performed a quality test to the U–Pb isotopic data in each sample, recalculating all the zircon ages following the procedure used in our samples (see the Supplement for complete description of the quality test and data selection). This dataset is used to fingerprint the zircon ages in the source areas using MDS and compare them with the extensive data collection of the synorogenic siliciclastic rock samples, which was expanded in this work from 13 to a total of 24 samples. We have also included new zircon age data on volcanic rocks from the UPa (two samples: P-381 and P-385) and from large olistoliths in the LPa (four samples: EC-PO-337; EC-PO-419; PET-1; RAB-1) to compare their age spectra with the detrital zircon samples, thus tracing the source areas for some of the large olistoliths and the flysch sequence.

The age data of the possible source areas were selected according to a conceptual paleogeographic model for the Upper Devonian–early Carboniferous (as provided in Dias da Silva et al., 2015 and Martínez Catalán et al., 2016). Following the reasoning explained in the Supplement, we define Source A samples as representative of the sources eroding from the peripheral bulge developed in the Autochthon (CIZ, WALZ-CZ and OMZ), but also including the UPa slice as one of the most distal section of north-central Gondwana passive margin. Source B reflects the NW Iberian Allochthon (GTMZ), defined as an accretionary complex built along the Devonian and emplaced onto the Autochthon during the early Carboniferous, forming the Parautochthon (UPa and LPa) as the lowermost tectonic sheet. The UPa is included in Source A because it was part of the Autochthon at the beginning of emplacement, although it was later incorporated into the Allochthon. So, Source A could have been at either side of the synorogenic basin margins, belonging to the peripheral bulge in early Variscan times (Late Devonian), and forming the GTMZ basal thrust sheet in the early Carboniferous, as the thrust front moved towards inland Gondwana. Source A samples were grouped by stratigraphic age, considering that the general stratigraphy of the Autochthon and UPa was not substantially imbricated by major Variscan thrust zones at time of its erosion. Opposingly, in Source B, the GTMZ allochthonous complexes were separated according to their tectonometamorphic unit or domain (UA, MA or LA) with all samples ranging from not fully constrained Ediacaran to Lower Ordovician stratigraphic ages and belonging to the Galician allochthonous units (Malpica-Tui, Órdenes and Cabo Ortegal complexes).

Using the MDS diagram with both potential sources as the base of our provenance interpretation, we have plotted the new and published zircon age populations of the synorogenic siliciclastic rocks, adding the new magmatic and inherited ages of the Middle Ordovician–Silurian volcanic rocks collected as olistoliths in the LPa and in the upper stratigraphic units of the UPa.

The final MDS plot (Fig. 11; Supplement for detailed information) shows two main age clusters: Cluster 1 – “Upper Parautochthon Middle Ordovician–Silurian volcanism” characterized by synorogenic sediments with abundant Cambrian and Ordovician zircon grains, a high concentration of Middle to Upper Ordovician ages, and minor amounts of Silurian and Devonian ages (Fig. 12); Cluster 2 – “Multiple Gondwana-derived sources” characterized by populations with a wide variety of sub-clusters (Groups 1 to 7, Figs. 11 and 13; Supplement), including the Autochthon and the UPa (Source A), and the allochthonous complexes (Source B). The defined groups show direct relation with the most probable sources (A and/or B), but they also represent different grades of sediment mixing and recycling (Fig. 11).

Figure 11MDS diagram for the samples studied and the reference populations in the autochthon and allochthon, with the differentiation between Cluster 1 (Fig. 12) and the seven groups of Cluster 2 (Fig. 13). Data to construct this plot are provided in the Tables S18 and S19. See more details on how this diagram was built in the Supplement.


Figure 12Age distribution plots of the Cluster 1 type populations: “Upper Parautochthon Middle Ordovician–Silurian volcanism”. See text and the Supplement for a detailed description.


Figure 13Age distribution plots of all groups belonging to the Cluster 2 populations: “multiple Gondwana-derived sources”. Legend is in Fig. 12. See text and the Supplement for a detailed description.


The MDS diagram (Fig. 11) and age distribution plots (Figs. 12 and 13) suggest that there is not a decipherable pattern in the provenance of sediments, both in time and space. It is possible to recognize provenance changes along and across the same stratigraphic units, sometimes sampled in different beds that are a few centimeters apart (samples MIR-41 and AD-PO-49). Our analysis evidence that zircon provenance varies, with sediments coming from both A and B sources at the same time and/or arriving to different sectors of the synorogenic basin. Mixing patterns between age groups can also be noted, with sediment recycling leading to dilution of sources towards those more typical of Gondwana as, for instance, Group 7 and Cluster 1 showing a trend to Group 5. Reversely, younger samples plot closer to sources attributable to the allochthonous complexes, as shown by the trend of Group 7, from samples SO-14, SO-1 and SO-2 towards the Upper Allochthon reference population, apparently marking the progradation of the allochthonous wedge onto the NW Iberian Autochthon from the Tournaisian to the upper Visean.

These fluctuations suggest variations on the topographic highs surrounding the synorogenic basin at both margins (accretionary complex and peripheral bulge) in the Upper Devonian–early Carboniferous (Fig. 14). The tectonic activity that controls the basin shape and sedimentation was the cause of highly erosive, large-scale mass-wasting processes leading to the large olistolith-bearing BIMF deposits. The synorogenic marine sediments (cohesive flysch and BIMFs) were gradually incorporated at the base of the accretionary wedge as a tectonic carpet, forming polygenic mélanges. The rapid frontal accretion of trench turbidites (Kusky et al., 2020) allowed their fast exhumation, leading to their recycling within the basin (wild flysch).

Figure 14Sketch of the orogenic collision at the Tournaisian–Visean displaying the trench-fill turbidites and the block-in-matrix deposits. The upper part represents the input of the zircon populations from different sources and zones.


6.3 Origin of Variscan detrital zircon grains

A contrasting aspect of NW Iberian synorogenic deposits with their equivalents in SW Iberia (e.g., Pereira et al., 2014, 2020a, b; Rodrigues et al., 2015; Pérez-Cáceres et al., 2017) is the scarcity of Variscan ages in the detrital zircon age populations in our study case. Only 14 detrital zircon samples out of 24 (11 new + 13 from previous works) have a minor population of Variscan zircons, compared with the predominant Variscan zircon populations in most of the synorogenic formations in SW Iberia.

Notwithstanding, synorogenic Variscan zircon grains are represented in all studied formations of the LPa, independently of the MDS age cluster they belong to (the clusters mostly define the “old” zircon age population patterns; Figs. 12 and 13). Because there is no evidence of volcanic activity associated with the synorogenic marine deposits of the LPa, the source of the rare Variscan detrital zircons must be located elsewhere. A possibility is the Upper Devonian–Carboniferous volcanism coeval with the development of the SW Iberian synorogenic basins (Oliveira et al., 2019a, b; Pereira et al., 2020b) or in the Visean ash deposits laying within the CZ condensed marine synorogenic sediments (Merino-Tomé et al., 2017). But because zircon age populations of synorogenic sediments include those identified in the volcanic and sedimentary olistoliths (Figs. 12 and 13), a more local derivation seems reasonable.

To identify the sources for the Variscan zircon grains we must check the main zircon forming events represented in the allochthonous complexes of the GTMZ, and in the underlying Autochthon (CIZ, WALZ, and less probably CZ and OMZ). The oldest Variscan zircon populations in the LPa range from ca. 400 to 380 Ma, which can be related to the Lower-Middle Devonian high-P and high-T metamorphism in the Upper Allochthon (Gómez-Barreiro et al., 2006, 2007; Gómez-Barreiro, 2007). An alternative contribution could be Emsian volcanism in the Autochthon, such as that identified in the CIZ to the south of Trás-os-Montes, in the core of the Tamames Syncline (ca. 395 Ma; Gutiérrez-Alonso et al., 2008). A younger group, between ca. 380–370 Ma may have its origin in metamorphic zircons growth during exhumation of the Upper Allochthon (Martínez Catalán et al., 2016, 2019) and to a small extent in prograde metamorphism of the Middle Allochthon (Pin et al., 2006; Arenas et al., 2007; Arenas and Sánchez-Martínez, 2015; Santos Zalduegui et al., 1996, and references therein). The age group in the range ca. 370–360 Ma may derive from the high-P and high-T metamorphic rocks of the Lower Allochthon (Abati et al., 2010; Díez Fernández et al., 2011; Santos Zalduegui et al., 1995). These early Variscan ages are found only in the allochthonous complexes and never in the CIZ (Source A). Younger zircon sources with ages of ca. 340–320 Ma are associated with low-P and high-T regional tectonometamorphic events and to magmatic pulses at ca. 340, 335 and 320 Ma (Martínez Catalán et al., 2003; Díez Férnandez and Pereira, 2016; Díez Fernández et al., 2017; López-Moro et al., 2017; Gutiérrez-Alonso et al., 2018; Dias da Silva et al., 2018). But zircons of this age interval have been found only in the Almendra Formation, the most external imbricate of the synorogenic deposits of the LPa, and in the San Clodio Series.

7 Conclusions

New results from field and geochronology studies on the Variscan (ca. 400–320 Ma) hinterland synorogenic marine deposits of the Parautochthon in NW Iberia are presented. They surround large parts of the Galicia–Trás-os-Montes Zone (GTMZ), separating the Allochthon and Parautochthon from the structurally underlying Autochthon of the Central Iberian Zone (CIZ). The new data are useful for better understanding the relationship of the NW Iberian Variscan hinterland marine basin with the structurally underlying and overlying units. We show that the two units defined in the Parautochthon at the eastern rim of the Morais and Bragança complexes cover an area larger than previously estimated, being exposed from Cabo Ortegal (NW Spain) to Trás-os-Montes (NE Portugal). The existence of a preorogenic highly folded Upper Parautochthon (UPa) and an imbricated Lower Parautochthon (LPa) composed of slices of Devono-Carboniferous turbidites and tectonically scrapped autochthonous Silurian carbonaceous–siliceous slates (SCSSs) is a general attribute of the whole NW Iberian Parautochthon. Regional tectonic and stratigraphic features show that the LPa represents a synorogenic basin that was gradually incorporated into the base of the allochthonous accretionary wedge while it was being emplaced onto the northern Gondwana margin, forming a continuous tectonic carpet at the base of the GTMZ.

Analysis of the stratigraphic and tectonic features of the synorogenic flysch in the LPa highlights the relevance of large-scale mass wasting deposits in the sequence. Block-in-matrix formations (BIMFs) represent sedimentary mélanges including large olistoliths derived from the accretionary wedge and probably the Autochthon exposed in a forebulge in front of the wedge. They appear surrounded by a chaotic matrix with slump folds and broken beds of the flysch sequence. The BIMFs formed due to gravitational instability triggered by tectonic activity within and at both margins of the basin. These deposits are frequently tectonized, forming imbricated thrust complexes with polygenic mélanges and tectonically scrapped autochthonous SCSSs at the base.

Zircon geochronology of the LPa siliciclastic rocks and of magmatic olistoliths derived from the UPa and the Autochthon constrains the provenance of the sediments and blocks in the LPa. Our study confirms the synorogenic nature of the LPa stratigraphic units, which include Emsian to Serpukhovian detrital zircon grains. These Variscan grains probably derive from nearby sources located in the allochthonous complexes (metamorphic ages of 400–360 Ma), and from Variscan granitoids and migmatites (ca. 340–320 Ma). A comparison of the detrital zircon populations of the NW Iberian synorogenic marine deposits, including the magmatic and inherited ages now obtained in the Middle Ordovician–Silurian volcanic rocks of the UPa (source) and LPa (olistoliths), with a compilation of reference samples from possible source areas, showed evidence of direct source-to-sink relationships of the Variscan hinterland basin with the accretionary complex (GTMZ) and the peripheral bulge in the Autochthon. Multidimensional scaling analysis showed evidence intrabasinal sediment recycling and mixing of sources in time and space. These are explained by (i) the tectonic instabilities within the basin and in its margins, (ii) the migration of the depocenter towards inland Gondwana, and (iii) the gradual incorporation of the synorogenic basin into the accretionary wedge leading to its denudation and recycling.

Data availability

Data are available in the Supplement, uploaded with this paper.


The supplement related to this article is available online at:

Author contributions

EGC, IDS, JRMC, JGB, GGA and ADM participated in the field work and the selection and preparation of the geochronology samples used in this work.

IDS, JMC, GGA, MH, AG and UL contributed to the U–Pb–Th isotopic analysis and age data processing using LA-ICP-MS in the Senckenberg geochronology lab in Dresden.

EGC, IDS, JMC and GGA were responsible for the preparation of the paper, supplements and illustrations, including field and microscope photos.

IDS was responsible for the geochronological analysis of the collected samples, the elaboration of the geochronological database and data quality test, and the statistical analysis and sink-to-source correlation using multidimensional scaling.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “The Iberian Massif in the frame of the European Variscan Belt”. It is not associated with a conference.


We appreciate the work of Kei Ogata (topical editor) and the valuable comments made by Manuel Francisco Pereira and John Wakabayashi (reviewers), which largely increased the quality of the text and figures. We also appreciate the constructive comments of Cristina Accotto that led to a deeper discussion of the data. This work is a contribution to IDL's (Portugal) Research Group 3 (Solid Earth dynamics, hazards, and resources).

Financial support

This research has been supported by the Fundação para a Ciência e a Tecnologia (grant nos. FCT/UIDB/50019/2020-IDL, SFRH/BPD/99550/2014, and Estímulo ao Emprego Científico – Norma Transitória); the Ministerio de Ciencia, Innovación y Universidades (grant no. CGL2016-78560-P); the Instituto Geológico y Minero de España (Plan Cartográfico del IGME); the European Commission Horizon 2020 Framework Programme (SYNTHESYS3 (grant no. DE-TAF-5798)); the Program XII grant of the University of Salamanca; the Spanish Ministry of Science, Innovation and Universities: project IBERCRUST (grant no. PGC2018–096534-B-100); and the Resolution of the Russian Federation Government (agreement no. 14.Y26.31.0012).

Review statement

This paper was edited by Kei Ogata and reviewed by Manuel Francisco Pereira and John Wakabayashi.


Abati, J., Castiñeiras, P., Arenas, R., Fernández-Suárez, J., Gómez Barreiro, J., and Wooden, J. L.: Using SHRIMP zircon dating to unravel tectonothermal events in arc environments. The early Palaeozoic arc of NW Iberia revisited, Terra Nova, 19, 432–439,, 2007. 

Abati, J., Gerdes, A., Fernández Suárez, J., Arenas, R., Whitehouse, M. J., and Díez Fernández, R.: Magmatism and early-Variscan continental subduction in the northern Gondwana margin recorded in zircons from the basal units of Galicia, NW Spain, Geol. Soc. Am. Bull., 122, 219–235,, 2010. 

Albert, R., Arenas, R., Gerdes, A., Sánchez Martínez, S., and Marko, L.: Provenance of the HP-HT subducted margin in the Variscan belt (Cabo Ortegal Complex, NW Iberian Massif), J. Metamorph. Geol., 33, 959–979,, 2015. 

Aldaya, F., Carls, P., Martínez-García, E., and Quiroga, J. L.: Nouvelles précisions sur la série de San Vitero (Zamora, nord-oueste de l'Espagne), C. R. Acad. Sci. Paris, 283, 881–883, 1976. 

Alonso, J. L., Marcos, A., Villa, E., Súarez, A., Merino-Tomé, O. A., and Fernández, L. P.: Mélanges and other types of block-in-matrix formations in the Cantabrian Zone (Variscan Orogen, northwest Spain): origin and significance, Int. Geol. Rev., 57, 563–580,, 2015. 

Alonso Alonso, J. L., Delgado Gutiérrez, C., Zubieta Freira, J. M., Pérez Rojas, A., Ferragne, A., and Rúiz Díaz, C.: Mapa Geológico de España, E. 1:50.000 y Memoria, Hoja 265-Laza, Segunda Serie, IGME, Madrid, Spain, 26 pp., 1981. 

Ancochea, E., Arenas, R., Brãdle, J. L., Peinado, M., and Sagredo, J.: Caracterización de las rocas metavolcánicas silúricas del Noroeste del Macizo Ibérico, Geosci. Aveiro, 3, 23–34, 1988. 

Antona, J. F. and Martínez Catalán, J. R.: Interpretación de la Formación San Vitero en relación con la Orogenia Hercínica, Cuad. Lab. Xeol. Laxe, 15, 257–269, 1990. 

Arce Duarte, J. M. and Fernández Tomás, J.: Mapa Geológico de España. Escala 1:50.000, Hoja 8 (7–3), Vivero. Segunda Serie, Mapa y Memoria, IGME, Madrid, Spain, 45 pp., 1976. 

Arce Duarte, J. M., Fernández Tomás, J., and Monteserín López, V.: Mapa Geológico de España. Escala 1:50.000, Hoja 7 (7–2), Cillero. Mapa y Memoria, Instituto Geológico Minero España, Madrid, Spain, 47 pp., 1977. 

Arenas, R.: Características y significado del volcanismo ordovícico-silúrico de la serie autóctona envolvente del Macizo de Cabo Ortegal (Galicia, NW España), Revista Materiales Procesos Geológicos, 11, 135–144, 1984. 

Arenas, R.: Evolución petrológica y geoquímica de la Unidad alóctona inferior del Complejo metamórfico básico-ultrabásico de Cabo Ortegal (Unidad de Moeche) y del Silúrico paraautóctono, Cadena Herciniana Ibérica (NW de España), Corpus Geologicum Gallaeciae, 4, 545 pp., 1988. 

Arenas, R. and Sanchez-Martinez, S.: Variscan ophiolites in NW Iberia: Tracking lost Paleozoic oceans and the assembly of Pangea, Episodes, 38, 315–333,, 2015. 

Arenas, R., Pascual, F. J. R., Garcia, F. D., and Catalan, J. R. M.: High-pressure micro-inclusions and development of an inverted metamorphic gradient in the Santiago Schists (Ordenes Complex, NW Iberian Massif, Spain): evidence of subduction and syncollisional decompression, J. Metamorph. Geol., 13, 141–164,, 1995. 

Arenas, R., Abati, J., Catalán, J. R. M., García, F. D., and Pascual, F. J. R.: P-T evolution of eclogites from the Agualada Unit (Ordenes Complex, northwest Iberian Massif, Spain): Implications for crustal subduction, Lithos, 40, 221–242,, 1997. 

Arenas, R., Martínez Catalán, J. R., Sánchez Martínez, S., Díaz García, F., Abati, J., Fernández-Suárez, J., Andonaegui, P., and Gómez Barreiro, J.: Paleozoic ophiolites in the Variscan suture of Galicia (northwest Spain): Distribution, characteristics and meaning, in: 4D Framework of Continental Crust, edited by: Hatcher Jr., R. D., Carlson, M. P., McBride, J. H., and Martínez Catalán, J. R., Geological Society of America, Boulder, Colorado, USA, 425–444,, 2007. 

Azor, A., Dias da Silva, Í., Gómez Barreiro, J., González-Clavijo, E., Martínez Catalán, J. R., Simancas, J. F., Martínez Poyatos, D., Pérez-Cáceres, I., González Lodeiro, F., Expósito, I., Casas, J. M., Clariana, P., García-Sansegundo, J., and Margalef, A.: Deformation and Structure, in: The Geology of Iberia: A Geodynamic Approach: Volume 2: The Variscan Cycle, edited by: Quesada, C. and Oliveira, J. T., Springer International Publishing, Cham, UK,, 2019. 

Ballèvre, M., Martínez Catalán, J. R., López Carmona, A., Pitra, P., Abati, J., Díez Fernández, R., Ducassou, C., Arenas, R., Bosse, V., Castiñeiras, P., Fernández-Suárez, J., Gómez Barreiro, J., Paquette, J. L., Peucat, J. J., Poujol, M., Ruffet, G., and Sánchez Martínez, S.: Correlation of the nappe stack in the Ibero-Armorican arc across the Bay of Biscay: a joint French-Spanish project, Geol. Soc. Spec. Publ., 405, 77–113,, 2014. 

Barrera Morate, J. L., González Lodeiro, F., Martín Parra, L. M., del Olmo Sanz, A., Farias Arquer, P., Marquínez García, J., Martínez Catalán, J. R., de Pablo Macía, J. G., and Rodríguez Fernández, L. R.: Mapa Geológico Nacional E. 1:200 000, 17–27, Orense–Verín, Instituto Geológico y Minero de España, Madrid, Spain, 1989. 

Braid, J. A., Murphy, J. B., Quesada, C., and Mortensen, J.: Tectonic escape of a crustal fragment during the closure of the Rheic Ocean: U-Pb detrital zircon data from the Late Palaeozoic Pulo do Lobo and South Portuguese zones, southern Iberia, J. Geol. Soc., 168, 383–392,, 2011. 

Bütler, E., Winkler, W., and Guillong, M.: Laser ablation U/Pb age patterns of detrital zircons in the Schlieren Flysch (Central Switzerland): new evidence on the detrital sources, Swiss J. Geosci., 104, 225,, 2011. 

Chiocci, F. L. and Casalbore, D.: Reprint of Unexpected fast rate of morphological evolution of geologically-active continental margins during Quaternary: Examples from selected areas in the Italian seas, Mar. Petrol. Geol., 87, 148–156,, 2017. 

Coleman, J. M. and Prior, D. B.: Mass Wasting on Continental Margins, Annu. Rev. Earth Pl. Sc., 16, 101–119,, 1988. 

Coke, C. and Gutiérrez Marco, J. C.: Braquiópodos Linguliformea del Ordovícico Inferior de la Serra do Marão (Zona Centroibérica, N de Portugal), Boletín Geológico y Minero, 112, 33–50, 2001. 

Dallmeyer, R. D., Martínez Catalán, J. R., Arenas, R., Gil Ibarguchi, J. I., Gutiérrez Alonso, G., Farias, P., Aller, J., and Bastida, F.: Diachronous Variscan tectonothermal activity in the NW Iberian Massif: Evidence from 40Ar/39Ar dating of regional fabrics, Tectonophysics, 277, 307–337,, 1997. 

DeCelles, P. G.: Foreland basin systems revisited: Variations in response to tectonic settings, in: Tectonics of sedimentary basins: Recent advances, edited by: Busby, C. and Azor, A., John Willey & Sons, Chichester, UK, 405–426,, 2012. 

Dias da Silva, Í.: Geología de las Zonas Centro Ibérica y Galicia – Trás-os-Montes en la parte oriental del Complejo de Morais, Portugal/España, Instituto Universitario de Geología “Isidro Parga Pondal” , Área de Xeoloxía e Minería do Seminario de Estudos Galegos, Coruña, Spain, 2014. 

Dias da Silva, Í., Valverde-Vaquero, P., González-Clavijo, E., Díez-Montes, A., and Martínez Catalán, J. R.: Structural and stratigraphical significance of U-Pb ages from the Mora and Saldanha volcanic complexes (NE Portugal, Iberian Variscides), Geol. Soc. Spec. Publ., 405, 115–135,, 2014. 

Dias da Silva, Í., Linnemann, U., Hofmann, M., González-Clavijo, E., Díez-Montes, A., and Martínez Catalán, J. R.: Detrital zircon and tectonostratigraphy of the Parautochthon under the Morais Complex (NE Portugal): implications for the Variscan accretionary history of the Iberian Massif, J. Geol. Soc., 172, 45–61,, 2015. 

Dias da Silva, Í., Díez Fernández, R., Díez-Montes, A., González Clavijo, E., and Foster, D. A.: Magmatic evolution in the N-Gondwana margin related to the opening of the Rheic Ocean – evidence from the Upper Parautochthon of the Galicia-Trás-os-Montes Zone and from the Central Iberian Zone (NW Iberian Massif), Int. J. Earth Sci., 105, 1127–1151,, 2016. 

Dias da Silva, Í., Pereira, M. F., Silva, J. B., and Gama, C.: Time-space distribution of silicic plutonism in a gneiss dome of the Iberian Variscan Belt: The Évora Massif (Ossa-Morena Zone, Portugal), Tectonophysics, 747/748, 298–317,, 2018. 

Dias da Silva, Í., González Clavijo, E., and Díez-Montes, A.: The collapse of the Variscan belt: a Variscan lateral extrusion thin-skinned structure in NW Iberia, Int. Geol. Rev., 63, 659–695,, 2021. 

Dickinson, W. R. and Valloni, R.: Plate settings and provenance of sands in modern ocean basins, Geology, 8, 82–86,<82:psapos>;2, 1980. 

Díez Fernández, R. and Pereira, M. F.: Extensional orogenic collapse captured by strike-slip tectonics: Constraints from structural geology and UPb geochronology of the Pinhel shear zone (Variscan orogen, Iberian Massif), Tectonophysics, 691, 290–310,, 2016. 

Díez Fernández, R., Martínez Catalán, J. R., Gerdes, A., Abati, J., Arenas, R., and Fernández-Suárez, J.: U-Pb ages of detrital zircons from the Basal allochthonous units of NW Iberia: Provenance and paleoposition on the northern margin of Gondwana during the Neoproterozoic and Paleozoic, Gondwana Res., 18, 385–399,, 2010. 

Díez Fernández, R., Martínez Catalán, J. R., Arenas, R., and Abati, J.: Tectonic evolution of a continental subduction-exhumation channel: Variscan structure of the basal allochthonous units in NW Spain, Tectonics, 30, 1–22,, 2011. 

Díez Fernández, R., Martínez Catalán, J. R., Arenas, R., Abati, J., Gerdes, A., and Fernández-Suárez, J.: U-Pb detrital zircon analysis of the lower allochthon of NW Iberia: age constraints, provenance and links with the Variscan mobile belt and Gondwanan cratons, J. Geol. Soc., 169, 655–665,, 2012. 

Díez Fernández, R., Foster, D., Gómez Barreiro, J., and Alonso-García, M.: Rheological control on the tectonic evolution of a continental suture zone: the Variscan example from NW Iberia (Spain), Int. J. Earth Sci., 102, 1305–1319,, 2013. 

Díez Fernández, R., Parra, L. M. M., and Rubio Pascual, F. J.: Extensional flow produces recumbent folds in syn-orogenic granitoids (Padrón migmatitic dome, NW Iberian Massif), Tectonophysics, 703/704, 69–84,, 2017. 

Díez-Montes, A., González Clavijo, E., Dias da Silva, Í., Gómez Barreiro, J., Martínez Catalán, J. R., and Castiñeiras, P.: Geochemical Evolution of Volcanism during the Upper Cambrian-Ordovician Extension of the North Gondwana Margin, X Congresso Ibérico de Geoquímica/XVIII Semana de Geoquímica, 17–19 October 2015, Alfragide, Portugal, 54–57, 2015. 

Dinis, P., Andersen, T., Machado, G., and Guimarães, F.: Detrital zircon U-Pb ages of a late-Variscan Carboniferous succession associated with the Porto-Tomar shear zone (West Portugal): Provenance implications, Sediment. Geol., 273/274, 19–29,, 2012. 

Ducassou, C., Poujol, M., Ruffet, G., Bruguier, O., and Ballèvre, M.: Relief variation and erosion of the Variscan belt: detrital geochronology of the Palaeozoic sediments from the Mauges Unit (Armorican Massif, France), Geol. Soc. Spec. Publ., 405, 137–167,, 2014. 

Eyles, N.: Marine debris flows: Late precambrian “tillites” of the Avalonian-Cadomian orogenic belt, Palaeogeogr. Palaeocl., 79, 73–98,, 1990. 

Farias, P.: La Geología de la región del sinforme de Verín (Cordillera Herciniana, NW de España), Instituto Universitario de Geología “Isidro Parga Pondal”, Área de Xeoloxía e Minería do Seminario de Estudos Galegos, La Coruña, Spain, 1990. 

Farias, P., Gallastegui, G., González-Lodeiro, F., Marquínez, J., Martín Parra, L. M., Martínez Catalán, J. R., de Pablo Maciá, J. G., and Rodríguez Fernández, L. R.: Aportaciones al conocimiento de la litoestratigrafía y estructura de Galicia Central, Memórias da Faculdade de Ciências da Universidade do Porto, 1, 411–431, 1987. 

Farias, P., Casado, B. O., Marcos, A., Ordóñez, A. R., and Fanning, C.: U-Pb zircon SHRIMP evidence for Cambrian volcanism in the Schistose Domain within the Galicia-Trás-os. Montes Zone (Variscan Orogen, NW Iberian Peninsula), Geol. Acta, 12, 209–218,, 2014. 

Fernández Pompa, F. and Monteserín López, V.: Mapa Geológico de España. Escala 1:50.000, Hoja 7 (6-3), Cedeira. Mapa y Memoria, Instituto Geológico Minero España, Madrid, Spain, 73 pp., 1976. 

Fernández Pompa, F. and Piera Rodríguez, T.: Mapa Geológico de España. Escala 1:50.000, Hoja 22 (6–4), Puentedeume, Mapa y Memoria, Instituto Geológico Minero España, Madrid, Spain, 45 pp., 1975. 

Fernández-Suárez, J., Gutiérrez-Alonso, G., Jenner, G. A., and Tubrett, M. N.: New ideas on the Proterozoic-Early Palaeozoic evolution of NW Iberia: insights from U-Pb detrital zircon ages, Precambrian Res., 102, 185–206,, 2000. 

Fernández-Suárez, J., Gutiérrez-Alonso, G., Cox, R., and Jenner, G. A.: Assembly of the Armorican microplate: a strike-slip terrane delivery? Evidence from U-Pb ages of detrital zircons, J. Geol., 110, 619–626,, 2002. 

Fernández-Suárez, J., Díaz García, F., Jeffries, T. E., Arenas, R., and Abati, J.: Constraints on the provenance of the uppermost allochthonous terrane of the NW Iberian Massif: inferences from detrital zircon U-Pb ages, Terra Nova, 15, 138–144,, 2003. 

Fernández-Suárez, J., Gutiérrez-Alonso, G., Pastor-Galán, D., Hofmann, M., Murphy, J. B., and Linnemann, U.: The Ediacaran – Early Cambrian detrital zircon record of NW Iberia: possible sources and paleogeographic constraints, Int. J. Earth Sci., 103, 1335–1357,, 2014. 

Festa, A., Ogata, K., Pini, G. A., Dilek, Y., and Alonso, J. L.: Origin and significance of olistostromes in the evolution of orogenic belts: A global synthesis, Gondwana Res., 39, 180–203,, 2016. 

Festa, A., Pini, G. A., Ogata, K., and Dilek, Y.: Diagnostic features and field-criteria in recognition of tectonic, sedimentary and diapiric mélanges in orogenic belts and exhumed subduction-accretion complexes, Gondwana Res., 74, 7–30,, 2019. 

Festa, A., Ogata, K., and Pini, G. A.: Polygenetic mélanges: a glimpse on tectonic, sedimentary and diapiric recycling in convergent margins, J. Geol. Soc., 177, 551–561,, 2020. 

Franke, W. and Engel, W.: Synorogenic sedimentation in the Variscan Belt of Europe, B. Soc. Géol. Fr., II, 25–33,, 1986. 

Garzanti, E., Doglioni, C., Vezzoli, G., and Andò, S.: Orogenic Belts and Orogenic Sediment Provenance, J. Geol., 115, 315–334,, 2007. 

Gil Ibarguchi, J. I.: Petrology of jadeite metagranite and associated orthogneiss from the Malpica-Tuy allochthon (Northwest Spain), Eur. J. Mineral., 7, 403–416,, 1995. 

Gil Ibarguchi, J. I. and Dallmeyer, R. D.: Hercynian blueschist metamorphism in North Portugal: tectonothermal implications, J. Metamorph. Geol., 9, 539–549,, 1991. 

Gómez Barreiro, J.: La Unidad de Fornás: Evolución tectonometamórfica del SO del Complejo de Órdenes, Instituto Universitario de Geología “Isidro Parga Pondal”, Área de Xeoloxía e Minería do Seminario de Estudos Galegos, Coruña, Spain, 2007. 

Gómez Barreiro, J., Wijbrans, J. R., Castiñeiras, P., Martínez Catalán, J. R., Arenas, R., Díaz García, F., and Abati, J.: 40Ar/39Ar laserprobe dating of mylonitic fabrics in a polyorogenic terrane of NW Iberia, J. Geol. Soc., 163, 61–73,, 2006. 

Gómez Barreiro, J., Martínez Catalán, J. R., Arenas, R., Castiñeiras, P., Abati, J., Díaz García, F., and Wijbrans, J. R.: Tectonic evolution of the upper allochtonon of the Órdenes complex (Northwestern Iberian Massif): Structural constraints to a poligenic peri-Gondwanan terrane, Geol. Soc. Am. Spec. Pap., 423, 315–332,, 2007. 

González Clavijo, E.: La Geología del sinforme de Alcañices, Oeste de Zamora, Instituto Universitario de Geología “Isidro Parga Pondal”, Área de Xeoloxía e Minería do Seminario de Estudos Galegos, La Coruña, Spain, 2006. 

González Clavijo, E. and Martínez Catalán, J. R.: Stratigraphic record of preorogenic to synorogenic sedimentation, and tectonic evolution of imbricate thrusts in the Alcañices synform (northwestern Iberian Massif), in: Variscan Appalachian Dynamics: The building of the Late Palaeozoic Basement, edited by: Martínez Catalán, J. R., Hatcher Jr., R. D., Arenas, R., and Días García, F., Geological Society of America, Boulder, USA,, 2002. 

González Clavijo, E., Gutiérrez-Marco, J. C., Jiménez Fuentes, E., Moro Benito, M. C., and Storch, P.: Graptolitos silúricos del Sinforme de Alcañices (prov. de Zamora, Zonas Centroibérica y Galaico-Trasmontana) XIII Reunión Internacional Proyecto 351 PICG Paleozoico Inferior del Noroeste de Gondwana, Libro de resúmenes y excursiones, Sociedad Española de Paleontología, Coruña, Spain, 71–74, 1997. 

González-Clavijo, E., Martínez-Catalán, J. R., Meireles, C., Belousovea, E., and Saeed, E.: Detrital zircon U-Pb ages in synorogenic deposits of the internal zones of the Variscan massif and their significance in the orogenic evolution, Géologie de la France, 1, 122–122, BRGM et la SGF, Orléans-Paris, France, 2012. 

González Clavijo, E., Dias da Silva, Í. F., Gutiérrez-Alonso, G., and Díez Montes, A.: U/Pb age of a large dacitic block locked in an Early Carboniferous synorogenic mélange in the Parautochthon of NW Iberia: New insights on the structure/sedimentation Variscan interplay, Tectonophysics, 681, 159–169,, 2016. 

Gutiérrez-Alonso, G., Fernández-Suárez, J., Jeffries, T. E., Jenner, G. A., Tubrett, M. N., Cox, R., and Jackson, S. E.: Terrane accretion and dispersal in the northern Gondwana margin, An Early Paleozoic analogue of a long-lived active margin, Tectonophysics, 365, 221–232,, 2003. 

Gutiérrez-Alonso, G., Murphy, J. B., Fernández-Suárez, J., and Hamilton, M. A.: Rifting along the northern Gondwana margin and the evolution of the Rheic Ocean: A Devonian age for the El Castillo volcanic rocks (Salamanca, Central Iberian Zone), Tectonophysics, 461, 157–165,, 2008. 

Gutiérrez-Alonso, G., Fernández-Suárez, J., Pastor-Galán, D., Johnston, S. T., Linnemann, U., Hofmann, M., Shaw, J., Colmenero, J. R., and Hernández, P.: Significance of detrital zircons in Siluro-Devonian rocks from Iberia, J. Geol. Soc., 172, 309–322,, 2015. 

Gutiérrez-Alonso, G., Fernández-Suárez, J., López-Carmona, A., and Gärtner, A.: Exhuming a cold case: The early granodiorites of the northwest Iberian Variscan belt – A Visean magmatic flare-up?, Lithosphere, 10, 194–216,, 2018. 

Gutiérrez-Alonso, G., López-Carmona, A., Núñez-Guerrero, E., Martínez García, A., Fernández-Suárez, J., Pastor-Galán, D., Gutiérrez-Marco, J. C., Bernárdez, E., Colmenero, J. R., Hofmann, M., and Linnemann, U.: Neoproterozoic-paleozoic detrital sources in the Variscan foreland of northern Iberia: primary v. recycled sediments, Geol. Soc. Spec. Publ., 503, 563–588,, 2020. 

Gutiérrez-Marco, J. C., Piçarra, J. M., Meireles, C. A., Cózar, P., García-Bellido, D. C., Pereira, Z., Vaz, N., Pereira, S., Lopes, G., Oliveira, J. T., Quesada, C., Zamora, S., Esteve, J., Colmenar, J., Bernárdez, E., Coronado, I., Lorenzo, S., Sá, A. A., Dias da Silva, Í., González-Clavijo, E., Díez-Montes, A., and Gómez-Barreiro, J.: Early Ordovician–Devonian Passive Margin Stage in the Gondwanan Units of the Iberian Massif, in: The Geology of Iberia: A Geodynamic Approach: Volume 2: The Variscan Cycle, edited by: Quesada, C. and Oliveira, J. T., Springer International Publishing, Cham, UK,, 2019. 

Hajná, J., Žák, J., Ackerman, L., Svojtka, M., and Pašava, J.: A giant late Precambrian chert-bearing olistostrome discovered in the Bohemian Massif: A record of Ocean Plate Stratigraphy (OPS) disrupted by mass-wasting along an outer trench slope, Gondwana Res., 74, 173–188,, 2019. 

Iglesias, M. and Robardet, M.: Silúrico de Galicia Media-Central. Su importancia en la paleogeografía Varisca, Cuadernos Laboratorio Xeolóxico Laxe, 1, 99–115, 1980. 

Jorge, R. C. G. S., Fernandes, P., Rodrigues, B., Pereira, Z., and Oliveira, J. T.: Geochemistry and provenance of the Carboniferous Baixo Alentejo Flysch Group, South Portuguese Zone, Sediment. Geol., 284/285, 133–148,, 2013. 

Keller, M., Bahlburg, H., and Reuther, C.-D.: The transition from passive to active margin sedimentation in the Cantabrian Mountains, Northern Spain: Devonian or Carboniferous?, Tectonophysics, 461, 414–427,, 2008. 

Krastel, S., Li, W., Urlaub, M., Georgiopoulou, A., Wynn, R. B., Schwenk, T., Stevenson, C., and Feldens, P.: Mass wasting along the NW African continental margin, Geol. Soc. Spec. Publ., 477, 151–167,, 2019. 

Kusky, T., Wang, J., Wang, L., Huang, B., Ning, W., Fu, D., Peng, H., Deng, H., Polat, A., Zhong, Y., and Shi, G.: Mélanges through time: Life cycle of the world's largest Archean mélange compared with Mesozoic and Paleozoic subduction-accretion-collision mélanges, Earth-Sci. Rev., 209, 103303,, 2020. 

Liang, X.-Q. and Li, X.-H.: Late Permian to Middle Triassic sedimentary records in Shiwandashan Basin: Implication for the Indosinian Yunkai Orogenic Belt, South China, Sediment. Geol., 177, 297–320,, 2005. 

Linnemann, U., Pereira, M. F., Jeffries, T. E., Drost, K., and Gerdes, A.: The Cadomian Orogeny and the opening of the Rheic Ocean: the diachrony of geotectonic processes constrained by LA-ICP-MS U-Pb zircon dating (Ossa-Morena and Saxo-Thuringian Zones, Iberian and Bohemian Massifs), Tectonophysics, 461, 21–43,, 2008. 

Linnemann, U., Ouzegane, K., Drareni, A., Hofmann, M., Becker, S., Gärtner, A., and Sagawe, A.: Sands of West Gondwana: An archive of secular magmatism and plate interactions: A case study from the Cambro-Ordovician section of the Tassili Ouan Ahaggar (Algerian Sahara) using U/Pb LA-ICP-MS detrital zircon ages, Lithos, 123, 188–203,, 2012. 

Linnemann, U., Gerdes, A., Hofmann, M., and Marko, L.: The Cadomian Orogen: Neoproterozoic to Early Cambrian crustal growth and orogenic zoning along the periphery of the West African Craton – Constraints from U-Pb zircon ages and Hf isotopes (Schwarzburg Antiform, Germany), Precambrian Res., 244, 236–278,, 2014. 

Linnemann, U., Pidal, A. P., Hofmann, M., Drost, K., Quesada, C., Gerdes, A., Marko, L., Gärtner, A., Zieger, J., Ulrich, J., Krause, R., Vickers-Rich, P., and Horak, J.: A  565 Ma old glaciation in the Ediacaran of peri-Gondwanan West Africa, Int. J. Earth Sci., 107, 885–911,, 2018. 

López-Carmona, A., Abati, J., and Reche, J.: Petrologic modeling of chloritoid-glaucophane schists from the NW Iberian Massif, Gondwana Res., 17, 377–391,, 2010. 

López-Carmona, A., Abati, J., Pitra, P., and Lee, J. W.: Retrogressed lawsonite blueschists from the NW Iberian Massif: P-T-t constraints from thermodynamic modelling and 40Ar/39Ar geochronology, Contrib. Mineral. Petr., 167, 1–20,, 2014. 

López-Moro, F. J., Romer, R., López-Plaza, M., and Gónzalez Sánchez, M.: Zircon and allanite U-Pb ID-TIMS ages of vaugnerites from the Calzadilla pluton, Salamanca (Spain): dating mantle-derived magmatism and post-magmatic subsolidus overprint, Geol. Acta, 15, 395–408,, 2017. 

Marcos, A. and Farias, P.: La estructura de las láminas inferiores del Complejo de Cabo Ortegal y su autóctono relativo (Galícia, NO España), Trabajos de Geología, 21, 53–72, 1999. 

Marcos, A. and Pulgar, J.: An approach to the tectonostratigraphic evolution of the Cantabrian foreland thrust and fold belt, Hercynian Cordillera of NW Spain, Neues Jahrbuch fuer Geologie und Palaeontologie, 163, 256–260, 1982. 

Marcos, A., Farias, P., Galán, G., Fernández, F., and Llana-Fúnez, S.: Tectonic framework of the Cabo Ortegal Complex: A slab of lower crust exhumed in the Variscan orogen (northwestern Iberian Peninsula), Geol. Soc. Am. Spec. Pap., 364, 143–162,, 2002. 

Marquínez García, J. L.: La geología del área esquistosa de Galicia Central (Cordillera Herciniana, NW de España), Instituto Geológico y Minero de España, Madrid, Spain, 1984. 

Martínez Catalán, J. R., Arenas, R., Díaz García, F., and Abati, J.: Variscan accretionary complex of northwest Iberia: Terrane correlation and succession of tectonothermal events, Geology, 25, 1103–1106,<1103:vaconi>;2, 1997. 

Martínez Catalán, J. R., Arenas, R., and Díez Balda, M. A.: Large extensional structures developed during the emplacement of a crystalline thrust sheet: the Mondoñedo nappe (NW Spain), J. Struct. Geol., 25, 1815–1839,, 2003. 

Martínez Catalán, J. R., Fernández-Suárez, J., Jenner, G. A., Belousova, E., and Díez Montes, A.: Provenance constraints from detrital zircon U-Pb ages in the NW Iberian Massif: implication for Palaeozoic plate configuration and Variscan evolution, J. Geol. Soc., 161, 463–476,, 2004. 

Martínez Catalán, J. R., Arenas, R., Díaz García, F., González Cuadra, P., Gómez Barreiro, J., Abati, J., Castiñeiras, P., Fernández-Suárez, J., Sánchez Martínez, S., Andonaegui, P., González Clavijo, E., Díez Montes, A., Rubio Pascual, F., and Valle Aguado, B.: Space and time in the tectonic evolution of the northwestern Iberian Massif: Implications for the Variscan belt, in: 4-D Framework of Continental Crust, edited by: Hatcher Jr., R. D., Carlson, M. P., McBride, J. H., and Martínez Catalán, J. R., Geologic Society of America, Boulder, USA,, 2007. 

Martínez Catalán, J. R., Fernández-Suárez, J., Meireles, C., González Clavijo, E., Belousova, E., and Saeed, A.: U-Pb detrital zircon ages in sinorogenic deposits of the NW Iberian Massif (Variscan belt): Interplay of Devonian-Carboniferous sedimentation and thrust tectonics, J. Geol. Soc., 165, 687–698,, 2008. 

Martínez Catalán, J. R., Arenas, R., Abati, J., Sánchez Martínez, S., Díaz García, F., Fernández-Suárez, J., González Cuadra, P., Castiñeiras, P., Gómez Barreiro, J., Díez Montes, A., González Clavijo, E., Rubio Pascual, F., Andonaegui, P., Jeffries, T. E., Alcock, J. E., Díez Fernández, R., and López Carmona, A.: A rootless suture and the loss of the roots of a mountain chain: The Variscan Belt of NW Iberia, C. R. Geosci., 341, 114–126,, 2009. 

Martínez Catalán, J. R., Rubio Pascual, F. J., Montes, A. D., Fernández, R. D., Barreiro, J. G., Dias Da Silva, Í., Clavijo, E. G., Ayarza, P., and Alcock, J. E.: The late Variscan HT/LP metamorphic event in NW and Central Iberia: relationships to crustal thickening, extension, orocline development and crustal evolution, Geol. Soc. Spec. Publ., 405, 225–247,, 2014. 

Martínez Catalán, J. R., González Clavijo, E., Meireles, C., Díez Fernández, R., and Bevis, J.: Relationships between syn-orogenic sedimentation and nappe emplacement in the hinterland of the Variscan belt in NW Iberia deduced from detrital zircons, Geol. Mag., 153, 38–60,, 2016. 

Martínez Catalán, J. R., Gómez Barreiro, J., Dias da Silva, Í., Chichorro, M., López-Carmona, A., Castiñeiras, P., Abati, J., Andonaegui, P., Fernández-Suárez, J., González Cuadra, P., and Benítez-Pérez, J. M.: Variscan Suture Zone and Suspect Terranes in the NW Iberian Massif: Allochthonous Complexes of the Galicia-Trás os Montes Zone (NW Iberia), in: The Geology of Iberia: A Geodynamic Approach: Volume 2: The Variscan Cycle, edited by: Quesada, C. and Oliveira, J. T., Springer International Publishing, Cham, UK,, 2019. 

Martínez García, E.: El Silúrico de San Vitero (Zamora), Comparación con las series vecinas e importancia orogénica, Act. Geol. Hisp., 4, 104–108, 1972. 

Matte, P.: La structure de la virgation hercynienne de Galice (Espagne), Rev. Géol. Alp., 44, 1–128, 1968. 

Meinhold, G., Morton, A. C., Fanning, C. M., Frei, D., Howard, J. D., Phillips, R. J., Strogen, D., and Whitham, A. G.: Evidence from detrital zircons for recycling of Mesoproterozoic and Neoproterozoic crust recorded in Paleozoic and Mesozoic sandstones of southern Libya, Earth Planet. Sci. Lett., 312, 164–175,, 2011. 

Meinhold, G., Morton, A. C., and Avigad, D.: New insights into peri-Gondwana paleogeography and the Gondwana super-fan system from detrital zircon U-Pb ages, Gondwana Res., 23, 661–665,, 2013. 

Meireles, C.: Litoestratigrafia do Paleozoico do Sector a Nordeste de Bragança (Trás-os-Montes), Instituto Universitario de Geología “Isidro Parga Pondal”,- Área de Xeoloxía e Minería do Seminario de Estudos Galegos, Coruña, Spain, 2013. 

Meireles, C., Sá, A., Piçarra, J. M., González Clavijo, E., and Ribeiro, A.: Novos avanços no conhecimento do limite Ordovícico-Silúrico na região de Trás-os-Montes (NE Portugal), Estremoz, Portugal, 2006. 

Meireles, C. A., Santos, J., Pereira, E., and Ribeiro, A.: Carta Geológica de Portugal à escala 1:50.000, Folha 3-D (Espinhosela), Instituto Geológico e Mineiro, Madrid, Spain, 1999a. 

Meireles, C. A., Santos, J., Pereira, E., and Ribeiro, A.: Carta Geológica de Portugal à escala 1:50.000, Folha 4-C (Deilão), Instituto Geológico e Mineiro, Madrid, Spain, 1999b. 

Merino-Tomé, O., Gutiérrez-Alonso, G., Villa, E., Fernández-Suárez, J., Llaneza, J. M., and Hofmann, M.: LA-ICP-MS U-Pb dating of Carboniferous ash layers in the Cantabrian Zone (N Spain): stratigraphic implications, J. Geol. Soc., 174, 836–849,, 2017. 

Mulder, J. A., Karlstrom, K. E., Fletcher, K., Heizler, M. T., Timmons, J. M., Crossey, L. J., Gehrels, G. E., and Pecha, M.: The syn-orogenic sedimentary record of the Grenville Orogeny in southwest Laurentia, Precambrian Res., 294, 33–52,, 2017. 

Munhá, J., Ribeiro, A., and Ribeiro, M. L.: Blueschists in the Iberian Variscan Chain (Trás-os-Montes, NE Portugal), Comunicações dos Serviços Geológicos de Portugal, 70, 31–53, 1984. 

Murphy, J. B., Gutiérrez-Alonso, G., Fernández-Suárez, J., and Braid, J. A.: Probing crustal and mantle lithosphere origin through Ordovician volcanic rocks along the Iberian passive margin of Gondwana, Tectonophysics, 461, 166–180,, 2008. 

Naidoo, T., Zimmermann, U., Vervoort, J., and Tait, J.: Evidence of early Archean crust in northwest Gondwana, from U-Pb and Hf isotope analysis of detrital zircon, in Ediacaran surpacrustal rocks of northern Spain, Int. J. Earth Sci., 107, 409–429,, 2018. 

Noronha, F., Ribeiro, M. A., Martins, H., and Lima, J.: Carta Geológica de Portugal à escala 1:50.000, Folha 6-D (Vila Pouca de Aguiar), Laboratório Nacional de Energia e Geologia, Amadora, Portugal, 1998. 

Nuño Ortea, C., López García, M. J., Ferragne, A., and Rúiz García, C.: Mapa Geológico de España, E. 1:50.000 y Memoria, Hoja 303-Verín, Segunda Serie, IGME, Madrid, Spain, 29 pp., 1981. 

Ogata, K., Festa, A., Pini, G. A., and Alonso, J. L.: Submarine Landslide Deposits in Orogenic Belts, in: Submarine Landslides, edited by: Ogata, K., Festa, A., and Pini, G. A., Geophysical Monograph Series, American Geophysical Union, Washington DC, USA, 1–26,, 2019. 

Oliveira, J. T., González-Clavijo, E., Alonso, J., Armendáriz, M., Bahamonde, J. R., Braid, J. A., Colmenero, J. R., Dias da Silva, Í., Fernandes, P., Fernández, L. P., Gabaldón, V., Jorge, R. S., Machado, G., Marcos, A., Merino-Tomé, Ó., Moreira, N., Murphy, J. B., Pinto de Jesus, A., Quesada, C., Rodrigues, B., Rosales, I., Sanz-López, J., Suárez, A., Villa, E., Piçarra, J. M., and Pereira, Z.: Synorogenic Basins, in: The Geology of Iberia: A Geodynamic Approach: Volume 2: The Variscan Cycle, edited by: Quesada, C. and Oliveira, J. T., Springer International Publishing, Cham, UK,, 2019a. 

Oliveira, J. T., Quesada, C., Pereira, Z., Matos, J. X., Solá, A. R., Rosa, D., Albardeiro, L., Díez-Montes, A., Morais, I., Inverno, C., Rosa, C., and Relvas, J.: South Portuguese Terrane: A Continental Affinity Exotic Unit, in: The Geology of Iberia: A Geodynamic Approach: Volume 2: The Variscan Cycle, edited by: Quesada, C. and Oliveira, J. T., Springer International Publishing, Cham, UK,, 2019b. 

Pastor-Galán, D., Gutiérrez-Alonso, G., Murphy, J. B., Fernández-Suárez, J., Hofmann, M., and Linnemann, U.: Provenance analysis of the Paleozoic sequences of the northern Gondwana margin in NW Iberia: Passive margin to Variscan collision and orocline development, Gondwana Res., 23, 1089–1103,, 2013. 

Pereira, E.: Estudo geológico-estructural da região de Celorico de Basto e sua intrepretação geodinâmica, PhD thesis, Universidade de Lisboa, Lisboa, Portugal, 274 pp., 1987. 

Pereira, E.: Carta Geológica de Portugal à escala 1:50.000, Folha 10-A (Celorico de Basto), Serviços Geológicos de Portugal, Lisbon, Portugal, 1989. 

Pereira, E., Ribeiro, A., and Silva, N.: Carta Geológica de Portugal à Escala 1:50.000, Folha 7-D (Macedo de Cavaleiros), Laboratório Nacional de Energia e Geologia, Lisboa, Portugal, 1998. 

Pereira, E., Ribeiro, A., and Castro, P.: Notícia Explicativa da Carta Geológica de Portugal à Escala 1:50.000, Folha 7-D (Macedo de Cavaleiros) Laboratório Nacional de Energía e Geologia, Lisboa, Portugal, 2000. 

Pereira, E., Pereira, D. Í., Rodrigues, J. F., Ribeiro, A., Noronha, F., Ferreira, N., Sá, C. M. D., Farinha Ramos, J., Moreira, A., and Oliveira, A. F.: Notícia Explicativa da Folha 2 da Carta Geológica de Portugal à Escala 1:200.000, Instituto Nacional de Engenharia, Tecnologia e Inovação, Lisboa, Portugal, 2006. 

Pereira, E., Ribeiro, A., Rebelo, J. A., and Castro, P.: Notícia Explicativa da Carta Geológica de Portugal à Escala 1:50.000, Folha 11-B (Mogadouro) Laboratório Nacional de Energía e Geologia, Lisboa, Portugal, 2008. 

Pereira, E., Ferreira da Silva, A., Rebelo, J. A., Ribeiro, A., and Dias, R.: Carta Geológica de Portugal à escala 1:50.000, Folha 11-D (Carviçais), Laboratório Nacional de Energía e Geologia, Lisboa, Portugal, 2009. 

Pereira, M. F., Chichorro, M., Solá, A. R., Silva, J. B., Sánchez-García, T., and Bellido, F.: Tracing the Cadomian magmatism with detrital/inherited zircon ages by in-situ U-Pb SHRIMP geochronology (Ossa-Morena Zone, SW Iberian Massif), Lithos, 123, 204–217,, 2011. 

Pereira, M. F., Chichorro, M., Johnston, S. T., Gutiérrez-Alonso, G., Silva, J. B., Linnemann, U., Hofmann, M., and Drost, K.: The missing Rheic Ocean magmatic arcs: Provenance analysis of Late Paleozoic sedimentary clastic rocks of SW Iberia, Gondwana Res., 22, 882–891,, 2012a. 

Pereira, M. F., Solá, A. R., Chichorro, M., Lopes, L., Gerdes, A., and Silva, J. B.: North-Gondwana assembly, break-up and paleogeography: U/Pb isotope evidence from detrital and igneous zircons of Ediacaran and Cambrian rocks of SW Iberia, Gondwana Res., 22, 866–881,, 2012b. 

Pereira, M. F., Ribeiro, C., Vilallonga, F., Chichorro, M., Drost, K., Silva, J. B., Albardeiro, L., Hofmann, M., and Linnemann, U.: Variability over time in the sources of South Portuguese Zone turbidites: evidence of denudation of different crustal blocks during the assembly of Pangaea, Int. J. Earth Sci., 103, 1453–1470,, 2014. 

Pereira, M. F., Gutíerrez-Alonso, G., Murphy, J. B., Drost, K., Gama, C., and Silva, J. B.: Birth and demise of the Rheic Ocean magmatic arc(s): Combined U–Pb and Hf isotope analyses in detrital zircon from SW Iberia siliciclastic strata, Lithos, 278, 383–399,, 2017. 

Pereira, M. F., Gama, C., Dias da Silva, Í., Fuenlabrada, J. M., Silva, J. B., and Medina, J.: Isotope geochemistry evidence for Laurussian-type sources of South Portuguese Zone Carboniferous turbidites (Variscan Orogeny), Geol. Soc. Spec. Publ., 503, SP503-2019-2163,, 2020a. 

Pereira, M. F., Gama, C., Dias da Silva, Í., Silva, J. B., Hofmann, M., Linnemann, U., and Gärtner, A.: Chronostratigraphic framework and provenance of the Ossa-Morena Zone Carboniferous basins (southwest Iberia), Solid Earth, 11, 1291–1312,, 2020b. 

Pereira, Z., Meireles, C., and Pereira, E.: Upper Devonian palynomorphs of NE sector of Trás-os-Montes (Central Iberian Zone), XV Reunión de Geología del Oeste Peninsular, in: Proceedings of the International Meeting on Cadomian Basement, 26 September–3 October 1999, Badajoz, Spain, 201–206, 1999. 

Pérez-Cáceres, I., Martínez Poyatos, D., Simancas, J. F., and Azor, A.: Testing the Avalonian affinity of the South Portuguese Zone and the Neoproterozoic evolution of SW Iberia through detrital zircon populations, Gondwana Res., 42, 177–192,, 2017. 

Pérez-Estaún, A.: Aportaciones al conocimiento del Carbonífero de San Clodio (Prov. de Lugo), Brev. Geol. Ast., 1, 3–5, 1974. 

Piçarra, J. M. and Meireles, C.: Identificação de graptólitos do Ludlow (Silúrico superior), na área de Guadramil (Bragança, Zona Centro Ibérica, Portugal): implicações na estratigrafia regional, Ciências da Terra (UNL), no. esp. V, A126–A129, 2003. 

Piçarra, J., Gutiérrez-Marco, J., Sarmiento, G. N., and Sá, A.: Novos dados de Conodontes e Graptólitos no Paleozóico parautóctone da Zona Galiza-Trás-os-Montes (Espanha e Portugal), in: Proceedings of the VII Congresso Nacional de Geologia, 2006, Évora, Portugal, 653–656, 2006a. 

Piçarra, J. M., Gutiérrez Marco, J. C., Sá, A. A., Meireles, C., and González Clavijo, E.: Silurian graptolite biostratigraphy of the Galicia-Trás-os-Montes Zone (Spain and Portugal), J. Geol. Soc. Sweden, 128, 185–188, 2006b. 

Pin, C., Paquette, J. L., Ábalos, B., Santos, F. J., and Gil Ibarguchi, J. I.: Composite origin of an early Variscan transported suture: Ophiolitic units of the Morais Nappe Complex (north Portugal), Tectonics, 25, TC5001,, 2006. 

Puetz, S. J.: A relational database of global U-Pb ages, Geosci. Front., 9, 877–891,, 2018. 

Quiroga de la Vega, J. L.: Estudio geológico del paleozoico del W de Zamora (Alba y Aliste), PhD thesis, University of Oviedo, Spain, 210 pp., 1981. 

Ribeiro, A.: Contribution à l'étude téctonique de Trás-os-Montes Oriental, Serviços Geológicos de Portugal, Lisboa, Portugal, 1974. 

Ribeiro, A.: A Evolução Geodinâmica de Portugal; os ciclos ante-mesozóicos, in: Geologia de Portugal, edited by: Dias, R., Araújo, A., Terrinha, P., and Kullberg, J. C., Escolar Editora, Lisboa, Portugal, 15–57, 2013. 

Ribeiro, A. and Sanderson, D.: SW-Iberia-Transpressional Orogeny in the Variscides, in: EUROPROBE- Lithosphere dynamics, Origin and evolution of continents, edited by: Gee, D. G. and Zeyen, H. J., Europrobe Secretriat, Uppsala, Sweden, 1996. 

Ribeiro, A., Pereira, E., Dias, R., Gil Ibarguchi, J. I., and Arenas, R.: Allochthonous Sequences, in: Pre-Mesozoic Geology of Iberia, edited by: Dallmeyer, R. D. and Garcia, E. M., Springer, Berlin and Heidelberg, Germany,, 1990. 

Ribeiro, A., Munhá, J., Dias, R., Mateus, A., Pereira, E., Ribeiro, M. L., Fonseca, P., Araújo, A., Oliveira, J. T., Romão, J., Chaminé, H., Coke, C., and Pedro, J. C.: Geodynamic Evolution of the SW Europe Variscides, Tectonics, 26, TC6009,, 2007. 

Ribeiro, M. A., Noronha, F., and Cuney, M.: Lithogeochemical study of the metassedimentary units of Vila Pouca de Aguiar Área (Westhern Trás-os-Montes, Northern Portugal), Univ. Porto. Mem. Mus. Lab. Mineral. Geol. Fac. Ciênc., 3, 299–303, 1993. 

Ribeiro, M. L.: Estudo litogeoquímico das formações metassedimentares encaixantes de mineralizações em Trás-os-Montes Ocidental, Implicações metalogenéticas, PhD thesis, University of Porto, Porto, Portugal, 231 pp., 1998. 

Ribeiro, M. L. and Ribeiro, A.: Signification paléogéographique et tectonique de la présence de galets de roches métamorphiques dans le flysch d'âge dévonien supérieur du Tras-Os-Montes oriental (Nord-Est du Portugal), C. R. Acad. Sc. Paris, 278, 3161–3163, 1974. 

Ribeiro, M. L. and Ribeiro, A.: Analise petrográfica e textural dos Gneisses de Saldanha (Trás-os-Montes oriental): elementos para nova intrepretação estratigráfica, Comunicações do Instituto Geológico e Mineiro, 91, 5–16, 2004. 

Riemer, W.: Datos para el conocimiento de la estratigrafía de Galicia, Notas y Comunicaciones, IGME, 81, 7–20, 1966. 

Rodrigues, B., Chew, D. M., Jorge, R. C. G. S., Fernandes, P., Veiga-Pires, C., and Oliveira, J. T.: Detrital zircon geochronology of the Carboniferous Baixo Alentejo Flysch Group (South Portugal); constraints on the provenance and geodynamic evolution of the South Portuguese Zone, J. Geol. Soc., 172, 294–308,, 2015. 

Rodrigues, J., Pereira, E., Ribeiro, A., and Meireles, C.: Organização tectonoestratigráfica do Complexo Parautóctone do NE de Portugal: uma proposta, in: Proceedings of the VI Congresso Nacional de Geologia, Lisbon, Portugal, 4–6 June 2003, 76–78, 2003. 

Rodrigues, J. F.: Estructuras do Arco de Santa Comba-Serra da Garraia, Parautoctone de Trás-os-Montes, PhD thesis, University of Lisboa, Lisboa, Portugal, 308 pp., 2008. 

Rodrigues, J. F., Pereira, E., and Ribeiro, A.: Estructura Interna do Complexo de Mantos Parautoctones, Sector de Murça-Mirandela (NE de Portugal), in: Geologia de Portugal no Contexto da Ibéria, edited by: Dias, R., Araújo, A., Terrinha, P., and Kullberg, J. C., Universidade de Évora, Évora, Spain, 2006. 

Rodrigues, J. F., Ribeiro, A., Pereira, E., Ribeiro, M. L., Ferreira, N., and Meireles, C. A.: Carta Geológica de Portugal à escala 1:50.000, Folha 7-C (Mirandela), Laboratório Nacional de Energia e Geologia, Amadora, Portugal, 2010. 

Rodrigues, J. F., Ribeiro, A., and Pereira, E.: Complexo de Mantos Parautóctones do NE de Portugal: estructura interna e tectonoestratigrafia, in: Geologia de Portugal, edited by: Dias, R., Araújo, A., Terrinha, P., and Kullberg, J. C., Escolar Editora, Lisboa, Portugal, 2013. 

Rodríguez, J., Cosca, M. A., Gil Ibarguchi, J. I., and Dallmeyer, R. D.: Strain partitioning and preservation of 40Ar/39Ar ages during Variscan exhumation of a subducted crust (Malpica-Tui complex, NW Spain), Lithos, 70, 111–139,, 2003. 

Romariz, C.: Graptolitos do Silúrico Português, Rev. Fac. Ciências, 10, 155–312, 1962. 

Romariz, C.: Graptolitos silúricos do Noroeste Peninsular, Comun. Serv. Geol. Port., 53, 107–156, 1969. 

Rubio Pascual, F., Arenas, R., García, F. D. A., Martínez Catalán, J. R., Abati, J., Catalán, J. R. M., Hatcher, R. D., Arenas, R., and García, F. D.: Contrasting high-pressure metabasites from the Santiago unit (Ordenes Complex, northwestern Iberian Massif, Spain), in: Variscan-Appalachian dynamics: The building of the late Paleozoic basement, edited by: Martínez Catalán, J. R., Hatcher Jr., R. D., Arenas, R., and Díaz García, F., Geological Society of America, Boulder, Colorado, USA,, 2002. 

Sá, A. A., Meireles, C., Castro, P., and Vaz, N.: Ordovician “olistoliths” from Casal do Rato Formation (Trás-os-Montes): contribution to their reappraisal, Comun. Inst. Geol. Min., 101, 307–311, 2014. 

Sánchez García, T., Chichorro, M., Solá, A. R., Álvaro, J. J., Díez-Montes, A., Bellido, F., Ribeiro, M. L., Quesada, C., Lopes, J. C., Dias da Silva, Í., González-Clavijo, E., Gómez Barreiro, J., and López-Carmona, A.: The Cambrian-Early Ordovician Rift Stage in the Gondwanan Units of the Iberian Massif, in: The Geology of Iberia: A Geodynamic Approach: Volume 2: The Variscan Cycle, edited by: Quesada, C. and Oliveira, J. T., Springer International Publishing, Cham, UK,, 2019. 

Santos Zalduegui, J. F., Schärer, U., and Gil Ibarguchi, J. I.: Isotope constraints on the age and origin of magmatism and metamorphism in the Malpica-Tuy allochthon, Galicia, NW Spain, Chem. Geol., 121, 91–103,, 1995. 

Santos Zalduegui, J. F., Schärer, U., Gil Ibarguchi, J. I., and Girardeau, J.: Origin and evolution of the Paleozoic Cabo Ortegal ultramafic-mafic complex (NW Spain): U-Pb, Rb-Sr and Pb-Pb isotope data, Chem. Geol., 129, 281–304,, 1996. 

Sarmiento, G. N. and García-López, S.: El Índice de Alteración de Color (CAI) de los conodontos: limitaciones y posibilidades. Ejemplos de su aplpicación en el Hercínico Ibérico, Rev. Soc. Geol. España, 9, 113–123, 1996. 

Sarmiento, G. N., Calvo, A. A., and González Clavijo, E.: Conodontos paleozoicos (Ashgill-Emsiense) del Sinforme de Alcañices (oeste de Zamora, España), in: Libro de Resúmenes y Excursiones, XIII Jornadas de Paleontología and V Reunión Internacional Proyecto edited by: Grandal d'Anglade, A., Gutiérrez-Marco, J. C., and Santos Fidalgo, L., Sociedad Española de Paleontología, Coruña, Spain, 108–111, 1997. 

Sarmiento, G. N., Piçarra, J. M., Rebelo, J. A., Robardet, M., Gutiérrez Marco, J. C., Storch, P., and Rábano, I.: Le Silurien du Synclinorium de Moncorvo (NE du Portugal): Biostratigraphie et importance paléogéographique, GEOBIOS, 32, 749–767,, 1998. 

Schermerhorn, L. J. G. and Kotsch, S.: First occurrence of lawsonite in Portugal and tectonic implications, Comunicações do Instituto Geológico e Mineiro, 70, 23–29, 1984. 

Shaw, J., Gutiérrez-Alonso, G., Johnston, S. T., and Pastor Galán, D.: Provenance variability along the Early Ordovician north Gondwana margin: Paleogeographic and tectonic implications of U-Pb detrital zircon ages from the Armorican Quartzite of the Iberian Variscan belt, Geol. Soc. Am. Bull., 126, 702–719,, 2014. 

Silva, J. B., Oliveira, J. T., and Ribeiro, A.: Structural Outline, in: Pre-Mesozoic Geology of Iberia, edited by: Dallmeyer, R. D. and Garcia, E. M., Springer, Berlin and Heidelberg, Germany,, 1990. 

Smeraglia, L., Aldega, L., Billi, A., Carminati, E., Di Fiore, F., Gerdes, A., Albert, R., Rossetti, F., and Vignaroli, G.: Development of an Intrawedge Tectonic Mélange by Out-of-Sequence Thrusting, Buttressing, and Intraformational Rheological Contrast, Mt. Massico Ridge, Apennines, Italy, Tectonics, 38, 1223–1249,, 2019. 

Sousa, M. B.: Considerações sobre a estratigrafia do Complexo Xisto-Grauváquico (CXG) e sua relação com o Paleozóico Inferior, Cuadernos Geología Ibérica, 9, 9–36, 1984. 

Stampfli, G. M., Hochard, C., Vérard, C., Wilhem, C., and von Raumer, J.: The formation of Pangea, Tectonophysics, 593, 1–19,, 2013. 

Stephan, T., Kroner, U. W. E., and Romer, R. L.: The pre-orogenic detrital zircon record of the Peri-Gondwanan crust, Geol. Mag., 156, 281–307,, 2018. 

Talavera, C., Montero, P., Martínez Poyatos, D., and Williams, I. S.: Ediacaran to Lower Ordovician age for rocks ascribed to the Schist-Graywacke Complex (Iberian Massif, Spain): Evidence from detrital zircon SHRIMP U-Pb geochronology, Gondwana Res., 22, 928–942,, 2012. 

Talavera, C., Martínez Poyatos, D., and González Lodeiro, F.: SHRIMP U-Pb geochronological constraints on the timing of the intra-Alcudian (Cadomian) angular unconformity in the Central Iberian Zone (Iberian Massif, Spain), Int. J. Earth Sci., 104, 1739–1757,, 2015. 

Teixeira, C. and Pais, J.: Sobre a presença de Devónico na região de Bragança (Guadramil e Mofreita) e de Alcañices (Zamora), Bol. Soc. Geol. Portugal, 18, 199–202, 1973. 

Teixeira, R. J. S., Neiva, A. M. R., Silva, P. B., Gomes, M. E. P., Andersen, T., and Ramos, J. M. F.: Combined U-Pb geochronology and Lu-Hf isotope systematics by LAM-ICPMS of zircons from granites and metasedimentary rocks of Carrazeda de Ansiães and Sabugal areas, Portugal, to constrain granite sources, Lithos, 125, 321–334,, 2011. 

Vacas, J. M. and Martínez Catalán, J. R.: El sinforme de Alcañices en la transversal de Manzanal del Barco, Stvdia Geológica Salmantina, 24, 151–175, 1987. 

Valladares, M. I., Barba, P., Ugidos, J. M., Colmenero, J. R., and Armenteros, I.: Upper Neoproterozoic-Lower Cambrian sedimentary successions in the Central Iberian Zone (Spain): sequence stratigraphy, petrology and chemostratigraphy. Implications for other European zones, Int. J. Earth Sci., 89, 2–20,, 2000. 

Valverde-Vaquero, P., Marcos, A., Farias, P., and Gallastegui, G.: U-Pb dating of Ordovician felsic vulcanics in the Schistose Domain of the Galicia-Trás-os-Montes Zone near Cabo Ortegal (NW Spain), Geol. Acta, 3, 27–37,, 2005.  

Valverde-Vaquero, P., Farias, P., Marcos, A., and Gallastegui, G.: U-Pb dating of Siluro-Ordovician volcanism in the Verín Synform (Ourense; Schistose Domain, Galicia-Trás-os-Montes Zone), Geogaceta, 41, 247–250, 2007. 

Vermeesch, P.: IsoplotR: A free and open toolbox for geochronology, Geosci. Front., 9, 1479–1493,, 2018. 

Wilmsen, M., Fürsich, F. T., Seyed-Emami, K., Majidifard, M. R., and Taheri, J.: The Cimmerian Orogeny in northern Iran: tectono-stratigraphic evidence from the foreland, Terra Nova, 21, 211–218,, 2009. 

Zimmermann, U., Andersen, T., Madland, M. V., and Larsen, I. S.: The role of U-Pb ages of detrital zircons in sedimentology – An alarming case study for the impact of sampling for provenance interpretation, Sediment. Geol., 320, 38–50,, 2015. 

Short summary
In northwest Iberia, a low-grade unit separates the unrooted accretionary prism from the autochthon. It was divided into Upper (preorogenic) and Lower (synorogenic) parautochthons around the Morais Complex. A review of the sedimentary sequence characteristics, supported by U–Pb ages in 17 new samples of synorogenic sediments and volcanic olistoliths, proposes this model to show the Galicia–Trás-os-Montes Zone forming a tectonic carpet exposed in most of the Parautochthon–Autochthon limit.