Coeval interaction between magmas of contrasting composition (Late Carboniferous-Early Permian Santa Eulália-Monforte massif, Ossa-Morena Zone): field relationships and geochronological constraints
Coeval interaction between magmas of contrasting composition (Late Carboniferous-Early Permian Santa Eulália-Monforte massif, Ossa-Morena Zone): field relationships and geochronological constraints
Geologica Acta: an international earth science journal, vol. 15, no. 4, pp. 409-428, 2017
Universitat de Barcelona
Received: 01 June 2017
Accepted: 01 November 2017
Published: 01 November 2017
Funding
Funding source: Fundação para a Ciência e Tecnologia
Contract number: “GOLD” (PTDC/GEO-GEO/2446/2012: COMPETE: FCOMP01-0124-FEDER-029192) and “LITHOS” (CGL2013-48408C3-1-P)
Award recipient: Coeval interaction between magmas of contrasting composition (Late Carboniferous-Early Permian Santa Eulália-Monforte massif, Ossa-Morena Zone): field relationships and geochronological constraints
Funding
Funding source: Universidad de Huelva
Award recipient: C. Rodríguez
Abstract: The Santa Eulália-Monforte massif is a post-kinematic Late Carboniferous-Early Permian calc-alkaline composite massif (LC-EP) located in the Ossa-Morena Zone (OMZ, Portugal). This paper examines the field relationships between pinkish granites and mafic-intermediate rocks from the external ring of this massif and presents new U-Pb zircon age determinations. The estimated 206Pb/238U ages, 297±4Ma for the pinkish granite and 303±3Ma for a gabbro-diorite point to a short time interval between the crystallization of both magmas. At outcrop scale, contacts of the mafic-intermediate rocks with the host pinkish granite are sharp and corroborate this age relationship, but do not justify why the host-enclave contacts often has curved and irregular shapes, indicating liquid-liquid interaction. A full analysis of the distribution of U-Pb zircon ages and respective Th/U ratios suggests that the compositionally distinct magmas were roughly contemporaneous. The obtained ages also fit the LC-EP Iberian calc-alkaline suite that was formed contemporaneously to the development of the Iberian-Armorican Arc and when the Paleotethyan oceanic Plate subducted under the Eurasian active margin. Taking this geodynamic setting as a reference, the LC-EP Iberian calc-alkaline magmatism can be interpreted as most probably related to the Cimmerian cycle instead of the traditionally accepted model that ascribe a closer connection of this magmatism with the Variscan cycle.
Keywords: iberian Massif, Late Carboniferous-Early Permian, Calc-alkaline bimodal plutonism, Zircon dating, Variscan and Cimmerian cycles.
INTRODUCTION
The Iberian Massif experienced a polymetamorphic evolution during the Carboniferous as a result of the collision between Laurussia and Gondwana (e.g. Variscan Orogeny; Martínez Catalán et al., 2009; Díez Fernández et al., 2016). The Laurussia-Gondwana convergence caused crustal thickening by the stacking of tectonic nappes and the rejuvenation of reliefs. The regional switch from contractional to extensional deformation was followed by an orogenic collapse and the generation of abundant crustal-derived magmatism (Escuder Viruete et al., 1994; Pereira et al., 2009; 2017a; Martínez Catalán et al., 2014; Alcock et al., 2015; Cambeses et al., 2015; Díez Fernández and Pereira, 2016). Mantle contributions were more abundant towards the SW Iberian Massif, as syn-kinematic magmatism related to intra-orogenic extension was slightly older in the Ossa-Morena and South Portuguese zones (Simancas et al., 2005; Pereira et al., 2009; 2015a) than in the Central Iberian Zone.
The age and compositional ranges of syn-kinematic (Carboniferous) and post-kinematic (Late CarboniferousEarly Permian, LC-EP) intrusions are wide in the NW (e.g. Central Iberian Zone; Castro et al., 2002 and references therein) and SW (Ossa-Morena Zone; Pereira et al., 2015a and references therein) Iberian Massif. Most syn- and postkinematic granitic rocks are S-type peraluminous, but there are also granitic rocks with transitional features between S- and I-type, I-type granites and hybrid granites, including metaluminous varieties (Dias et al., 1998, 2002; Castro et al., 1999, 2002; Bea et al., 1999, 2003; Fernández-Suárez et al., 2000, 2011; Neiva and Gomes 2001; González Menéndez et al., 2006; Moita et al., 2009; Neiva et al., 2009; Solá et al., 2009; Lima et al., 2011, 2012; Teixeira et al., 2011; Villaseca et al., 2012; Merino Martínez et al., 2014; Pereira et al., 2015a, 2017a; López-Moro et al., 2017). These Carboniferous and Permian calc-alkaline granitic rocks are often closely and spatially associated with calc-alkaline mafic-intermediate rocks (e.g.Castro et al., 2003; Bea, 2004; Larrea et al., 2004; González Menéndez et al., 2006; Jesus et al., 2007; Pin et al., 2008; Solá et al., 2009; Moita et al., 2015) and are characterized by many structures and textures formed through the complex interplay of distinct petrogenetic processes. Such composite plutons typically preserve an integrated record of prolonged magmatic evolution, including pluton accumulation, crystallization and emplacement within the crust that evidence the coeval interaction between two or more magmas of contrasting composition and physical properties (Turnbull et al., 2010). Mixing of magmas within a magma chamber to produce a homogeneous hybrid, followed by mingling and partial mixing of the hybrid magma with one of the remaining endmembers is commonly identified in the successive stages of the magmatic history of composite plutons (Barbarin, 2005 and references therein). By studying the temporal and spatial interaction of magmas of contrasting compositions that remained in the contact in liquid state, insights can be gained into the duration of physical and thermal processes responsible for the formation of plutons.
Our study presents original field data and SHRIMP zircon U-Pb dating for a post-kinematic intrusion in the Ossa-Morena Zone with the aim of studying the temporal and spatial interaction of a felsic and a mafic-intermediate magma forming the external ring of the Santa EuláliaMonforte massif. The obtained results will shed light on the timing of crystallization and emplacement of the Santa Eulália-Monforte intrusion, and on the global geodynamic context of the post-kinematic Iberian calc-alkaline magmatism. The results achieved broaden the knowledge on the spatial and temporal relationships between distinct LC-EP post-kinematic intrusions of the Iberian Massif. Our data also show that the crystallization ages of the Santa Eulália-Monforte I-type granite and associated mafic rocks are younger than the Nisa-Albuquerque batholith. Both intrusions cut across the Variscan structures, establishing a maximum range of about 13Ma, from the Late Moscovian to Asselian, for the LC-EP post-kinematic magmatism in the SW Iberian Massif. The subject of this study was chosen because L.G. Corretgé, who receives an honorable mention in this volume, was one of the petrologists who mostly contributed to gather data on the magmatism of the Iberian Massif, thus, fostering scientific cooperation between Portuguese and Spanish geologists.
LP-EC POST-KINEMATIC INTRUSIONS OF THE IBERIAN MASSIF
The LC-EP granitic rocks of the Iberian Massif have been traditionally classified as Late-Variscan because they crosscut the main Variscan structures and syn-kinematic intrusive rocks (ca. 306-287Ma; Ferreira et al., 1987; Dias et al., 1998; Fernández-Suárez et al., 2011; Neiva et al., 2012; Orejana et al., 2012; Villaseca et al., 2012). However, these post-kinematic intrusions are also recently considered as unrelated to the Gondwana-Laurussia collision sensu stricto .e.g.Gutiérrez-Alonso et al., 2011; Pereira et al., 2015b).
The LC-EP post-kinematic intrusions are abundant in the Central Iberian Zone, but are relatively less represented in the Ossa-Morena Zone. The Nisa-Albuquerque massif is a composite elongated body oriented parallel to, and crossing, the boundary between the Central Iberian and Ossa-Morena zones (González Menéndez, 2002) (Fig. 1). This post-kinematic intrusion consists of a dominant S-type monzogranite-syenogranite (ca. 309-307Ma) that surrounds a discontinuous core of I-type tonalitegranodiorite (ca. 306Ma) (SHRIMP zircon U-Pb dating; Solá et al., 2009). The Santa Eulália-Monforte massif, located in the northern domains of the Ossa-Morena Zone, is a composite intrusion (Gonçalves and Coelho, 1969-1970; Gonçalves, 1971; Gonçalves et al., 1972, 1975; Carrilho Lopes et al., 1998; Sant’Ovaia et al., 2015) consisting of S- and I-type granites and mafic-intermediate bodies (González Menéndez et al., 2006). Rb-Sr wholerock ages of ca. 281Ma (Mendes, 1967-1968) and of ca. 290Ma (Pinto et al., 1987) have been, for decades, the only available isotopic ages performed on granites from the Santa Eulália-Monforte massif.
THE SANTA EULÁLIA-MONFORTE MASSIF
Host rocks of the Santa Eulália-Monforte massif
The Santa Eulália-Monforte massif intrudes two important domains of the Ossa-Morena Zone (Figs. 1 and 2): i) the strongly deformed and variable metamorphosed Ediacaran to Ordovician sedimentary and igneous rocks of the Coimbra-Córdoba shear zone (Burg et al., 1981; Pereira et al., 2008, 2010); and ii) the Cambrian-Ordovician sedimentary and igneous rocks of the Alter do Chão-Elvas sector (Apalategui et al., 1990; Oliveira et al., 1991).


The original contacts and textures of the EdiacaranCambrian sedimentary and volcanic rocks (Pereira et al., 2006; Linnemann et al., 2008) and Cambrian-Ordovician plutonic rocks (Pereira et al., 2011; Sánchez-García et al., 2013; Díez Fernández et al., 2015) are best exposed in the northern and southern margins (Fig. 2). In the northern and southern margins of the Coimbra-Córdoba shear zone, as well as in the Alter do Chão-Elvas sector, despite the Variscan ductile deformation and greenschist facies metamorphism, the Ediacaran succession-consisting of metapelite, metagreywacke, black metachert, metamafic rocks, and rare carbonate rocks is unconformably overlain by a Cambrian sequence (Gonçalves, 1971; Gonçalves et al., 1975; Oliveira et al., 1991; Pereira et al., 2006). From the base to the top, the Cambrian sequence consists of: i) felsic metavolcanic rocks, quartzite and marble, and ii) schists, metagreywackes, quartzites, metaconglomerates and metamafics rocks. Cambrian magmatism includes calc-alkaline and MORB compositions (Pereira and Quesada, 2006; Sánchez-García et al., 2010). The Ediacaran-Cambrian succession is intruded by Cambrian-Ordovician peralkaline and alkaline igneous rocks (Gonçalves and Assunção, 1970; Gonçalves et al., 1975: Díez Fernández et al., 2015).
In the Central Unit of the Coimbra-Córdoba shear zone, the Carboniferous (ca. 340-333Ma) ductile deformation and associated metamorphism, reaching amphibolite to granulite and ecologite facies conditions, transformed the original contact and textures of the Ediacaran-Ordovician rocks (Pereira et al., 2008, 2010, 2012) (Fig. 2). CambrianOrdovician plutons occur as elongated bodies of weakly- to moderately-foliated rocks, or as narrow bands of orthogneiss alternating with schists and paragneisses affected by strong Variscan deformation and metamorphism (Gonçalves et al., 1972; Pereira et al., 2008, 2010; Díez Fernández et al., 2015) (Fig. 1).
Field relationships, petrography and geochemistry
The Santa Eulália-Monforte massif (Fig. 2) is a ringshaped composite intrusion 18km wide and 32km long elongated body in E-W direction (Gonçalves and Coelho, 1969-1970; Gonçalves, 1971; Gonçalves et al., 1972, 1975; Oliveira, 1975; Carrilho Lopes et al., 1998; González Menéndez et al., 2002, 2006). It includes a 10km wide and 16km long core of greyish fine- to medium-grained biotite granite (Fig. 3A) surrounded by a ring consisting of a complex association between pinkish coarse-grained biotite granite and mafic-intermediate rocks (Figs. 3, 4, 5 and 6).




The core of this ring-shaped intrusion consists of S-(prevailing) and I-type high-K calc-alkaline peraluminous granite ranging in composition from monzogranite to granite (Fig. 7A; B). Quartz, plagioclase, microcline, biotite and muscovite are the main minerals, being also present cordierite (Fig. 3A). The geochemistry of the greyish granite suggests the derivation from a mixed source of metasedimentary (dominant) with intermediate to acid igneous rocks extracted from the mid-lower crust (González Menéndez et al., 2006). The external pinkish granite is composed of quartz, biotite, and K-feldspar prevails in respect to plagioclase. It also contains chlorite, amphibole, allanite and zircon. The available geochemical data (Carrilho Lopes et al., 1998; González Menéndez et al., 2006) indicate that it is an I-type high-K calc-alkaline peraluminous granite (Fig. 7A; C). The pinkish granite is interpreted to derive from partial melting of intermediate composition gneissic rocks from the mid-lower crust (González Menéndez et al., 2006).

The pinkish granite of the external ring (Fig. 6A) includes variable size and irregular, 0.5-1.5km wide and 0.8-6km long gabbro-diorite to granodioritequartzdiorite-tonalite enclaves, and kilometer-scale xenoliths consisting of Ediacaran and Cambrian host rocks (Fig. 6B). Mafic-intermediate rocks range compositionally (Fig. 7) from gabbro-diorite, which are mainly composed by plagioclase, amphibole, pyroxene, biotite, apatite and zircon, to granodiorite-quartzdioritetonalite consisting of plagioclase, K-feldspar, quartz, biotite, amphibole, apatite and zircon. The maficintermediate associations have been interpreted to derive from distinct sources (González Menéndez et al., 2006): an enriched mantle (gabbro-diorite) and from partial melting of a lower mafic crust (granodioritequartzdiorite-tonalite).
Field relationships show that the greyish granite from the core (Fig. 3A) intrudes both the pinkish granite (Fig. 6A) and the associated mafic-intermediate rocks (Fig. 4A). Maficintermediate bodies at the eastern margin of the Santa Eulália-Monforte massif are intruded by the pinkish granite (Figs. 3B and 6C, D, E). Angular mafic enclaves (Fig. 5B) seem to represent early cooled pieces of gabbro caught up in the granodiorite-quartzdiorite-tonalite magma, but the variation from sharp though irregular contacts to more diffuse contacts suggests that they are contemporaneous (Fig. 4E). It is notable the presence of composite enclaves (Fig. 3E) and irregular patches of granodioritequartzdiorite-tonalite (Fig. 3D) hosted by gabbro-diorite, together with a channel-like enclave swarm of maficintermediate rocks (Fig. 4C; D) suggest intermingling of two magmas. Gabbro-dioritic rocks include leucocratic veins and lenses of quartzdiorite and/or tonalite (Fig. 3C; D). The occurrence of K-feldspar within the gabbro-dioritic enclaves, which may be the same as that of the host granitic rock (Fig. 5), suggests that granitic magma could have interacted with the mafic-intermediate magmas, probably at a time when both were not completely crystallized. The existence of mafic-intermediate enclaves with irregular, rounded and elliptical shapes, suggests mingling with the host granitic magma (Fig. 5).
U-PB GEOCHRONOLOGY
Sample preparation and analytical methods
Zircon separation from a sample of granite (MFT-1) and a sample of gabbro-diorite (BCA-2) of the external ring of the Santa Eulália-Monforte massif (sampling sites on Fig. 2) was performed by means of traditional techniques using dense liquids and magnetic separation at the University of Évora (Portugal). Selected crystals were hand-picked using a binocular lens. Zircon concentrates were cast on a 3.5cm diameter epoxy mount, together with zircon standards (TEMORA zircon, SL13 zircon and GAL zircon), then polished and documented using SEM-CL (Fig. 8), at the IBERSIMS (University of Granada, Spain). Mounts were coated with gold (80-microns thick) and inserted into the SHRIMP for analysis.

Each selected spot was rastered with the primary beam during 120s prior to analysis, and then analyzed over 6 scans following the isotope peak sequence 196Zr2O, 204Pb, 204.1 background, 206Pb, 207Pb, 208Pb, 238U, 248ThO, 254UO. Every peak of each scan was measured sequentially 10 times with the following total counts per scan: 2s for mass 196; 5s for masses 238, 248, and 254; 15s for masses 204, 206, and 208; and 20s for mass 207. The primary beam, composed of 16.16.2+, was set to an intensity of 4 to 5nA, using a 120-micron Kohler aperture, which generates 17 x 20micron elliptical spots on the target. The secondary beam exit slit was set at 80microns, achieving a resolution of about 5,000 at 1% peak height. All calibration procedures were performed on the standards included on the same mount. Mass calibration was performed on the GAL zircon (ca. 480Ma, very high U, Th and common lead content; Montero et al., 2008). Analytical sessions initially involved the measurement of SL13 zircon (Claoué-Long et al., 1995), which was used as a concentration standard (238ppm U). TEMORA zircon (ca. 417Ma, Black et al., 2003), used as isotope ratios standard, was then measured every four unknowns.
Data reduction was performed using the SHRIMPTOOLS software specifically developed for IBERSIMS (available at www.ugr.es/~fbea). The intensity of each measured isotope was calculated in two steps using the software: first, the STATA lettervalue display algorithm was used to find outliers in the ten replicates measured at each peak during each scan, discarding them and averaging the rest once they had been normalized to SBM measurements; then, for each isotope, a robust regression of each scan average versus time, if measured, was performed. The final result for each isotope was calculated as the value at the midtime of the analysis resulting from a robust regression line. Errors (95% confidence level) were calculated as the standard error of the linear prediction at the midpoint of the analysis. 206Pb/238U was calculated from the measured 206Pb+/238U+ and UO+/U+ following the method described by Williams (1998). For high-U zircons (U>2,500ppm) 206Pb/238U was further corrected using the algorithm of Williams and Hergt (2000). Plotted and tabulated analytical uncertainties are 1σ precision estimates. Uncertainties in calculated mean ages are 95% confidence limits (tσ, where t is the Student’s t multiplier) and, for mean 206Pb/238U ages, include the uncertainty in Pb/U calibration (0.3-0.5%). Common Pb corrections assumed a model common Pb composition appropriate to the age of each spot (Cumming and Richards, 1975). U-Th-Pb data are presented in Table 1 (electronic appendix, available at www.geologica-acta.com) and plotted in Figure 9, using Isoplot 4 (Ludwig, 2009).
RESULTS
Sample MFt-1 (granite)
Most zircons are stubby and elongated euhedral prisms (150-300μm diameter in the long-axis) (Fig. 8). Morphologically complex zircons analyzed show cores of variable size (dark-cathodoluminescence (CL) subeuhedral cores with concentric zoning and wavy extinction, subeuhedral cores of a patchy to irregular texture with linear or wavy dark and bright-CL bands, sub-rounded dark-CL and unzoned cores). Cores are surrounded by concentric zoned to unzoned rims (Fig. 8). Few grains are simple, showing banded or concentric zoning and some show discordant zircon overgrowths or recrystallization textures. Twentynine U-Th-Pb isotopic analyses of 25 representative zircon grains are listed in Table 1. Twenty-five 206Pb/238U ages in the range ca. 536 to 273Ma, were obtained using 204-lead correction and show a discordance ≤5%. The oldest age of ca. 537Ma (Cambrian), obtained in analysis 8.2, is interpreted as an inherited old core and, indeed, close inspection showed that it is surrounded by a young overgrowth (8.1). The remaining 19 206Pb/238U ages were obtained in zircons with a wide range of U (295-2,277ppm) and Th (118-1,763ppm) contents, and an average Th/U ratio of 0.37 (error standard deviation=0.03; Fig. 9A). These 19 U-Th-Pb isotopic compositions yielded a weighted mean 206Pb/238U age of 296±4Ma (MSWD=0.62). However, the omission of three analyses with high-U (10.1, 14.1 and 16.1) gave a weighted mean 206Pb/238U age of 297±4Ma (MSWD=0.43; Fig. 9B), which is probably the best estimate of the crystallization age of the granite.

Sample BcA-2 (gabbro-diorite)
Zircons of this sample present morphological differences with respect to those of the granite described above. CL images of stubby and elongated euhedral prisms (150-350μm diameter in long-axis) show simpler internal textures. Most zircon grains have banded zoning (Fig. 8), but a few show oscillatory concentric zoning or are unzoned.
Twenty-four U-Th-Pb isotopic analyses of representative zircon grains are listed in Table 1. Common lead 204-lead corrected isotope ratios gave 206Pb/238U ages of ca. 319248Ma, having a cluster between ca. 310Ma and 290Ma. Sixteen U-Th-Pb analyses with discordance of ≤5% were performed on zircons characterized by a wide range of U (89-629ppm) and Th (54-343ppm) content and average Th/U ratio of 0.64 (error standard deviation=0.02; Fig. 9A). This group of analyses yielded a concordia age of 303±3Ma (MSWD=0.079; Fig. 9D) overlapping the weighted mean 206Pb/238U age of 303±4Ma (MSWD=1.01; Fig. 9E), which is probably the best estimate for the crystallization age of the gabbro-diorite.
DISCUSSION
Significance of the obtained U-Th-Pb data
The new obtained U-Pb ages of the Santa EuláliaMonforte massif suggest that the gabbro-diorite (ca. 303Ma) crystallized somewhat 6Ma earlier than the external granite (ca. 297Ma) (Figs. 9 and 10). However, the overlapping ages of zircons shown in Figure 10 suggest that the maficintermediate rocks and the pinkish granite can also be roughly contemporaneous, as documented by field data. In some cases, it has been demonstrated that zircon populations in a magma chamber are a mix of autocrysts, antecrysts, and xenocrysts, and that the resultant plutonic rocks are, thus, a mechanical mixture of crystals (e.g.Charlier et al., 2005; Miller et al., 2007) that may have different individual ages. The obtained data suggest that the external ring of the Santa Eulália-Monforte massif formed by the crystallization of injections of compositionally distinct magmas. These magmas were close enough in time so that the earlier gabbro-dioritic magma had not cooled yet and had differentiated significantly before the addition of the later granitic magma.

The Th/U ratios obtained from analyzed zircons of both samples (Fig. 10) prove the coexistence of two compositionally distinct magmas from which zircon crystallized during the evolutionary history of the bimodal magmatic chamber. The average Th/U ratio of the granite (0.40) is significantly lower than that of the gabbro-diorite (0.64), indicating felsic-intermediate metaluminous sources for the two magmas (Heaman et al., 1990; Hanchar and Miller, 1993, Hoskin and Schaltegger, 2003). In Figure 6, most zircon analyses of granite (sample MFT-1) are plotted in the range of 0.3<Th/U ratio<0.5, closer to the field of felsic peraluminous sources, in which probably there is abundant monazite (Th>U), and lower than those from the gabbro-diorite (sample BCA-2), in the range of 0.5<Th/U ratio<0.7, in which there is no other U-Th accessory.
Age relationship of the Santa Eulália-Monforte massif with other LC-EP Iberian calc-alkaline plutonic rocks
The reported weighted mean 206Pb/238U ages (this study) for the plutonic rocks of the external ring of the Santa Eulália-Monforte massif is of 303±3Ma (Late Carboniferous) for the gabbro-diorite and of 297±4Ma (Early Permian) for the granite. T-test indicates that the difference between these two zircon age populations is statistically significant (p-value is 0.02, i.e. lower than 0.05). Field observations of the consistent intrusive relationships between the greyish granite and the external ring (Gonçalves and Coelho, 1969-1970; Gonçalves, 1971; Oliveira, 1975; Gonçalves et al., 1972, 1975; Carrillho Lopes et al., 1998; González Ménendez et al., 2006; Sant’Ovaia et al., 2015) also suggest that these intrusions are coeval. Crystallization ages of the calc-alkaline magmatic rocks of the Santa Eulália-Monforte massif are within the wider range of ages obtained for other post-kinematic magmatic rocks of the Iberian Massif: i) the neighboring calcalkaline Nisa-Alpalhão and Los Pedroches batholiths (U-Pb zircon ages of ca. 314-304Ma; Carracedo et al., 2008; Solá et al., 2009) (Fig. 1); ii) the calc-alkaline plutonic rocks of the Central Iberian and West-Asturian Leonese zones (U-Pb zircon ages of ca. 312-286Ma; Valverde-Vaquero, 1997; Dias et al., 1998; Montero et al., 2004; Valle Aguado et al., 2005; Bea et al., 2006; Neiva et al., 2009; Fernández-Suárez et al., 2011; Gutiérrez-Alonso et al., 2011; Valverde-Vaquero et al., 2011; Díaz Alvarado et al., 2013; Díez Fernández and Pereira, 2017); and iii) the volcanic and plutonic rocks of the Central-eastern Pyrenees (U-Pb zircon ages of ca. 314-283Ma; Roberts et al. 2000; Aguilar et al., 2014; Druguet et al., 2014; Pereira et al., 2014; and 40Ar/39Ar ages in the range of ca. 291-285Ma; Solé et al., 2002). Our new U-Pb ages are ca. 20Ma younger with regard to the published age estimates for the Variscan plutonism and high-medium grade metamorphism recorded in the Ossa-Morena Zone (ca. 353-324Ma; Fig. 1). This age relationship indicates that the undeformed calc-alkaline magmatic rocks of the Santa Eulália-Monforte massif are consistent with a post-Variscan intrusion (i.e. postkinematic), representing so far the youngest LC-EP calc-alkaline magmatic event of the Ossa-Morena Zone. The emplacement of the Santa Eulália-Monforte massif immediately followed the neighboring calc-alkaline granite intrusion of the Nisa-Albuquerque massif (ca. 309-306Ma; Solá et al., 2009) and, in part, was almost synchronous to the latest magmatic pulses of the Los Pedroches massif (ca. 309-304Ma; Carracedo et al., 2008) (Fig. 1).
Geodynamic model
There is a general acceptance of a single, consensual interpretation that supports the growth of the Iberian LC-EP batholiths within the framework of the Variscan orogenic cycle (“late-Variscan batholiths”). This interpretation is considered valid even if the origin of these calc-alkaline magmas is about 60-80Ma younger than the development of the Rheic Ocean active margin (Martínez Catalán et al., 2009; Díez Fernández et al., 2016; Pereira et al., 2017b). This paradigm disregards the wider context of the amalgamation of the Pangaea and the spatial proximity of Iberia relative to the Eurasian active margin during the LC-EP closing of the Paleotethys Ocean (Pereira et al., 2015b).
The generation of LC-EP post-kinematic I-type calcalkaline Iberian batholiths in a tectonic setting other than subduction conflicts with essential thermal requirements and phase equilibria constraints on metaluminous magma generation and fractionation (Pereira et al., 2015b). The emplacement and thermal effects of the LC-EP calcalkaline batholiths into a fertile middle crustal level, dominated by several kilometer-thick fertile EdiacaranCambrian metasedimentary successions, explains the late generation of S-type peraluminous granitic rocks and local hybridization of the intrusive magma, which became more peraluminous and potassic (Díaz Alvarado et al., 2011; 2013). The lower crust is usually mentioned as the source of the Iberian calc-alkaline batholiths based on the assumptions that they are intrusive into the middle crust and that melts produced from the mantle underlying the continental crust must be basaltic (Dias et al., 1998; Neiva et al., 2009, 2012; Villaseca et al., 2009, 2012). The new conception of subduction-related magmatism ownig to plume-assisted relamination (Gerya et al., 2004; Castro et al., 2010; Gerya and Meilick, 2011; Vogt et al., 2012) assumes that the generation of magmas of intermediate composition (andesite-diorite) occurs by means of the melting of subducted materials in silicic composite plumes, i.e. composed of oceanic crust and sediments. These materials are later relaminated to a level below the lower crust, where they split by melt segregation into liquids that are emplaced in the middle and upper crust (batholiths) and residues that remain in the lower crust as mafic granulites (Castro, 2013, 2014). In the case of the LC-EP calc-alkaline Iberian magmatism, several factors indicate an extra-crustal origin (Castro et al., 2010, 2013; Castro, 2013, 2014; Pereira et al., 2015b). The LC-EP Iberian batholiths have close compositional similarities with the calc-alkaline silicic magmatism of active continental margins as the Andes (Patagonian batholith) and the North American Cordillera (Sierra Nevada batholith). They also include arc-related appinitic rocks, suggesting a link to a subduction setting (Murphy, 2013) similar to that of the post-Paleozoic Cordilleran batholiths formed at the Pacific active margin of the Americas (Castro et al., 2011). The period during which LC-EP calc-alkaline magmatism was generated (ca. 315-280Ma) coincides with the development of the Iberian-Armorican Arc (Fig. 11) (i.e. Cantabrian orocline; Gutiérrez-Alonso et al., 2008, 2011; Fernández-Suárez et al., 2011; Weil et al., 2012; Johnston et al., 2013), which was probably related to the Pangaea self-subduction (Gutiérrez-Alonso et al., 2008, 2011). In a global context, during the Late Carboniferous the western zones of the Iberian Massif were affected by the late stages of the Laurussia-Gondwana convergence, whereas the subduction of the Paleotethyan Ocean prevailed during and after LC-EP times (Cimmerian cycle) along the Eurasian eastern domains (Stampfli and Borel, 2002; Stampfli et al., 2013) (Fig. 11). The progression of the Paleotethys subduction led to the collision between Cimmerian terranes (detached from the non-collisional northern margin of Gondwana) and the Eurasian active margin (Cocks and Torsvik, 2006; Stampfli and Kozur, 2006), which can explain the inception of a LC-EP Paleotethyan magmatic arc in Iberia (Fig. 11; Pereira et al., 2015b).

CONCLUSIONS
This paper presents the first SHRIMP U-Pb zircon dating of the Santa Eulália-Monforte massif. Concordia and weighted average zircon ages from granite and gabbrodiorite samples indicate magmatic crystallization at ca. 297Ma and 303Ma, respectively. However, field evidence and full analysis of zircon ages do not definitively rule out that they may have been contemporaneous during a certain stage of magma emplacement.
The Santa Eulália-Monforte massif represents so far the youngest LC-EP calc-alkaline magmatic event of the Ossa-Morena Zone (ca. 303-297Ma), which immediately followed the neighboring calc-alkaline granite production of the Nisa-Albuquerque massif (ca. 309-306Ma) and in part, was almost synchronous to the latest magmatic pulses of the Los Pedroches massif (ca. 309-304Ma). The LC-EP Iberian magmatism was coeval with the formation of the Iberian-Armorican arc (Cantabrian orocline). The LC-EP Iberian magmatism can be framed into a plausible model of a transient continental magmatic arc developed in the Eurasian convergent Plate margin above the subducted Paleotethys oceanic Plate.
Acknowledgments
Financial support has been provided by Fundação para a Ciência e Tecnologia (Portugal) through research projects “GOLD” (PTDC/GEO-GEO/2446/2012: COMPETE: FCOMP01-0124-FEDER-029192) and “LITHOS” (CGL2013-48408C3-1-P). C. Rodríguez appreciates the financial support from Universidad de Huelva (Spain) through its postdoctoral program (under the “Estrategia Politíca Científica de la UHU 2016/2017 in collaboration with the Universidade de Évora (Portugal)). This is the IBERSIMS publication Nº 44. The authors are grateful for the critical spirit of the comments made by the editors of this volume, by J.B. Silva, R. Díez Fernández and an anonymous reviewer, who helped to greatly improve this work.
REFERENCES
Aguilar, C., Liesa, M., Castiñeiras, P., Navidad, M., 2014. Late Variscan metamorphic and magmatic evolution in the eastern Pyrenees revealed by U-Pb age zircon dating. Journal of the Geological Society, London, 171(2), 181-192.
Alcock, J.E., Martínez Catalán, J.R., Rubio Pascual, F.J., Montes, A.D., Díez Fernández, R., Gómez Barreiro, J., Arenas, R., Dias da Silva, Í., González Clavijo, E., 2015. 2-D thermal modeling of HT-LP metamorphism in NW and Central Iberia: Implications for Variscan magmatism, rheology of the lithosphere and orogenic evolution. Tectonophysics, 657, 21-37.
Apalategui, O., Eguiluz, L., Quesada, C., 1990. Ossa-Morena zone: Structure. In: Dallmeyer, R.D., Martínez García, E., (eds.). Pre-Mesozoic Geology of Iberia. Berlin, Germany, Springer-Verlag, 280-292.
Barbarin, B., 2005. Mafic magmatic enclaves and mafic rocks associated with some granitoids of the central Sierra Nevada batholith, California: nature, origin, and relations with the hosts. Lithos, 80, 155-177.
Bea, F., 2004. La naturaleza del magmatismo de la Zona Centro Ibérica: consideraciones generales y ensayo de correlación. In: Vera, J.A. (ed.). Geología de España. Ed. Sociedad Geológica de España (SGE) - Instituto Geológico y Minero de España (IGME), Madrid, 128-133.
Bea, F., Montero, P., Molina, J.F., 1999. Mafic precursors, Peraluminous Granitoids, and Late Lamprophyres in the Avila Batholith: A Model for the Generation of Variscan Batholiths in Iberia. Journal of Geology, 107, 399-419.
Bea, F., Montero, P., Zinger, T., 2003. The nature, origin, and thermal influence of the granite source layer of Central Iberia. Journal of Geology, 111, 579-595.
Bea, F., Montero, P., González Lodeiro, F., Talavera, C., Molina, J.F., Scarrow, J.H., Whitehouse, M.J., Zinger, T., 2006. Zircon thermometry and U-Pb ion-microprobe dating of the gabbros and associated migmatites of the Variscan Toledo Anatectic Complex, Central Iberia. Journal of the Geological Society, 163, 847-855.
Black, L.P., Kamo, S.L., Williams, I.S., Mundil, R., Davis, D.W., Korsch, R.J., Foudoulis, C., 2003. The application of SHRIMP to Phanerozoic geochronology, a critical appraisal of four zircon standards. Chemical Geology, 200, 171-188. DOI:10.1016/S0009-2541 (03) 00166-9.
Burg, J., Iglesias, M., Laurent, P., Ribeiro, A., 1981. Variscan intracontinental deformation: The Coimbra-Córdoba shear zone (SW Iberian Peninsula). Tectonophysics, 78, 161-177.
Cambeses, A., Scarrow, J.H., Montero, P., Molina, J.F., Moreno, J.A., 2015. SHRIMP U-Pb zircon dating of the Valencia del Ventoso plutonic complex, Ossa-Morena Zone, SW Iberia: Early Carboniferous intra-orogenic extension related ‘calc-alkaline’ magmatism. Gondwana Research, 8, 735-756. DOI:10.1016/j.gr.2014.05.013
Carracedo, M., Paquette, J.L., Alonso Olazabal, A., Santos Zalduegui, J.F., García de Madinabeitia, S., Tiepolo, M., Gil Ibarguchi, J.I., 2008. U-Pb dating of granodiorite and granite units of the Los Pedroches batholith. Implications for geodynamic models of the southern Central Iberian Zone (Iberian Massif). International Journal of Earth Sciences (Geologishe Rundschaw), 98: 1659. DOI: 10.1007/s00531-008-0317-0.
Carrilho Lopes, C.J.M., Munha, J.M., Wu, C.T., Oliveira, V.M.J., 1998. O complexo plutónico de Monforte-Santa Eulália (Alentejo-NE, Portugal Central): caracterização geoquímica e considerações petrogenéticas. [In Portuguese] Comunicações do Instituto Geológico e Mineiro, 83, 127-142.
Castro, A., 2013. Tonalite-granodiorite suites as cotectic systems: A review of experimental studies with applications to granitoid petrogenesis. Earth Science Reviews, 124, 68-95.
Castro, A., 2014. The off-crust origin of granite batholiths. Geoscience Frontiers, 5, 63-75.
Castro A., Patiño Douce, A.E., Corretgé, L.G., de la Rosa, J.D., El-Biad, M., El-Hmidi, H., 1999. Origin of peraluminous granites and granodiorites, Iberian Massif, Spain: an experimental test of granite petrogenesis. Contributions to Mineralogy and Petrology, 135, 255-276.
Castro, A., Corretgé, L.G., de la Rosa, J., Enrique, P., Martínez, F.J., Pascual, E., Lago, M., Arranz, E., Galé, C., Fernández, C., Donaire, T., López, S., 2002. Paleozoic magmatism.In: Gibbons, W., Moreno, M. T. (eds.). The Geology of Spain. London. Geological Society, 649pp.
Castro, A., Corretgé, L.G., de la Rosa, J.D., Fernández, C., López, S., Chacón, H., 2003. The appinite–migmatite complex of Sanabria, NW Iberian massif, Spain. Journal of Petrology, 44, 1309-1334.
Castro, A., Gerya, T., García-Casco, A., Fernández, C., Díaz Alvarado, J., Moreno-Ventas, I., Loew, I., 2010. Melting relations of MORB-sediment mélanges in underplated mantle wedge plumes. Implications for the origin of cordilleran-type batholiths. Journal of Petrology, 51, 1267-1295.
Castro, A., Moreno-Ventas, I., Fernández, C., Vujovich, G., Gallastegui, G., Heredia, N., Martino, R.D., Becchio, R., Corretgé, L.G., Díaz-Alvarado, J., García-Arias, M., Liu, D.Y., 2011. Petrology and SHRIMP U-Pb zircon geochronology of Cordilleran granitoids of the Bariloche area, Argentina. Journal of South American Earth Sciences, 32, 508-530.
Castro, A., Vogt, K., Gerya, T.V., 2013. Generation of new continental crust by sublithospheric silicic-magma relamination in arcs: A test of Taylor’s andesite model. Gondwana Research, 23, 68-95.
Charlier, B.L.A., Wilson, C.J.N., Lowenstern, J.B., Blake, S., Van Calsteren, P.W., Davidson, J.P., 2005. Magma generation at a large, hyperactive silicic volcano (Taupo, New Zealand) revealed by U-Th and U-Pb systematics in zircons. Journal of Petrology, 46, 3-32.
Claoué-Long, J.C., Compston, W., Roberts, J., Fanning, C.M., 1995. Two Carboniferous ages: a comparison of SHRIMP zircon dating with conventional zircon ages and 40Ar/39Ar analysis. In: Berggren, W. A., Kent, D. V., Aubry, M. P., Hardenbol, J., (eds.). Geochronology Time Scales and Global Stratigraphic Correlation. Society for Sedimentary Geology (SEPM), Special Publication No. 4, 3-21.
Cocks, L.R.M., Torsvik, T.H., 2006. European geography in a global context from the Vendian to the end of the Palaeozoic. In: Gee, D.G., Stephenson, R.A., (eds.). European Lithosphere Dynamics. Geological Society of London, 83-95.
Cumming, G.L., Richards, J.R., 1975. Ore lead isotope ratios in a continuously changing earth. Earth and Planetary Science Letters, 28, 155-171.
Dias, G., Leterrier, J., Mendes, A., Simoes, P.P., Bertrand, J.M., 1998. U–Pb zircon and monazite geochronology of postcollisional Hercynian granitoids from the Central Iberian Zone (Northern Portugal). Lithos, 45, 349-369.
Dias, G., Simões, P.P., Ferreira, N., Leterrier, J., 2002. Mantle and crustal sources in the genesis of late-Hercynian granitoids (NW Portugal): geochemical and Sr-Nd isotopic constraints. Gondwana Research, 5, 287-305.
Díaz Alvarado, J., Castro, A., Fernández, C., Moreno-Ventas, I., 2011. Assessing bulk assimilation in cordierite-bearing granitoids from the Central System batholith, Spain; experimental, geochemical and geochronological constraints. Journal of Petrology, 52, 223-256.
Díaz Alvarado, J., Fernandez, C., Castro, A., Moreno-Ventas, I., 2013. SHRIMP U-Pb zircon geochronology and thermal modeling of multilayer granitoid intrusions. Implications for the building and thermal evolution of the Central System batholith, Iberian Massif, Spain. Lithos, 175-176, 104-123.
Díez Fernández, R., Pereira, M.F., 2016. Extensional orogenic collapse captured by strike-slip tectonics: Constraints fromstructural geology and U–Pb geochronology of the Pinhel shear zone (Variscan orogen, Iberian Massif). Tectonophysics, 691, 290-310.
Díez Fernández, R., Pereira, M.F., 2017. Strike-slip shear zones of the Iberian Massif: are they coeval? Lithosphere, 9 (5), 726-744.
Díez Fernández, R., Pereira, M.F., Foster, D.A., 2015. Peralkaline and alkaline magmatism of the Ossa-Morena zone (SW Iberia): age, source, and implications for the Paleozoic evolution of Gondwanan lithosphere. Lithosphere, 7, 73-92.
Díez Fernández, R., Arenas, R., Pereira, M.F., Sánchez Martínez, S., Albert, R., Martín Parra, L.M., Rubio Pascual, F.J., Matas, J., 2016. Tectonic evolution of Variscan Iberia: GondwanaLaurussia collision revisited. Earth Science Reviews, 162, 269-292.
Druguet, E., Castro, A., Chichorro, M., Pereira, M.F., Fernández, C., 2014. Zircon geochronology of intrusive rocks from Cap de Creus, Eastern Pyrenees. Geological Magazine, 151(6), 1095-1114.
Escuder Viruete, J., Arenas, R., Martínez Catalán, J.R., 1994. Tectonothermal evolution associated with Variscan crustal extension in the Tormes Gneiss Dome (NW Salamanca, Iberian Massif, Spain). Tectonophysics, 238, 1-22.
Fernández-Suárez, J., Dunning, G.R., Jenner, G.A., Gutiérrez-Alonso, G., 2000. Variscan collisional magmatism and deformation in NW Iberia: constraints from U-Pb geochronology of granitoids. Journal of the Geological Society, 157, 565-576.
Fernández-Suárez, J., Gutiérre-Alonso, G., Johnston, S.T., Jeffries, T.E., Pastor-Galan, D., Jenner, G.A., Murphy, J.B., 2011. Iberian late-Variscan granitoids: Some considerations on crustal sources and the significance of “mantle extraction ages”. Lithos, 123, 121-132.
Ferreira, N., Iglesias Ponce de León, M., Noronha, F., Ribeiro, A., Ribeiro, M.L., 1987. Granitóides da Zona Centro Ibérica e seu enquadramento geodinâmico. In: Bea, F., Carnicero, A., Gonzalo, J.C., López-Plaza, M., Rodríguez Alonso, M.D. (eds.). Geología de los granitoides y rocas asociadas del Macizo Hespérico. Libro Homenaje a L.C. García de Figuerola. Ediciones Rueda, Madrid, 542pp.
Frost, B.R., Barnes, C.G., Collins, W.J., Arculus, R.J., Ellis, D.J. Frost, C.D., 2001. A Geochemical Classification for Granitic Rocks. Journal of Petrology, 42(11), 2033-2048.
Gerya, T.V., Meilick, F.I., 2011. Geodynamic regimes of subduction under an active margin: Effects of rheological weakening by fluids and melts. Journal of Metamorphic Geology, 29, 7-31.
Gerya, T.V., Yuen, D.A., Sevre, E.O.D., 2004. Dynamical causes for incipient magma chambers above slabs. Geology, 32, 8992.
Gonçalves, F., 1971. Subsídios para o conhecimento geológico do Nordeste Alentejano. Memórias dos Serviços Geológicos de Portugal, 18, 62pp. [In Portuguese]
Gonçalves, F., Coelho, A.V.P., 1969-1970. Nota prévia sobre o provável carácter subvulcânico do maciço granítico de Santa Eulália (Alto Alentejo). Boletim do Museu e Laboratório de Mineralogia e Geologia da Faculdade de Ciências, 11(2), 251-263. [In Portuguese]
Gonçalves, F., Assunção, C.F.T., 1970. Carta Geológica de Portugal na escala 1:50,000. Notícia Explicativa da Folha 37A, Elvas: Lisboa, Portugal, Serviços Geológicos de Portugal, 50pp. [In Portuguese]
Gonçalves, F., Assunção, C.F.T., Coelho, A.V.P., 1972. Carta Geológica de Portugal na escala 1:50,000. Notícia Explicativa da Folha 33-C, Campo Maior: Lisboa, Portugal, Serviços Geológicos de Portugal, 41pp.
Gonçalves, F., Fernandes, A.P., 1973. Carta Geológica de Portugal na escala 1:50000. Notícia Explicativa da Folha 32-B, Portalegre: Lisboa, Portugal, Serviços Geológicos de Portugal, 45pp.
Gonçalves, F., Zbyszewsky, G., Coelho, A.V.P., 1975. Carta Geológica de Portugal na escala 1:50000. Notícia Explicativa da Folha 32-D, Sousel: Lisboa, Portugal, Serviços Geológicos de Portugal, 49pp.
González Menéndez, L., 2002. Petrología del batolito granítico de Nisa-Alburquerque. Revista de la Sociedad Geológica de España, 15(3-4), 233-246.
González Menéndez, L., Azor, A., Acosta, A., 2002. Estudio petrológico del batolito de Santa Eulália-Monforte (Alto Alentejo, Portugal). Geogaceta, 32, 147-150.
González Menéndez, L., Azor, A., Pereira, M.D., Acosta, A., 2006. Petrogénesis del plutón de Santa Eulália (Alto Alentejo, Portugal). Revista de la Sociedad Geológica de España, 19(1-2), 69-86.
Gutiérrez-Alonso, G., Fernández-Suárez, J., Weil, A.B., Murphy, J.B., Nance, R.D., Corfu, F., Johnston, S.T., 2008. Selfsubduction of the Pangaean global plate. Nature Geoscience, 1(8), 549–553.
Gutiérrez-Alonso, G., Fernández-Suárez, J., Jeffries, T.E., Johnston, S.T., Pastor-Galán, D., Murphy, J.B., Piedad Franco, M., Gonzalo, J.C., 2011. Diachronous post-orogenic magmatism within a developing orocline in Iberia, European Variscides. Tectonics, 30(5), TC5008, DOI:10.1029/2010TC002845
Hanchar, J.M., Miller, C.F., 1993. Zircon zonation patterns as revealed by cathodoluminescence and backscattered electron images: implications for interpretation of complex crustal histories. Chemical Geology, 110, 1-13.
Heaman, L.M., Bowins, R., Crocket, J., 1990. The chemical composition of igneous zircon studies: implications for geochemical tracer studies. Geochimica et Cosmochimica Acta, 64, 1905-1923.
Hoskin, P.W.O., Schaltegger, U., 2003. The composition of zircon and igneous and metamorphic petrogenesis. Reviews in Mineralogy and Geochemistry, 53, 27-62.
Jesus, A., Munhá, J., Mateus, A., Tassinari, C., Nutman, A., 2007. The Beja layered gabbroic sequence (Ossa-Morena Zone, Southern Portugal): geochronology and geodynamic implications. Geodinamica Acta, 20, 139-157.
Johnston, S.T., Weil, A.B., Gutiérrez-Alonso, G., 2013. Oroclines, Thick and thin. Gelogical Society of America(GSA) Bulletin 125(5-6), 643-663.
Larrea, F.J., Carracedo, M., Alonso Olazabal, A., Donaire, T., Pascual, E., 2004. El batolito de Los Pedroches. In: Vera, J.A. (ed.). Geología de España. Ediciones de la. Sociedad Geológica de España (SGA) - Instituto Geológico y Minero de España (IGME), Madrid, 122-124.
Lee, C.-T.A., Morton, D.M., Kistler, R.W., Baird, A.K., 2007. Petrology and tectonics of phanerozoic continent formation: From island arcs to accretion and continental arc magmatism. Earth and Planetary Science Letters, 263, 370-387.
Lima, S.M., Corfu, F., Neiva, A.M.R., Ramos, M.F., 2011. Dissecting Complex Magmatic Processes: an in-depth U-Pb Study of the Pavia Pluton, Ossa-Morena Zone, Portugal. Journal of Petrology, 53, 1887-1911.
Lima, S.M., Corfu, F., Neiva, A.M.R., Ramos, M.F., 2012. Dissecting complex magmatic processes: an in-depth U-Pb Study of the Pavia Pluton, Ossa-Morena Zone, Portugal. Journal of Petrology, 53(9), 1887-1911.
Linnemann, U., Pereira, M.F., Jeffreis, T., Drost, K., Gerdes, A., 2008. Cadomian Orogeny and the opening of the Rheic Ocean: New insights in the diacrony 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.
López Moro, F.J., López-Plaza, M., Gutiérrez-Alonso, G., Fernández-Suárez, F., López-Carmona, A., Hofmann, M., Romer, R.L., 2017. Crustal melting and recycling: geochronology and sources of Variscan syn-kinematic anatectic granitoids of the Tormes Dome (Central Iberian Zone). A U-Pb LA-ICP-MS study. International Journal of Earth Sciences (Geologische Runbschaw), 1-20, DOI: 10.1007/s00531-017-1483-8.
Ludwig, K.R., 2009. Squid 2.50: a User’s Manual Berkeley Geochronology Center, Berkeley, California, USA, 95pp.
Martínez Catalán, J.R., Arenas, R., Abati, J., Sanchéz Martínez, S., Díaz Garcia, F., Fernández-Suárez, J., González Cuadra, P., Castineiras, P., Gomez Barreiro, J., Díez Montes,A., Clavijo, E., Pascual, F.J., Andonaegui, P., Jeffries, T.E., Alcock, J.E., Díez Fernández, R., López Carmona, A., 2009. A rootless suture and the loss of the roots of a mountain chain: the Variscan Belt of NW Iberia. Comptes Rendus Geoscience, 341, 114-126.
Martínez Catalán, J.R., Rubio Pascual, F.J., Díez Montes, A., Díez Fernández, R., Gómez Barreiro, J., Dias da Silva, Í., González Clavijo, E., Ayarza, P., Alcock, J.E., 2014. The late Variscan HT/LP metamorphic event in NW and Central Iberia: relationships to crustal thickening, extension, orocline development and crustal evolution. Geological Society of London, Special Publication, 405, 225-247.
Mendes, F., 1967-1968. Contribution à l’étude géochronologique, par la méthode au strontium, des formations cristallines du Portugal. Boletim Museu e. Laboratório de Mineralogia e Geologia da Faculdade de Ciências da Universidade de Lisboa, 11(1), 3-155.
Merino Martínez, E., Villaseca, C., Orejana, D., Pérez-Soba, C., Belousova, E., Andersen, T., 2014. Tracing magma sources of three different S-type peraluminous granitoid series by in situ U-Pb geochronology and Hf isotope zircon composition: the Variscan Montes de Toledo batholith (central Spain). Lithos, 200-201, 273-298.
Miller, J.S., Matzel, J.E.P., Miller, C.F., Burgess, S.D., Miller, R.B., 2007. Zircon growth and recycling during assembly of large composite arc plutons. Journal of Volcanology and Geothermal Research, 167, 282-299.
Moita, P., Santos, J.F., Pereira, M.F., 2009. Layered granitoids: interaction between continental crust recycling processes and mantle-derived magmatism. Examples from the Évora Massif (Ossa-Morena Zone, southwest Iberia, Portugal). Lithos, 111(3-4), 125-141.
Moita, P., Santos, J.F., Pereira, M.F., Costa, M.M., Corfu, F., 2015. The quartz-dioritic Hospitais intrusion (SW Iberian Massif) and its mafic microgranular enclaves - evidence for mineral clustering. Lithos, 224-225, 78-100.
Montero, P., Bea, F., Zinger, T.F., Scarrow, J.H., Molina, J.F., Whitehouse, M., 2004. 55 million years of continuous antexis in Central Iberia: single-zircon dating of the Peña Negra Complex. Journal of the Geological Society of London, 161, 255-263.
Montero, P., Bea, F., Corretgé, L.G., Floor, P., Whitehouse, M.J., 2008. U-Pb ion microprobe dating and Sr-Nd isotope geology of the Galiñeiro Igneous Complex. A model for the peraluminous/peralkaline duality of the Cambro-Ordovician magmatism of Iberia. Lithos, 107, 227-238.
Murphy, J.B., 2013. Appinite suites: A record of the role of water in the genesis, transport, emplacement and crystallization of magma. Earth Science Reviews, 119, 35-59.
Neiva, A.M.R., Gomes M.E.P., 2001. Diferentes tipos de granitos e seus processos petrogenéticos: Granitos Hercínicos Portugueses. Memórias da Academia das Ciências de Lisboa XXXIX, 53-95.
Neiva, A. M. R., Williams, I.S., Ramos, J.M.F., Gomes, M.E.P., Silva, M.M.V.G., Antunes, I.M.H.R., 2009. Geochemical and isotopic constraints on the petrogenesis of Early Ordovician granodiorite and Variscan two-mica granites from the Gouveia area, central Portugal, Lithos, 111, 186202.
Neiva, A.M.R., Williams, I.S., Lima, S.M., Teixeira, R.J.S., 2012. U-Pb and 39Ar/40Ar data constraining the ages of the source, emplacement and recrystallization/cooling events from late- to post-D. Variscan granites of the Gouveia area, central Portugal. Lithos, 153, 72-83.
Oliveira, V.M.J., 1975. Cartografia inédita do granito de Santa Eulália (anéis G1, G2 e G3). Unpublished map. Serviço de Fomento Mineiro da Direcção Geral de Geologia e Minas, Portugal.
Oliveira, J.T., Oliveira, V.M.J., Piçarra, J.M., 1991. Traços gerais da evolução tectonoestratigráfica da zona de Ossa Morena, em Portugal: Síntese crítica do estado actual dos conhecimentos. Comunicações dos Serviços Geológicos de Portugal, 77, 3-26.
Orejana, D., Villaseca, C., Valverde-Vaquero, P., Belousova, E.A., Armstrong, R.A., 2012. U-Pb geochronology and zircon composition of late Variscan S- and I-type granitoids from the Spanish Central System. International Journal of Earth Sciences. DOI:10.1007/s00531-012-0750-y
Pereira, M.F., 1999. Caracterização da estrutura dos domínios setentrionais da Zona de Ossa-Morena e seu limite com a Zona Centro-Ibérica, no nordeste Alentejano. Unpublished PhD Thesis, Universidade de Évora, 115pp.
Pereira, M.F., Quesada, C., 2006. Ediacaran to Visean crustal growth processes in the Ossa-Morena zone (SW Iberia), International Geoscience Programme 497 Évora Meeting 2006. Conference Abstracts and Field Trip Guide: Madrid, Spain, Publicaciones IGME, 115pp.
Pereira, M.F., Chichorro, M., Linnemann, U., Eguiluz, L., and Silva, J.B., 2006. Inherited arc signature in Ediacaran and Early Cambrian basins of the Ossa-Morena zone (Iberian Massif, Portugal): Paleogeographic link with European and North African Cadomian correlatives. Precambrian Research, 144, 297-315.
Pereira, M.F., Apraiz, A., Silva, J.B., and Chichorro, M., 2008. Tectonothermal analysis of high temperature mylonitization in the Coimbra-Cordoba shear zone (SW Iberian Massif, Ouguela tectonic unit, Portugal) Evidence of intracontinental transcurrent transport during the amalgamation of Pangea. Tectonophysics, 461, 378-394.
Pereira, M.F., Chichorro, M., Williams, I.S., Silva, J.B., Fernandez, C., Diaz-Azpiroz, M., Apraiz, A., Castro, A., 2009. Variscan intra-orogenic extensional tectonics in the Ossa-Morena Zone (Évora-Aracena-Lora del Rio metamorphic belt, SW Iberian Massif): SHRIMP zircon U–Th–Pb geochronology. In: Murphy, J.B., Keppie, J.D., Hynes, A.J. (eds.). Ancient Orogens and Modern Analogues. Geological Society of London, Special Publications 327, 215-237.
Pereira, M.F., Apraiz, A., Chichorro, M., Silva, J.B., Armstrong, R.A., 2010. Exhumation of high-pressure rocks in northern Gondwana during the early Carboniferous (Coimbra-Cordoba shear zone, SW Iberian Massif): Tectonothermal analysis and U-Th-Pb SHRIMP in-situ zircon geochronology. Gondwana Research, 17, 440-460.
Pereira, M.F., Chichorro, M., Solá, A.R., Silva, J.B., Sánchez-García, T., Bellido, F., 2011. 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.
Pereira, M.F., Chichorro, M., Silva, J.B., Ordonez-Casado, B., Lee, J.K.W., Williams, I.S., 2012. Early Carboniferous wrenching, exhumation of high-grade metamorphic rocks and basin instability in SW Iberia: Constraints derived from structural geology and U–Pb and 40Ar/39Ar geochronology. Tectonophysics, 558-559, 28-44.
Pereira, M.F., Castro, A., Chichorro, C., Fernández, C., Díaz-Alvarado, J. Martí, J., Rodríguez, C., 2014. Chronological link between deep-seated processes in magma chambers and eruptions: Permo-Carboniferous magmatism in the core of Pangaea (Southern Pyrenees). Gondwana Research, 25, 290-308.
Pereira, M.F., Chichorro, M., Moita, P., Santos, J.F., Solá, A.M.R., Williams, I.S., Silva, J.B., Armstrong, R.A., 2015a. The multistage crystallization of zircon in calc-alkaline granitoids: U-Pb age constraints on the timing of Variscan tectonic activity in SW Iberia. International Journal of Earth Sciences. http://dx.DOI.org/10.1007/s00531-015-1149-3.
Pereira, M.F., Castro, A., Fernández, C., 2015b. The inception of a Paleotethyan magmatic arc in Iberia. Geosciences Frontiers, 6, 297-306.
Pereira, M.F., Díez Fernández, R., Gama, C., Hofmann, M., Gärtner, A., Linnemann, U., 2017a. S-type granite generation and emplacement during a regional switch from extensional to contractional deformation (Central Iberian Zone, Iberian autochthonous domain, Variscan Orogeny). International Journal of Earth Sciences, DOI: 10.1007/s00531-017-1488-3.
Pereira, M. F., Gutíerrez-Alonso, G., Murphy, J. B., Drost, K., Gama, C., Silva, J. B., 2017b. 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-281, 383-399.
Pin, C., Fonseca, P.E., Paquette, J.L., Castro, P., Matte, Ph., 2008. The ca. 350Ma Beja igneous complex: a record of transcurrent slab break-off in the southern Iberia Variscan Belt? Tectonophysics, 461, 356-377.
Pinto, M.S., Casquet, C., Ibarrola, E., Corretgé, L.G., Ferreira, M.P., 1987. Síntese geocronológica dos granitóides do Maciço Hespérico. In: Bea, F., Carnicero, A., Gonzalo, J.C., López Plaza, M., Rodríguez Alonso, M.D. (eds.). Geología de los granitoides y rocas associadas del Macizo Hespérico, Libro de Homenaje a L. C. García de Figuerola. Editorial Rueda, Madrid, 69-86pp.
Roberts, M.P., Pin, C., Clemens, J.D., Paquette, J.L., 2000. Petrogenesis of mafic to felsic plutonic rock associations: the calc-alkaline Quérigut complex, French Pyrenees. Journal of Petrology, 41, 809-844.
Sánchez-García, T., Bellido, F., Pereira, M.F., Chichorro, M., Quesada, C., Pin, C., Silva, J.B., 2010. Rift-related volcanism predating the birth of the Rheic Ocean (Ossa-Morena zone, SW Iberia). Gondwana Research, 17, 392-407.
Sánchez-García, T., Pereira, M.F., Bellido, F., Chichorro, M., Silva, J.B., Valverde-Vaquero, P., Pin, C., Solá, A.R., 2013. Early Cambrian granitoids of North Gondwana margin in the transition from a convergent setting to intra-continental rifting (Ossa-Morena zone, SW Iberia). International Journal of Earth Sciences, 103, 1203-1218.
Sant’Ovaia, H., Nogueira, P., Carrilho Lopes, J., Gomes, C., Ribeiro, M.A., Martins, H.C.B., Dória, A., Cruz, C., Lopes, L., Sardinha, R., Rocha, A., Noronha, F., 2015. Building up of a nested granite intrusion: magnetic fabric, gravity modelling and fluid inclusion planes studies in Santa Eulália Plutonic Complex (Ossa Morena Zone, Portugal). Geological Magazine, 152 (4), 648-667.
Simancas, J.F., Tahiri, A., Azor, A., Lodeiro, F.G., Poyatos, D.J.M., Hadi, H.E., 2005. The tectonic frame of the Variscan– Alleghanian orogen in southern Europe and northern Africa. Tectonophysics, 98, 181-198.
Solá, A.R., Williams, I.S., Neiva, A.M.R., Ribeiro, M.L., 2009. U–Th–Pb SHRIMP ages and oxygen isotope composition of zircon from two contrasting late Variscan granitoids, NisaAlbuquerque batholith, SW Iberian Massif: Petrologic and regional implications. Lithos, 111, 156-167.
Solé, J., Cosca, M., Sharp, Z., Enrique, P., 2002. 40Ar/39Ar geochronology and stable isotope geochemistry of lateHercynian intrusions from north-eastern Iberia with implications for argon loss in K-feldspar. International Journal of Earth Sciences, 91, 865-881.
Stampfli, G.M., Borel, G.D., 2002. A plate tectonic model for the Paleozoic and Mesozoic constrained by dynamic plate boundaries and restored synthetic oceanic isochrones: Earth and Planetary Science Letters, 196, 17-33
Stampfli, G.M., Kozur, H., 2006. Europe from the Variscan to the Alpine cycles. In: Gee, D.G., Stephenson, R., (eds.). European Lithosphere Dynamics. Geological Society of London Memoir, 32, 57-82.
Stampfli, G.M., Hochard, C., Vérard, C., Wilhem, C., von Raumer, J., 2013. The formation of Pangea. Tectonophysics, 593, 1-19.
Teixeira, R.J.S., Neiva, A.M.R., Silva, P.B., Gomes, M.E.P., Andersen, T., Ramos, J.M.F., 2011. 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.
Turnbull, R., Weaver, S., Tulloch, A., Cole, J., Handler, M., Ireland, T., 2010. Field and geochemical constraints on maficfelsic interactions, and processes in high-level arc magma chambers: an example from the Halfmoon Pluton, New Zealand. Journal of Petrology, 51(7), 1477-1505.
Valle Aguado, B., Azevedo, M.R., Schaltegger, U., Martínez Catalán, J.R., Nolan, J., 2005. U-Pb zircon and monazite geochronology of Variscan magmatism related to synconvergence extension in Central Northern Portugal. Lithos, 82, 169-184.
Valverde-Vaquero, P., 1997. An integrated field, geochemical and U-Pb geochronological study of the southwest Hermitage Flexure (Newfoundland Appalachians, Canada) and the Sierra De Guadarrama (Iberian Massif, Central Spain): a contribution to the understanding of the geological evolution of circum- Atlantic Peri-Gondwana. Unpublished PhD Thesis. Memorial University of Newfoundland, 312pp.
Valverde-Vaquero, P., Bento dos Santos, T., Clavijo, E.G., Díez Montes, A., Ribeiro, M.L., Solá, A.R., Dias da Silva, Í., 2011. The Berlengas Archipelago granitoids within the frame of the Variscan Orogeny, W-Portugal: new data and insights. VII Hutton Symposium on Granites and Related Rocks. Abstracts Book, Avila, Spain, 131-132.
Villaseca, C., Bellido, F., Pérez-Soba, C., Billström, K., 2009. Multiple crustal sources for post-tectonic I-type granites in the Hercynian Iberian Belt. Mineralogy and Petrology, 96, 197-211.
Villaseca, C., Orejana, D., Belousova, E.A., 2012. Recycled metaigneous crustal sources for S- and I-type Variscan granitoids from the Spanish Central System batholith: constraints from Hf isotope zircon composition. Lithos, 153, 84-93.
Vogt, K., Gerya, T.V., Castro, A., 2012. Crustal growth at active continental margins: Numerical modeling. Physics of the Earth and Planetary Interiors, 192-193. 1-20.
Weil, A., Gutiérrez-Alonso, G., Johnston, S.T., Pastor-Galán, D., 2012. Kinematic constraints on buckling a lithospheric-scale orocline along the northern margin of Gondwana: A geologic synthesis. Tectonophysics, 582, 25-49.
Williams, I.S., 1998. U-Th-Pb Geochronology by Ion Microprobe. In: McKibben, M. A., Shanks III, W. C., Ridley, W. I. (eds.). Applications of microanalytical techniques to understanding mineralizing processes. Reviews in Economic Geology, 7, 1-35.
Williams, I.S., Hergt, J.M., 2000. U-Pb dating of Tasmanian dolerites: a cautionary tale of SHRIMP analysis of high-U zircon. In: J.D. Woodhead, J.M. Hergt, W.P. Noble (eds.). Beyond 2000: new frontiers in isotope science. Lorne, Abstracts and Proceedings, 185-188.
ELECTRONIC APPENDIX I

SHRIMPTOOLS data processing; IBERSIMS laboratory, University of Granada, Spain
Errors are at one sigma level. The error in 206/238 averaging the standard has been already propagated. Th/U ratios are atomic (232/238), not in weight
Point-to point errors, calculated on replicates of the TEMORA standard, are: 0.36% for 206Pb/238U, and 0.32% for 207Pb/206Pb

SHRIMPTOOLS data processing; IBERSIMS laboratory, University of Granada, Spain
Errors are at one sigma level. The error in 206/238 averaging the standard has been already propagated. Th/U ratios are atomic (232/238), not in weight
Point-to point errors, calculated on replicates of the TEMORA standard, are: 0.36% for 206Pb/238U, and 0.32% for 207Pb/206Pb

SHRIMPTOOLS data processing; IBERSIMS laboratory, University of Granada, Spain
Errors are at one sigma level. The error in 206/238 averaging the standard has been already propagated. Th/U ratios are atomic (232/238), not in weight
Point-to point errors, calculated on replicates of the TEMORA standard, are: 0.36% for 206Pb/238U, and 0.32% for 207Pb/206Pb

SHRIMPTOOLS data processing; IBERSIMS laboratory, University of Granada, Spain
Errors are at one sigma level. The error in 206/238 averaging the standard has been already propagated. Th/U ratios are atomic (232/238), not in weight
Point-to point errors, calculated on replicates of the TEMORA standard, are: 0.36% for 206Pb/238U, and 0.32% for 207Pb/206Pb

SHRIMPTOOLS data processing; IBERSIMS laboratory, University of Granada, Spain
Errors are at one sigma level. The error in 206/238 averaging the standard has been already propagated. Th/U ratios are atomic (232/238), not in weight
Point-to point errors, calculated on replicates of the TEMORA standard, are: 0.36% for 206Pb/238U, and 0.32% for 207Pb/206Pb

SHRIMPTOOLS data processing; IBERSIMS laboratory, University of Granada, Spain
Errors are at one sigma level. The error in 206/238 averaging the standard has been already propagated. Th/U ratios are atomic (232/238), not in weight
Point-to point errors, calculated on replicates of the TEMORA standard, are: 0.36% for 206Pb/238U, and 0.32% for 207Pb/206Pb

SHRIMPTOOLS data processing; IBERSIMS laboratory, University of Granada, Spain
Errors are at one sigma level. The error in 206/238 averaging the standard has been already propagated. Th/U ratios are atomic (232/238), not in weight
Point-to point errors, calculated on replicates of the TEMORA standard, are: 0.36% for 206Pb/238U, and 0.32% for 207Pb/206Pb

SHRIMPTOOLS data processing; IBERSIMS laboratory, University of Granada, Spain
Errors are at one sigma level. The error in 206/238 averaging the standard has been already propagated. Th/U ratios are atomic (232/238), not in weight
Point-to point errors, calculated on replicates of the TEMORA standard, are: 0.36% for 206Pb/238U, and 0.32% for 207Pb/206Pb