Petrology and Geochemistry of Highly Differentiated Tholeiitic Magmas: Granophyres in the Messejana–Plasencia Great Dyke (Central Iberia)

Abstract

The Messejana–Plasencia great dyke (MPGD) is a Late Triassic tholeiitic gabbro intrusion related to the Central Atlantic Magmatic Province. Its large outcrop extent (~530 km), combined with its prolongation below the Duero basin (additional 100 km), makes it one of the world’s largest dykes known. We have studied felsic granophyric bodies appearing in its northernmost segment at different scales, from mm-sized (interstitial micrographic pockets) to felsic dykes of up to 10 m thick and 1.5 km long, intruding within the gabbros. Significant differences exist in the mineral and whole-rock composition of gabbros and granophyres, including the Sr–Nd isotopic ratios. The chemical variation in the gabbros is coherent with fractionation of olivine, clinopyroxene and plagioclase at depth. However, the presence of a compositional gap between gabbros and granophyres (absence of intermediate compositions) and the formation of these late-stage intergranular felsic melts within the gabbro mesostasis suggest that they could be derived by liquid immiscibility. The Sr–Nd isotopic heterogeneity in the MPGD gabbros and the presence of zircons with Variscan ages (~286 Ma), inherited from granulitic rocks, indicate that the mafic magmas experienced some degree of lower crust assimilation during fractionation close to the Moho depth. On the contrary, the scarce xenocrystic Variscan zircon crystals found in a granophyric dyke within the MPGD gabbro display similar textures and ages (~299 Ma) to those of the country rock granites and point to contamination at a different crustal level.

1. Introduction

The Messejana–Plasencia great dyke (MPGD) intruded the Variscan basement of Western Iberia at circa 200 Ma [1]. It has an outcropping length of 530 km and a SW–NE orientation, from the Alentejo (southwestern Portugal) to Avila (central Spain), plus more than ~100 km below the Duero Cenozoic Basin [2]. Such great length, together with its variable width, ranging from 25 to 200 m (Figure 1), makes it the largest dyke and the longest simple tectonic structure crosscutting the Iberian Peninsula [2] and even one of the world’s largest known dykes (slightly longer than the Great Dyke of Zimbabwe, although much thinner). It belongs to the so-called Central Atlantic Magmatic Province (CAMP) and is related to the radiating dyke swarms generated during the opening of the Central Atlantic Ocean in the Triassic–Jurassic transition [3,4].

The MPGD has been studied from the early 1950s, mostly in the Portuguese sector [5,6], to recent years [7,8,9]. Most researchers agree that the dolerites of the MPGD might be the result of partial melting of lithospheric mantle sources enriched by subduction-related components and that they also reflect some degree of crustal contamination, either at the mantle source [9,10] or en route due to lower crustal rocks assimilation [7].

Figure 1

Geological sketch maps showing (A) the Messejana–Plasencia great dyke (MPGD) in the Iberian Peninsula modified after [11] and (B) the northeasternmost segment of the MPGD. Sample locations are shown with stars including the following: (1) Cid quarry (103309, 103310, 103311); (2) NE Sanchorreja (106285, 106286, 106289, 113105, 113106); (3) SW Sanchorreja (106284); (4) Balbarda (103301); (5) Villatoro (103318); (6) Villatoro mountain pass (103317, 104367, 106292, 113104) and (7) Villafranca quarry (103326).

[Figure omitted. See PDF]

In its northern sector, the MPGD shows dispersed granophyric veinlets, dykelets and dykes intruding the main gabbroic dolerite dyke [12]. This is a common feature to other gabbroic dykes of the Central Atlantic Magmatic Province (CAMP), as granophyres have been described in the Shelburne dyke of Nova Scotia (Canada) [13], and in the doleritic dyke swarms of the Guyana craton and northern Brazil [14], among others. Rhyolitic interstitial melts have also been described in tholeiitic continental basalts from Nova Scotia (CAM Province) [15]. Two accepted origins for these granophyres are based on open-system magmatic processes, either combining crystal fractionation and variable crustal assimilation [14] or by partial melting of the crust induced by the high thermal energy of gabbroic intrusions [13]. Nevertheless, the formation of granophyres in tholeiitic magmas linked to a closed-system magmatic evolution involving extreme crystal fractionation (e.g., [16]) or melt immiscibility has also been proposed in tholeiitic bodies (e.g., [15,17,18]).

The present work is focused on the easternmost segment of the MPGD, between Villafranca de la Sierra and Sanchorreja villages (45 km long and with an average width of 100 m) (Figure 1). We analysed the mineral chemistry and the major, trace element and isotopic (Sr–Nd) bulk rock geochemistry of the small-volume felsic rocks associated with the dyke, together with data on the zircon populations (zircon petrography and U–Pb geochronology) of two samples (gabbro and granophyre). This study represents the first attempt to constrain the origin of these SiO2-rich felsic melts and their genetic relationship with the host gabbros of the MPGD.

2. Geological Setting

The Iberian Massif is a basement area representing the westernmost part of the European Variscan Belt, a consequence of the continental collision of Gondwana and Laurussia. It outcrops in central and western Iberia and has been subdivided into several zones according to their stratigraphic, structural and petrologic characteristics (e.g., [19]). These zones may be considered independent microterranes originally close to the Gondwana margin or internal–external domains assembled during the Variscan collision in a single terrane (e.g., [20]). The Central Iberian Zone (CIZ) of the Iberian Massif is representative of the innermost region of the Iberian Variscan Belt and is mainly composed of Neoproterozoic to Paleozoic metasedimentary and metaigneous rocks, together with abundant granitic plutons, which intruded during the last stages of the orogeny.

Several post-Variscan magmatic events have been recognised in the CIZ, which consist mainly of small subvolcanic intrusions, dykes and volcanic products of alkaline, tholeiitic or transitional composition. The oldest one corresponds to a suite of volatile-rich mafic–ultramafic alkaline dykes (diabases and alkaline lamprophyres) of Permian age (274–264 Ma) ([21] and references therein) accompanied by monzogabbroic to syenitic porphyries. They are found in central Iberia [21,22] and northern Portugal [23], and have been related to a rifting context after the Permo-Carboniferous extension in western Europe, resulting in the rise of the asthenospheric mantle and melting of lithospheric and asthenospheric mantle sources [8,21].

The main post-Variscan tholeiitic magmatism is represented by the large Messejana–Plasencia dyke–fault system, which is a conspicuous geological intrusion within the Iberian Massif due to its length (>530 km) and width (25–200 m; locally > 300 m). It is mainly made of gabbros outcropping discontinuously from its SW endpoint (Alentejo) to the NE until it is covered by Cenozoic sediments in the Duero basin. According to the most accurate geochronology studies based on ⁴⁰Ar/³⁹Ar dating, the dyke was likely emplaced as a result of several magma injections during the Jurassic–Triassic transition (~206–200 Ma) [1,10,24]. However, paleomagnetic data on the MPGD, as well as other CAMP mafic dykes, point to magma intrusion taking no longer than 1 Ma [25]. The main petrographic features and preliminary data on the whole-rock and mineral geochemistry have been described in several studies focused on different sectors along the dyke [7,9,11,12,26,27,28,29]. These works have highlighted the presence of variable textures in the gabbros from coarse-grained to fine-grained (sometimes associated with core–margin dyke zonation) and scarce inner pegmatoid bodies. The whole-rock composition of the MPGD gabbros reveals an affinity to low-Ti continental tholeiites and a moderate enrichment in some incompatible trace elements (LILE), and the Sr–Nd–Pb isotope signatures are in accordance with a significant degree of crystal fractionation and small percentages of assimilation of lower-crustal rocks [7].

The easternmost segment of the MPGD intrudes mainly in Variscan peraluminous granitic rocks. Between Villafranca de la Sierra and Sanchorreja (Figure 1), the dyke includes diverse granophyric bodies cutting the mafic rocks in its central zone. The presence of occasional quartz-bearing leucocratic bodies of granophyric texture inside the MPGD has already been described [12,26], mainly in this NE segment, but also in its central part, where these felsic structures are less abundant. However, very scarce geochemical data exist on these rocks, and their origin has not been discussed yet.

3. Materials and Methods

The major element mineral composition was determined at the Centro de Microscopía Electrónica “Luis Bru” (Complutense University of Madrid) using a JEOL JZA-8900 M electron microprobe with four wavelength dispersive spectrometers. The accelerating voltage during the analyses was 15 kV and the electron beam current was 20 nA, with a beam diameter of 5 μm. The counting time was 10 s on the peak and 5 s on each background position. The following mineral standards were used: sillimanite for Al; albite for Si and Na; almandine for Fe and Mn; kaersutite for Mg; microcline for K; ilmenite for Ti; Ni alloy for Ni and chromite for Cr. Corrections were made using the ZAF method. Analytical precision, expressed as a relative percentage, is 0.5–6% for oxides with concentration > 1.5 wt.% and <10% for oxides with concentration < 1.5 wt.%.

The whole-rock composition of 8 samples (7 granophyres and 1 gabbro) was determined at the Actlabs laboratories (Ancaster, ON, Canada) for major and trace elements. The fusion technique was used by employing lithium metaborate/tetraborate. The molten bead was then digested in a weak nitric acid solution. The samples were analysed by inductively coupled plasma optic emission spectrometry (ICP–OES) for major elements, while trace elements were determined by ICP mass spectrometry (ICP–MS), using a Thermo Jarrell Ash ENVIRO II simultaneous and sequential ICP or a Perkin Elmer Optima 3000 ICP. Uncertainties in major elements were always below 0.01%. All elements were analysed under the control of certified international standards (SY-3, NIST 694, W-2a, DNC-1, BIR-1 and GBW 07113). Further details on calibration, detection limits of trace elements, etc., can be found at http://www.actlabs.com/ (accessed on 2 February 2024).

Five granophyre samples and one gabbro were selected to determine their Sr–Nd isotopic signatures. The analyses were performed at the CAI de Geocronología y Geoquímica Isotópica of the Complutense University in Madrid, using an automated VG Sector 54 multicollector thermal ionisation mass spectrometer with data acquired in multidynamic mode. Isotopic ratios of Sr and Nd were measured on an aliquot of whole-rock powder. The Sr–Nd analytical procedures used in this laboratory have been described elsewhere [30]. Repeated analysis of NBS 987 gave ⁸⁷Sr/⁸⁶Sr = 0.710249 ± 30 (2σ, n = 15) and for the JM Nd standard, ¹⁴³Nd/¹⁴⁴Nd = 0.511809 ± 20 (2σ, n = 13). The 2σ error on calculated ε(Nd) values was ±0.4.

Two samples were selected for U–Pb geochronology including one granophyre and one gabbro. Zircons were separated and concentrated using standard techniques, and a representative selection was handpicked from each sample. They were mounted in epoxy resin for microanalytical analysis at the IBERSIM laboratories (University of Granada, Spain), together with the reference zircons TEMORA 1, SL13 and REG. The mount was polished, and the zircon central portions were imaged with transmitted and reflected light and with cathodoluminescence on a scanning electron microscope. U, Th and Pb isotopes were analysed on the Sensitive High-Resolution Ion Microprobe (SHRIMP II) with each analysis consisting of six scans through the relevant mass range. The primary beam, composed of ¹⁶O¹⁶O⁺, was set to an intensity of about 5 nA, with a 120 μm Kohler aperture, which generates 17 μm × 20 μm elliptical spots on the target. All calibration procedures were performed on the standards included on the same mount. One TEMORA zircon standard was analysed for every four unknown analyses. Data were reduced using SHRIMPTOOLS software, specifically developed for IBERSIMS by Fernando Bea. More details can be found at http://www.ugr.es/~ibersims/ibersims/Sample_analysis_and_methods.html (accessed on 2 February 2024). Pb/U ratios were normalised relative to a value of 0.06683 for the ²⁰⁶Pb/²³⁸U ratio of the TEMORA reference zircon, equivalent to an age of 417 ± 1 Ma (IDTIMS in zircon) [31]. Concentration data were normalised against zircon standard SL13 (U = 210 ppm) [32]. Uncertainties given for individual analyses (ratios and ages) are at the 1σ level, and all ages are 204-corrected ²⁰⁶Pb/²³⁸U ages. Concordia and probability density plots were drawn using Isoplot 3.0 software [33].

4. Results

This section describes the outcrops, petrography, mineral chemistry and whole-rock geochemistry of granophyre samples from the MPGD. The gabbros, which constitute the main lithology of this dyke, have been studied in detail in previous works [7,10,29]. Nevertheless, a brief summary of their petrology and mineral chemistry is also presented in this section, together with the new whole-rock analysis of a gabbro sample.

4.1. Outcrop Features and Petrography

The MPGD granophyres appear on variable scales. They can be major leucogranite bodies forming veinlets, irregular dykes of small to moderate size (3 cm to 1 m) or a large dyke (up to 10 m thickness and more than 1.5 km length) within the main gabbroic intrusion (Figure 2A,B); but also interstitial pockets with micropegmatitic quartz–alkali feldspar intergrowths, usually radiating from plagioclase crystals of the diabase (Figure 2C–F).

These leucocratic rocks are found in Sanchorreja, El Cid quarry (north of Sanchorreja), Villatoro pass and the Villafranca de la Sierra quarry (Figure 1). The 15 m wide leucogranitic dyke outcrops in Sanchorreja, show small miarolitic cavities, scarce biotite and small (2–3 cm) fine-grained mafic enclaves. In the rest of the outcrops, granophyres mainly form up to 5–30 cm wide dykelets with sharp contacts and scarce miarolitic cavities (Figure 2A,B). Moreover, abundant interstitial micrographic quartz–alkali feldspar intergrowths appear within the crystal framework of the host diabase (Figure 2C–F) in this NE MPGD segment. These features are consistent with the more evolved composition of this dyke sector when compared with the younger SW segment (Alentejo, Portugal), which also displays less radiogenic BSE-like Sr–Nd isotopic ratios (εNd > 0) [8]. Similar intergrowth textures are also characteristic of the granophyre dykes (Figure 3).

The gabbros from the MPGD display mainly a coarse-grained diabasic texture in its core, which varies to fine-grained subofitic towards the contact. The mineralogy is typical of continental tholeiites with plagioclase, clinopyroxene, Ca-poor pigeonite and orthopyroxene as major phases, as well as minor amphibole, biotite, ilmenite, ulvöspinel, apatite, quartz and alkali feldspar. Accessory olivine may also be present, although usually as a pseudomorph (e.g., [7]). Interstitial quartz between plagioclase laths is common and it sometimes forms micrographic intergrowth with alkali feldspar or albite (Figure 2C–F). Very coarse-grained amphibole-rich pegmatoids are occasional, with the amphibole growing interstitially to plagioclase laths.

The pyroxene composition in the gabbros evolves from Mg-rich to Fe–Ti-rich [29]. This variation is observed in several pyroxene types of the MPGD gabbros. Thus, enstatites (Mg# = 0.80–0.72) evolve towards Ca-poor pigeonite (Mg# = 0.75–0.57) and ferroan enstatite (Mg# = 0.68–0.56). Augite is more abundant than low-Ca pyroxenes (orthopyroxene and pigeonite) and exhibits a broad composition of Wo_39–34, En_44–35 and Fs_22–25 (towards Wo_38–46, En_21–35 and Fs_23–36 in pegmatitic facies). Plagioclase displays a wide compositional range from An₇₅ to An₂₀ (sometimes as core-to-rim zoning), although it frequently shows narrower values typical of labradorite composition (An_55–65). Values of An_18–03 are observed in interstitial plagioclases associated with quartz–alkali feldspar microgranophyric intergrowths (Figure 4A). H₂O-rich mafic minerals (amphibole and biotite) appear mostly as accessory phases rimming previous augite or pigeonite crystals, except in the amphibole-rich pegmatitic gabbros. Amphibole ranges from Fe-edenite to Mg-hornblende, with Mg# in the range of 0.3–0.6 and a negative correlation between Mg and Ti [29]. Biotite composition is typical of subalkaline mafic rocks and is characterised by high Fe (Mg# = 0.2–0.45) and TiO₂ (>3 wt%) contents. It is worth mentioning that biotite and apatite crystals of the MPGD gabbros are mostly concentrated in the felsic interstitial granophyric intergrowths, and their chemical composition is very similar and overlaps that of associated granitic granophyres (Figure 4B–F).

Interstitial leucocratic granophyric intergrowth of quartz and albite or alkali feldspar can be found within the diabases (Figure 2C,D), mostly in the MPGD NE segment (Figure 1). The gabbros showing these felsic micro-aggregates are also associated with small felsic granophyre dykes of variable thickness (from a few cm to up to 2 m) and sharp contacts, along with subvertical irregular veins permeating the main gabbroic body. Granophyre mineralogy is simple: quartz, alkali feldspar and scarce albite as major phases, and minor biotite, titanite, acicular apatite, tourmaline and ilmenite, as accessories. The small veins or dykes display an anhedral microgranular texture, also with some subhedral plagioclases, whereas larger dykes develop granophyric texture, with quartz–alkali feldspar intergrowth among plagioclase and quartz crystals (Figure 3). Miarolitic cavities are common in granophyric dykes, filled with prehnite, epidote, Fe-Mg chlorite and fibrous clay minerals.

4.2. Granophyre Mineral Chemistry

4.2.1. Feldspars

Euhedral and subhedral plagioclase crystals with simple twinning, as well as those from micrographic intergrowths, were analysed (Table 1). In most cases, the composition is that of pure albite (Ab₉₉–Ab₉₄). Alkali feldspars are orthoclase (Or₉₉ Or₈₁) or sanidine (~Or₈₀, rarely reaching Or₄₈) (Figure 4A). The composition of both feldspars is similar to that of those present in the micrographic intergrowths of gabbros (interstitial granophyres of Figure 4A and Table 1).

4.2.2. Biotite

Biotite composition varies significantly depending on the outcrop considered (Figure 4B,C). Two broad biotite types can be distinguished as follows: (1) higher Mg–Ti biotite (mainly from Sanchorreja outcrop 2; Figure 1) and (2) lower Mg–Ti biotite (from other outcrops). These elements are negatively correlated (Figure 4B).

As a whole, biotite shows a negative correlation between Al–Ti and Mg#, and it displays lower TiO₂ and higher Al₂O₃ and Na₂O contents when compared with biotite from the host gabbros (Figure 4B–D). It also shows fairly variable K₂O values. Biotite from the Villafranca and Villatoro samples has a composition closer to that of the gabbro biotites, mainly with those from the pegmatoid facies (Figure 4B,C).

When plotted in the Mg–Al_total diagram of [36] (not shown), granophyre biotite falls within the field of calc-alkaline series. Most granophyre biotites also plot in calc-alkaline fields according to the classification of [35] (Figure 4E). In these diagrams, gabbro biotites plot in Al-poor subalkaline fields, although there is some overlap with the granophyre compositional field.

The high F content of biotite from these felsic granophyres is remarkable (up to 3.1 wt%, Table 2), which is much higher than that of micas from the MPGD gabbros (mostly F < 0.36 wt%) or from Variscan granites (mostly F < 1.27 wt%).

4.2.3. Apatite

Slight chemical differences can also be observed in apatite from the granophyres. The F/Cl ratio of apatite from the Sanchorreja samples (large dyke, outcrop 2 in Figure 1) is much lower than that of the Villatoro dykelets, and also from some gabbro apatites, although granophyre apatite can show higher Na₂O + K₂O contents (Figure 4F, Table 2).

4.3. Whole-Rock Composition

4.3.1. Major and Trace Elements

The MPGD granophyres display a silica-rich composition (SiO₂ = 70.55–74.85 wt%) and draw a large compositional gap with the associated gabbros (Figure 5). Their composition is slightly peraluminous and shows high K contents (K/Na ratios mostly > 1.5), plotting in high-K calc-alkaline to shoshonitic fields, in contrast to the host gabbros (Figure 5C,D).

Granophyre major element contents are somewhat homogeneous (Table 3) and characterised by MgO < 1.8 wt%, Fe₂O₃ < 2.9 wt% and low TiO₂ (<0.33 wt%) and CaO (<0.71 wt%) values when compared with the MPGD gabbros (Figure 6). Alkalies are more variable in composition, displaying low Na₂O concentrations (mostly below 2.5 wt%) and relatively high K₂O (3.2–8.6 wt%). A subtle correlation can be observed within granophyre samples with increasing SiO₂ and Al₂O₃ and decreasing Fe₂O₃ and TiO₂ towards lower MgO (Figure 6). Their trace element composition is characterised by relatively high LILE (Ba = 217–741 ppm; Rb = 105–278 ppm), Th (3.2–10.2 ppm), U (2–3.5 ppm) and Pb (up to 36 ppm) values and moderate to low contents of Sr (22–117 ppm), other HFSE (Zr = 38–107 ppm; Hf = 1.5–3 ppm; Nb = 6.1–10.4 ppm; Ta = 1.1–1.9 ppm) and REE (e.g., Ce = 12.3–38.6 ppm; Eu = 0.29–0.57 ppm; Yb = 0.88–2.19 ppm) (Figure 7). Such composition yields a slightly fractionated chondrite-normalised LREE pattern (flat for the HREE) with a marked negative Eu anomaly and a Primitive Mantle-normalised spectrum with high relative abundance of Rb, Ba, Th, U and K, a strong positive Pb anomaly and relatively deep throughs at Nb, Sr, Eu and Ti (Figure 8). A slight positive Ta peak in the composition of the studied granophyres is also remarkable. The Cr concentrations are variable, reaching relatively high values (200 ppm). The fact that the samples with higher Cr contents (113105 and 113106) are the least evolved granophyres, and the lack of correlation between Cr and the other metallic elements (e.g., Mg, Ti, Ni, Co, V, Sc), suggests that the highest Cr values are actually not reliable but an artefact likely caused by Cr addition during sample crushing using a Cr-rich disc mill.

The only gabbroic sample analysed for this work (sample 113104) plots towards the differentiated pole of the MPGD gabbros compositional field (Figure 6). It displays low MgO (3.41 wt%) and relatively high TiO₂ (1.68 wt%) and Fe₂O₃ (12.47 wt%). Its trace element contents and the Chondrite- and Primitive Mantle-normalised patterns also overlap those of the literature MPGD gabbros, plotting towards the higher values for most incompatible elements (Figure 7 and Figure 8). The overall Cr contents of the gabbros are slightly heterogeneous. It should be noted that the gabbros collected by our research team [8] in the northeastern segment of the MPGD display a tendency towards higher Cr concentrations when compared with the samples analysed by other authors (Figure 7). As mentioned above for the granophyres, this discrepancy is only shown by Cr and points to contamination during sample crushing.

A comparison between the MPGD granophyre composition and the main felsic SCS Variscan magmatism can be found in Figure 6 and Figure 7. It can be appreciated that the granophyre compositional field broadly overlaps the major and trace element field of the SCS granites and also reproduces similar elemental correlations between MgO and SiO₂, Fe₂O₃, TiO₂, P₂O₅, Ba, Rb and Cr. However, the MgO–Al₂O₃ correlation of the granophyres does not match that of the granites, and the alkalies display a more heterogeneous and distinct composition: the Na₂O values are clearly lower and the K₂O concentrations are broadly higher than those in the Variscan granites (Figure 6). Several differences are also apparent in some trace elements, mainly the low Sr, Zr, Nb, Th, U and REE contents (Figure 7).

When comparing the composition of the granophyres with that of the host mafic rocks (mainly gabbros), there is a marked gap in major and trace element compositions between the Mg-poor pole of these mafic rocks and the granophyres (e.g., SiO₂ ~52–71 wt%, Fe₂O₃ ~3–9 wt%, Rb ~49–105 ppm, Sr ~61–180; Figure 6 and Figure 7). Moreover, the chemical trend drawn by the mafic rocks undergoes an abrupt change if the granophyre pole is considered as a hypothetical cogenetic line of evolution, implying a sudden increase in Si, K, Ba, Rb and U and decrease in Ti, Al, Fe, Mn, Ca, Na, Sr, Zr, Hf, Nb and REE (Figure 6 and Figure 7).

4.3.2. Sr–Nd Isotopes

The Sr–Nd isotopic ratios of the analysed samples were calculated to 200 Ma, as well as the rest of magmatic rocks plotted in Figure 9 (MPGD gabbros, SCS mafic and felsic Variscan magmatic rocks and SCS orthogneisses). The granophyre data represent the first attempt to constrain the isotopic composition of these evolved felsic rocks from the Messejana–Plasencia dyke. They display fairly variable radiogenic Sr signatures [(⁸⁷Sr/⁸⁶Sr)₀ = 0.708292–0.717215] but restricted εNd values (−3.7 to −5.8) (Figure 9 and Table 4). The most Sr radiogenic values are found in miarolitic granophyres from the large felsic dyke intruding the MPGD in the NE area of Sanchorreja (outcrop 2 in Figure 1). The granophyre Sr–Nd isotopic composition is intermediate between the less radiogenic values of the MPGD gabbros and the field of the SCS granites and orthogneisses wall rocks (Figure 9). The gabbro analysed for this study (sample 113104) shows a BSE-like composition [(⁸⁷Sr/⁸⁶Sr)₀ = 0.706401; εNd = −1.0] which plots within the isotopic field previously defined for the MPGD gabbros, which, in turn, matches the composition of the SCS Variscan mafic magmatic series, mainly of calc-alkaline affinity [8] (Figure 9).

4.3.3. U–Pb Zircon Geochronology

Zircons from two samples of the MPGD (gabbro 103324 and granophyre 113106) were separated for U–Pb geochronology. The gabbro provided a higher number of zircon grains, from which 60 were selected, whereas only 11 zircons were extracted from the granophyre. Their size, external morphology and inner texture were studied by transmitted and reflected microscopy and CL imaging, and significant differences were found between the two samples.

The gabbro zircons are generally small: most of them are below 100 μm (91%), and only a few (9%) can reach up to 200 μm. These latter larger grains are irregular in shape, likely representing prism fragments, and develop marked oscillatory zoning (A-2.1 in Figure 10). However, the bulk of zircons in this sample are small rounded grains or stubby prisms, half of them displaying subhedral, variably resorbed inherited cores. These cores show blurred or oscillatory zoning and a dark thin rim, and they are surrounded by a dark or bright regrowth, sometimes formed by several concentric bands. Zircon grains without inherited cores and displaying a homogeneous dark or bright CL response are also relatively abundant (B-14.1 in Figure 10).

The 11 zircon grains from the granophyre sample have variable sizes (50–170 μm) and shapes, predominating the euhedral stubby or elongated prisms, with aspect ratios from 1:2 to 1:4 and a tendency to develop bipyramidal ends (Figure 10). CL images reveal the frequent presence of anhedral to subhedral inherited cores with blurred or oscillatory zoning. These cores appear in 40% of the zircons (stubby prisms) and represent more than 50% of their volume. The surrounding rims and the elongated prisms without cores are characterised by magmatic oscillatory zoning (Figure 10).

The U–Pb isotopic signatures determined by SHRIMP in zircons from both samples are shown in Table 5. A total of 35 crystals were analysed (25 from sample 103324 and 10 from sample 113106), a few of them in several places (core and rim), resulting in 39 analyses. The isotopic data were plotted in Wetherill Concordia diagrams, and the calculated ages were represented in density probability plots (Figure 11). Only data with discordance below 10% and common Pb values lower than 1.1 were plotted and considered in the discussion.

The ²⁰⁶Pb/²³⁸U ages obtained in zircons from the gabbro (103324) are characterised by an abundant group restricted to a narrow range (272–301 Ma), some scattering from 358 to 506 Ma, and a single younger age of 209 Ma, which is the only age close to the age of the MPGD intrusion (~206–200 Ma) [1,10,24]. Thus, the bulk of zircons have to be considered of xenocrystic origin. The oldest ages (358–506 Ma) were obtained from the larger irregular zircon grains and stubby prisms with oscillatory zoning. The most abundant ages clearly record the Variscan cycle, yielding a peak at ~280 Ma in the probability density plot (Figure 11B), and are restricted to the overgrowths that surround the cores or to rounded homogenous grains (Figure 10). No ages were determined in the inherited cores because the initial objective was the U–Pb dating of the MPGD mafic magmatism, and age determination in cores was left aside. However, it is expected that they would have provided very old ages (see discussion). Based on the similitude in size, shape and internal structure of the analysed zircons with zircons from central Spain lower crustal granulites [44,46], the density plot of Figure 11B shows a comparison with results from these SCS granulites. It can be observed that the main U–Pb age peak produced by the gabbro is very close to that given by the granulitic zircons (~286 Ma) and that the oldest ages in the gabbro sample somewhat reproduce the pre-Variscan Palaeozoic age range of the granulites.

The scarce ²⁰⁶Pb/²³⁸U ages of zircons from the granophyre (113106) are scattered in a range from 298 to 689 Ma (Figure 11C). Thus, all zircons from this sample have to be considered xenocrystic. The distribution of ages yields several peaks of similar height (Figure 11D), but they have no statistical relevance given the low number of analyses considered. However, it is interesting to note that all pre-Variscan ages (>315 Ma) were obtained from inherited cores, whereas the Variscan ages (290–315 Ma) come from elongated zircons without inheritances or rims around inherited cores (Figure 10). Given that Variscan magmatic rocks are likely related to the origin of the analysed zircons, we also show in Figure 11D the distribution of zircon ages from the central Spain granites (see references in the caption of Figure 11). The broad similitude between the main peak of granite intrusion for central Spain (~299 Ma) and that provided by the youngest zircons from the analysed granophyre (~309 Ma) are noteworthy. The distribution of the pre-Variscan zircon ages yields several peaks in the same interval from the Ordovician to the Neoproterozoic if compared with data from the SCS country-rock granites (Figure 11D). Zircons with similar petrographic features and U–Pb age distribution to those of the studied granophyres are also present (although scarce) in rocks from the lower crustal levels of central Spain [47], and they are more abundant in lower crustal rocks further south of the MPGD outcropping section (i.e., Calatrava volcanic field) [44].

5. Discussion

5.1. Mafic Magma Evolution

The subalkaline composition of gabbroic and diabasic samples from the MPGD and the trend of increasing Fe in the AFM diagram (Figure 5A,B) are features typical of tholeiitic suites. The TiO₂ value under 1.5 wt% in the most primitive samples relates the composition of the studied mafic rocks with that of the low-Ti continental flood basalts (CFB) (e.g., [48]).

The major element geochemistry of MPGD gabbros describes a moderate to strong compositional heterogeneity, characterised by a positive correlation between Mg and Ca, whereas Ti, Na and K display a moderate enrichment towards lower Mg contents (Figure 6). Such behaviour is coherent with some degree of crystal fractionation. Al and Fe depict a different evolution pattern, defined by an initial Al enrichment at constant Fe, followed by an Fe increase at constant Al (starting at MgO ~6 wt%; Figure 6). These latter features, together with the overall early depletion of Ca and Mg, likely imply that fractionation was initially dominated by Ca-rich clinopyroxene (augite). The low Cr and Ni concentrations of the most primitive samples (Cr < 300 ppm, Ni < 140 ppm; Figure 7) imply that they cannot be considered primary magmas (Cr = 500–1000 ppm; Ni = 200–500 ppm) and that some degree of early olivine fractionation also took place. On the contrary, plagioclase likely acquired a more relevant role afterwards. The chemical variation in trace elements agrees with this possibility, as all incompatible elements increase in composition in the more evolved samples, with the exception of Sr, which reproduces the same pattern of Al (Figure 7) and supports that the change in behaviour is controlled by an increase in the modal abundance of plagioclase during fractionation. The lack of a final depletion in Fe and Ti towards the most evolved gabbroic samples indicates that the degree of differentiation did not achieve fractionation of Fe-Ti oxides, typical of the tholeiitic evolution.

An equivalent fractionated mineral assemblage has been proposed for several mafic magmas of the CAMP (including samples from the Messejana–Plasencia dyke) based on chemical modelling using MELTS software [10]. Their results indicate that plagioclase is an early crystallising phase, but their data also present an evolution pattern with increasing Al towards lower Mg, which seems to contradict the typical evolution of tholeiitic magmas. The stabilisation of plagioclase in crystallising tholeiitic melts at low pressures prior to augite has been described in experimental studies (e.g., Refs. [49,50]). However, the plagioclase/pyroxene ratio, which can change with small variations in parental melt composition, has a strong effect on the magma chemical evolution [49] and can lead to kink patterns with an initial increase in the Al content followed by depletion [50].

5.2. Origin of the Granophyric Melt

Two broad possible scenarios can be considered in the origin of these granophyres: they might represent felsic melts unrelated to the MPGD or they may be genetically linked with the mafic magma via extensive differentiation, as shown by experimental works focused on the formation of silicic melt due to crystal fractionation of basalt (e.g., [51]). The following two main reasons argue against the possibility of MPGD granophyres being related to granitic magmas of crustal origin totally disconnected from the mafic melt: first, felsic magmatic events in central Spain ended about 85 Ma before the intrusion of the MPGD (e.g., [45]), and second, the granophyres composition is not in accordance with that of minimum temperature anatectic granitic melts (Figure 12). The granophyres are found exclusively as intrusive bodies within the gabbros at several scales including the following: small interstitial pockets within the gabbro minerals, small veins of irregular contact permeating the mafic rocks and, occasionally, as a large dyke of 1 km in length (Figure 2). In all cases, these felsic bodies are only found within the mafic dyke, not in the wall-rock granites. This layout of the granophyric melts points to a direct connection between the felsic pockets within the gabbros and the small veinlets and dykes, which might be considered segregates formed due to the squeezing and extraction of a residual felsic melt during the last stages of crystallisation. The presence of similar small-volume granophyres within tholeiitic gabbros has already been described (e.g., [52,53,54]), including the CAMP magmas [11,16], and explained as the result of mafic magma evolution. Accordingly, it is likely that the granophyres formed during the last stages of the mafic magma differentiation within the MPGD intrusion.

The major and trace element geochemistry of the felsic rocks associated with the gabbroic intrusion is heterogeneous and clearly separated from that of the mafic rocks (Figure 6 and Figure 7), independently of the size of these bodies (veinlets, small dykes or larger dykes). This is manifested by a compositional gap that is noticeable for SiO₂, Fe₂O₃, CaO, Rb, Ba and Sr, among other elements (Figure 6 and Figure 7). The absence of volcanic and plutonic rocks of intermediate composition, identified as the Daly gap, is not uncommon in highly differentiated mafic magmas, but their origin is controversial (e.g., [57]). Some studies suggest that melts of intermediate composition may form during crystal fractionation of basaltic magmas and that their absence or scarcity would be the result of viscosity–density barriers [58], fast evolution of intermediate liquids [59] or other mechanical factors related to magma chamber dynamics (e.g., low crystallinity percentage and magma convection) [60]. However, the presence of interstitial quartz–alkali feldspar-bearing micro-aggregates within the gabbro mineral groundmass (Figure 2C–F) implies that the differentiation process leading to granophyre formation can occur in the dyke itself after its intrusion, rather than within a deeper magma chamber. The fact that the micro-aggregates, veinlets and dykes represent felsic compositions, with no intermediate terms, suggests that simple crystal fractionation would not be the main differentiation mechanism in the generation of these granophyres.

An alternative hypothesis proposed to explain the Daly compositional gap is liquid immiscibility (e.g., [57]). This process has been associated with the presence of glass droplets in the mesostasis of basalts (e.g., [17]), but it has also been identified in mafic plutonic intrusions, such as Skaergaard, Sept Iles, Bushveld and Duluth (e.g., [61,62,63,64,65,66]). These studies describe the coeval segregation of silica and iron-rich melts at the mm- to cm-scale.

Immiscibility along tholeiitic trends involves the unmixing of ferrobasaltic and rhyolitic melts (e.g., [66]). The compositional space of such an immiscibility gap has been represented in triangular diagrams such as those of Figure 13. In all these plots, the granophyres from the MPGD display a composition very close to that of the silica-rich immiscible liquids, represented either by glass droplets in basalts (Figure 13A) or melts from immiscibility experiments in a tholeiitic system (Figure 13B,C).

The bulk of the gabbroic samples plot in an intermediate position between the Si-rich and Fe-rich poles, within the wide one-liquid field represented by lavas from some major tholeiitic provinces (Figure 13B,C). The mineral chemistry of the granophyres is slightly different from that of the gabbros (e.g., higher Na-Al in biotites, higher Na + K in apatite; Figure 4), but mostly overlaps with the mineral geochemical evolution trends (Figure 4). Experimental and empirical studies ([66] and references therein) have established that the onset of immiscibility in tholeiitic mafic systems is associated with plagioclase with a composition in the range of An_36–60 and clinopyroxenes with Mg# < 64. Such values are reached in the studied MPGD gabbros (Figure 4A), which agree with the evolved contents usually associated with liquid immiscibility. The evolution depicted by the MPGD mafic rocks is also in accordance with the liquid line of descent described by tholeiitic suites, which have experienced immiscibility, such as the slight increase in Ti, Fe and P in the gabbros (Figure 6). The relative depletion in Eu, Sr and Ti in the granophyres (Figure 8) is coherent with their generation during the last stages of differentiation of the mafic magma, after fractionation of plagioclase and Fe–Ti oxides. Moreover, other chemical features of the granophyres, such as the relative enrichment in Rb and elevated Ba/La (20–94) and Ta/Nb (0.12–0.31) ratios (Figure 8), can be considered good indicators of liquid immiscibility in silicic rocks. Liquid immiscibility in tholeiitic systems can partition trace elements, e.g., favouring the incorporation of LILE to the felsic melt over REE, thus imposing these distinct geochemical imprints in the two immiscible melts [68]. Such partitioning could explain the similar or lower REE contents of the granophyres when compared to the gabbros.

A possible weakness of this magma immiscibility model is the low ratio of iron-rich melts to silica-rich melts, the latter being represented by the minor granophyric bodies within the MPGD. Only one sample (103329) shows a particularly high Fe–Mn–Ti-rich and Si–Al–Na-poor composition (Figure 6), which plots towards the Fe-rich pole of the experimental immiscibility fields (Figure 13B,C). This sample might be representative of cumulate rocks derived from this Fe-rich immiscible counterpart and may have been undervalued during sampling due to their similar appearance to any other MPGD gabbro. In any case, the experiments suggest that the relative amount of Fe-rich melt produced during magma immiscibility is low [66]. In summary, these data, together with the lack of intermediate compositions and the close textural relationship between gabbros and granophyres, support the likely origin of the MPGD granophyres as immiscible Si-rich melts generated during the last stages of differentiation of the mafic magma.

5.3. Magma Ascent and Contamination

The slightly evolved composition of the MPGD gabbros and their isotopic heterogeneity points to an open-system differentiation process, although it is not straightforward at which lithospheric level it took place. As stated above, the low Mg and Cr values of these gabbros are not characteristic of primary magmas and point to a certain degree of olivine–clinopyroxene–plagioclase fractionation. The fact that both plagioclase and clinopyroxene (instead of only plagioclase) are relevant early-crystallising minerals suggests that fractionation took place at moderate pressures. Fractional crystallisation experiments with tholeiitic melts at 0.7 GPa [50], equivalent to the lower crust (~23–25 km), yielded trends similar to the MPGD gabbros (e.g., initial Al₂O₃ increment).

The entrainment of zircon xenocrysts are also relevant data. The presence of such zircons in mafic magmas is usually interpreted as the result of crustal assimilation. The gabbro sample 103324 includes relatively abundant zircon xenocrysts with an internal CL structure similar to that of metamorphic zircons of the SCS granulitic xenoliths (Figure 10): small, rounded grains with a homogeneous dark CL response or dark mantles surrounding older inherited zones [44,46]. Moreover, the distribution of U–Pb zircon ages both in gabbro xenocrysts and zircons from granulites overlaps and displays a remarkable peak at about 280–290 Ma (Figure 11). The study of these xenoliths [44,46,69], which are included in Permian mafic magmas, has shown that (1) they were extracted from the SCS lower crust (~25–30 km), (2) this deep crustal level is abundant in felsic granulites and (3) the above ages represent the last stages of granulitisation in the SCS lower crust at the end of the Variscan cycle. These textural features and geochronology data of zircon xenocrysts suggest that the mafic magma likely interacted with the lower crust. Another possibility is that these zircons were incorporated from the mantle source. Examples of mantle zircons have been described in orogenic peridotites and mantle xenoliths, either associated with mantle metasomatism (e.g., [70]) or recycling of crustal components (e.g., [71,72]). Such recycling is possible via lower crustal delamination into the mantle, which has been proposed for the central zone of the Iberian Variscan Belt [73]. Unfortunately, we do not have solid arguments to support this possibility.

Accordingly, we think that a process of mafic magma stagnation, fractionation and contamination with lower crustal granulites is likely behind the evolution of the MPGD magmas. Such a process has already been proposed in previous studies [7,10], although the conclusions differ with respect to the nature of the assimilating component. The possibility of crustal contamination could also account for the variability in the Sr–Nd radiogenic ratios (Figure 9). We tested this possibility with a crystal fractionation and assimilation (AFC) model (Figure 14) based on the equations of [74]. The model (see details on the parameters used in Table 6) illustrates the isotopic and trace element variation in an initial mafic magma with a composition similar to the most depleted MPGD gabbros. The contaminant used in the model is the average of SCS lower crustal granulite xenoliths (grey area) [69]. The fractionating phases are cpx, pl and ol, as interpreted at the beginning of the discussion. The model fits well with the chemical variation in the gabbros when using an assimilation/fractionation ratio (r) of 0.25, which implies assimilation of lower crustal rocks in the approximate range of 10–12%. A model considering the SCS granites as the contaminating component does not fit the isotopic variation in the gabbros due to their high radiogenic Sr values. Taking into account that the analysed rocks do not represent a primary composition, these crustal contamination rates could be considered minimum values.

The more radiogenic isotope composition of the granophyres with respect to that of the MPGD gabbros (Figure 9), along with the inclusion of typically igneous Variscan zircons, points also to crustal contamination, although involving a different crustal protolith. These zircon xenocrysts display CL textures similar to those of Variscan granites or pre-Variscan meta-igneous rocks: large euhedral bipyramidal prisms with oscillatory zoning (Figure 10). Their U–Pb ages also reproduce the distribution observed in zircons from Variscan granites (Figure 11D), with a large peak at ~300 Ma (crystallisation age) and older (usually Cadomian) inheritances. Such features are characteristic of the SCS granites [45] and some lower crustal granulites [44,47] and likely reflect a contamination process occurring at a different crustal level than that described above. As the MPGD granophyres are formed and segregated within the gabbros, the contamination process has to occur initially in the mafic magma. We applied an AFC model for Sr–Nd isotopes (Figure 15) starting with the composition of the least primitive gabbros (see details on the parameters used in Table 6) because the granophyres were derived from the more differentiated mafic magmas. The fractionating minerals used in the model are the most abundant phases in the gabbros (cpx, pl and opx), and the contaminant is a representative SCS granite taken from [40]. An assimilation/fractionation ratio (r) of 0.25 is used. The resulting model only fits well for the least radiogenic granophyres, yielding assimilation rates up to 21%. The evolution line does not reach the more radiogenic samples from NE Sanchorreja, which would require a much more radiogenic contaminant and/or unrealistic assimilation rates.

Another process, apart from crustal assimilation, is necessary to explain the composition of the two granophyres with ⁸⁷Sr/⁸⁶Sr values higher than 0.714. This increase in radiogenic Sr could be associated with the alteration of these rocks, caused by their interaction with meteoric or hydrothermal fluids and manifested in the variable transformation of the main igneous minerals to secondary phases. The principle behind it is based on the washing of radiogenic ⁸⁷Sr from K-feldspar lattice sites due to an extensive and prolonged rock–fluid interaction, which could lead to preferential mobilisation of Sr vs. Nd. A link between hydrothermal fluid circulation and the disturbance of Rb–Sr isotope systematics has already been proposed for granitic rocks [75,76]. Other studies focused on the effects of weathering or hydrothermal alteration on the Rb–Sr isotopic system [77,78,79] also suggest that an increase in the alteration degree can lead to higher ⁸⁷Sr/⁸⁶Sr ratios. The presence of miarolitic cavities filled with low-T minerals (prehnite, chlorite, clay minerals) in these samples from the NE Sanchorreja granophyre dyke, is indicative of the involvement of late-stage fluids. Also, the region where the samples were collected has been subject to several hydrothermal events after the MPGD intrusion [80]. Anyhow, the nature of this fluid and the way it may have interacted with the MPGD felsic rocks is not straightforward and is beyond the scope of this work. However, it should be noted that the higher K₂O contents and general depletion in HREE and most HFSE of the samples from NE Sanchorreja, when compared to other granophyre dykes, are compatible with potassic alteration [81].

Table 6

Parameters used in the AFC models for the Messejana–Plasencia Great Dyke gabbros and granophyres.

AFC model for the MPGD gabbros
	87Sr/86Sr	143Nd/144Nd	Sr	Nd	Rb	Ba
Initial melt composition	0.70560	0.512443	150	10.5	14.5	122
Contaminant composition	0.71216	0.511983	248	32	104	1037
Fractionating phases: Cpx (67%), Pl (27%) and Ol (6%)
AFC model for the MPGD granophyres
	87Sr/86Sr	143Nd/144Nd	Sr	Nd
Initial melt composition	0.70636	0.512271	214	20
Contaminant composition	0.71915	0.512016	96	32
Fractionating phases: Cpx (50%), Pl (40%) and Opx (10%)
Mineral/mafic magma partition coefficients
Mineral	Element	Kd	Reference
Clinopyroxene	Sr	0.1	[82]
Clinopyroxene	Nd	0.38	[82]
Clinopyroxene	Rb	0.015	[83]
Clinopyroxene	Ba	0.005	[82]
Plagioclase	Sr	2	[84]
Plagioclase	Nd	0.3	[84]
Plagioclase	Rb	0.08	[84]
Plagioclase	Ba	0.2	[84]
Olivine	Sr	0.0094	[83]
Olivine	Nd	0.008	[85]
Olivine	Rb	0.0133	[83]
Olivine	Ba	0.01	[83]
Orthopyroxene	Sr	0.0026	[82]
Orthopyroxene	Nd	0.013	[82]

The bulk of the analysed granophyres constitute a relatively heterogeneous group, as highlighted by the variable concentrations of some elements (e.g., Na, K, Rb, Ba; Figure 6 and Figure 7). This variability, which is also reflected in the mineral chemistry (Figure 4), suggests that granophyres from each sector followed independent differentiation processes as disconnected magma bodies after their segregation from the mafic intrusion. The larger size of the NE Sanchorreja dyke, along with its highly radiogenic Sr ratios and the entrainment of zircon xenocrysts, calls for a more complex differentiation scenario (a larger magma body). Overall, the above data point to the joint participation of crystal fractionation, crustal rocks assimilation and posterior fluid interaction.

6. Conclusions

The Messejana–Plasencia great dyke (MPGD) is a complex Late Triassic tholeiitic gabbro intrusion which, with a minimum confirmed length of ~630 km, can be considered one of the world’s largest dykes known.

The major and trace element composition of the mafic rocks does not represent that of primary mantle-derived melts. Nevertheless, primary tholeiitic (low-Ti) magmas would ascend from their lithospheric mantle source and stagnate close to the Moho, where they would differentiate via crystal fractionation of olivine, clinopyroxene and plagioclase. Their Sr–Nd radiogenic ratios, close to the BSE values, display a significant heterogeneity and correlation with the SiO₂ content. These geochemical features, together with the inclusion of xenocrystic zircons with similar texture and age to those of SCS lower crustal granulites, indicate that these magmas assimilated rocks from the deep crust. These evolved mafic magmas would rise and intrude the upper crust in successive pulses.

The granophyric bodies found in the NE segment of this gabbroic intrusion represent highly evolved magmas with a composition clearly separated from that of the gabbros by a compositional gap (lack of intermediate values). Several features support their formation via liquid immiscibility including the following: (1) the presence in the same outcrop of granophyre melts at very different scales (interstitial mm-sized pockets within the gabbros, small veins and m-sized dykes), (2) the compositional gap with the host gabbros and (3) the similitude of their major element composition with that of experimental immiscible felsic melts formed from mafic tholeiitic magmas.

The more radiogenic Sr–Nd ratios of the MPGD granophyres with respect to the gabbros, and the presence of xenocrystic zircons with internal textures and ages similar to those of the SCS granites, suggests that the ascending magmas were also contaminated by crustal rocks equivalent to these latter intrusions. These felsic melts would segregate and evolve in independent magma bodies of different volumes. Hydrothermal alteration can explain the markedly high Sr radiogenic values of samples from the larger Sanchorreja granophyre dyke.

Author Contributions

Petrography and major element mineral chemistry, C.V. and P.A.; geochemical models, E.L. and D.O. All authors collaborated in the creation of figures and tables. D.O. conceived this study and wrote an original draft, which was later improved with inputs from C.V., E.L. and P.A. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

The original contributions presented in this study are included in this article. Further inquiries can be directed to the corresponding author.

Acknowledgments

We acknowledge Alfredo Fernández Larios for his assistance with the electron microprobe analyses in the CAI of Microscopía Electrónica (UCM). We also thank Pilar Montero for the work carried out during zircon preparation and analysis in the IBERSIMS laboratory (University of Granada). This work has also benefited from the constructive comments of four anonymous reviewers, which have helped to improve the final manuscript.

Conflicts of Interest

The authors declare that they have no personal circumstances or interests that may be perceived as inappropriately influencing the representation or interpretation of the reported research results. The funders had no role in the design of this study; in the collection, analyses, or interpretation of data; in the writing of this manuscript; or in the decision to publish the results.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables

Figure 2. (A) Example of a field outcrop of a granophyric dyke intruding the gabbroic Messejana–Plasencia dyke (hammer for scale) (stop 1 of Figure 1). (B) Granophyric cm-sized dykelete intruding diabasic gabbro in the Villatoro pass (stop 6 of Figure 1). See a one-euro coin in the left corner for scale. (C–F) Plane- and cross-polarised light images of thin sections illustrating interstitial quartz-alkali feldspar granophyric intergrowth within the diabasic gabbro: (C,D): sample 103324 50× PPL and XPL; (E,F): sample 103301 50× PPl and XPL, respectively. The white bar in images (C–F) represents 550 μm. Mineral abbreviations after [34].

Figure 3. Cross-polarised light photographs (A,B) and Back-Scattered Electron (BSE) images (C–F) of selected granophyres of the Messejana–Plasencia dyke. (A,B) Thin section of granophyre 113106 (50× PPL) showing granophyric intergrowths of quartz and alkali feldspar within a medium-grained matrix of alkali feldspar and quartz crystals. (C) Gabbro 103324 showing interstitial granitic granophyre intergrowths of quartz and alkali feldspar. (D,E) Granophyre 113106 showing quartz-K-feldspar intergrowths and accessory biotite and apatite crystals. (F) Granophyre 113105 showing biotite partially altered and apatite crystals. The white bar represents 550 μm and 200 μm in images A–B and C–F, respectively. Mineral abbreviations after [34].

Figure 4. Selected major oxide mineral chemistry of MPGD gabbros and granophyres. (A) Feldspar composition in Ab-Or-An plots. (B–E) Biotite composition of selected oxides (wt%) and AlIV (apfu) vs. Mg# or Fe# (in atomic ratios). FeO-MgO-Al2O3 plot with compositional fields after [35]. Fields I, II and III are aluminium-potassic (peraluminous), calc-alkaline and monzonitic granitic series, respectively. (F) Apatite composition in a Na2O + K2O vs. F/Cl plot (in wt%). In some diagrams, granophyre fields from different locations (see key to Figure 1B) are marked.

Figure 5. Whole-rock major element composition (wt%) of MPGD rocks: (A) TAS diagram: total alkalis vs. silica. Classification fields taken from [37] and references therein. (B) Triangular (Na2O + K2O)–FeOt–MgO diagram for tholeiitic–calc-alkaline series discrimination. (C) K2O vs. SiO2 diagram. (D) A/NK vs. A/CNK diagram. Literature data from [7,8].

Figure 6. Major elements vs. MgO (wt%) diagrams for the MPGD samples. Literature data from [7,8]. The composition of SCS granites is also plotted for comparison (grey squares) [38,39,40].

Figure 7. Trace elements (ppm) vs. MgO (wt%) diagrams for the MPGD samples. Literature data from [7,8]. The composition of SCS granites is also plotted for comparison (grey squares) [38,39,40].

Figure 8. (A) Chondrite-normalised REE and (B) Primitive Mantle-normalised multi-trace element composition of the MPGD samples. Literature data from [7,8]. Chondrite and Primitive Mantle normalisation values after [41,42], respectively.

Figure 9. Initial 87Sr/86Sr and 143Nd/144Nd ratios of the MPGD samples. Literature data from [7,8]. The fields of SCS Variscan mafic magmatism, SCS orthogneisses and SCS Variscan felsic magmatism were taken from [21,43,40], respectively.

Figure 10. CL images of representative zircon grains from samples 103324 (gabbro) and 113106 (granophyre) from the MPGD, showing the analysed spots (red ellipse) and age obtained.

Figure 11. SHRIMP U–Pb data for zircons from the MPGD samples. (A,B) Concordia and probability plot diagrams of sample 103324 (gabbro). (C,D) Concordia and probability plot diagrams of sample 113106 (granophyre). Error ellipses are given at the 2σ level. The distribution of zircon ages from the SCS lower crust granulites and SCS Variscan granites, which were taken from [44,45], respectively, is plotted for comparison in diagrams (B,D).

Figure 12. Normative composition of the MPGD granophyres in the Ab–Or–Qtz system. The ternary cotectic curves and eutectic minima for 0.2, 0.5 and 1 GPa [55], and the composition of typical anatectic granites are also shown. Modified after [56].

Figure 13. Triangular diagrams showing the major element composition of the MPGD samples together with immiscibility fields from experiments in tholeiitic systems. (A): Plot modified after [17], including the composition of Fe-rich and Si-rich glasses formed by liquid immiscibility. The immiscibility field of fayalite–leucite–silica was taken from [67]. (B,C): Plots modified after [66] showing tholeiitic basalts and their products of differentiation and the Fe-rich and Si-rich experimental and natural immiscible melts. The thick grey line separates the two-liquid (up) from the one-liquid (down) fields. Literature data from [7,8].

Figure 14. Assimilation and Fractional Crystallisation (AFC) models [74] with Sr–Nd isotopic ratios and trace element composition of the MPGD gabbros. Literature data from [7,8]. The black arrow represents the chemical evolution predicted by the model with an indication of the percentage of fractionation. The best fit is obtained for an assimilation/fractionation ratio (r) of 0.25. The initial melt and contaminant compositions are represented with a white circle and a black star, respectively. Note the change in scale in the isotopic diagram (left).

Figure 15. Assimilation and Fractional Crystallisation (AFC) model [74] for the Sr–Nd isotopic ratios of the MPGD granophyres. Literature data from [7,8]. The black line represents the evolution predicted by the model with an indication of the percentage of fractionation. The initial melt and contaminant compositions are represented with a white circle and a black star, respectively. The compositional field of SCS granites is also shown for comparison (grey field).

Table 1

Major element chemistry (wt%) of representative feldspars from the gabbros and granophyres of the Messejana–Plasencia dyke.

Sample	103324	103324	103324	103324	103317	103317	103317	103327	103327	103327	103327	103309	104367
Analysis	48	37	39	43	28	42	45	33	35	2	19	74	81
Rock	Gb	Gb	Gb	Gb	Gb	Gb	Gb	Gb peg	Gb peg	Gb peg	Gb peg	Gb	Grn
Mineral	Kfs	Kfs	Pl	Pl	Pl	Pl	Pl	Pl	Pl	Pl	Kfs	Pl	Pl
Location	interst	interst	rim	core			interst	interst			interst
Outcrop	7	7	7	7	6	6	6	7	7	7	7	1	6
SiO2	64.99	62.94	59.96	52.63	49.46	50.19	66.27	66.83	55.51	61.61	63.92	50.52	67.52
TiO2	0.01	0.08	0.02	0.11	0.01	0.04	0.00	0.00	0.06	0.04	0.08	0.03	0.03
Al2O3	18.73	18.15	23.44	28.75	30.18	29.56	20.19	19.67	27.29	23.72	18.97	30.90	19.77
FeOt	0.08	0.03	0.19	0.47	0.74	0.78	0.12	0.00	0.71	0.35	0.13	0.65	0.04
MnO	0.04	0.00	0.00	0.00	0.00	0.06	0.03	0.00	0.00	0.00	0.00	0.00	0.01
MgO	0.02	0.00	0.00	0.04	0.13	0.20	0.06	0.00	0.02	0.01	0.00	0.02	0.00
CaO	0.24	0.05	5.74	11.53	16.74	14.92	1.48	0.54	10.43	5.65	0.36	14.73	0.09
Na2O	3.25	0.68	8.32	4.85	3.00	3.83	11.01	11.42	5.75	8.00	2.54	3.08	12.02
K2O	12.20	15.95	0.77	0.31	0.12	0.22	0.42	0.04	0.43	0.82	12.56	0.16	0.03
Total	99.57	97.88	98.45	98.73	100.56	100.16	99.63	98.50	100.27	100.25	98.57	100.15	99.75
Cations calculated based on 8 O
Si	2.980	2.979	2.721	2.420	2.273	2.312	2.929	2.968	2.509	2.740	2.953	2.307	2.968
Ti	0.000	0.003	0.001	0.004	0.000	0.001	0.000	0.000	0.002	0.001	0.003	0.001	0.001
Al	1.013	1.013	1.254	1.559	1.635	1.605	1.052	1.030	1.454	1.244	1.050	1.664	1.024
Fe3+	0.003	0.001	0.006	0.016	0.026	0.027	0.004	0.000	0.024	0.012	0.004	0.022	0.001
Mg	0.001	0.000	0.000	0.002	0.009	0.014	0.004	0.000	0.002	0.001	0.000	0.001	0.000
Mn	0.002	0.000	0.000	0.000	0.000	0.002	0.001	0.000	0.000	0.000	0.000	0.000	0.001
Ca	0.012	0.003	0.279	0.568	0.824	0.736	0.070	0.026	0.505	0.269	0.018	0.721	0.004
Na	0.289	0.063	0.732	0.432	0.267	0.342	0.943	0.983	0.504	0.690	0.231	0.273	1.024
K	0.714	0.963	0.044	0.018	0.007	0.013	0.024	0.002	0.025	0.047	0.752	0.009	0.002
Suma cat.	5.013	5.024	5.038	5.020	5.041	5.053	5.027	5.009	5.024	5.003	5.010	4.998	5.025
An	1.18	0.27	26.44	55.78	74.90	67.12	6.77	2.52	48.85	26.79	1.80	71.88	0.41
Ab	28.48	6.10	69.36	42.46	24.26	31.22	90.86	97.26	48.75	68.59	23.09	27.22	99.43
Or	70.34	93.64	4.20	1.77	0.66	1.19	2.29	0.22	2.39	4.63	75.11	0.90	0.16
Sample	103311	103311	103326	103326	106284	106284	106292	106292	113105	113105	113106	113106
Analysis	32	60	11	13	95	97	77	81	30	31	50	57
Rock	Grn	Grn	Grn	Grn	Grn	Grn	Grn	Grn	Grn	Grn	Grn	Grn
Mineral	Pl	Kfs	Pl	Kfs	Kfs	Kfs	Kfs	Kfs	Kfs	Kfs	Kfs	Kfs
Location										simpl
Outcrop	1	1	7	7	3	3	6	6	2	2	2	2
SiO2	66.75	64.58	68.21	64.01	63.26	63.15	63.18	63.35	63.97	65.33	63.24	63.51
TiO2	0.06	0.02	0.03	0.04	0.04	0.12	0.03	0.00	0.02	0.04	0.00	0.04
Al2O3	20.14	18.39	19.50	18.35	18.25	18.13	18.03	18.16	18.11	18.32	20.03	19.09
FeOt	0.08	0.26	0.03	0.02	0.09	0.81	0.22	0.10	0.00	0.03	0.05	0.03
MnO	0.01	0.02	0.04	0.00	0.00	0.00	0.01	0.05	0.00	0.03	0.00	0.00
MgO	0.01	0.00	0.00	0.01	0.04	0.17	0.01	0.00	0.02	0.01	0.00	0.03
CaO	0.77	0.08	0.34	0.05	0.01	0.01	0.22	0.02	0.18	0.08	0.40	0.17
Na2O	11.17	0.42	11.62	0.23	0.14	0.12	0.47	0.27	2.30	4.15	4.88	2.44
K2O	0.04	15.53	0.04	16.49	16.81	16.68	15.85	16.65	13.62	10.77	9.76	13.29
Total	99.03	99.33	99.81	99.41	98.65	99.26	98.04	98.61	98.21	98.75	98.37	98.59
Cations calculated based on 8 O
Si	2.951	2.997	2.987	2.986	2.979	2.968	2.985	2.984	2.990	3.001	2.918	2.955
Ti	0.002	0.001	0.001	0.001	0.001	0.004	0.001	0.000	0.001	0.001	0.000	0.001
Al	1.049	1.006	1.007	1.010	1.013	1.004	1.004	1.008	0.998	0.992	1.090	1.047
Fe3+	0.003	0.009	0.001	0.001	0.003	0.029	0.008	0.004	0.000	0.001	0.002	0.001
Mg	0.001	0.000	0.000	0.001	0.003	0.012	0.000	0.000	0.001	0.000	0.000	0.002
Mn	0.000	0.001	0.002	0.000	0.000	0.000	0.000	0.002	0.000	0.001	0.000	0.000
Ca	0.036	0.004	0.016	0.003	0.000	0.000	0.011	0.001	0.009	0.004	0.020	0.009
Na	0.957	0.038	0.987	0.020	0.013	0.011	0.043	0.024	0.208	0.370	0.437	0.220
K	0.002	0.919	0.002	0.982	1.010	1.000	0.956	1.000	0.812	0.631	0.575	0.789
Suma cat.	5.001	4.976	5.002	5.004	5.023	5.028	5.009	5.023	5.019	5.001	5.042	5.024
An	3.66	0.42	1.61	0.27	0.04	0.02	1.09	0.10	0.90	0.38	1.90	0.84
Ab	96.13	3.97	98.16	2.02	1.22	1.08	4.23	2.39	20.22	36.82	42.36	21.67
Or	0.20	95.61	0.24	97.36	98.73	98.89	94.68	97.51	78.89	62.81	55.73	77.50

Mineral abbreviations after [34]. Gb = gabbro, Gb peg = pegmatitic gabbro, Grn = felsic granophyre, interst = interstitial granophyric intergrowth. Outcrop numbers from Figure 1.

Table 2

Major element chemistry (wt%) of representative biotite and apatite from the gabbros and granophyres of the Messejana–Plasencia dyke.

Sample	103322	103322	103325	103325	103325	103327	103327	103301	103309	103309	103311	103311	113105
Analysis	112	121	13	16	19	50	53	62	1	3	28	34	34
Rock	Gb	Gb	Gb	Gb	Gb	Gb peg	Gb peg	Gb	Gb	Gb	Grn	Grn	Grn
Mineral	Bt	Bt	Bt	Bt	Bt	Bt	Bt	Bt	Bt	Bt	Bt	Bt	Bt
Outcrop	7	7	7	7	7	7	7	4	1	1	1	1	2
SiO2	36.40	37.57	36.02	34.66	35.80	36.72	36.25	34.99	36.01	36.20	34.12	34.91	35.40
TiO2	4.81	3.90	5.08	5.51	4.12	3.74	5.14	3.93	4.89	4.45	1.95	1.65	7.49
Al2O3	12.97	13.04	13.03	13.07	13.66	12.60	13.05	13.05	12.51	12.31	16.83	17.19	14.56
FeOt	25.97	23.03	27.48	28.06	26.58	22.26	26.14	25.76	25.55	23.41	25.64	25.67	16.15
MnO	0.17	0.07	0.07	0.15	0.09	0.17	0.05	0.10	0.04	0.13	0.25	0.16	0.19
MgO	7.57	10.10	5.10	4.47	7.01	9.97	7.16	9.18	6.43	8.78	6.02	5.85	11.89
CaO	0.01	0.04	0.11	0.07	0.04	0.00	0.02	0.01	0.00	0.00	0.07	0.10	0.01
Na2O	0.24	0.21	0.21	0.17	0.13	0.38	0.34	0.15	0.28	0.24	0.50	0.45	0.45
K2O	9.61	9.51	9.23	9.25	9.26	9.47	8.95	8.03	9.08	9.20	8.18	8.08	9.11
P2O5	nd.	nd.	nd.	nd.	nd.	nd.	nd.	nd.	nd.	nd.	nd.	nd.	nd.
F	0.26	0.36	nd.	nd.	nd.	nd.	nd.	0.15	0.36	0.53	nd.	nd.	0.20
Cl	nd.	nd.	nd.	nd.	nd.	nd.	nd.	0.42	0.36	0.27	nd.	nd.	0.31
Total	97.99	97.82	96.33	95.41	96.68	95.31	97.12	95.87	95.49	95.51	93.57	94.05	95.96
Mg#	0.34	0.44	0.25	0.22	0.32	0.44	0.33	0.39	0.31	0.40	0.29	0.29	0.56
Cations calculated based on 22 O (OH, Cl, F) for biotite and 25 O (OH, Cl, F) for apatite
Si	5.588	5.666	5.679	5.561	5.595	5.715	5.618	5.470	5.647	5.614	5.466	5.537	5.311
AlIV	2.346	2.317	2.321	2.439	2.405	2.285	2.382	2.404	2.313	2.251	2.534	2.463	2.576
AlVI	0.000	0.000	0.099	0.032	0.111	0.026	0.002	0.000	0.000	0.000	0.644	0.751	0.000
Ti	0.490	0.425	0.602	0.665	0.484	0.438	0.599	0.336	0.537	0.384	0.235	0.196	0.732
Fe3+	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
Fe2+	3.334	2.904	3.622	3.765	3.473	2.897	3.388	3.368	3.350	3.037	3.435	3.404	2.027
Mn	0.022	0.009	0.010	0.021	0.012	0.022	0.007	0.013	0.005	0.016	0.034	0.021	0.025
Mg	1.732	2.271	1.199	1.068	1.634	2.312	1.655	2.139	1.503	2.029	1.437	1.384	2.661
Ca	0.001	0.007	0.018	0.013	0.007	0.000	0.004	0.002	0.000	0.000	0.011	0.017	0.002
Na	0.071	0.061	0.063	0.054	0.038	0.114	0.102	0.045	0.084	0.072	0.154	0.137	0.131
K	1.882	1.829	1.857	1.893	1.845	1.880	1.770	1.601	1.815	1.820	1.672	1.636	1.744
P	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
F	0.128	0.173	0.000	0.000	0.000	0.000	0.000	0.075	0.180	0.257	0.000	0.000	0.094
Cl	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.111	0.095	0.072	0.000	0.000	0.078
Suma cat	15.594	15.661	15.470	15.511	15.605	15.689	15.526	15.563	15.529	15.552	15.623	15.546	15.380
Sample	113106	106285	106285	106286	106286	106289	106292	106292	103301	103324	113106	104367
Analysis	54	34	38	42	44	13	76	84	65	42	56	8
Rock	Grn	Grn	Grn	Grn	Grn	Grn	Grn	Grn	Gb	Gb	Grn	Grn
Mineral	Bt	Bt	Bt	Bt	Bt	Bt	Bt	Bt	Ap	Ap	Ap	Ap
Outcrop	2	2	2	2	2	2	6	6	4	7	2	6
SiO2	36.70	33.21	33.93	34.75	33.33	33.61	33.46	35.28	0.46	0.14	0.14	0.21
TiO2	5.55	3.57	3.60	2.33	2.29	3.14	4.43	3.57	0.00	0.00	0.00	0.00
Al2O3	15.88	17.74	16.62	17.63	18.27	18.51	14.39	13.89	0.00	0.00	0.00	0.00
FeOt	15.01	25.94	26.02	25.57	26.12	26.59	24.15	22.18	0.39	0.56	1.97	0.19
MnO	0.07	0.33	0.27	0.29	0.34	0.38	0.12	0.13	0.05	0.00	0.69	0.00
MgO	14.32	3.68	5.03	5.81	4.39	2.90	7.40	8.72	0.07	0.00	0.39	0.00
CaO	0.02	0.02	0.05	0.01	0.02	0.02	0.05	0.08	55.59	55.94	50.77	55.08
Na2O	0.63	0.46	0.62	0.61	0.62	0.45	0.53	0.52	0.05	0.06	0.38	0.08
K2O	8.81	9.09	8.88	8.98	8.18	9.19	8.93	8.61	0.01	0.02	0.06	0.04
P2O5	nd.	nd.	nd.	nd.	nd.	nd.	nd.	nd.	41.20	41.63	40.86	42.30
F	0.88	2.79	2.73	3.10	3.01	2.67	1.37	1.56	1.95	2.59	1.73	2.26
Cl	0.27	0.69	0.40	0.43	0.50	1.29	1.13	1.07	0.68	0.13	1.94	0.09
Total	98.22	97.51	98.20	99.76	97.12	98.76	95.96	95.61	99.52	100.06	99.03	99.43
Mg#	0.63	0.20	0.25	0.29	0.23	0.16	0.35	0.41
Cations calculated based on 22 O (OH, Cl, F) for biotite and 25 O (OH, Cl, F) for apatite
Si	5.267	4.976	5.052	5.052	4.980	4.979	5.152	5.354	0.077	0.023	0.025	0.036
AlIV	2.687	3.024	2.916	2.948	3.020	3.021	2.611	2.484	0.000	0.000	0.000	0.000
AlVI	0.000	0.108	0.000	0.072	0.196	0.210	0.000	0.000	0.000	0.000	0.000	0.000
Ti	0.554	0.402	0.371	0.255	0.257	0.349	0.277	0.246	0.000	0.000	0.000	0.000
Fe3+	0.000	0.000	0.053	0.279	0.067	0.000	0.000	0.000	0.000	0.000	0.000	0.000
Fe2+	1.801	3.250	3.186	2.829	3.195	3.294	3.110	2.815	0.054	0.077	0.284	0.026
Mn	0.009	0.042	0.034	0.035	0.043	0.048	0.016	0.016	0.008	0.000	0.100	0.000
Mg	3.064	0.821	1.117	1.260	0.978	0.641	1.699	1.974	0.017	0.000	0.100	0.000
Ca	0.004	0.003	0.009	0.002	0.003	0.003	0.009	0.013	9.904	9.899	9.376	9.939
Na	0.174	0.135	0.178	0.172	0.180	0.130	0.159	0.152	0.015	0.020	0.126	0.026
K	1.613	1.737	1.686	1.666	1.559	1.737	1.754	1.666	0.002	0.004	0.014	0.008
P	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	5.898	5.965	5.935	5.948
F	0.398	1.319	1.284	1.425	1.420	1.253	0.666	0.747	0.993	1.302	0.903	1.141
Cl	0.067	0.175	0.101	0.106	0.125	0.325	0.295	0.276	0.186	0.036	0.539	0.023
Suma cat	15.638	15.991	15.987	16.102	16.025	15.989	15.749	15.744	17.154	17.326	17.401	17.148

Mineral abbreviations after [34]. Gb = gabbro, Gb peg = pegmatitic gabbro, Grn = felsic granophyre. Outcrop numbers from Figure 1.

Table 3

Whole-rock composition (major and trace elements) of samples from the Messejana–Plasencia dyke ¹.

Sample	103310	103326	106284	106285	106292	113105	113106	113104
Type 2	Grn	Grn	Grn	Grn	Grn	Grn	Grn	Gb
Location	1	7	3	2	6	2	2	6
SiO2	72.11	74.85	73.62	70.55	71.82	73.86	72.95	51.77
TiO2	0.24	0.20	0.28	0.32	0.20	0.07	0.10	1.68
Al2O3	12.42	12.84	11.88	12.93	12.20	14.00	14.22	16.31
Fe2O3 (*)	2.90	1.31	2.62	2.70	2.88	0.68	0.63	12.47
MnO	0.03	0.02	0.01	0.01	0.03	0.01	0.01	0.15
MgO	0.80	0.78	1.78	1.14	1.60	0.33	0.42	3.41
CaO	0.71	0.64	0.46	0.31	0.56	0.39	0.38	10.29
Na2O	1.87	4.75	0.06	0.65	2.52	1.55	1.64	2.60
K2O	6.36	3.17	4.70	7.61	5.60	8.57	8.36	0.72
P2O5	0.18	0.15	0.18	0.20	0.18	0.22	0.20	0.17
LOI	2.42	1.31	3.52	2.29	1.68	0.93	1.09	1.18
Total	100.05	100.02	99.12	98.71	99.28	100.61	100.00	100.75
ACNK 3	1.6	1.0	2.0	1.3	1.1	1.1	1.1	0.7
Rb	221	105	172	278	155	252	237	24
Ba	741	217	428	370	653	438	451	200
Cs	3.90	1.20	19.4	12.1	2.58	8.20	4.50	2.10
Sr	117	61.0	32.4	28.8	42.1	49.0	45.0	208
Pb	bdl	bdl	bdl	32.9	bdl	36.0	35.0	bdl
Th	5.93	10.2	6.67	6.82	8.76	3.24	3.85	2.72
U	2.74	2.63	2.19	3.42	3.33	2.06	2.12	0.73
Zr	87.0	107	84.5	93.6	89.9	38.0	54.0	127
Nb	8.50	10.3	9.26	10.4	10.1	6.10	6.90	8.60
Y	20.4	21.5	16.8	23.0	21.4	10.8	11.6	28.8
Co	4.00	3.00	3.68	4.49	4.91	5.00	7.00	37.0
V	16.0	16.0	21.5	22.2	13.6	7.0	9.0	374
Ni	bdl	bdl	bdl	bdl	bdl	bdl	bdl	40.0
Cr	22.0	87.0	39.4	39.3	143	180	200	150
Ga	14.0	15.0	15.4	16.5	13.4	13.0	14.0	20.0
Ta	1.64	1.70	1.18	1.75	1.59	1.91	1.43	0.59
Hf	2.50	3.00	2.36	2.87	2.62	1.50	1.80	3.70
La	7.91	10.0	12.3	18.3	10.1	6.62	7.77	14.0
Ce	20.2	25.1	26.8	38.6	25.5	12.3	15.2	30.8
Pr	2.75	3.24	3.33	4.70	3.27	1.52	1.85	3.96
Nd	11.9	13.6	13.1	17.9	13.3	5.88	6.82	17.7
Sm	3.05	3.28	3.00	4.22	3.46	1.47	1.79	4.54
Eu	0.30	0.31	0.39	0.57	0.33	0.36	0.41	1.44
Gd	3.17	3.39	2.90	4.31	3.57	1.59	1.91	4.99
Tb	0.58	0.63	0.50	0.73	0.62	0.29	0.35	0.90
Dy	3.66	3.89	3.06	4.20	3.90	1.68	2.04	5.37
Ho	0.73	0.77	0.64	0.84	0.78	0.33	0.36	1.15
Er	2.12	2.26	2.04	2.49	2.31	0.95	0.97	3.32
Tm	0.32	0.34	0.30	0.36	0.34	0.14	0.14	0.48
Yb	1.93	2.06	1.87	2.19	2.06	0.91	0.88	2.94
Lu	0.28	0.29	0.27	0.32	0.28	0.14	0.14	0.46

¹ Major and trace element concentrations in wt.% and ppm, respectively. ² Grn: granophyre; Gb: gabbro. ³ ACNK = [Al₂O₃/(CaO + Na₂O + K₂O)] in molecular values. (*) Total Fe expressed as Fe₂O₃.

Table 4

Sr–Nd isotopic composition of the Messejana–Plasencia dyke granophyres and gabbro.

Sample	103310	103326	106285	113105	113106	113104
Type 1	Grn	Grn	Grn	Grn	Grn	Gb
Location	1	7	2	2	2	6
Rb (ppm)	221	105	278	252	237	24
Sr (ppm)	117	61	29	49	45	208
87Rb/86Sr	5.48	4.99	27.96	14.96	15.32	0.33
87Sr/86Sr ± 2σ	0.723865 ± 04	0.723491 ± 05	0.786676 ± 07	0.759765 ± 03	0.758034 ± 03	0.707350 ± 02
87Sr/86Sr200 Ma	0.708292	0.709300	0.707157	0.717215	0.714467	0.706401
Sm (ppm)	3.1	3.3	4.2	1.5	1.8	4.5
Nd (ppm)	11.9	13.6	17.9	5.9	6.8	17.7
147Sm/144Nd	0.155	0.146	0.143	0.151	0.159	0.155
143Nd/144Nd ±2σ	0.512395 ± 03	0.512365 ± 03	0.512269 ± 05	0.512350 ± 02	0.512360 ± 02	0.512533 ± 02
143Nd/144Nd 200 Ma	0.512192	0.512174	0.512083	0.512153	0.512152	0.512330
ε(Nd)200 Ma	–3.7	–4.0	–5.8	–4.4	–4.5	–1.0

Rb, Sr, Sm and Nd concentrations determined by ICP–MS. ¹ Grn: granophyre; Gb: gabbro.

Table 5

Ion microprobe (SHRIMP) U–Pb data of zircon from the analysed Messejana–Plasencia dyke samples.

Spot Number	Common 206Pb (%)	U (ppm)	Th (ppm)	Th/U	Radiogenic Ratios				Age (Ma)			Disc. (%) 1
					207Pb/235U ± σ (%)	206Pb/238U ± σ (%)	ρ	207Pb/206Pb ± σ (%)	207Pb/206Pb ± σ (Ma)	206Pb/238U ± σ (Ma)	207Pb/235U ± σ (Ma)
103324 (Gabbro)
A-1.1	1.4	400	114	0.3	0.6135 ± 3.2	0.0805 ± 1.8	0.40	0.0553 ± 2.7	424.5 ± 58.2	498.9 ± 8.5	485.8 ± 12.5	–2.6
A-1.2	0.1	760	256	0.3	0.6341 ± 2.0	0.0816 ± 1.7	0.62	0.0564 ± 1.0	466.5 ± 21.6	505.6 ± 8.4	498.6 ± 8.0	–1.4
A-2.1	1.2	75	107	1.5	0.2610 ± 9.8	0.0433 ± 1.7	0.13	0.0437 ± 9.6	0.0 ± 0.0	273.3 ± 4.6	235.5 ± 20.8	–16.0
A-3.1	0.4	90	100	1.1	0.3039 ± 8.6	0.0446 ± 2.0	0.17	0.0494 ± 8.3	166.1 ± 183.8	281.5 ± 5.4	269.4 ± 20.4	–4.4
A-4.1	0.1	557	323	0.6	0.5067 ± 2.2	0.0685 ± 1.9	0.61	0.0536 ± 1.1	354.7 ± 25.2	427.4 ± 7.8	416.2 ± 7.6	–2.6
B-1.1	0.3	1558	2128	1.4	0.2209 ± 3.6	0.0331 ± 2.9	0.58	0.0484 ± 2.0	119.9 ± 47.6	209.8 ± 5.9	202.7 ± 6.6	–3.6
B-2.1	1.1	129	201	1.6	0.2507 ± 9.7	0.0413 ± 3.6	0.27	0.0441 ± 8.9	0.0 ± 99.9	260.6 ± 9.3	227.1 ± 19.8	–14.8
B-3.1	0.2	1063	794	0.8	0.5957 ± 2.4	0.0702 ± 2.0	0.62	0.0616 ± 1.2	658.9 ± 24.6	437.2 ± 8.6	474.5 ± 9.1	7.8
B-4.1	0.1	2275	46	0.0	0.3978 ± 4.4	0.0572 ± 2.8	0.47	0.0504 ± 3.3	214.9 ± 75.0	358.6 ± 9.9	340.0 ± 12.7	–5.4
B-5.1	0.3	126	114	0.9	0.3373 ± 4.2	0.0437 ± 2.2	0.37	0.0560 ± 3.6	451.5 ± 77.4	275.7 ± 5.9	295.1 ± 10.8	6.6
B-6.1	0.7	120	76	0.7	0.2871 ± 4.7	0.0447 ± 2.0	0.30	0.0466 ± 4.3	0.0 ± 126.7	282.1 ± 5.5	256.3 ± 10.8	–10.0
B-7.1	0.0	230	131	0.6	0.3276 ± 2.8	0.0458 ± 1.5	0.40	0.0519 ± 2.3	281.5 ± 51.0	288.5 ± 4.4	287.7 ± 6.9	–0.2
B-8.1	0.1	315	115	0.4	0.3147 ± 4.1	0.0436 ± 2.0	0.35	0.0524 ± 3.5	302.3 ± 78.4	274.9 ± 5.4	277.8 ± 9.9	1.0
B-9.1	0.1	139	60	0.4	0.3514 ± 5.3	0.0478 ± 3.1	0.42	0.0533 ± 4.3	340.3 ± 94.8	301.2 ± 9.1	305.8 ± 14.2	1.6
B-10.1	0.1	218	71	0.3	0.3511 ± 4.7	0.0436 ± 4.0	0.61	0.0584 ± 2.5	543.7 ± 54.0	275.3 ± 10.8	305.6 ± 12.6	10.0
B-11.1	0.3	326	62	0.2	0.3053 ± 3.7	0.0444 ± 1.7	0.34	0.0498 ± 3.3	187.3 ± 74.2	280.2 ± 4.7	270.5 ± 8.9	–3.6
B-12.1	0.2	299	72	0.2	0.2959 ± 4.0	0.0438 ± 1.9	0.35	0.0490 ± 3.5	146.9 ± 80.0	276.4 ± 5.2	263.2 ± 9.4	–5.0
B-13.1	0.7	65	37	0.6	0.3300 ± 4.5	0.0462 ± 2.9	0.47	0.0518 ± 3.4	276.9 ± 75.8	291.1 ± 8.2	289.6 ± 11.4	–0.6
B-14.1	0.0	449	151	0.3	0.3168 ± 2.3	0.0455 ± 2.2	0.67	0.0505 ± 0.8	220.1 ± 17.8	286.6 ± 6.1	279.4 ± 5.7	–2.6
B-15.1	0.3	114	0	0.0	0.3227 ± 3.5	0.0432 ± 2.7	0.55	0.0542 ± 2.3	378.9 ± 50.0	272.6 ± 7.2	284.0 ± 8.8	4.0
B-16.1	0.4	104	44	0.4	0.3143 ± 4.3	0.0451 ± 2.6	0.43	0.0506 ± 3.4	221.9 ± 77.4	284.1 ± 7.2	277.5 ± 10.5	–2.4
B-17.1	0.0	333	8	0.0	0.3119 ± 2.7	0.0435 ± 2.0	0.54	0.0520 ± 1.8	283.9 ± 40.0	274.7 ± 5.5	275.7 ± 6.6	0.4
B-18.1	0.7	82	117	1.5	0.3262 ± 4.8	0.0459 ± 2.7	0.41	0.0516 ± 3.9	265.9 ± 87.4	289.2 ± 7.6	286.7 ± 12.0	–0.8
B-19.1	0.0	260	40	0.2	0.3939 ± 8.6	0.0579 ± 6.9	0.58	0.0493 ± 5.1	162.5 ± 115.4	363.0 ± 24.2	337.2 ± 24.9	–7.6
B-20.1	0.1	142	76	0.5	0.3211 ± 3.4	0.0466 ± 1.1	0.22	0.0500 ± 3.2	195.1 ± 73.4	293.5 ± 3.1	282.8 ± 8.5	–3.8
B-21.1	0.6	107	1	0.0	0.3188 ± 3.8	0.0441 ± 2.3	0.43	0.0524 ± 3.0	302.9 ± 66.6	278.3 ± 6.2	281.0 ± 9.3	1.0
113106 (granophyre)
1.1	0.2	261	141	0.6	0.7556 ± 2.4	0.0937 ± 1.3	0.38	0.0585 ± 2.0	547 ± 43	578 ± 7	572 ± 11	–1.0
2.1	0.2	1109	37	0.0	0.4345 ± 2.3	0.0584 ± 2.0	0.64	0.0539 ± 0.9	368 ± 21	366 ± 7	366 ± 7	0.0
3.1	10.2	817	239	0.3	0.3660 ± 4.8	0.0501 ± 1.6	0.24	0.0530 ± 4.5	327 ± 100	315 ± 5	317 ± 13	0.4
4.1	0.1	1163	340	0.3	0.3483 ± 2.0	0.0491 ± 1.8	0.63	0.0514 ± 0.9	259 ± 20	309 ± 5	303 ± 5	–2.0
5.1	0.7	1146	30	0.0	0.6012 ± 3.7	0.0781 ± 3.5	0.69	0.0559 ± 1.0	447 ± 21	485 ± 17	478 ± 14	–1.4
6.1	1.4	1195	157	0.1	0.3785 ± 3.8	0.0463 ± 1.3	0.24	0.0593 ± 3.6	579 ± 76	292 ± 4	326 ± 11	10.6
6.2	0.2	926	231	0.3	0.6811 ± 2.1	0.0842 ± 1.9	0.65	0.0587 ± 0.8	555 ± 18	521 ± 9	527 ± 9	1.2
7.1	0.0	870	139	0.2	0.5906 ± 1.6	0.0751 ± 1.3	0.58	0.0570 ± 0.9	492 ± 19	467 ± 6	471 ± 6	0.8
8.1	3.1	289	115	0.4	0.6443 ± 3.6	0.0741 ± 2.5	0.50	0.0631 ± 2.5	710 ± 53	461 ± 11	505 ± 14	8.8
8.2	18.9	297	192	0.7	0.9724 ± 20.1	0.0942 ± 3.6	0.13	0.0749 ± 19.8	1065 ± 354	580 ± 20	690 ± 106	15.8
9.1	0.1	1109	16	0.0	0.3342 ± 3.2	0.0474 ± 2.3	0.50	0.0512 ± 2.3	249 ± 52	298 ± 7	293 ± 8	–1.8
9.2	1.1	87	15	0.2	0.7106 ± 6.1	0.0971 ± 1.3	0.15	0.0531 ± 6.0	332 ± 130	598 ± 8	545 ± 26	–9.6
10.2	1.0	826	575	0.7	1.0089 ± 4.2	0.1127 ± 1.8	0.31	0.0649 ± 3.8	771 ± 78	689 ± 12	708 ± 22	2.8

Uncertainties are given at the one-sigma level. Common Pb corrected using measured ²⁰⁴Pb. ρ: error correlation between ²⁰⁶Pb/²³⁸U and ²⁰⁷Pb/²³⁵U ratios. ¹ Disc. = U−Pb discordance; difference between the ²⁰⁴Pb-corrected ²⁰⁶Pb/²³⁸U and ²⁰⁷Pb/²³⁵U ages.

References

1. Dunn, A.M.; Reynolds, P.H.; Clarke, D.B.; Ugidos, J.M. A comparison of the age and composition of the Sherburne Dyke, Nova Scotia, and the Messejana Dyke, Spain. Can. J. Earth Sci.; 1998; 35, pp. 1110-1115. [DOI: https://dx.doi.org/10.1139/e98-058]

2. de Vicente, G.; Olaiz, A.; Muñoz-Martín, A.; Cunha, P.P. Longest and still longer: The Messejana-Plasencia dyke and its links with later Alpine deformation belt in Iberia. Tectonophysics; 2021; 815, 229009. [DOI: https://dx.doi.org/10.1016/j.tecto.2021.229009]

3. Greenough, J.D.; Hodych, J.D. Evidence for lateral magma injection in the Early Mesozoic dykes of eastern North America. Mafic Dykes and Emplacement Mechanisms; Parker, R.J.; Rickwood, P.C.; Tucker, D.H. Balkema: Rotterdam, The Netherlands, 1990; pp. 35-46.

4. Dalziel, I.W.D.; Lawver, L.A.; Murphy, J.B. Plumes, orogenesis, and supercontinental fragmentation. Earth Planet. Sci. Lett.; 2000; 78, pp. 1-11. [DOI: https://dx.doi.org/10.1016/S0012-821X(00)00061-3]

5. Torre de Assunçao, C. Sobre uma intrusão doleritica no Antracolitico do Baixo-Alentejo. Bol. Soc. Portu. Cienc. Naturais; 1949; 2, pp. 66-74.

6. Torre de Assunçao, C. Os doleritos pigeoníticos de quimismo pacífico do sul de Portugal. Boletim do Museu e Laboratório Mineralógico e Geológico da Faculdade de Ciências; Universidade de Lisboa: Lisboa, Portugal, 1951; Volume 19, pp. 109-117.

7. Cebriá, J.M.; López-Ruiz, J.; Doblas, M.; Martins, L.T.; Munha, J. Geochemistry of the Early Jurassic Messejana-Plasencia dyke (Portugal-Spain); Implications on the Origin of the Central Atlantic Magmatic Province. J. Petrol.; 2003; 44, pp. 547-568. [DOI: https://dx.doi.org/10.1093/petrology/44.3.547]

8. Villaseca, C.; Orejana, D.; Higueras, P.; Pérez-Soba, C.; García Serrano, J.; Lorenzo, S. The evolution of the subcontinental mantle beneath the Central Iberian Zone: Geochemical tracking of its mafic magmatism from the Neoproterozoic to the Cenozoic. Earth Sci. Rev.; 2022; 228, 103997. [DOI: https://dx.doi.org/10.1016/j.earscirev.2022.103997]

9. Alibert, C. A Sr-Nd isotope and REE study of Late Triassic dolerites from the Pyrenees (France) and the Messejana dyke (Spain and Portugal). Earth Planet. Sci. Lett.; 1985; 73, pp. 81-90. [DOI: https://dx.doi.org/10.1016/0012-821X(85)90036-6]

10. Callegaro, S.; Rapaille, C.; Marzoli, A.; Bertrand, H.; Chiaradia, M.; Reisberg, L.; Bellieni, G.; Martins, L.; Madeira, J.; Mata, J. et al. Enriched mantle source for the Central Atlantic magmatic province: New supporting evidence from southwestern Europe. Lithos; 2014; 188, pp. 15-32. [DOI: https://dx.doi.org/10.1016/j.lithos.2013.10.021]

11. Schermerhorn, L.J.G.; Priem, H.N.A.; Boelrijk, N.A.I.M.; Hebeda, E.H.; Verdurmen, E.A.T.; Verschure, R.H. Age and origin of the Messejana dolerite fault–dike system (Portugal and Spain) in the light of the opening of the North Atlantic Ocean. J. Geol.; 1978; 86, pp. 299-309. [DOI: https://dx.doi.org/10.1086/649692]

12. García de Figuerola, L.C.; Carnicero, A. El extremo noreste del gran dique del Alentejo-Plasencia. Stud. Geol. Salmant.; 1973; 6, pp. 73-84.

13. Dostal, J.; Dunning, M. Geochemical constraints on the origin and evolution of early Mesozoic dikes in Atlantic Canada. Eur. J. Mineral.; 1998; 10, pp. 79-93. [DOI: https://dx.doi.org/10.1127/ejm/10/1/0079]

14. Nomade, S.; Pouclet, A.; Chen, Y. The French Guyana doleritic dykes: Geochemical evidence of three populations and new data for the Jurassic Central Atlantic Magmatic Province. J. Geodyn.; 2002; 34, pp. 595-614. [DOI: https://dx.doi.org/10.1016/S0264-3707(02)00034-0]

15. Kontak, D.J.; De Young, M.Y.D.; Dostal, J. Late-stage crystallization history of the Jurassic North Mountain Basalt, Nova Scotia, Canada. I. Textural and chemical evidence for pervasive development of silicate-liquid immiscibility. Can. Mineral.; 2002; 40, pp. 1287-1311. [DOI: https://dx.doi.org/10.2113/gscanmin.40.5.1287]

16. Halder, M.; Paul, D.; Yang, S.Y. Origin of silicic rocks of the Deccan Traps continental flood basalt province: Inferences from field observations, petrography, and geochemistry. Geochem. Ger.; 2023; 83, 125958. [DOI: https://dx.doi.org/10.1016/j.chemer.2023.125958]

17. Philpotts, A.R. Composition of immiscible liquids in volcanic rocks. Contrib. Mineral. Petrol.; 1982; 80, pp. 201-218. [DOI: https://dx.doi.org/10.1007/BF00371350]

18. Andreeva, O.A.; Andreeva, I.A.; Yarmolyuk, V.V.; Ji, J.Q.; Zhou, X.; Borisovskii, S.E. Silicate Liquid Immiscibility as a Result of Fenner-Type Crystal Fractionation of Wangtian’e Tholeiitic Melts, Northeast China. Petrology; 2020; 28, pp. 357-373. [DOI: https://dx.doi.org/10.1134/S0869591120040025]

19. Martínez-Catalán, J.R. The Central Iberian arc, an orocline centered in the Iberian Massif and some implications for the Variscan belt. Int. J. Earth Sci.; 2012; 101, pp. 1299-1314. [DOI: https://dx.doi.org/10.1007/s00531-011-0715-6]

20. Kroner, U.; Romer, R.L. Two plates—Many subduction zones: The Variscan orogeny reconsidered. Gondwana Res.; 2013; 24, pp. 298-329. [DOI: https://dx.doi.org/10.1016/j.gr.2013.03.001]

21. Orejana, D.; Villaseca, C.; Kristoffersen, M. Geochemistry and geochronology of mafic rocks from the Spanish Central System: Constraints on the mantle evolution beneath central Spain. Geosci. Front.; 2020; 11, pp. 1651-1667. [DOI: https://dx.doi.org/10.1016/j.gsf.2020.01.002]

22. Scarrow, J.H.; Molina, J.F.; Bea, F.; Montero, P.; Vaughan, A.P.M. Lamprophyre dikes as tectonic markers of late orogenic transtension timing and kinematics: A case study from the Central Iberian Zone. Tectonics; 2011; 30, TC4007. [DOI: https://dx.doi.org/10.1029/2010TC002755]

23. Oliveira, A.J.T.; Martins, H.C.B.; da Silva, H.M.S.M. Permo-Carboniferous hypabyssal magmatism in northern Portugal: The case of the Lamas de Olo microgranite and lamprophyre dykes. Iber. J. Geol.; 2022; 48, pp. 1-28. [DOI: https://dx.doi.org/10.1007/s41513-021-00179-8]

24. Sebai, A.; Féraud, G.; Bertrand, H.; Hanes, J. 40Ar/39Ar dating and geochemistry of tholeiitic magmatism related to the early opening of the Central Atlantic rift. Earth Planet. Sci. Lett.; 1991; 104, pp. 455-472. [DOI: https://dx.doi.org/10.1016/0012-821X(91)90222-4]

25. Osete, M.L.; Palencia Ortas, A. Iberian palaeomagnetic poles for the last 300 Ma. Fis. Tierra; 2006; 18, pp. 157-181.

26. García de Figuerola, L.C.; Corretge, L.G.; Bea, F. El dique de Alentejo-Plasencia y haces de diques básicos de Extremadura. Bol. Geol. Min. Esp.; 1974; 85, pp. 308-337.

27. Bertrand, H.; Millot, G. Le magmatisme tholéitique continental de la marge Ibérique, précurseur de l’overture de l’Atlantique Central: Les dolérites du dyke de Messejana–Plasencia (Portugal–Espagne). C. R. Geosci.; 1987; 304, pp. 215-220.

28. Martins, L.T.; Cebriá, J.M.; Doblas, M.; López-Ruiz, J.; Munhá, J. Petrologia e geoquímica do dique de Messejana. Mem. Mus. Lab. Geol. Mineral. Porto; 1995; 4, pp. 747-750.

29. Andonaegui, P.; Villaseca, C.; López García, J.A. Caracterización de los minerales del dique gabroideo de Messejana-Plasencia en su sector nor-oriental. Geogaceta; 2005; 38, pp. 219-222.

30. Dahlquist, J.A.; Galindo, C.; Morales Camera, M.M.; Moreno, J.A.; Alasino, P.H.; Basei, M.A.S.; Macchioli Grande, M. A combined zircon Hf isotope and whole rock Nd and Sr isotopes study of Carboniferous A-type granites, Sierras Pampeanas of Argentina. J. S. Am. Earth Sci.; 2020; 100, 102545. [DOI: https://dx.doi.org/10.1016/j.jsames.2020.102545]

31. Black, L.P.; Kamo, S.L.; Allen, C.M.; Aleinikoff, J.N.; Davis, D.W.; Korsch, R.J.; Foudoulis, C. TEMORA 1: A new zircon standard for Phanerozoic U–Pb geochronology. Chem. Geol.; 2003; 200, pp. 155-170. [DOI: https://dx.doi.org/10.1016/S0009-2541(03)00165-7]

32. Black, L.P.; Kamo, S.L.; Allen, C.M.; Davis, D.W.; Aleinikoff, J.N.; Valley, J.W.; Mundil, R.; Campbell, I.H.; Korsch, R.J.; Williams, I.S. et al. Improved Pb/U microprobe geochronology by the monitoring of a trace-element-related matrix effect; SHRIMP, ID-TIMS, ELA-ICP-MS and oxygen isotope documentation for a series of zircon standards. Chem. Geol.; 2004; 205, pp. 115-140. [DOI: https://dx.doi.org/10.1016/j.chemgeo.2004.01.003]

33. Ludwig, K.R. ISOPLOT/Ex, version 3. A geochronological toolkit for Microsoft Excel. Berkeley Geochr. Cent. Spec. Pub.; 2003; 4, pp. 1-71.

34. Whitney, D.L.; Evans, B.W. Abbreviations for names of rock-forming minerals. Am. Mineral.; 2010; 95, pp. 185-187. [DOI: https://dx.doi.org/10.2138/am.2010.3371]

35. Rossi, P.; Chevremont, P. Classification des associations magmatiques granitoides. Geochronique; 1987; 21, pp. 14-18.

36. Nachit, H.; Razafimahefa, N.; Stussi, J.M.; Carron, J.P. Biotite chemical-compositions and magmatic typology of the granitoids. C. R. Geosci.; 1985; 301, pp. 813-817.

37. Rollinson, H.R. Using Geochemical Data: Evaluation, Presentation, Interpretation; Longman Scientific & Technical: New York, NY, USA, 1993; 352p.

38. Brandebourger, E. Les Granitoides Hercyniens Tardifs de la Sierra de Guadarrama (Système Central, Espagne): Pétrographie et Geochimie. Ph.D. Thesis; Université de Lorraine: Metz, France, 1984; 209p.

39. Casillas, R. Las Asociaciones Plutónicas Tardihercínicas del Sector Occidental de la Sierra de Guadarrama, Sistema Central Español, (Las Navas del Marqués–San Martín de Valdeiglesias). Petrología, Geoquímica, Génesis y Evolución. Ph.D. Thesis; Universidad Complutense de Madrid: Madrid, Spain, 1989; 316p.

40. Villaseca, C.; Barbero, L.; Rogers, G. Crustal origin of Hercynian peraluminous granitic batholiths of central Spain: Petrological, geochemical and isotopic (Sr,Nd) arguments. Lithos; 1998; 43, pp. 55-79. [DOI: https://dx.doi.org/10.1016/S0024-4937(98)00002-4]

41. Sun, S.S.; McDonough, W.F. Chemical and isotopic systematics of oceanic basalts; implications for mantle composition and processes. Magmatism in Ocean Basins; Saunders, A.D.; Norrey, M.J. Geological Society Special Publication: Oxford, UK, 1989; Volume 42, pp. 313-345.

42. McDonough, W.F.; Sun, S.S. The composition of the Earth. Chem. Geol.; 1995; 120, pp. 223-253. [DOI: https://dx.doi.org/10.1016/0009-2541(94)00140-4]

43. Villaseca, C.; Merino Martínez, E.; Orejana, D.; Andersen, T.; Belousova, E.A. Zircon Hf signatures from granitic orthogneisses of the Spanish Central System: Significance and sources of the Cambro-Ordovician magmatism in the Iberian Variscan Belt. Gondwana Res.; 2016; 34, pp. 60-83. [DOI: https://dx.doi.org/10.1016/j.gr.2016.03.004]

44. Orejana, D.; Villaseca, C.; Elburg, M.A.; Merino-Martínez, E.; Serrano, J.G. Nature and evolution of the lower crust under central Spain: Inferences from granulite xenoliths (Calatrava Volcanic Field-Spanish central system). Geosci. Front.; 2023; 14, 101525. [DOI: https://dx.doi.org/10.1016/j.gsf.2022.101525]

45. Orejana, D.; Villaseca, C.; Valverde-Vaquero, P.; Belousova, E.A.; Armstrong, R.A. U–Pb geochronology and zircon composition of late Variscan S- and I-type granitoids from the Spanish Central System batholith. Int. J. Earth Sci.; 2012; 101, pp. 1789-1815. [DOI: https://dx.doi.org/10.1007/s00531-012-0750-y]

46. Orejana, D.; Villaseca, C.; Armstrong, R.A.; Jeffries, T.E. Geochronology and trace element chemistry of zircon and garnet from granulite xenoliths: Constraints on the tectonothermal evolution of the lower crust under central Spain. Lithos; 2011; 124, pp. 103-116. [DOI: https://dx.doi.org/10.1016/j.lithos.2010.10.011]

47. Fernández Suárez, J.; Arenas, R.; Jeffries, T.E.; Whitehouse, M.J.; Villaseca, C. A U–Pb study of zircons from a lower crustal granulite xenolith of the Spanish Central system: A record of Iberian lithospheric evolution from the Neoproterozoic to the Triassic. J. Geol.; 2006; 114, pp. 471-483. [DOI: https://dx.doi.org/10.1086/504180]

48. Cox, K.G.; MacDonald, R.; Hornung, G. Geochemical and petrographic provinces in the Karroo basalts of southern Africa. Am. Mineral.; 1967; 52, pp. 1451-1474.

49. Thy, P.; Lesher, C.E.; Nielsen, T.F.D.; Brooks, C.K. Experimental constraints on the Skaergaard liquid line of descent. Lithos; 2006; 92, pp. 154-180. [DOI: https://dx.doi.org/10.1016/j.lithos.2006.03.031]

50. Villiger, S.; Ulmer, P.; Müntener, O. Equilibrium and fractional crystallization experiments at 0.7 GPa: The effect of pressure on phase relations and liquid compositions of tholeiitic magmas. J. Petrol.; 2007; 48, pp. 159-184. [DOI: https://dx.doi.org/10.1093/petrology/egl058]

51. Sisson, T.W.; Ratajeski, K.; Hankins, W.B.; Glazner, A.F. Voluminous granitic magmas from common basaltic sources. Contrib. Mineral. Petrol.; 2005; 148, pp. 635-661. [DOI: https://dx.doi.org/10.1007/s00410-004-0632-9]

52. Hotz, P.E. Petrology of granophyre in diabase near Dillsburg, Pennsylvania. Geol. Soc. Am. Bull.; 1953; 64, pp. 675-704. [DOI: https://dx.doi.org/10.1130/0016-7606(1953)64[675:POGIDN]2.0.CO;2]

53. Marsh, B.D. Solidification fronts and magmatic evolution. Mineral. Mag.; 1996; 60, pp. 5-40. [DOI: https://dx.doi.org/10.1180/minmag.1996.060.398.03]

54. White, C.M. The graveyard point intrusion: An example of extreme differentiation of Snake River Plain basalt in a shallow crustal pluton. J. Petrol.; 2007; 48, pp. 303-325. [DOI: https://dx.doi.org/10.1093/petrology/egl062]

55. Johannes, E.; Holtz, F. Petrogenesis and Experimental Petrology of Granitic Rocks. Minerals, Rocks and Mountains; Springer: Berlin, Germany, 1996; 335p.

56. Winter, J.D. An Introduction to Igneous and Metamorphic Petrology; 1st ed. Prentince Hall: Upper Saddle River, NJ, USA, 2001; 697p.

57. Charlier, B.; Namur, O.; Toplis, M.J.; Schiano, P.; Cluzel, N.; Higgins, M.D.; Vander Auwera, J. Large-scale silicate liquid immiscibility during differentiation of tholeiitic basalt to granite and the origin of the Daly gap. Geology; 2011; 39, pp. 907-910. [DOI: https://dx.doi.org/10.1130/G32091.1]

58. Geist, D.; Howard, K.A.; Larson, P. The Generation of Oceanic Rhyolites by Crystal Fractionation—The Basalt-Rhyolite Association at Volcan-Alcedo, Galapagos-Archipelago. J. Petrol.; 1995; 36, pp. 965-982. [DOI: https://dx.doi.org/10.1093/petrology/36.4.965]

59. Peccerillo, A. Relationships between Mafic and Peralkaline Silicic Magmatism in Continental Rift Settings: A Petrological, Geochemical and Isotopic Study of the Gedemsa Volcano, Central Ethiopian Rift. J. Petrol.; 2003; 44, pp. 2003-2032. [DOI: https://dx.doi.org/10.1093/petrology/egg068]

60. Dufek, J.; Bachmann, O. Quantum magmatism: Magmatic compositional gaps generated by melt-crystal dynamics. Geology; 2010; 38, pp. 687-690. [DOI: https://dx.doi.org/10.1130/G30831.1]

61. McBirney, A.R. Differentiation of the Skaergaard intrusion. Nature; 1975; 253, pp. 691-694. [DOI: https://dx.doi.org/10.1038/253691a0]

62. Jakobsen, J.K.; Veksler, I.V.; Tegner, C.; Brooks, C.K. Immiscible iron- and silica-rich melts in basalt petrogenesis documented in the Skaergaard intrusion. Geology; 2005; 33, pp. 885-888. [DOI: https://dx.doi.org/10.1130/G21724.1]

63. Ripley, E.M.; Severson, M.J.; Hauck, S.A. Evidence for sulfide and Fe-Ti-P-rich liquid immiscibility in the Duluth Complex, Minnesota. Econ. Geol. Bull. Soc.; 1998; 93, pp. 1052-1062. [DOI: https://dx.doi.org/10.2113/gsecongeo.93.7.1052]

64. Namur, O.; Charlier, B.; Holness, M.B. Dual origin of Fe-Ti-P gabbros by immiscibility and fractional crystallization of evolved tholeiitic basalts in the Sept Iles layered intrusion. Lithos; 2012; 154, pp. 100-114. [DOI: https://dx.doi.org/10.1016/j.lithos.2012.06.034]

65. VanTongeren, J.A.; Mathez, E.A. Large-scale liquid immiscibility at the top of the Bushveld Complex, South Africa. Geology; 2012; 40, pp. 491-494. [DOI: https://dx.doi.org/10.1130/G32980.1]

66. Charlier, B.; Namur, O.; Grove, T.L. Compositional and kinetic controls on liquid immiscibility in ferrobasalt-rhyolite volcanic and plutonic series. Geochim. Cosmochim. Acta; 2013; 113, pp. 79-93. [DOI: https://dx.doi.org/10.1016/j.gca.2013.03.017]

67. Roedder, E.W. Low temperature liquid immiscibility field in the system K₂O–FeO–Al₂O₃–SiO₂. Am. Mineral.; 1951; 36, pp. 282-286.

68. Veksler, I.V.; Dorfman, A.M.; Danyushevsky, L.V.; Jakobsen, J.K.; Dingwell, D.B. Immiscible silicate liquid partition coefficients: Implications for crystal-melt element partitioning and basalt petrogenesis. Contrib. Mineral. Petrol.; 2006; 152, pp. 685-702. [DOI: https://dx.doi.org/10.1007/s00410-006-0127-y]

69. Villaseca, C.; Downes, H.; Pin, C.; Barbero, L. Nature and composition of the lower continental crust in central Spain and the granulite-granite linkage: Inferences from granulitic xenoliths. J. Petrol.; 1999; 40, pp. 1465-1496. [DOI: https://dx.doi.org/10.1093/petroj/40.10.1465]

70. Grieco, G.; Ferrario, A.; von Quadt, A.; Koeppel, V.; Mathez, E.A. The zircon-bearing chromitites of the phlogopite peridotite of Finero (Ivrea zone, southern Alps): Evidence and geochronology of a metasomatized mantle slab. J. Petrol.; 2001; 42, pp. 89-101. [DOI: https://dx.doi.org/10.1093/petrology/42.1.89]

71. Liu, Y.; Gao, S.; Hu, Z.; Gao, C.; Zong, K.; Wang, D. Continental and oceanic crust recycling-induced melt-peridotite interactions in the Trans-North China Orogen: U–Pb dating, Hf isotopes and trace elements in zircons from mantle xenoliths. J. Petrol.; 2010; 51, pp. 537-571. [DOI: https://dx.doi.org/10.1093/petrology/egp082]

72. González-Jiménez, J.M.; Marchesi, C.; Griffin, W.L.; Gervilla, F.; Belousova, E.A.; Garrido, C.J.; Romero, R.; Talavera, C.; Leisen, M.; O’Reilly, S.Y. et al. Zircon recycling and crystallization during formation of chromite- and Ni-arsenide ores in the subcontinental lithospheric mantle (Serranía de Ronda, Spain). Ore Geol. Rev.; 2017; 90, pp. 193-209. [DOI: https://dx.doi.org/10.1016/j.oregeorev.2017.02.012]

73. Gutiérrez-Alonso, G.; Murphy, J.B.; Fernández-Suárez, J.; Weil, A.B.; Franco, M.P.; Gonzalo, J.C. Lithospheric delamination in the core of Pangea: Sm–Nd insights from the Iberian mantle. Geology; 2011; 39, pp. 155-158. [DOI: https://dx.doi.org/10.1130/G31468.1]

74. DePaolo, D.J. Trace element and isotopic effects of combined wallrock assimilation and fractional crystallization. Earth Planet. Sci. Lett.; 1981; 53, pp. 189-202. [DOI: https://dx.doi.org/10.1016/0012-821X(81)90153-9]

75. Wenner, D.; Taylor, H. Oxygen and hydrogen isotope studies of a precambrian granite-rhyolite terrane, St-Francois mountains, southeastern Missouri. Geol. Soc. Am. Bull.; 1976; 87, pp. 1587-1598. [DOI: https://dx.doi.org/10.1130/0016-7606(1976)87<1587:OAHISO>2.0.CO;2]

76. Walraven, F.; Strydom, J.H.; Strydom, N. Rb–Sr open-system behaviour and its application as a pathfinder for Sn mineralisation in granites of the Bushveld Complex, South Africa. J. Geochem. Expl.; 1990; 37, pp. 333-350. [DOI: https://dx.doi.org/10.1016/0375-6742(90)90107-L]

77. Brantley, S.L.; Chesley, J.T.; Stillings, L.L. Isotopic ratios and release rates of strontium from weathering feldspars. Geochim. Cosmochim. Acta; 1998; 62, pp. 1493-1500. [DOI: https://dx.doi.org/10.1016/S0016-7037(98)00082-9]

78. Ma, Y.; Liu, C. Sr isotope evolution during chemical weathering of granites—Impact of relative weathering rates of minerals. Sci. China (Earth Sci.); 2001; 44, pp. 726-734. [DOI: https://dx.doi.org/10.1007/BF02907202]

79. Kaur, P.; Chaudhri, N.; Hofmann, A.W.; Raczek, I.; Okrusch, M.; Skora, S.; Baumgartner, L.P. Two-Stage, extreme albitization of A-type granites from Rajasthan, NW India. J. Petrol.; 2012; 53, pp. 919-948. [DOI: https://dx.doi.org/10.1093/petrology/egs003]

80. Ubanell, A.G.; Garzón, G.; de la Peña, J.A.; Bustillo, A.; Martil, R. Estudio de procesos de alteración hidrotermal en rocas graníticas y sedimentarias (provincia de Ávila). Est. Geol.; 1978; 34, pp. 151-160.

81. Parsapoor, A.; Khalili, M.; Maghami, M. Discrimination between mineralized and unmineralized alteration zones using primary geochemical haloes in the Darreh-Zar porphyry copper deposit in Kerman, southeastern Iran. J. Afr. Earth Sci.; 2017; 132, pp. 109-126. [DOI: https://dx.doi.org/10.1016/j.jafrearsci.2017.04.027]

82. Irving, A.J.; Frey, F.A. Trace element abundances in megacrysts and their host basalts: Constraints on partition coefficients and megacryst genesis. Geochim. Cosmochim. Acta; 1984; 48, pp. 1201-1221. [DOI: https://dx.doi.org/10.1016/0016-7037(84)90056-5]

83. Philpotts, J.A.; Schnetzler, C.C. Phenocryst-matrix partition coefficients for K, Rb, Sr and Ba, with applications to anorthosite and basalt genesis. Geochim. Cosmochim. Acta; 1970; 34, pp. 307-322. [DOI: https://dx.doi.org/10.1016/0016-7037(70)90108-0]

84. Bindeman, I.N.; Davis, A.M.; Drake, M.J. Ion Microprobe Study of Plagioclase-Basalt Partition Experiments at Natural Con-centration Levels of Trace Elements. Geochim. Cosmochim. Acta; 1998; 62, pp. 1175-1193. [DOI: https://dx.doi.org/10.1016/S0016-7037(98)00047-7]

85. Paster, T.P.; Schauwecker, D.S.; Haskin, L.A. The behavior of some trace elements during solidification of the Skaergaard layered series. Geochim. Cosmochim. Acta; 1974; 38, pp. 1549-1577. [DOI: https://dx.doi.org/10.1016/0016-7037(74)90174-4]