Pyroxenites and Megacrysts From Alkaline Melts of the Calatrava Volcanic Field (Central Spain): Inferences From Trace Element Geochemistry and Sr-Nd Isotope Composition

Alkaline volcanic rocks from explosive monogenetic centers often carry an unusual cargo of crystals and rock fragments, which may provide valuable constraints on magma source, ascent and eruption. One of such examples is the Cenozoic Calatrava Volcanic Field in central Spain, a still poorly explored area to address these issues. Clinopyroxene, amphibole and phlogopite appear either as megacryst/phenocrysts or forming fine-grained cumulates (pyroxenite enclaves s.l.) in some eruptive centers of this volcanic field. They have previously been interpreted as cogenetic high-P minerals formed within the upper lithospheric mantle. The presence of Fe-Na-rich green and Mg-Cr-rich colorless clinopyroxene types as phenocryst cores or as oscillatory zoned crystals in pyroxenite enclaves points to a complex evolution of mineral fractionates from petrogenetically related magmas. In trace element chemistry all studied clinopyroxene types show parallel rare earth element patterns irrespective of whether they are megacrysts, colorless or green core phenocrysts, or zoned crystals within pyroxenite cumulates. This similarity indicates a genetic relationship between all the fractionated minerals. This is in agreement with the overlapping of initial 143Nd/144Nd and 87Sr/86Sr ratios of pyroxenite enclaves (0.512793–0.512885 and 0.703268–0.703778) that is within the chemical field of the host magmas and the Calatrava volcanics. The initial 143Nd/144Nd and 87Sr/86Sr ratios of megacrystic clinopyroxene, amphibole and phlogopite show a more restricted range (0.512832–0.512890 and 0.703217–0.703466), also falling within the isotopic composition of the Calatrava volcanic rocks. Deep magmatic systems beneath monogenetic volcanic fields involve several stages of melt accumulation, fractionation and contamination at variable depths. Trace element and isotope mineral chemistry are powerful tools to understand the history of ascent and stagnation of alkaline basaltic magmas and discriminate between magma mixing, wall-rock contamination and closed magmatic system evolution. In our study, we establish a cogenetic origin for green and colorless clinopyroxene as high-pressure precipitates from liquids of different fractionation degrees (up to 80%, for the highly evolved melts equilibrated with the green clinopyroxene), originated from a highly solidified front of silica-undersaturated alkaline magmas at mantle reservoirs.


INTRODUCTION
The study of monogenetic volcanic fields provide important information regarding primary magma genesis, deep-seated differentiation processes and the architecture of the subvolcanic plumbing systems (e.g. Mattsson et al., 2013;Re et al., 2017). These magmas commonly carry a heterogeneous crystal cargo with complex compositional and textural features, relevant to interpreting the depths and rates of magma storage, transport, fractionation and mixing (e.g. Jankovics et al., 2016). Such data are also useful to understand the mechanisms which trigger volcanic eruption and control its characteristics (Mattsson et al., 2013).
Many monogenetic volcanoes are formed by alkaline rocks which transport mafic megacrysts and ultramafic xenoliths (e.g. Irving and Frey, 1984;Praegel, 1981;Jankovics et al., 2016), whereas pyroxenite enclaves are not so common (Downes, 2007). These enclaves have been interpreted as recycled subducted lithosphere (Allègre and Turcotte, 1986), metasomatic products (Garrido and Bodinier, 1999) or high pressure cumulates (e.g. Bodinier et al., 1987;Wilkinson and Stolz, 1997). When formed by deep-seated fractionation, clinopyroxenite enclaves contain volatile-rich minerals, such as amphibole and phlogopite, and are accompanied by equivalent megacrysts (Orejana et al., 2006 and references therein). Green and colorless anhedral to subhedral phenocryst/megacryst cores can also be found in alkaline volcanic rocks and have been studied in order to unravel their enigmatic origin (e.g. Duda and Schmincke, 1985;Ubide et al., 2014;Jankovics et al., 2016). In situ trace element data from such a variety of relic clinopyroxene cores are still sparse and completely lacking in most of the Cenozoic volcanic fields of the Iberian Peninsula, where mafic megacrysts and phenocrysts are common within the pyroclastic deposits of these mafic volcanic centers.
Mafic minerals are important petrogenetic indicators of the evolution of basic-ultrabasic magmas (e.g. Green and Ringwood, 1964;Irving and Frey, 1984;Liotard et al., 1988). Their geochemical composition helps to constrain whether they derive from a common parental melt or from diverse melt batches, and if any genetic affinity is possible between the host volcanic magma and any other melt formed in the context of an active magma plumbing system. A study on the major element geochemistry of mafic megacrysts (clinopyroxene, amphibole and phlogopite) and associated pyroxenite enclaves from two volcanoes of the Calatrava volcanic field was recently published by our research team (Villaseca et al., 2019a). The two host magmas are highly porphyritic, showing a complex crystal population of mafic megacrysts (>5 mm), macrocrysts (5-0.5 mm) and phenocrysts with significant geochemical similarities. Moreover, their major element composition resembles that of equivalent phases in clinopyroxenite enclaves, suggesting that all of them constitute a cogenetic suite of igneous origin. These fractionates were interpreted as high-P cumulates from mantle magma reservoirs whose rigid chilled margins (solid crusts) were fragmented by subsequent deep CO 2 boiling of the fractionated host magmas. The variable crystal cargo dragged by the volcanic magmas (crystals and pyroxenite cumulates) added to variable amounts of peridotite xenoliths entrained from the mantle wall-rocks (Villaseca et al., 2019a).
This work adds new trace element compositional data of minerals (megacrysts, phenocrysts and clinopyroxenite enclaves) within the host alkaline magmas of previously sampled Calatrava volcanoes. We have also determined the whole-rock composition of host rocks and selected clinopyroxenite enclaves along with the Sr-Nd isotope composition of host magmas, pyroxenites and megacrysts. The main objective of this study is to constrain the degree of consanguinity between megacrysts, enclaves and their host magmas and check the previously established cumulate origin. With this objective, our work aims to deepen the understanding of the origin of pyroxenitic and related (mica-and amphibole-rich) cumulates, as well as the generation of diverse types of colored (green to colorless) clinopyroxene pheno-tomegacrysts, in silica-undersaturated alkaline magmas. This work also contributes to a better characterisation of the mantle beneath Central Iberia and to the understanding of differentiation processes operating beneath monogenetic volcanic fields.

GEOLOGICAL SETTING
The Cenozoic Calatrava Volcanic Field comprises more than 200 monogenetic volcanic centers in an area of around 5500 km 2 (Ancochea, 1982). Volcanism in the Calatrava Field took place in two different stages (Ancochea, 1982): (1) a minor ultrapotassic event around 8.7-6.4 Ma and (2) alkaline basalts, basanites and olivine nephelinites and melilitites from 3.7 to 0.7 Ma. The two studied volcanic centers have been recently dated by Ar-Ar on mafic megacrysts yielding ages of 2.2 ± 0.04 Ma (phlogopite from the Cerro Pelado) and 2.8 ± 0.1 Ma (amphibole from the El Aprisco) (Villaseca et al., 2019b). The Pliocene ages of mafic megacrysts from both volcanoes are in agreement with their derivation from the host magmatic systems and represent ages much younger than mantle metasomatic events dated in peridotite xenoliths from the El Aprisco volcano (Villaseca et al., 2019b).
Two volcanoes were sampled due to the common presence of mafic megacrysts and pyroxenite enclaves in their pyroclastic deposits: the Cerro Pelado scoria cone (olivine nephelinite) and the El Aprisco maar (olivine melilitite) (Supplementary Figure S1). The variety and size of these mafic megacrysts are higher in the Cerro Pelado pyroclastic fall deposits than in other Calatrava centers (Villaseca et al., 2019a). Moreover, mantle xenoliths are also present in both volcanoes and they have been described and interpreted as peridotites that have undergone small-to-moderate degrees of partial melting (≤10%) (Villaseca et al., 2010). These peridotite xenoliths have been overprinted by different metasomatic agents, mostly alkaline silica-undersaturated melts (Villaseca et al., 2010;Lierenfeld and Mattsson, 2015) but also carbonatite fluids (e.g. González-Jiménez et al., 2014;Villaseca et al., 2019b).
The trace element composition and Sr-Nd isotopic data of primary alkaline magmas suggest that most of the Calatrava volcanics are derived from enriched asthenospheric sources, similar to those defined for the European asthenospheric mantle (Cebriá and López Ruiz, 1995;Granet et al., 1995). Clinopyroxene, amphibole and phlogopite megacrysts are found in different Calatrava volcanoes, but the hydrous mafic megacrysts (amphibole, phlogopite) recorded in the pyroclastic fall deposits of the Cerro Pelado scoria cone outstand by their large size (up to 8.5 cm; Villaseca et al., 2019a). Marginal pyroclastic surge layers of olivine melilitites from the El Aprisco maar also show fragments of clinopyroxene and amphibole megacrysts with a heterogeneous size. Ultramafic clinopyroxenerich enclaves with variable amounts of hydrous minerals (amphibole and phlogopite), and scarce olivine, can be found associated with these megacrysts in both volcanoes. A rare variety of phlogopite-rich clinopyroxenite (glimmerite, mica > 90 vol%) occurs in the El Aprisco maar (e.g. Villaseca et al., 2019a).

SAMPLES AND ANALYTICAL METHODS
Four samples were collected from the host mafic rocks. It is important to note that, due to their pyroclastic nature, the bulk of these samples may include contributions from xenoliths (peridotites of mantle derivation and shallow crustal fragments), crystal cargo and hydrothermal alteration. Accordingly, it has been difficult to obtain geochemical data representative of the erupting melts (see also Lierenfeld and Mattsson, 2015). Only one sample from the Cerro Pelado volcanic center (CAL-73) plotted within the main Calatrava volcanic rock data of Cebriá and López Ruiz (1995), whereas the other three (one from El Aprisco and other two from Cerro Pelado) must be considered crystal-contaminated samples. On the other hand, three clinopyroxenite enclaves were also sampled, and there are five additional samples from the Calatrava volcanoes (eruptive centers in Table 1).
Whole-rock samples were fused using a lithium metaboratetetraborate mixture. The melt produced by this process was completely dissolved with 5% HNO 3 . Major oxide analyses were carried out using a Variant Vista 735 ICP while trace elements were diluted and analyzed by a Perkin Elmer Sciex ELAN 6000 mass spectrometer (ICP-MS) following code 4 Lithoresearch at Activation Laboratories (ACTLABS, Canada).
Relative uncertainties for major elements are bracketed between 1 and 3%, except for MnO (5-10%). The precision for Rb, Sr, Zr, Y, V, Hf and most of the rare earth elements (REE) range from 1 to 5%, and between 5 and 10% for the rest of trace elements. Some samples have concentrations in Pb below detection limits (5 ppm). The precision was evaluated from repeated analyses of the international standards NIST 694, DNC-1, GBW071113, W-2, SY-4, JR-1 and BIR-1a (Supplementary Table S1). More information on the procedure, precision and accuracy of ACTLABS ICP-MS analyses can be found at www.actlabs.com. Whole-rock chemical compositions are shown in Table 1.
Sr-Nd isotopic ratios were measured on four volcanic rocks (including samples from the Cerro Pelado and El Aprisco host rocks, together with volcanics from other eruptive centers of the Calatrava Volcanic Field), three clinopyroxenites and four mafic megacrysts ( Table 2) at the Centro de Asistencia a la Investigación (CAI) of Geochronology and Isotope Geochemistry (Complutense University of Madrid, Spain). Whole-rock and megacryst samples were dissolved in ultrapure reagents and the isotopes were subsequently isolated by exchange chromatography. Isotope analyses were carried out using a Phoenix-IsotopX Multicollector Thermal Ionization Mass Spectrometer with data acquired in multidynamic mode. The analytical procedures used in this laboratory have been described elsewhere (Reyes et al., 1997). Repeated analyses on the NBS-987 standard gave 87 Sr/ 86 Sr = 0.710240 ± 0.00005 (2σ, n = 8) and for the La Jolla JNdi-1 standard, values of 143 Nd/ 144 Nd = 0.512109 ± 0.000005 (2σ, n = 16) were obtained (Supplementary Table S1). The BHVO-2 standard was also used as internal reference for precision. The 2σ analytical errors are 0.01% for 87 Sr/ 86 Sr and 0.006% for 143 Nd/ 144 Nd.
Trace element analysis of mafic minerals (clinopyroxene, amphibole, phlogopite) was performed on the mineral grains previously analyzed by EMP (Villaseca et al., 2019a). Two clinopyroxene and amphibole megacrysts (from each volcanic center) and one phlogopite megacryst (from the Cerro Pelado scoria cone), along with two phenocryst-rich samples (one from each volcano), were selected for trace element analyses (Supplementary Tables S2, S3). Two phlogopite-amphibole pyroxenites, one amphibole pyroxenite and a phlogopite-rich (glimmerite) pyroxenite were also used for trace element studies. Trace element mineral compositions were obtained by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) at the Instituto Andaluz de Ciencias de la Tierra (CSIC, UGR, Granada) using a CETAC-Photon Machines 193 nm laser attached to an Agilent 8800 ICP-MS. The diameter of the laser beam was 40 to 60 µm associated to repetition rates of 10 Hz and laser fluence at the target of ca. 8 J/cm 2 . A 30 s gas blank was analyzed first to establish the background, followed by 60 s measurements for the remainder of the analysis. Relative element sensitivities were calibrated with a NIST SRM 611 glass standard. Then, each analysis was normalized to Si or Ca using concentrations determined by electron microprobe. Precision and accuracy were assessed from repeated analyses of the USGS-BIR 1G standard and were estimated to be between 2 and 10% for most of the analyzed trace elements.

PETROGRAPHY
Below we summarize the textural features of the studied mafic crystals and clinopyroxe-rich enclaves that appear as crystal cargo in alkaline basaltic rocks of the Calatrava Volcanic Field. A detailed petrographic description of the studied samples can be found in Villaseca et al. (2019a). Clinopyroxene, amphibole and phlogopite megacrysts are single pyroclasts (lapillus) or small fragments in the Cerro Pelado and El Aprisco pyroclastic rocks. Clinopyroxene megacrysts are colorless and mostly unzoned, except for a thin reaction rim with the host magma. In contrast, the Cerro Pelado volcanic rocks show two types of relic clinopyroxene phenocryst cores: (1) anhedral green and (2) anhedral to euhedral colorless. In the El Aprisco maar only the second type of relic clinopyroxene core (colorless) was observed. Both phenocryst cores develop pale purple rims petrographically similar to melilitite/nephelinite groundmass clinopyroxene (Figures 1A,C). Amphibole makes up the largest megacrysts (up to 8.5 cm in the Cerro Pelado volcano), and phlogopite is usually unzoned and no larger than 4.5 cm. Most megacrysts display reaction rims or corrosion gulfs.
The studied pyroclastic rocks also include small black and elongated enclaves that can be classified as amphibole clinopyroxenite, amphibole-phlogopite clinopyroxenite, phlogopite clinopyroxenite and glimmerite (>90% phlogopite) (Villaseca et al., 2019a). Their size ranges from 1 to 7 cm. They show mostly medium-to fined-grained equigranular texture ( Figure 1B). Pyroxenites lack typical metamorphic textures, which is indicative of a recent magmatic origin before volcanic entrainment. Clinopyroxene can be colorless or pale green and appear mutually intergrowing as complex oscillatory zoned crystals ( Figure 1D), mantled by an external discordant thin purple rim, suggesting a genetic link with the two types of phenocryst cores described in the nephelinite from the Cerro Pelado cone. Amphibole shows slight zoning in the form of irregular patches. Occasionally, phlogopite may show kinked textures in the glimmeritic enclaves, likely associated to flowage and collisions with the conduit walls during ascent (Peterson and LeCheminant, 1993). The studied volcanic rocks also include a varied xenocryst population from the fragmentation of mantle xenoliths (e.g. Fo-rich olivine, orthopyroxene and Al-Mg-rich spinel).

WHOLE-ROCK GEOCHEMISTRY OF VOLCANIC HOST ROCKS AND PYROXENITE ENCLAVES
Although studied volcanic rocks have been classified as melilitites (El Aprisco) or nephelinites (Cerro Pelado) by petrographic means (Ancochea, 1982), their variable silica and alkali contents causes them to plot in the foidite, basanite and basalt fields when using the TAS diagram (Figure 2A). The complexity of sampling pyroclastic rocks representative of the transporting melt is evident, as they appear as mingled rocks comprising a mixture of crystal cargo (phenocrysts, fragmented antecrysts, xenoliths), juvenile material and accidental crustal lithics (see also Lierenfeld and Mattsson, 2015). Previous attempts to record the original composition of host magma that transported such a great quantity of solid cargo failed to reproduce the melilitic or nephelinitic original character of the volcanic host (Lierenfeld and Mattsson, 2015). In this work, only sample CAL-73 from the Cerro Pelado nephelinite volcano (sample provided by J.M. Cebriá) yields a consistent composition when compared to other nephelinites from the Calatrava Field (Figure 2; Table 1). Most crystal-rich samples from these volcanoes have low alkali contents and plot outside the typical compositional field of Calatrava volcanic rocks (Figure 2A). On the basis of the silica-undersaturated diagram of Le Bas (1989), studied rocks of the two volcanoes (including pyroxenites) plot in the nephelinite field ( Figure 2B). The higher MgO contents (>13.8 wt%, Table 1) of phenocryst-rich rocks of these two volcanic centers, with respect to other Calatrava volcanic rocks (mostly in the range of 8.4 to 12.8 wt%) are likely related to the relative abundance of mafic crystal cargo (clinopyroxene, olivine, amphibole and phlogopite). When using binary major element diagrams, most of the Calatrava volcanic rocks define a broad trend that might be induced by operation of minor fractional crystallization processes (Cebriá, 1992; Figure 2). In these diagrams, pyroxenites usually plot defining contrasted evolution lines or they are clearly outside of the main Calatrava compositional field (Figures 2C,D, respectively). Pyroxenite composition seems to be mainly controlled by mafic phases, especially by clinopyroxene or amphibole chemistry.
Major and trace element variation diagrams for pyroxenites and host-rocks display a marked data scatter, but in some diagrams, a clear distinction is drawn between pyroxenites and host magma compositional fields (Figures 2, 3). Pyroxenites have markedly low contents of Na, P, Nb, Ta, Th, U, Zr, Hf, LILE (except for K, Rb and Ba in phogopite-bearing types: 115601, 116598), Cr, Ni, and REE-Y (Figure 3). On the contrary, phlogopite-bearing pyroxenites are markedly rich in Al, K, Ti, Rb, Ba and Sc-V. Phlogopite-rich pyroxenites also display higher SiO 2 , Al 2 O 3 contents than the variety without mica (Table 1 and Figure 2C). Most of the chemical features of pyroxenites seem to be controlled by the combination of their mafic phases and Fe-Ti-rich minerals (Ti-magnetite, ilmenite, titanite). Pyroxenites in most trace element diagrams plot outside the compositional field of Calatrava volcanic rocks (Figure 3).
In chondrite normalized REE patterns all the pyroxenites show a similar pattern, with a mostly flat LREE pattern and a steeper shape from Eu to Lu (Figure 4). On the contrary, melilitites/nephelinites of the Calatrava Field (including sample CAL-73 from El Aprisco) exhibit a significant homogeneity in composition, displaying a strongly fractionated REE pattern (Figure 4). The analyzed pyroxenites show marked negative anomalies in Th, Nb (Ta) and P (and positive in Ti), contrasting with the patterns of the Calatrava volcanic rocks. Phlogopitebearing pyroxenites show slightly positive Rb, Ba and K anomalies (Figure 4). The CAL-73 sample and other Calatrava rocks show a similar pattern when compared to the Laterift volcanic rocks of the Neogene Bohemian Massif, with characteristic negative anomalies in K and Zr and positive anomalies in P (Ulrych et al., 2011).
The Sr and Nd isotopic composition of pyroxenites and the Calatrava volcanic rocks fall in the same range of 0.703270-0.703890 and 0.512790-0.512890, plotting close to the HIMU and FOZO components, and overlapping the EAR, common asthenospheric mantle reservoir of the circum-Mediterranean area ( Figure 5 and Table 2). This good fitting in isotopic composition suggests that the pyroxenite cumulates likely crystallized directly from magmas cogenetic with the host volcanic melts.  (Villaseca et al., 2019a). Two types of phenocryst cores appear in the Cerro Pelado volcanics: colorless and green clinopyroxenes. They show a contrasted composition, with the green cores exhibiting lower Mg#, TiO 2 and Al 2 O 3 , and higher FeO and Na 2 O contents than the colorless clinopyroxene (Villaseca et al., 2019a).

TRACE ELEMENT AND
The trace element composition of clinopyroxene megacrysts, plotted on chondrite-and primitive mantle-normalized diagrams, displays a relatively uniform pattern in both volcanic centers, although slightly higher contents can be found in megacrysts from the El Aprisco melilitite (Figure 6). They all present characteristic convex-upward LREE patterns, with strong depletion in the HREEs and negative anomalies in Rb-Ba, K, Pb and P (Figures 6A,B). The colorless cores of clinopyroxene phenocrysts exhibit an overlapping composition with respect to megacrysts (Supplementary Table S2 and Figure 6C). This feature suggests a genetic connection between both types of crystals. The trace element chemistry of these Calatrava clinopyroxene megacrysts-phenocrysts is very similar to that described in types II and III of megacrysts from the Eifel (Shaw and Eyzaguirre, 2000) but is clearly different to the xenocrystic cores of phenocrysts in alkaline basalts from the Carpathian-Pannonian region (Jankovics et al., 2016).
The green core phenocrysts show REE patterns parallel to those from colorless megacrysts and phenocrysts (Figure 6C), reinforcing the cogenetic character of all studied clinopyroxene types. Nevertheless, the green cores have higher REE, Zr-Hf and Nb-Ta and lower Ti-Cr-Ni contents than colorless clinopyroxenes (Figures 6C,D, 7). The chemistry of purple rims around green and colorless phenocryst cores, showing intermediate values in Mg# and Ni (Cr) (Figure 7), suggests that the host magma that transported the crystal cargo might have an intermediate chemical composition between those deep magmas that originated both phenocryst core types. Nevertheless, the wide REE and other incompatible trace element contents shown by purple clinopyroxene rims (Figures 6G,H) suggests a complex melt evolution of the host magma during volcanic emplacement.

Amphibole
The amphibole megacrysts found in the Cerro Pelado and El Aprisco volcanoes constitute a homogeneous group displaying a moderate chemical variation (Villaseca et al., 2019a). Amphibole Mg# may reach high values (up to 0.83), which is indicative of a fairly primitive composition and is positively correlated with Al 2 O 3 and K 2 O (Villaseca et al., 2019a). Patchy zoning and smaller amphibole crystals can be observed occasionally within FIGURE 3 | Trace element variation diagrams for the Calatrava Volcanic Field (CVF) rocks and associated pyroxenite enclaves. Same symbols as in Figure 2. The Calatrava compositional field does not overlap with most of the crystal-and xenolith-rich (contaminated) magmatic rocks from the El Aprisco and the Cerro Pelado volcanic centers.
megacrysts. Amphibole phenocrysts have never been found in the studied pyroclastic rocks of these two volcanoes.
The amphibole megacrysts exhibit a REE pattern similar to that of clinopyroxene megacrysts, characterized by a convexupward shape for the LREE and depleted HREE ( Figure 8A). Nevertheless, amphibole primitive mantle-normalized patterns are very different to those of clinopyroxene due to its higher LILE (K, Rb, Ba > Sr), Nb-Ta and Pb contents ( Figure 8B). The trace element chemistry of the Calatrava amphibole megacrysts (Supplementary Table S3) is very similar to that described in megacrysts from the Eifel (Shaw and Eyzaguirre, 2000) or from the Bohemian Massif (Mayer et al., 2014;Ulrych et al., 2018), as was the case of clinopyroxenes. Moreover, the D amp/cpx (REE,Th,U, Zr and Hf) estimated for both mineral megacrysts are in the range of 0.6-2, which is very similar to the values obtained from natural rocks and experimental data (Chazot et al., 1996;Grégoire et al., 2000;Raffone et al., 2009;Ubide et al., 2014), indicating the approach of chemical equilibrium between them.
The phlogopite megacryst has very low (subchondritic) REE-Y contents (Figure 8C). Its LILE (K, Rb, Ba) and Cr-Ni-(Ti) contents are higher than those of amphibole, as can be seen in multitrace primitive mantle-normalized diagrams ( Figure 8D). To our knowledge, there are no published trace element contents of phlogopite megacrysts from any Cenozoic volcano of the European Western circum-Mediterranean area, although this hydrous mafic mineral is an occasional megacryst in some of these volcanic centers (Hegner et al., 1995;Shaw and Eyzaguirre, 2000;Mayer et al., 2014). Only phlogopite megacryst of Cenozoic melilitites from the Ahaggar Swell (Algeria) has similar trace element patterns with slightly higher contents than the studied crystals (Kaczmarek et al., 2016). It is interesting to note that D Phl/Amp for Nb-Ta (0.5-0.6), Zr-Hf (0.06-0.12), Rb (16-20), Ba (8-9) and Sr (0.3-0.4) are broadly similar to those reported in other studies (Grégoire et al., 2000;Kaczmarek et al., 2016).
Regarding the isotopic composition, the four analyzed megacrysts (two amphiboles from both volcanic centers and the phlogopite and clinopyroxene from the Cerro Pelado volcano) ( Table 2) yield similar 87 Sr/ 86 Sr (0.703217-0.703480) and 143 Nd/ 144 Nd ratios (0.512832-0.512890) and plot within the wider isotopic field of the Calatrava volcanic rocks (Figure 5). These isotopic signatures suggest that the studied mafic megacrysts crystallized from melilititic-nephelinitic melts are similar in composition to their host magmas.

TRACE ELEMENT GEOCHEMISTRY OF PYROXENITE MINERALS
The composition of clinopyroxene from the pyroxenite enclaves overlaps the field of megacrysts, although defining a wider range for the REE contents, plotting toward either more enriched but also less differentiated compositions (Figures 6E-F). The two clinopyroxenites included in pyroclasts from the El Aprisco maar, which show colorless-greenish zoned crystals, have the lowest (sample 115591) and highest (sample 115601) REE contents of all the studied high-P clinopyroxenes ( Figure 6E).
Zoning profiles in these clinopyroxenes were previously characterized by the alternation of Si-Mg-and Fe-Al-Ti-Narich bands (Villaseca et al., 2019a). With respect to trace element contents, the green bands ( Figure 1D) have slightly higher P, REE, Y, Th, U, Nb and Ta contents than the colorless ones (Figure 9). This chemical contrast is smaller than that shown by green and colorless phenocryst cores, at least for some elements (e.g. Zr, Sc, Cr and Ni) (Figures 6, 7). In any case, the similarity in composition of the colorless phenocryst cores, megacrysts and most clinopyroxenes from the pyroxenite enclaves is remarkable. Similarly, some greenish clinopyroxene  Boynton (1985), whereas primitive mantle (PM) normalizing values are from Sun and McDonough (1989). The compositional field of clinopyroxene megacrysts (yellow field) is plotted in diagrams (C-H). from pyroxenite cumulates shows positive Zr (Hf) and Nb (Ta) anomalies in multi-trace element patterns (Figure 6F), although less marked than green phenocryst cores of sample 116596.
Amphiboles from the pyroxenite enclaves yield a more heterogeneous composition than megacrysts, although certain overlapping between both fields exists (Figure 8). Amphiboles from pyroxenite 115591 (El Aprisco) display the lowest Mg# values (up to 0.59), and something similar can be observed for REE, Y, Th, P and Ti (Figures 8A,B). On the contrary, amphiboles from this sample exhibit higher Pb and Zr (Hf) contents, similar to the accompanying greenish clinopyroxene, and consequently yield positive Pb and Zr (Hf) peaks in the multi-trace element pattern ( Figure 8B).
The trace element composition of phlogopite is highly variable (Supplementary Table S3) reaching extreme values in the clinopyroxenite and glimmerite enclaves, whereas phlogopite megacrysts and phenocrysts show an intermediate composition within that broader range but display similar normalized trace element patterns (Figures 8C,D). The mica from the glimmerite enclave has Mg# values similar to megacrysts and phenocrysts (averaging 0.83) and the highest REE contents, in contrast to phlogopite from pyroxenite enclave 115601, which shows the lowest REE and Th-U contents and the lowest Mg# (averaging 0.71) of all the studied micas (Figures 8C,D).

Origin of Megacrysts and Phenocrysts
It is not unusual that large mafic minerals (clinopyroxene, amphibole, phlogopite) were found within alkaline mafic to ultramafic rocks. The origin of these megacrysts is controversial and could be either xenocrystic (i.e. mantle-derived) or phenocrystic/antecrystic (i.e. melt-derived). In the former case, megacrysts could be considered fragments of mantle wall-rocks, whereas melt-derived megacrysts could be either (i) pieces of mantellic cumulates unrelated to the host melt or (ii) minerals crystallized at high pressure conditions from their host volcanic magma or a cogenetic melt fraction (e.g. Shaw and Eyzaguirre, 2000;Gernon et al., 2016). The textural features and major element chemistry of the studied mafic megacryst-phenocrysts strongly support that the most primitive minerals were close to equilibrium and crystallized from primitive melts with a composition similar to that of the El Aprisco and Cerro Pelado volcanic rocks (Villaseca et al., 2019a). Nevertheless, the existence of different types of phenocryst cores with a variable composition, and the presence of reaction textures and mineral zoning indicates a more complex scenario.
A xenocrystic origin of the megacrysts is not favored on the basis of several petrological and geochemical features. The crystal size (a few millimeters) and the interstitial character of the metasomatic hydrous phases in peridotites (amphibole and phlogopite) do not support the possibility of the large mafic crystals being fragments from these mantle rocks. The major element mineral composition is also different in the peridotites and the megacrysts and phenocryst cores: mainly, the lower Mg# and Cr and higher Ti values shown by the latter minerals (Villaseca et al., 2019a). Mineral trace element data reinforces this conclusion, since a clear disagreement exists between the geochemistry of the lherzolite minerals and that of megacrysts. These discrepancies are apparent for most elements in the El Aprisco volcano (clinopyroxene and amphibole composition;  Table S3. Chondrite and primitive mantle normalizing values are from Boynton (1985) and Sun and McDonough (1989), respectively. Figures 10C,D, 11E,F) and also for REE, Th, U, Nb and Ti of lherzolite clinopyroxene from the Cerro Pelado scoria cone (Figures 10A,B). Nevertheless, some similarities can be found for both major and trace elements only when comparing with the heterogeneous composition of clinopyroxene from wehrlite xenoliths from the Cerro Pelado volcano. However, this similitude is likely associated with the metasomatic origin of these wehrlites, since its metasomatic signature is genetically related to reaction with the host nephelinitic magma (Villaseca et al., 2010). Accordingly, the metasomatic amphibole and phlogopite of the Cerro Pelado wehrlite xenoliths display similar trace element patterns to those shown by the associated megacrysts, which are high-pressure fractionates derived from similar melts (Figures 11A-D). Nevertheless, this geochemical resemblance does not necessarily imply that megacrysts represent wehrlite fragments, which is not supported by petrographic features. Moreover, green clinopyroxene phenocryst cores have chemical features clearly different to wehrlitic clinopyroxene, including their notorious Zr positive peak in multitrace element patterns (Figure 6), also indicating a non-xenocrystic (mantle fragmentation) origin for these deep-seated crystals.
It has been long recognized that convex-upward REE patterns, such as those depicted by the Calatrava clinopyroxene and amphibole megacrysts, are typical of minerals crystallized from alkaline mafic magmas at high-pressure (Irving and Frey, 1984). Colorless clinopyroxene phenocryst cores display a chemical composition overlapping that of the equivalent megacrysts. This feature, which has also been observed for the major elements (Villaseca et al., 2019a), is apparent for the trace elements (Figures 6C,D) and implies that they all must share a common origin as deep mineral fractionates. Studied clinopyroxene and amphibole megacrysts are in chemical equilibrium, as their intermineral partition coefficients (see comments above) are close to those previously published for basic-ultrabasic lithotypes (Vannucci et al., 1995;Raffone et al., 2009). This suggests that they are cogenetic phases derived from the same magma. This is confirmed by the similar composition of the liquids in equilibrium with studied clinopyroxene and amphibole calculated using mineral/melt distribution coefficients of high-T and high-P basaltic systems (Hart and Dunn, 1993;La Tourrette et al., 1995;Ionov et al., 1997; Figure 12). The results compare well with the whole-rock data of the Calatrava host magmas, with REE patterns of similar slope and equivalent positive and negative trace element anomalies (Figure 12). Nevertheless, the higher chemical heterogeneity shown by the calculated melts with megacryst/phenocryst clinopyroxene, mainly for the LIL FIGURE 9 | Chemical profiles showing compositional variations in a greenish-to-colorless zoned clinopyroxene crystal of pyroxenite 115601. The Na 2 O vs. distance plot is taken from Villaseca et al. (2019a). Green clinopyroxene area is marked, and laser spots are depicted in Figure 1D of this work. Note the high Na 2 O, REE and HFSE contents in the green bands relative to the surrounding (colorless) clinopyroxene. elements (Rb, Ba, Th, U), is likely associated with the variable composition of these mafic minerals and the evolution of magma at deep stagnation levels. Even so, the broad similarity between the calculated melts and the Calatrava volcanics indicates that mafic megacrysts derive from high-pressure fractional crystallization of melts with a composition similar to that of the host alkaline magmas.
On the other hand, the Fe-rich green clinopyroxene phenocryst cores are slightly trace element-enriched (remarkably for Zr), and some purple clinopyroxene rims are clearly different from other clinopyroxenes due to their elevated trace element concentrations (Figures 6G,H). The calculated geochemical heterogeneity of equilibrated melts with these clinopyroxene types (Figures 12A-D) supports the involvement of magmas with a variable degree of evolution for their origin, within a complex scenario of crystallization at different levels in the upper lithosphere (see discussion below).

Origin of Pyroxenite and Glimmerite Enclaves
The Calatrava clinopyroxenite and glimmerite enclaves do not have characteristics of xenolithic mantle rocks. While minerals in peridotites and metasomatic veins are usually unzoned and very homogeneous in major element composition (Villaseca et al., 2010;Lustrino et al., 2016), clinopyroxene and amphibole of the studied pyroxenites display oscillatory or patchy zoning, (respectively), commonly interpreted as igneous textures. The absence of exsolution lamellae in clinopyroxene from these enclaves also contradicts their origin as wall-rock mantle fragments entrapped by the ascending basic magma (Villaseca et al., 2019a).
The mineral geochemistry also favors the magmatic origin of the Calatrava pyroxenite and glimmerite enclaves, due to their low Mg# values when compared to those of the mantle minerals (<0.85; Villaseca et al., 2019a). The trace element composition of the main mafic minerals of these enclaves (colorless clinopyroxene, amphibole and phlogopite) is coincident with that of the megacrysts and colorless phenocryst cores dispersed within the host volcanic magma (Figures 6,  8), which implies that they crystallized from similar cogenetic basic melts. Nevertheless, liquids in equilibrium with greenish clinopyroxene from these enclaves show REE patterns parallel to those from the Calatrava volcanics but showing different REE contents ( Figure 12C) or having spikier trace element patterns with a marked Zr positive anomaly that is absent in Calatrava volcanic rocks (Figure 12D), as was the case with the data from green clinopyroxene phenocryst cores.

Relationships Between Magmatism and Metasomatism Within Peridotite Wall-Rocks
As shown above, we have observed similarities among the trace element compositions of clinopyroxene, amphibole and phlogopite megacrysts in the Cerro Pelado center and the equivalent metasomatic phases in associated peridotite xenoliths, mostly in the case of the "Fe-Ti metasomatised" wehrlites. These  Boynton (1985) and Sun and McDonough (1989), respectively. latter enclaves have neoblastic spongy clinopyroxene (samples 72674 and 116648) and interstitial phlogopite, amphibole and alkaline glass. Taking into account that the studied megacrystsphenocrysts crystallized from magmas similar in composition to their present host melts, the compositional similarities might indicate that the metasomatic agents which reacted with the Cerro Pelado peridotites and the host alkaline mafic magmas were genetically related. In fact, previous studies in those mantle xenoliths already established that they had undergone mafic silicate metasomatism, and that the most metasomatised wehrlites were probably connected with reactive infiltration of alkaline melts similar to their host magmas (Villaseca et al., 2010). Thus, although some interstitial amphibole and phlogopite in lherzolites record reaction with the Calatrava alkaline magma, these hydrous minerals would rarely give rise to euhedral cm-sized megacrysts, and their much higher Mg and Cr contents make it easy to discriminate magmatic megacrysts of melt-mantle metasomatic minerals from strongly transformed mantle xenoliths (Yu et al., 2019). Neither amphibole-(or phlogopite-) rich pegmatite-like veins nor large idiomorphic crystals (megacrysts) have been recorded in the studied peridotite suites.
Mantle xenoliths from the El Aprisco maar have metasomatic phases remarkably different in composition to the studied megacrysts and phenocrysts. The study of these mantle xenoliths has enabled the identification of at least three old metasomatic events occurring during the Cretaceous (∼118 Ma), the Oligocene (∼29 Ma) and the Miocene (∼16-4 Ma) (Villaseca et al., 2019b). As in the case of the Cerro Pelado peridotites, the high Cr, Ni and Sc contents of metasomatic amphiboles from the El Aprisco peridotite suites also agree with a longer time of equilibration within the mantle in comparison to the much younger age of mafic minerals formed during the Calatrava volcanic event (2.8 to 2.2 Ma ages for megacrysts formed in both volcanic centers, after Villaseca et al., 2019b), even if their formation occurred at similar P-T conditions at mantle depths.

Melting, Storage and Melt Evolution Within a Magmatic Plumbing System
The main alkaline silica-undersaturated magmas of the Calatrava Volcanic Field have been interpreted as being generated from a mantle asthenosphere source with trace element and isotopic ratios similar to the HIMU (or EAR) OIB-reservoir FIGURE 11 | Chondrite-normalized REE and primitive mantle-normalized trace element patterns of amphibole (A,B,E,F) and phlogopite (C,D) of peridotite xenoliths from studied volcanic centers. Amphibole trace element composition of mantle xenoliths (G1 to G3 groups) from the El Aprisco maar is taken from Villaseca et al. (2010Villaseca et al. ( , 2019b and González-Jiménez et al. (2014). Those from the Cerro Pelado scoria cone are unpublished data. Amphibole megacryst field (in yellow) and the average (black line) is taken from Figure 8. Chondrite values are from Boynton (1985) whereas primitive mantle normalizing values are from Sun and McDonough (1989). (Cebriá and López Ruiz, 1995). Phlogopite and garnet are considered residual phases in the mantle sources, consistent with the fractionated HREE patterns and negative K anomalies shown by Calatrava volcanic rocks (Figure 4). The determination of accurate P-T conditions of mantle melting is challenging. We have considered the geothermobarometry approach of Lee et al. (2009), which estimates the P-T condition of basaltic melt extraction in equilibrium with an averaged lherzolitic peridotite source. The high MgO content (>10 wt%) of most of the Calatrava magmas (   Boynton (1985) and Sun and McDonough (1989), respectively. pyroxenite enclaves have been excluded due to their high crystal cargo or cumulate origin, respectively. The resulting pressure and temperature of melt segregation from the mantle source are in the range of 2.3 to 5.6 GPa and 1,300 to 1,540 • C, respectively ( Table 1). These values have to be considered as rough estimates of partial melting conditions for Calatrava magmas from a lherzolite protolith. It is interesting to note that these estimations imply melt extraction within the stability field of garnet, in agreement with the chemical features of these alkaline magmas.
During their ascent, some of the Calatrava magmas were arrested in shallow lithospheric mantle reservoirs, where fractional crystallization was triggered. This is illustrated by the formation of a suite of mafic megacryst-phenocrysts and pyroxenitic cumulates, such as those described in this study. The conditions of crystallization indicate pressures of 1.2-1.6 GPa and temperatures in the range 1,275-1,350 • C, corresponding to the uppermost mantle (Villaseca et al., 2019a). The relative wide depth range for magma stagnation derived from the pressure data could imply either a polybaric crystallization during melt ascent or uncertainties in the P-T calculations. In any case, the much lower pressure provided by clinopyroxene phenocryst purple rims (∼0.3-0.8 GPa, Villaseca et al., 2019a) points to the formation of these overgrowths during the magma transport through the crust.
The presence of clinopyroxene phenocrysts of contrasted composition (colorless and green cores), sharing similar purple rims, reveals a complex magmatic story. The green cores are enriched in Fe and incompatible trace elements when compared to the colorless cores, and the corresponding melts in equilibrium must represent a more evolved composition (Figure 13). This feature is indicative of a magma differentiation process occurring at depth and responsible for the formation of a heterogeneous suite of mafic minerals with a variable degree of evolution. Clinopyroxene and amphibole might be the main liquidus minerals during magma crystallization at depth due to their modal abundance as megacrysts and mineral phases in pyroxenite enclaves. Phlogopite and olivine (and accessory Fe-Ti-rich minerals) were also involved in minor proportions. As an approximation of this process, we have applied a model of fractional crystallization starting from a melt C 0 representing a primitive host basaltic composition (see caption to Figure 13) and assuming clinopyroxene + amphibole + minor olivine FIGURE 13 | Gd/Yb vs Zr/Hf ratios of melts equilibrated with clinopyroxenes. The compositional field of Calatrava Volcanic Field (CVF) magmas is also shown for comparison. The clinopyroxene compositions used for calculations are averaged values for selected colorless megacrysts (115603,114404,111639) and phenocryst cores (115605,116596), together with green phenocryst cores (116596) and pyroxenites (115601,115591). Purple rims overgrow phenocrysts (115605, 116596) and pyroxenite crystals (115591) (see Supplementary Table S2 for averaged clinopyroxene values). The C 0 for the fractional crystallization model is the CAL-73 nephelinite (from the Cerro Pelado volcano, CP in figure). The arrow represents the trend of chemical evolution of melts based on the fractionation model (up to 80% of fractional crystallization) considering extraction of 45% clinopyroxene + 45% amphibole + 10% olivine.
fractionation. The distribution coefficients used are those of La Tourrette et al. (1995) (Zr-Hf) and Irving and Frey (1984) (REE) for amphibole; Hart and Dunn (1993) for clinopyroxene and Adam and Green (2006) for olivine. The results indicate that a high degree of fractionation is needed (up to 80%) in order to obtain the high Zr/Hf and the slightly low Gd/Yb ratios of the evolved melts equilibrated with green clinopyroxene (Figure 13). The presence of mica-rich (glimmerite) enclaves is an evidence in favor of phlogopite crystallization during differentiation, but its participation does not change significantly the above estimated evolution pattern.
A fractional crystallization process also explains the clinopyroxene zoning in pyroxenite enclaves, with variable Mg# (around 0.65 in the green bands and 0.82 in the colorless core) and trace element contents within the same crystal (Figure 9). The fact that this zoning is frequently oscillatory requires a connection with a more primitive magma batch acting as a recharge source. The geochemical similarities shared by the evolved green phenocrysts cores and the pale green clinopyroxenes in pyroxenites supports their formation in a similar context: magma pools stagnated within the mantle. We think that magmas slightly more primitive than the volcanic host melts evolved by fractional crystallization within mantle reservoirs, generating a wide range of liquid compositions (either more primitive or significantly more fractionated) at local scale. This large chemical heterogeneity contrasts with the coincident isotopic composition of megacrysts, pyroxenite enclaves and host rocks (Figure 5), which is indicative of the cogenetic nature of the involved magmas. This model contrast with previous works studying clinopyroxene green phenocryst cores in alkaline volcanoes of the European circum-Mediterranean region, where their origin from significantly more evolved melts than host lavas was interpreted as xenocrystic fragments from noncogenetic rocks (Aulinas et al., 2013;Jankovics et al., 2013;Matusiak-Malek et al., 2018).
A typical feature of chilled margins in cooling magmatic bodies is the presence of interstitial highly fractionated liquids within crystal mush carapaces (e.g. Marsh, 1996). The pyroxenitic cumulates found in the Calatrava volcanic rocks can be considered pieces of this crystalline margin, and the complex zoned crystals they include serve as evidence of the involvement of highly evolved melts. These latter liquids could give rise to single microcrysts (green phenocryst cores), subsequently incorporated as crystal cargo in ascending and less fractionated magmas from the same magmatic system. Nevertheless, there are several differences between the trace element composition of melts in equilibrium with green clinopyroxene and the most evolved Calatrava volcanic rocks, such as the marked Zr content of green phenocryst cores (Figures 12, 13). The highly evolved liquids which gave rise to the low-Mg green clinopyroxene have not been erupted in Calatrava volcanoes, and were probably trapped within the cumulate chilled margins of deep mantle reservoirs. The small volume of highly fractionated liquids generated hinders to collect them to be transported toward the Earth's surface; they are relatively immobile.
The major and trace element composition of the clinopyroxene purple rims is broadly intermediate between primitive colorless and evolved green cores (Figures 8, 13; see also Villaseca et al., 2019a). Purple rims can be found overgrowing both types of phenocryst cores and also forming the most external margin of clinopyroxene in some marginal areas of pyroxenite enclaves. It is thus likely that rims were formed from the magma that entrapped megacrysts, phenocrysts and enclaves, and erupted at the Earth's surface. This conclusion is supported by the broadly similar composition of calculated melts in equilibrium with thin purple rims and the host melts (Figure 13). Other rims provide equilibrated melts with trace element contents significantly higher than the Calatrava magmas (Figures 12, 13). This feature is in accordance with an increase in the cooling and decompression rates, during the rise of magma toward the surface, which favors the enrichment in incompatible trace elements as a consequence of modification of the partitioning of elements and crystallization being governed by kinetic effects (e.g. Re et al., 2017 and references therein).
The crystal fractionation process described in this study adds novel information regarding the magmatic differentiation associated with deep plumbing systems beneath monogenetic volcanos. Mafic mineral crystallization is one of the main factors controlling magma differentiation inducing changes in magma buoyancy, volatile composition and, eventually, the onset of volcanic eruption (Mattsson et al., 2013;Re et al., 2017). The cumulate origin of clinopyroxene, amphibole and phlogopite, assessed via trace element modeling, implies that the primary magmas were volatile-rich and that a significant proportion of minerals was removed from the melt due to crystal fractionation at mantle reservoirs. It is commonly accepted that volcanic eruptions of carbonate-rich nephelinitic and melilititic magmas generating large maar-diatreme volcanoes can be triggered by pressure oversaturation associated with deep exsolution of CO 2 due to silicate crystallization (e.g. Mattsson et al., 2013). This process is favored by decompression during magma ascent or due to recharge of deep juvenile volatile-rich melts into the magma reservoirs (Ubide et al., 2019 and references therein). Although such discussion is beyond the scope of this work, our data leaves open the possibility of the CO 2 increase being related with deep-seated crystal fractionation (Villaseca et al., 2019a). The combination of CO 2 increase in the environs of the solidification front and upward movement of magma within the upper mantle could explain magma boiling, fragmentation and a rapid ascent to the surface, carrying a heterogeneous suite of crystals, pyroxenitic enclaves and mantle xenoliths.

CONCLUSION
This is the first study focused on the characterisation and interpretation of the trace element mineral chemistry of phenocrysts, megacrysts, and clinopyroxenite and glimmerite enclaves carried by alkaline rocks of the Calatrava volcanic field (El Aprisco and Cerro Pelado volcanoes), central Spain. Mineral data were complemented with whole-rock and isotopic (Sr-Nd) data of megacrysts, enclaves and volcanic rocks. Mineral analyses focused on clinopyroxene, amphibole and phlogopite, which appear as megacrysts and phenocrysts, and as main constituents of the pyroxenite enclaves. The mafic megacrysts display a chemical composition more akin to an origin as deep segregates rather than fragments of a peridotitic mantle. They show clear differences in composition when compared to minerals from the mantle xenoliths carried by the Calatrava volcanics but are similar to the same phases in clinopyroxenite and glimmerite enclaves.
The presence of disequilibrium textures, complex zoning and contrasting compositions in clinopyroxene crystals indicates that they did not crystallize from a single melt. Clinopyroxene is the most common and widespread phase, occurring as megacrysts, macrocrysts and phenocrysts in these volcanic rocks. Two types of high-P crystals of different colors have been observed: green and colorless phenocryst cores. The former are less abundant and have low Mg# (Cr, Ni, Sc) and high Zr-Sr contents, which is indicative of their derivation from melts much more fractionated than those related to the primitive high Mg# colorless cores. The contrasting compositions also coexists in complexly zoned crystals of clinopyroxenites, suggesting that the colorless and green types might be cognate and derive from mixing of primitive and evolved fractionated liquids within the same magmatic system. The geochemistry of clinopyroxene supports a fractionation process mainly controlled by clinopyroxene and amphibole crystallization.
The similarity in mineral chemistry and the homogeneous whole-rock Sr-Nd ratios of all the analyzed minerals and rocks suggest that the deep fractionating magma batches and the melts transporting this complex crystal cargo were likely generated from similar sources. The possibility that this suite of fractionated magmas derive from a single magma batch is supported by the overlapping composition of the melts in equilibrium with the mafic megacrysts/phenocrysts and the volcanic host magmas and also may explain the presence of oscillatory zoning at micronscale in some cumulates, linked to crystallization fronts in deep magma chambers, where repeated infiltration of recharging external liquids is probably limited. This model contrasts with previous studies that suggested that clinopyroxene green phenocryst cores in alkaline melts are xenocrystic fragments from non-cogenetic rocks.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/Supplementary Material.

AUTHOR CONTRIBUTIONS
All authors collaborated in the sampling, petrographic study and the performing of analytical data and contributed to the final version of the manuscript, figures, and tables. CV conceived the study and wrote an original draft later improved with inputs from JG and DO.

FUNDING
This work is included in the objectives of, and supported by, the CGL2016-78796 project of the Ministerio de Economía y Competitividad de España, and the 910492 UCM project.

ACKNOWLEDGMENTS
We acknowledge Alfredo Fernández Larios for his assistance with the electron microprobe in the Centro Nacional de Microscopía Electrónica Luis Bru (UCM). Manuel Alpiste is thanked for his assistance with the LA-ICP-MS of the Instituto Andaluz de Ciencias de la Tierra (Granada). We also greatly appreciate the comments made by PS, MG, and the associated editor Teresa Ubide on a previous version of the manuscript.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart. 2020.00132/full#supplementary-material FIGURE S1 | (A) The Calatrava Volcanic Field and location of the two studied volcanoes: El Aprisco and Cerro Pelado (Ancochea, 1982;Cebriá et al., 2011). (B) Sketch map of the Iberian Peninsula showing the location of the main Cenozoic volcanic fields (in black; Lustrino and Wilson, 2007).
TABLE S1 | Compositional data on certified standards analyzed with the samples in this study.