Mononuclear η⁶-arene ruthenium(II) complexes with pyrazolyl-pyridazine ligands: synthesis, CT-DNA binding, reactivity towards glutathione, and cytotoxicity.

Contributions to Mineralogy and Petrology (2023) 178:63 https://doi.org/10.1007/s00410-023-02043-7 ORIGINAL PAPER Ultra‑crystalline pyroclastic deposits and rhyolitic lavas controlled by crystal mushes: insights from the Acoculco Caldera Complex, México Mario E. Boijseauneau‑López1 Benjamin J. Andrews3 · Giovanni Sosa‑Ceballos2 · Gabriela A. Farfán3 · José Luis Macías2 · Received: 2 June 2022 / Accepted: 31 July 2023 / Published online: 20 August 2023 © The Author(s) 2023 Abstract The Acoculco Caldera Complex (ACC), located in eastern Mexico, began its activity during the Pleistocene ~ 2.7 Ma. One of the most relevant and largest rhyolitic eruption in the complex, the Piedras Encimadas Ignimbrite (PEI), occurred during the late post-caldera phase at ~ 1.2 Ma. This ignimbrite is unique with respect to the other caldera products and other contemporaneous ignimbrites in the Trans-Mexican Volcanic Belt (TMVB) because of its ultra-high crystallinity and the absence of pumice fragments. The PEI is made almost entirely of crystals where the main constituents are k-feldspars and silica polymorphs that range from ≤ 5 µm to tens of centimeters in size. XRD on bulk rock, geochemical modeling, FTIR, Raman, and EPMA analyses were carried out in all mineral phases to assess the origin and the causes of high crystallinity within the PEI. We interpret the high crystallinity on the basis of magmatic crystallization of a magma body that was remobilized and altered by post-depositional hydrothermal alteration processes. We suggest that ACC rhyolites are geochemically influenced by at least one crystal mush established during the Pleistocene. We suggest that the PEI could be the result of an erupted crystal mush (melt + crystals), or a cumulate, or an ancient and crystallized reservoir generated after the first ACC collapse due to intrusion or underplating of mafic hot magmas. Extensional episodes within the ACC facilitated the ascent of mafic magmas. This interaction increased the liquid fraction of the mush through partial melting/crystal dissolution, generating a drop in density and viscosity in the mush, thus triggering eruption. The PEI provides evidences for an association between the geochemically-diverse ACC rhyolites with the complex interaction between mafic transitional alkaline magmas and a crustal mush system, promoted by continuous changes in the stress field during the Pleistocene. Keywords Crystal-Mush · Acoculco Caldera Complex · Piedras Encimadas Ignimbrite · Trydimite · Cristobalite · Hydrothermal Alteration · Geothermal Exploration Introduction Communicated by Gordon Moore. * Mario E. Boijseauneau‑López cswasa@gmail.com 1 Posgrado en Ciencias de la Tierra, Escuela Nacional de Estudios Superiores, UNAM, Campus Morelia, 58190 Morelia Michoacán, México 2 Institute of Geophysics, UNAM, Campus Morelia, IGUM, 58190 Morelia Michoacán, Mexico 3 National Museum of Natural History, Smithsonian Institution, Washington D.C., USA The origin and evolution of rhyolitic magmas have been studied for decades (e.g. Whitney et al. 1989; Anderson et al. 2000; Bindeman and Valley 2003; Hildreth and Wilson 2007; Bindeman 2008; Scaillet et al. 2016; Jolles and Lange 2021; Clemens et al. 2022; Troch et al. 2022). Some of the largest (VEI > 6) ignimbrites are typically associated with rhyolitic magmas, but generation of its evolved magmas (regardless of the process, fractional crystallization or partial melting) involves segregation of meltrich, crystal-poor magma from a crystal-rich source—the so-called crystal mush (Bachmann and Bergantz 2003, 2008a, b; Bachmann and Bergantz 2004a; Cashman et al. 2017). However, the depths of magma production, 13 Vol.:(0123456789) 63 Page 2 of 24 extraction, final storage, and trigger mechanism of these eruptions remain controversial (Bachmann et al. 2002; Gottsmann et al. 2009; Malfait et al. 2014). It is widely accepted that the segregation of melt-rich forms aphyric rhyolitic magmas, crystal-poor magma from a crystal-rich source (Gualda et al. 2019). However, previous studies rarely have investigated if crystal-rich rhyolitic ignimbrites could be associated with the remobilization and reactivation of melt-dominant bodies that had cooled and crystallized close to the point of critical crystallinity or rheological lock-up at >  ~ 50 vol.% of crystals (Marsh 1981; Brophy 1991; Vigneresse et al. 1996; Huber et al. 2010; Petford 2003); that is, the mobilization of a crystal mush including melt + crystals. Only, crystal-rich monotonous intermediates and aphyric rhyolitic ignimbrites have been stated as result of the activity of a crystal mushes (e.g., Hildreth 1981; Bachmann and Bergantz 2004a and 2008b; Hildreth and Wilson 2007; Blundy and Cashman 2008; Cashman and Giordano 2014; Zhang et al. 2018; Giordano and Caricchi 2022). Even more difficult to interpret are the cases of volcanic complexes with different ignimbrite types (crystal rich and crystal poor) and a series of rhyolitic lavas with no clear relationship to a crystal mush nor to a classical ignimbrite generation model. Crystalrich, pumice-free ignimbrites with > 70 vol% crystals are sparsely described in the literature. This type of rhyolitic ignimbrites already has been interpreted as part of caldera eruptions related with a crystal mush; however, their origins are not fully understood (e.g., Sohn et al. 2009; Zhang et al. 2018). The Acoculco Caldera Complex (ACC) in eastern Mexico (Fig. 1a) is an excellent case study because it hosts a great variety of rhyolites, comprising a syn-collapse ignimbrite, a series of annular rhyolitic domes (located in the rim of the caldera), and two post-collapse ignimbrites related to large eruptions triggered by local tectonics (Avellán et al. 2019, 2020). Furthermore, the ACC has pervasive hydrothermal alteration features distributed inside and outside the caldera, fumaroles, and hot spring waters. It is well documented how hydrothermalism generates alteration of the rocks in which the fluid circulation modifies the natural paragenesis and crystallinity through mineral chemical and structural modifications, as well as the precipitation of new mineral phases or the replacement of others (i.e., Schwartz 1959; Browne 1978; Giggenbach 1984; Kimball 1990; Mathieu 2018). The Piedras Encimadas Ignimbrite (PEI), the biggest felsic eruption reported in the ACC, is the first crystal-rich rhyolitic volcanic deposit described in Mexico. Due to its high crystal contents, we will examine if the hydrothermal alteration reported in the ACC changes the crystallinity or if the minerals are product of magmatic crystallization. Moreover, this study has broader applications for the understanding 13 Contributions to Mineralogy and Petrology (2023) 178:63 of how crystal-rich deposits are produced and may help to determine the conditions in which a crystal mush model is applicable. Here we characterize the geochemical, mineralogical, and trace element record of the PEI post-collapse ignimbrite and other felsic products of the ACC (Fig. 1b). We also propose a new model to understand the evolution of the ACC plumbing system and magmatism. This model links the different erupted compositions to a coherent differentiation scheme and is compatible with a process of crystal mush remobilization that controls the compositional variation within the ACC. Geologic setting The ACC is located on the eastern part of the Trans-Mexican Volcanic Belt (TMBV) (Fig. 1a). This volcanic arc is the product of the subduction of the Cocos and Rivera oceanic plates underneath the North American continental plate (Pardo and Suarez 1995; Gómez-Tuena et al. 2018). The local basement beneath the ACC region is formed by Cretaceous limestones (Avellán et al. 2019; López-Hernández et al. 2009), metamorphosed by a series of granitic intrusions of unknown age. Our knowledge of the local basement is the result of two exploratory wells performed by the Comisión Federal de Electricidad (CFE, the National Electricity Company). In addition to the metamorphic and sedimentary rocks, there are a series of aplitic dikes-sills (183 ± 36 ka, Avellán et al. 2019) and lava flows with basaltic andesitic and andesitic compositions (García-Palomo et al. 2018). Regional changes in stress orientations during the Miocene created a normal fault system with a NW–SE trend (García-Palomo et al. 2018; Gómez-Alvarez et al. 2021). This deformation in the southernmost part of the ACC zone developed to form the Tlaxco NW–SE fault system (GómezAlvarez et al. 2021). During the Pliocene–Pleistocene, the local stress changed, generating new normal fault systems in the Apan-Acoculco area with NW–SE orientations (GarcíaPalomo et al. 2018). This structural system dominated until the late Pleistocene when the fault system orientation changed to NE–SW, producing graben-horst structures such as the Apan, Tlaxco-Chignauapan grabens, and the RosarioAcoculco horst (García-Palomo et al. 2002, 2018). Volcanic activity at the ACC started at 2.7 Ma and has been continuous until 0.016 Ma (Avellán et al. 2020). The volcanic evolution of the ACC was characterized by Avellán et al. (2019) into the categories syn-caldera, early postcaldera, late post-caldera, and extra-caldera units. The precaldera volcanism of the ACC comprises andesitic-dacitic domes; the syn-caldera stage generated the Acoculco andesitic ignimbrite (~ 2.7 Ma) and the first collapse of the region (Avellán et al. 2019). The collapse modified the local Contributions to Mineralogy and Petrology (2023) 178:63 Page 3 of 24 63 Fig. 1  a Location of the Acoculco Caldera Complex (ACC) within the Trans-Mexican Volcanic Belt (TMVB). b Simplified geological map of the ACC, after Sosa-Ceballos et al. (2018) and Avellán et al. (2019). The map contains the distribution of the basement rocks as limestones (Ksl), Peñuela and Quexnol (Mv), and the volcanic units; pre-caldera volcanism (Pc), syn-caldera volcanism (Sc), early post-caldera volcanism (Epc), late post-caldera volcanism (Lpc), and Extra-caldera volcanism (ATVF). The red lines represent the faults in the area after Calcagno et al. (2019) that form normal faults (NE– SW) and graben–horst structures. The localization of PEI samples is showed by a black mark on the map (modified from Pérez-Orózco et al. 2021) stress field and promoted the ascent of peralkaline rhyolites mixed with calc-alkaline magmas (Sosa-Ceballos et al. 2018; Avellán et al. 2019). The post-caldera stage is characterized by bimodal volcanism: first with voluminous mafic rocks of the early post-caldera volcanism (2.6–2.2 Ma) followed by the late post-caldera volcanism (2–1 Ma), which comprised a series of rhyolitic domes, the PEI (Fig. 1b), Tecoloquillo ignimbrite (Avellán et al. 2020), and, finally, with a series of minor Quaternary monogenetic eruptions of mafic composition. The most voluminous episodes of rhyolitic volcanism within the ACC area are represented by the PEI (Fig. 1b) and the Tecoloquillo ignimbrite. PEI is a welded ignimbrite with several flow units that appear as massive, light-gray to white, beds (Avellán et al. 2019, 2020). Each bed consists of microcrystalline matrix-supported particles with feldspar and quartz phenocrysts. It has an approximated volume of 26 ­km3 (Avellán et al. 2020). The PEI represents caldera reactivation at 1.2 Ma, and its volcanic center was probably at the eastern portion of the ACC ring structure (López-Hernández et al. 2009; Avellán et al. 2020). The Later, another explosive 13 63 Page 4 of 24 eruption occurred at ∼ 0.8 Ma at the south part of the caldera producing the Tecoloquillo rhyolitic ignimbrite and its summit dome (Avellán et al. 2019). Pérez-Orozco et al. (2021) reveal that Tecoloquillo products show evidence of magma mixing between calc-alkaline and peralkaline melts. Furthermore, melt inclusions suggest that two magma bodies, located approximately 6–10 km in-depth and, isolated from one another, sourced the eruption (Pérez-Orozco et al. 2021). Older felsic products in the region are represented by the Peñuela Dacitic Dome Complex (~ 13–10 Ma) and Terrerillos Domes (~ 3 Ma) (García-Palomo et al. 2002; Avellán et al. 2019, 2020). It is not clear if magmatism is continuous between the Peñuelas and Terrerillos Dome complexes; however, felsic products are present in the late post-caldera units (~ 2–0.016 Ma) forming lava flows and domes with predominant rhyolitic compositions (Avellán et al. 2019, 2020). In addition, during the formation of the ACC, mafic volcanism (2.4–0.19 Ma) related to the Apan-Tezontepec Volcanic Field (ATVF) is interbedded with the products of the ACC and correspond to the extra-caldera stage (GarcíaPalomo et al. 2002, 2018; López-Hernández et al. 2009; Avellán et al. 2019). The interest in exploring the ACC region arose in the 1980s with prospective works carried out by the Mexican Federal Electricity Commission (CFE, for its acronym in Spanish) to determine its geothermal potential. Despite the hydrothermal alteration being the most outstanding feature of most ACC deposits (Canet et al. 2015), there are few spots of superficial geothermal manifestations such as springs and fumaroles (Lopez-Hernandez and Castillo-Hernandez 1997). The alteration features are reported for an area of several square km in the center of the caldera complex, and outside, in the eastern border, affecting pyroclastic deposits, and, to a lesser degree, lavas (Canet et al. 2015; Sánchez-Córdova et al. 2020; Pandarinath et al. 2020). Inside the caldera, the most common hydrothermal alteration is the silicic type. There are some regions where acid-sulfate alteration types can be found (vuggy silica and advanced argillic alterations); however, outside of the caldera, the hydrothermal alteration has no studies. Methods We collected bulk rock samples (comprising crystals and matrix) from the bottom, middle, and top portions. We study four sections of the PEI in proximal and distant deposits in gullies of the western caldera borders (Fig. 1b). The field strategy consisted of sampling the deposit based on vertical variations. Samples were crushed and hand-picked to obtain lithic-free fragments for analysis. Whole-rock analysis of major and trace element concentrations were determined by fusion inductively coupled plasma emission spectroscopy 13 Contributions to Mineralogy and Petrology (2023) 178:63 (ICP-ES) and inductively coupled plasma mass spectrometry (ICP-MS) at the Activation Laboratories, Ancaster, Canada (https://​actla​bs.​com/). Fifteen rock samples were cut and processed into thin polished sections for petrographic and electron microprobe analysis (EPMA). To investigate how hydrothermal alteration could modify the crystallinity of the PEI deposits, Fourier-transform infrared spectroscopy (FTIR), Raman spectroscopy, and X-ray diffraction (XRD) were carried out on mineral samples and on the matrix of different samples. On each sample (bottom, middle top), we realize 1 XRD analysis, more than 35 FTIR analysis and more than 30 Raman analysis per each thin section. FTIR analyses were performed with a Thermo Scientific apparatus at the Laboratorio de Microanálisis, CEMIEGeo UNAM. Micro-FTIR analyses were measured with 512 scans at a resolution of 4 ­cm−1 and were collected using a KBr beam splitter. The sample was placed on a water-free KBr base for support, and the background scan and sample analyses were taken through the KBr base under the same analytical conditions. Micro-FTIR analyses were obtained over a wavelength range of 5000–600 ­cm−1, with a liquidnitrogen cooled HgCdTe2 (MCT) detector and an aperture set at 10 μm square. This provided a strong intensity and detector response for the mid-IR water absorbance bands. Raman analyses were performed with a Thermo Scientific DXR apparatus at the Laboratorio de Microanálisis, CEMIEGeo UNAM. Raman measurements were carried out in the matrix and minerals of all the PEI deposits. Raman analyses were measured with a 532 nm laser over a wavelength range of 3568–20 ­cm−1, at a resolution of 1.7–2.7 ­cm−1, with a 50 µm aperture and 8.8 mV of laser power. Mineral compositions were determined for a representative subset of the PEI samples. Major-element compositions of plagioclase and Fe-Ti oxides were acquired using a JEOL JXA-8230 electron microprobe housed at the Laboratorio de Microanálisis, CEMIE-Geo UNAM. Analytical conditions were 15 keV accelerating voltage and 10 nA current for major elements and 15 keV and 20 nA for selected trace elements; a focused beam was used in both cases. A defocused electron beam of 10 µm was used to analyze glass groundmass to avoid Na migration. Analytic uncertainties in EPMA were determined with secondary standards with a variety of minerals and glasses. Standards were analyzed throughout each analytical session to monitor instrument drift and correct compositional offsets. Mineral phase identification was achieved using powder X-ray diffraction (XRD) on a Rigaku D/MAX-Rapid microX-ray diffractometer at the Smithsonian National Museum of Natural History Department of Mineral Sciences, EUA. The instrument uses Mo Kα radiation (λ = 0.71069 Å) and collects data on a 2D imaging plate detector with Debye–Scherrer geometry. Sample powders were gently re-ground in an agate mortar and pestle and were mounted in 1 µm-wide Page 5 of 24 63 Contributions to Mineralogy and Petrology (2023) 178:63 Kapton capillaries. The samples were run for 10 min each with omega fixed at 0° and phi rotating at 1°/second. Three sample replicates were collected for each sample by focusing the X-ray beam on three different regions of the powderfilled capillary. XRD image files were background-corrected (manual setting = 4) and integrated into intensity vs 2θ patterns from 4.0 to 45.0 2θ degrees and 81.6337–430.00 β using 2DP software. Mineral phases in our samples were initially identified using the search/match function of Jade (2021) (Materials Data, Inc.) with the PDF-4 2021 ICDD database. Contributions of mineral phases in each of the different samples (represented as phase percentages and weight percentages) were estimated via whole pattern Rietveld refinement analyses using GSASII software (Toby and Von Dreele 2013; Post and Bish 1989). The background was modeled using a Chebyschev polynomial. Anorthoclase [(Al, ­Si3) ­(Na0.667K0.333) ­O8; Harlow 1982], quartz (Antao et al. 2008), cristobalite (Downs and Palmer 1994), tridymite with a monoclinic structure (Hirose et al. 2005), and tridymite with an orthorhombic structure (Lee and Xu 2019) were used as model structures. Although there are likely other feldspar phases present, our focus was to identify and quantify the tridymite and cristobalite phases in these samples. We selected anorthoclase as a single representative feldspar in this system to minimize the number of phases in the Rietveld refinement. Errors in phase estimates represent standard deviations in measurements between the three replicates. Given the mush-like nature of the PEI deposits, we tested with a geochemical model (Fig. 13) if rhyolites have a relationship with melt extracts produced in a crystal mush, probably peralkaline in composition, or if these rocks are not affected by the crystal mush. We chose to perform partial melt models because most authors, regardless of the tectonic context in which volcanism occurs, show that partial melt episodes can occur in the cumulate part of crystal mushes (Foley et al. 2020 and cites therein). This can be promoted by the flow of volatiles and heat coming, in general, from some mafic-intrusive melts in the base of the cumulate portion of the mush (Bachmann et al. 2002; Deering et al. 2011; Foley et al. 2020; Lubbers et al. 2020; Pamukcu et al. 2013; Sliwinski et al. 2017; Szymanowski et al. 2017; Tavazzani et al. 2020). The process above described eventually triggers the mobilization of the crystal mush. For the model, we assume the following: (i) the mush extracts behave like a partial melt, (ii) PEI deposits represent the eruption of a portion of a crystal mush (partial melts + crystals), (iii) we use different granitic rock compositions for the model because a crystal mush can spawn from the crystallization of a magma body or by the partial melting of a pre-existing crystallized body in any tectonic setting. However, we choose the chemical composition of a xenolith of an intrusive granitic rock reported by Sosa-Ceballos et al. (2018) and PEI as a proxy of the partial melting source that produce the melt extracts; (iv) the injection of hot mafic magma promotes the partial melting but lead to produce a minimal chemical interaction with the mush and unevenly distributed thermal flux (Foley et al. 2020); (v) we consider the multiple crystal mush scenario, hosted in different portions of the upper crust. Due to the geochemical modeling complexity, we decided to explore partial melting as a process that could occur in one crystal mush. To run the model, we consider the partial melting of the granitic proxy in a range of F = 0.05–0.45 (5–45%) for the xenolith sample, F = 0.1–1(0.1,0.15,0.2, …,1) for the PEI sample; we used a paragenesis of 0.45 K-Fls + 0.35 Plg + 0.10 Qz + 0.05Opx + 0.05Cpx for the xenolith sample and 0.45 K-Fls + 0.35 Qz + 0.10 Plg + 0.05Opx + 0.05Cpx for the PEI model. We used distribution coefficients of some REE, Rb and Zr in mineral phases published by Rollinson (2014) with the equation Cl = C0/ (F + (D(1-F)). Results The PEI deposits, like most of the ACC rhyolitic deposits, are dispersed in the outer zone of the caldera (Fig. 1b). The PEI is massive, gray white, and crystal rich (Fig. 2). The maximum observed thickness is ~ 100 m, and deposits are generally welded at the bottom and more friable toward the top (Fig. 2b, c). The contact between the local basement and the base of the deposit, or the lag breccia, does not crop out around the ACC, although Avellán et al. (2019) and López-Hernández et al. (2009) suggest that the eastern border of the caldera was the emission point of the PEI. The deposits are composed almost entirely of microphenocrysts and phenocrysts (~ 1% vol of crystal with sizes > 15 cm) of alkali feldspar, quartz, plagioclase, and Fe-Ti oxides (Fig. 3 a, b, c, d) with a crystal content up to ~ 90 wt. % in welded portions and ~ 60 wt. % in friable deposits (Fig. 3). Neither pumice, glass chips, nor lithic fragments are found within the deposits. Portions of the deposit appear to contain an ashy matrix (Fig. 2 d), but the matrix comprises aggregates of quartz and feldspar with sizes > 5 µm, as observed in the welded portion (Fig. 3). The friable portion of the PEI has portions with argillic hydrothermal alteration (Fig. 2a). This alteration was reported by Canet et al. (2015) in the inner part of the caldera. The PEI is rhyolitic in bulk rock composition and has slight internal vertical chemical variations (Fig. 4, Table 1) The compositions of the PEI sample become more silicic and exhibit a decrease in ­Fe2O3 and MgO toward the top of the deposit. We analyzed unaltered samples with LOI values below 1 wt. %. Trace element concentrations normalized to the primitive mantle (Sun and McDonough 1989) show negative anomalies of Ba, Nb, Ta, U, Th, K, and Zr, 13 63 Page 6 of 24 Contributions to Mineralogy and Petrology (2023) 178:63 Fig. 2  Photographs showing the PEI field features of the deposit: a The non-welded top of the deposit, b The welded middle part of the deposit, c welded bottom part of the deposit, and d a cm-scale feldspar phenocryst and a positive anomaly of Pb, Sr, Eu, P, and Ti. Furthermore, heavy rare earth elements (HREE) present a horizontal trend and light rare earth elements (LREE) are more enriched than HREE (Fig. 5). These anomalies and trace element concentrations are also present in the other felsic members of ACC such as Tecoloquillo ignimbrite (TQ) and xenolith fragments of granitic composition reported by SosaCeballos et al. (2018). The xenolith fragment may come from felsic intrusions associated with the ACC magmatic system. Felsic intrusions are reported by CFE exploratory drilling in the region, where they are described as aplitic dikes. ACC's pyroclastic and lava felsic products show low Eu/Eu* ≤ 0.3 and low Ba < 500 ppm (Fig. 6b, e). Moreover, these rocks show a positive trend in Ba/Rb vs Rb, Rb vs La, low #MgO (1–16), and a wide range in the La/Sm vs La diagrams (Fig. 6d). The composition of rhyolites in the ACC is diverse but can be distinguished between A-type and calc-alkaline I-Type rhyolites as a result of the analysis of Zr + Nb + Y + Ce and trace elements (Fig. 6a). Calc-alkaline rhyolites tend to have higher Ba/Sr and lower La/Yb ratios and higher abundances of HFSE and HREE than A-type rhyolites (Fig. 6a). The larger phenocrysts of feldspar are average 1 cm in size (but can exceed 15 cm in length; Fig. 2d) and have 13 subhedral broken shapes (Fig. 3a). Microphenocrysts and phenocrysts show disequilibrium features such as exsolution textures, rounded borders, and resorbed rims (Fig. 3c, d). Furthermore, all the PEI deposits show relics of mafic minerals. Feldspar phenocrysts have anorthoclase and sanidine compositions with 2–8 mol% An and 20–78 mol% Or, whereas the matrix is formed by orthoclase and microcline (Fig. 9). XRD and Raman analyses, however, show that phenocrysts in the samples have an anorthoclase structure (Figs. 10, 11). Oxides in the PEI deposits are ferropseudobrookite–ilmenite and ulvospinel–wustite. Furthermore, oxides display variations of FeO, ­Fe2O3, and ­TiO2 from the center to the rim of the phenocrysts (Supplementary material). Core to rim transects in oxides show Fe enrichment with Ti depletion (Fig. 8). The composition of the PEI matrix samples exhibits a mixture of quartz and felspar (Fig. 10a). EPMA imaging and compositional analysis of the PEI samples reveals that the matrix is composed of microcrystalline texture (< 5 µm) of quartz and feldspar with, 6–13 wt. % of ­Na2O + ­K2O (Fig. 7, 8, 9). FTIR spectroscopy shows primarily anorthoclase and quartz crystals with no hydroxy or water structural defects in the 3000–3800 ­cm−1 bands. Page 7 of 24 63 Contributions to Mineralogy and Petrology (2023) 178:63 Fig. 3  Photomicrographs of four PEI samples in cross-polarized light. a The mineral assemblage formed by phenocrysts of feldspar and quartz and a groundmass formed by tridymite and cristobalite crystals; b The few plagioclase phenocrysts often have disequilib- rium textures, such as rounded borders and sieve texture; c a feldspar phenocryst with exsolution texture and rounded borders; d a feldspar phenocryst with rounded borders and twining Rietveld refinements of X-ray diffraction patterns (Fig. 10a) estimate the relative contribution of different mineral phases outlined in Table 2. Anorthoclase and quartz dominate, as expected, with smaller contributions of cristobalite and tridymite as monoclinic and orthorhombic phases. Cristobalite is not present in PE1904. Orthorhombic tridymite is present in sample AC100. Due to the absence of clay minerals and glass in the XRD patterns, we use the results of Rietveld refinements as a proxy of %vol for each sample (Table 2). Discussion Matrix‑free deposits: a hydrothermal fingerprint The PEI deposits contain microcrysts of quartz, cristobalite, and tridymite as confirmed in FTIR and Raman analysis (Fig. 11a, b). Likewise, XRD measurements show crystal structures of polymorphs of silica (Fig. 10). These mineral phases are related to hydrothermal alteration and 13 63 Page 8 of 24 Fig. 4  Total alkali vs silica diagram of the Acoculco Caldera Complex (ACC, red field) rocks, PE = Piedras Encimadas Ignimbrite simples, Tc = Tecoloquillo ignimbrite, Xen = xenolith, blue cross symbols represent G1, and diamond symbols represents G2 rhyolite ACC lavas of the early and late post-caldera events are reported in other deposits in the inner portion of the ACC (Canet et al. 2015; Sánchez-Córdova et al. 2020). In the less crystalline portions of the PEI, the content of quartz, cristobalite, and tridymite aggregates can reach 40–50 vol% in the welded portions at the bottom part of the deposit. Additionally, to the top of the deposit, crystals are supported by a cryptocrystalline matrix formed by silica polymorphs (Fig. 2a). Thus, we suggest that late hydrothermal alteration affected the deposits even in the external parts of the caldera. Syn- and post-depositional alteration of primary volcanic deposits with the formation of clay minerals have been observed in pyroclastic surge deposits, e.g., at Vulcano, Italy (Capaccioni and Coniglio 1995) and El Chichón, México (Macías et al. 1997). However, we do not have clay minerals. This evidence could suggest that the matrix-free nature of the welded portion of the PEI deposits could be produced by hydrothermal alteration. One of the pending questions to be resolved is how the ­SiO2 polymorphs of the PEI were formed and in which pressure and temperature conditions. Cristobalite and tridymite crystals commonly precipitate at low pressures up to ∼0.4 GPa and temperatures between 870 °C and 1470 °C, but these minerals also exist as metastable phases at lower temperature conditions (Kihara et al. 1986; Heaney et al. 1994); Deer et al. 2004). Tridymite crystals tend to exhibit a hexagonal structure above ~ 380 °C with slight distortions of the Si–O–Si bond angles and Si–O bond length (Kihara et al. 1986). Between 380 °C and 110 °C, tridymite recrystallizes to form an orthorhombic structure, and below 110 °C, monoclinic tridymite is ubiquitous (Heaney et al. 1994). All the PEI deposits contain sharp, well-defined XRD peaks of monoclinic (Fig. 10) tridymite and cristobalite. However, orthorhombic 13 Contributions to Mineralogy and Petrology (2023) 178:63 tridymite is only clearly present in samples AC100 and PE-1902 (Fig. 11). AC100 is located in the SW border of ACC near the topographic rim described by Sosa-Ceballos et al. (2018). Therefore, orthorhombic tridymite found in the PEI deposits could suggest that hydrothermal fluids circulated and reached temperatures up to 110 °C and that these fluids promoted the devitrification process of the deposits. Furthermore, the wüstite composition of the Fe-Ti oxides (Fig. 8) suggest hydrothermal alteration process. Wüstite minerals are reported in hydrothermal veins formed in low fugacity conditions or during the exsolution of pre-existing oxides (Nadoll and Mauk 2011), as observed in S-type granites where the exsolution of oxides should occur during the partial melting of metapelitic rocks (Seifert et al. 2010). The devitrification of rhyolitic pyroclastic rocks is a relatively common process, resulting in the production of microcrystalline silica and feldspars from the fast emplaced and slow cooling deposits (Rowe et al. 2012) and by postdepositional hydrothermal circulation. These chemical alteration processes lead to mass fluxes of elements, both from the altered glass into solution and vice versa, as it has been also reported in mafic rocks from the ocean floor (Furnes 1978); Zhou and Fyfe 1989; Crovisier et al. 1992; Stroncik and Schmincke 2001; Walton et al. 2003; Pauly et al. 2011). In ACC, this process was studied and quantified by Sánchez-Córdova et al. (2020). As a result of the hydrothermal alteration process, the supersaturation of fluids causes the precipitation of neoformation minerals (Hay and Iijima 1968; Ibrahim and Hall 1996; Stroncik and Schmincke 2001; Kousehlar et al. 2012) such as the tridymite and cristobalite observed in the PEI samples (Fig. 7a). Furthermore, the silicification was well-studied in the central ACC (SánchezCórdova et al. 2020). The alteration processes are, thus, pervasive and likely affect the matrix to a greater degree because of the devitrification process after emplacement of the pyroclastic density currents. Feldspars do not show signs of alteration such as albitization (despite of low %Or), a distinct reddening, relicts of its ancient feldspar/plagioclase compositions and/or the absence of micrometric/nanometric pores containing white mica (Engvik et al. 2008), alteration to sericite, zeolites, or other typical alteration minerals (Figs. 10, 11). Thus, the pre-existing ash matrix of the PEI could be transformed to tridymite + cristobalite + K-feldspars (Rowe et al. 2012), as it was observed in the XRD, Raman, FTIR, and EPMA maps (Fig. 7, 10, 11). If this is true, the hypothetical glass (pumice fragments) or matrix (formed by fine ash particles) content could reach less than 40% vol of deposits, as is observed in petrography and quantified by Rietveld refinement of XRD patterns (Table 2). Thus, given the pervasive hydrothermal alteration described in the ACC, we propose the matrix was replaced by secondary minerals. Such replacement processes have been described in other 13.73 Pyroclast 76.91 Type SiO2 5.19 20 2 2.1 11.8 Cs 29.8 31.3 7 Nb Mo Sn 37.3 439 47.4 181 Y Zr Sb 15 207 30 Rb Sr 2.5 0.8 3 3 128 <5 1.8 2.3 18 Ge 22 As 40 24 S Ga 4 21 7 55 Cu Zn 3 1 5 Co 30 Cr 5 99.4 Ni 6 V 2 2 6 Sc 6 #MgO Be 3 101 Total 0.03 0.26 0.02 0.25 P2O5 0.23 0.013 LOI 4.78 0.179 K2O 0.19 4.6 0.44 3.92 CaO Na2O TiO2 0.03 0.04 0.11 MnO MgO 0.10 0.87 0.18 1.60 Fe2O3 0.97 596,510 1278 1278 ± 14 FeO 12.53 2,192,035 North 1.78 586,067 East Al2O3 Pyroclast 74.15 611 Age (Ma) Fe2O3(T) 2,204,497 611 ± 72 Age Ma Lpc PEIAc100 Ec TQ-Ac12 Units Sample 1.4 0.8 3 2 35.5 425 34.9 17 118 1.4 22 0.006 18 2 1 9 4 4 2 99.21 0.82 0.245 5.15 4.49 0.22 0.02 0.01 0.83 0.09 0.92 14.46 72.86 Pyroclast 2,204,605 599,601 1278 1278 ± 14 PEI-1902 Lpc 1.8 1 4 2 36.8 403 40.4 17 126 1.6 20 0.029 54 2 6 5 5 1 100.4 0.91 0.05 0.227 5.04 4.34 0.23 0.03 0.017 1.93 0.22 2.15 13.38 73.99 Pyroclast 2,203,900 597,037 1278 1278 ± 14 PEI-1904 Lpc 6.2 1.2 6 7 30.7 232 60.7 27 202 14 1.9 22 20 56 6 2 38 6 7 3 7 99.76 2.83 0.02 0.178 5.18 3.48 0.55 0.12 0.042 1.45 0.16 1.61 12.56 73.2 Lava 2,194,317 582,915 Ac7 Ec 5.1 1.2 5 8 29.2 160 42.7 8 164 13 2.2 20 40 49 3 1 27 6 3 4 100.2 0.51 0.01 0.135 5.06 3.69 0.38 0.08 0.049 1.60 0.18 1.78 11.86 76.62 Lava 2,208,919 574,826 2553 2553 ± 110 Ac90 Epc 2.3 0.5 2 4 26.4 302 27 75 143 6 1.7 22 60 46 2 12 12 4 4 11 98.86 1.14 0.337 4.56 3.98 0.67 0.24 0.059 1.89 0.21 2.1 13.81 71.98 Lava 2,213,875 590,387 1438 1438 ± 24 Ac107 Lpc 4.1 0.8 4 7 28.3 196 51.1 21 162 9 1.8 21 60 53 22 5 6 3 6 100.8 0.62 0.178 5.07 3.92 0.41 0.09 0.052 1.21 0.13 1.34 12.86 76.27 Lava 2,212,926 588,634 Ac108 Lpc 2.7 0.5 4 2 31.9 500 52.6 2 171 6 1.9 22 50 73 2 31 5 5 2 99.77 0.48 0.21 4.73 4.19 0.08 0.05 0.039 1.93 0.21 2.14 11.96 75.88 Lava 2,196,999 584,295 1870 1870 ± 36 Ac113 Lpc Lpc 100.1 4 4 0.6 3.9 3.9 3 7 27.5 324 40.5 88 138 7 1.7 21 40 69 7 3 18 12 5 0.9 3 4 26.8 324 43.1 91 133 7 1.8 22 30 60 1 70 12 5 16 99.32 12 2.97 0.07 0.362 4.62 4.03 1.04 0.38 0.076 1.99 0.22 2.21 13.21 71.13 Lava 2,200,734 576,764 1700 1700 ± 400 Ac35 0.37 0.07 0.355 4.42 4.24 0.95 0.29 0.09 2.02 0.23 2.25 13.58 72.7 Lava 2,210,789 589,651 1300 1300 ± 600 Ac115 Lpc 5.1 1.3 6 9 30.7 185 44.6 7 174 13 2.2 21 20 52 10 4 2 6 3 5 100.9 3.58 0.01 0.139 5.71 2.87 0.36 0.07 0.046 1.29 0.14 1.43 11.55 75.08 Pyroclast 2,199,284 575,914 1084 1084 ± 22 Ac37 Lpc 1.3 0.4 3 4 33.1 499 42.6 122 121 1.8 24 30 72 3 2 7 11 4 6 5 99.15 1.18 0.05 0.471 4.44 4.73 0.66 0.12 0.038 2.28 0.25 2.53 14.88 70.04 Lava 2,205,382 588,654 1600 1600 ± 100 Ac51 Lpc Table 1  Major and trace geochemistry of ACC samples. Error for major elements is ~ 0.01 wt% and trace elements ~ 1 ppm in average Lpc 3.5 0.6 3 4 18.5 227 27.2 52 150 6 1.6 20 20 46 2 2 6 6 3 3 2 99.72 1.44 0.02 0.24 4.67 4.01 0.47 0.04 0.028 1.42 0.16 1.58 14.22 73.01 Lava 2,204,474 581,401 1400 1400 ± 200 Ac56 Lpc 4.4 0.7 4 7 25.4 168 44.6 11 164 11 2 21 20 51 2 2 14 6 2 3 100.6 2.82 0.01 0.143 5.46 3.2 0.34 0.04 0.048 1.13 0.13 1.26 12.58 74.7 Lava 2,199,502 578,889 1283 1283 ± 88 Ac72 Lpc 5.2 1.3 6 9 34.3 271 53.4 4 165 15 2.1 21 20 60 5 3 1 7 3 3 100.7 3.67 0.135 5.27 3.32 0.31 0.05 0.048 1.55 0.17 1.72 11.46 74.71 Lava 2,190,307 595,629 1145 1145 ± 14 Ac82 Lpc 2.5 0.2 3 5 28.5 446 87.9 74 145 7 1.5 22 40 58 4 2 22 19 5 5 3 100.1 1.93 0.22 0.385 4.44 4.04 0.56 0.08 0.045 2.63 0.29 2.92 13.24 72.19 Lava 2,208,434 572,356 1066 1066 ± 42 Ac89 Lpc 3.4 1 3 3 27.2 181 54.9 21 151 9 1.9 21 30 47 3 2 20 5 2 5 99.72 0.67 0.02 0.182 4.96 3.67 0.3 0.08 0.051 1.30 0.14 1.44 12.99 75.36 Lava 2,213,879 583,057 1360 1,360 ± 15 Ac92 Epc 3.1 0.2 20.4 136 38.6 19 136 1.2 15 30 42 <1 2 54 6 4 2 7 99.66 1.14 0.01 0.126 3.09 2.52 0.3 0.09 0.023 1.10 0.12 1.22 8.02 83.11 Xen 2,197,201 584,862 Ac112 Contributions to Mineralogy and Petrology (2023) 178:63 Page 9 of 24 63 13 13 11.4 3.97 0.33 16 13.2 4.42 6.77 1.27 3.5 0.37 16 19.1 6.23 7.18 W Tl Pb Th U Eu* 9.3 1.9524 Ta/U PM 0.2885 0.2787 8.4881 Ta* Nb/Th PM 0.2660 0.5441 0.1931 Ta/U Nb* 1.9524 8.4881 0.5633 2.2576 0.64 1.6387 Ba/Rb Nb/Th 2.1 2.49 5.5 3.39 Hf 0.66 4.59 0.69 6.43 Ta 5.59 0.79 Yb Lu 4.89 0.845 Er Tm 1.35 7.55 1.59 Dy Ho 5.82 1.03 6.1 1.2 Gd 0.44 7.6 40.3 Tb 8.12 0.374 Sm Eu 13.6 45 Pr 52 108 65.3 125 La Ce Nd 163 133 Ba Lpc PEIAc100 Ec TQ-Ac12 Sample Table 1  (continued) Units Lpc 1.9524 8.4881 0.3319 0.3009 0.6480 2.5540 1.69 5.36 3.92 13.9 14 0.27 1.1 2.54 10.7 0.655 3.95 0.585 3.49 1.12 5.48 0.85 4.88 0.543 5.8 28.5 8.07 68.7 35.8 200 PEI-1902 Lpc 1.9524 8.4881 0.3359 0.3097 0.6559 2.6286 1.21 7.59 4.01 14 16 0.36 1.2 2.63 9.8 0.726 4.39 0.643 4.09 1.44 7.14 1.14 6.82 0.434 8.29 43.2 12.2 100 52.8 152 PEI-1904 Ec 1.9524 8.4881 0.3272 0.1739 0.6388 1.4760 0.65 8.59 6.09 20.8 17 1.23 373 3.89 7.1 1.02 6.7 0.991 6.25 1.98 9.72 1.5 8.37 0.352 8.8 45.5 13.3 121 64.1 131 Ac7 Epc 1.9524 8.4881 0.2984 0.2098 0.5826 1.7805 0.45 6.68 5.51 16.4 15 0.72 252 3.21 5 0.69 4.88 0.738 4.2 1.39 6.71 1.08 5.73 0.229 7.51 40.7 11.8 110 58.3 74 Ac90 Lpc 1.9524 8.4881 0.3430 0.1932 0.6696 1.6398 3.01 4.67 4.48 16.1 17 0.97 84 3 7.7 0.507 3.21 0.475 2.94 0.89 4.71 0.75 4.25 1.02 5.06 25.8 7.41 89.5 37.7 430 Ac107 Lpc 1.9524 8.4881 0.3475 0.2045 0.6785 1.7362 1.44 8.00 4.79 16.3 14 0.91 233 3.25 6.1 0.797 5.17 0.795 5.06 1.63 8.4 1.27 7.51 0.642 8.47 43.4 12.3 110 62 233 Ac108 Lpc 1.9524 8.4881 0.3589 0.2306 0.7006 1.9571 0.23 9.38 4.81 16.3 25 0.71 254 3.37 12.5 0.956 6.29 0.874 5.8 1.87 9.9 1.56 8.91 0.375 9.83 52.3 14.8 116 66.7 39 Ac113 Lpc 1.9524 8.4881 0.3088 0.2133 0.6029 1.8108 3.78 7.69 4.81 14.8 17 0.63 12.1 2.9 7.9 0.717 4.66 0.694 4.42 1.46 7.16 1.16 7.12 1.2 8.22 42.3 12.1 125 54.5 503 Ac115 Lpc 1.9524 8.4881 0.3390 0.2282 0.6619 1.9366 3.31 6.89 4.2 14.2 16 0.89 156 2.78 8.6 0.674 4.76 0.645 4.2 1.36 6.99 1.08 6.59 1.14 7.17 37.3 10.4 98.3 48.9 457 Ac35 Lpc 1.9524 8.4881 0.2921 0.2219 0.5703 1.8834 0.39 7.28 5.26 16.3 12 0.8 2.9 3 5.5 0.732 5.12 0.748 4.34 1.42 7.02 1.15 6.43 0.248 8.04 41.9 12.6 112 59.2 67 Ac37 Lpc 1.9524 8.4881 0.4008 0.3333 0.7825 2.8291 5.22 8.41 3.54 11.7 15 0.39 58.3 2.77 11.7 0.721 4.83 0.69 4.64 1.5 8.11 1.36 7.79 2.22 8.98 42.3 11.9 110 51.1 632 Ac51 Lpc 1.9524 8.4881 0.2767 0.1627 0.5403 1.3806 4.07 5.17 3.85 13.4 16 0.57 38.7 2.08 6.1 0.437 3.16 0.452 2.88 0.92 5.05 0.81 4.68 0.761 5.61 32.8 9.86 62.1 50.9 611 Ac56 Lpc 1.9524 8.4881 0.3344 0.1691 0.6529 1.4350 0.82 7.99 4.84 17.7 17 1.2 122 3.16 5.5 0.735 5.21 0.7 4.82 1.55 7.77 1.26 7.45 0.327 8.49 44.7 13.2 112 63.8 134 Ac72 Lpc 1.9524 8.4881 0.3110 0.2377 0.6072 2.0176 0.13 8.60 5.55 17 15 0.88 2.7 3.37 7.6 0.856 6.04 0.895 5.23 1.7 8.57 1.35 7.34 0.186 9.69 51.5 14.8 124 67 22 Ac82 Lpc 1.9524 8.4881 0.2959 0.2269 0.5777 1.9257 2.08 13.80 5.21 14.8 12 0.68 196 3.01 10.3 1.35 8.75 1.27 8.49 2.81 14.3 2.29 13.7 0.855 13.9 70.2 19 119 89.9 302 Ac89 Lpc 1.9524 8.4881 0.2884 0.2094 0.5631 1.7778 1.40 8.59 4.99 15.3 13 0.51 5.4 2.81 5.1 0.852 5.63 0.849 5.06 1.63 8.03 1.35 7.52 0.608 9.54 49.2 14 101 61.8 212 Ac92 Epc 1.9524 8.4881 0.2157 0.1892 0.4212 1.6063 0.63 5.61 6.41 12.7 10 0.22 489 2.7 4.7 0.674 4.39 0.61 3.91 1.26 6.04 0.92 5.3 0.238 5.91 29.5 8.69 83.7 42.6 85 Ac112 63 Page 10 of 24 Contributions to Mineralogy and Petrology (2023) 178:63 Contributions to Mineralogy and Petrology (2023) 178:63 Page 11 of 24 63 Fig. 5  Trace elements, including rare earth elements (REE), normalized to the chondrite (right) and the primitive mantle (left; Sun and McDonough 1989). Symbols as in Fig. 4. The relative enrichment of LILE (large-ion lithophile elements) when compared to the HFSE (high-field-strength elements) and the Nb–Ta–Pb anomalies reveal the calc-alkaline nature of the felsic ACC. The xenolith sample corresponds to a sample reported by Sosa-Ceballos et al. (2018). Tc compositions were obtained from Pérez-Orozco et al. (2021) volcanic systems to have occurred during or after the depositional process. deposits represent explosive, volcanic eruptions. One hypothesis that simplifies the origin of high crystallinity is that the PEI formed as part of an intrusive body that was exhumed and exposed by erosion. We think this hypothesis is unlikely because the deposits are welded at the bottom, but not the top, the deposit has a tabular shape, and there is a systematic thinning of the deposits that fill the valleys toward the east (Fig. 1b). Field evidence suggests these rocks are ignimbrites that erupted from the eastern border of the caldera, where the deposits are thicker and thin toward the northeast (López-Hernández et al. 2009; Origin of the deposits The PEI deposits were originally described as ignimbrites (López-Hernández et al. 2009; Avellán et al. 2019, 2020). Nonetheless, its high crystallinity (almost totally formed by crystals) and absence of pumice, ash, and low lithic contents have not been described in the literature. Given its high crystal content, we first investigate if these 13 63 Page 12 of 24 Contributions to Mineralogy and Petrology (2023) 178:63 Fig. 6  a Zr + Nb + Y + Ce vs ­SiO2 diagram of the ACC rhyolites, a dotted line separates A- from I-type rhyolites based on Whalen et al. (1987). b Rb/Ba vs Rb diagram displays the relation between crys- tallization vs partial melting of the ACC. c #Mg vs S ­ iO2 of the PEI rhyolites shows the low #Mg contents of the most of samples. d La/ Sm vs La diagram, e Eu/Eu* vs S ­ iO2 diagram. Symbols as in Fig. 4 Avellán et al. 2019). Due to hydrothermal alteration and the massive structure of the deposits (Fig. 2a, b, c), it is hard to tell if PEI was formed by the aggradation of one eruption or by multi-PDC. However, Avellán et al. (2019, 2020) describe several flow units that suggest a multi-PDC origin for the PEI. The PEI deposits are highly crystalline: the top of the deposit contains more than 70% phenocrysts with sizes ≥ 2 mm, whereas the middle and bottom portions contain up to 40–60% of phenocrysts (Fig. 3a). Moreover, the phenocrysts are embedded in a groundmass formed exclusively by microcrysts of polymorphs of silica and K-feldspar, devoid of glass. The presence of K-feldspar megacrysts suggests that the PEI magmas crystallized near its rheologic lock-up threshold, as K-feldspar phenocrystals do not precipitate significantly while they are mobile (e.g., Glazner and Johnson (2013). Arzilli et al. (2020) suggest that changes in crystal fraction of alkali feldspar can range from hours to several 13 Contributions to Mineralogy and Petrology (2023) 178:63 Page 13 of 24 63 Fig. 7  Matrix composition shown as a X-ray maps of the matrix and feldspar. Semiquantitative X-ray maps show the homogeneous composition of the anorthoclase phenocryst and the composition of the matrix that is formed entirely of quartz polymorphs and k-feldspars. The lateral scale shows the relative concentrations of Si, K, and Na in relative intensities [arbitrary units]. b Total alkali vs silica diagram of the matrix made with quantitative EPMA analysis to have a compositional reference days as a function of undercooling and ­H 2O content in the melt generating the alkali feldspar nucleation delay and growth. Therefore, we suggest that the low changes in Or content in the PEI phenocrysts of anorthoclases suggest large periods of storage and growth under same P–T conditions. Thus, the size and abundance of crystals (some of them with exsolution textures) and the absence of pumice strongly suggest that the PEI deposits represent a crystal mush (crystals + melt) mobilized to the surface. 13 63 Page 14 of 24 Fig. 8  a and b Micrographs of ferropseudobrookite and ulvospinel crystal whit a–a´ and b–b´ transects. c and d ­FeOtotal (blue line and rectangles) and ­TiO2 (green line and rectangles) concentrations vs distance (µm). e and f Oxides in rutile– wustite–hematite classification diagrams. Ferropseudobrookite has a relatively consistent composition across the transect a—a’. Disequilibrium in a and b is evidenced by the rounded form and bay rims of the crystals. The compositional changes in ulvospinel tend to be more extreme than those observed in ferropseudobrookite. g and h Oxide exsolution 13 Contributions to Mineralogy and Petrology (2023) 178:63 Contributions to Mineralogy and Petrology (2023) 178:63 Fig. 9  Feldspar compositional transect, a secondary electron image of an anorthoclase crystal, b compositional transect (A–A’) of anorthoclase. Note the slight compositional variation from center to rim Page 15 of 24 63 where the center tends to be more potassic. c) Felspar classification diagram showing low felspar compositional variations Fig. 10  Representative X-ray diffraction patterns of volcanic ash samples (black) with Rietveld refinement fits (grey) and background fits (red-dashed) averaged across three sample replicates. Peak positions of reference mineral phases are displayed above with the most dominant peaks outlined in black (Harlow 1982; Levien et al. 1980; Downs and Palmer 1994; Hirose et al. 2005; Lee and Xu 2019) 13 63 Page 16 of 24 Contributions to Mineralogy and Petrology (2023) 178:63 Table 2  Mineral phases determined by XRD on PEI samples Sample Phase Phase % s.d Wt % s.d ACC100 Anorthoclase Quartz Cristobalite Tridymite–Monoclinic Tridymite–Orthorhombic Total Anorthoclase Quartz Cristobalite Tridymite–Monoclinic Total Anorthoclase Quartz Cristobalite Tridymite–Monoclinic Total Anorthoclase Quartz Cristobalite Tridymite–Monoclinic Total 51.6 16.2 9.6 6.2 16.4 100 51.1 34.2 13.8 0.9 100 42.2 52 4.3 1.5 100 76 18.4 4 1.6 100 0.8 3.4 2.6 0.5 1.6 63.1 4.9 2.6 20.3 9 100 78.4 13 4.7 3.8 100 70.3 21.5 1.6 6.6 100 88.5 5.3 1.1 5.1 100 1 0.6 0.4 0.9 0.6 PE1902 PE1904 PE1908 7.5 5.5 0.2 0.3 1.6 2.9 5.7 0.5 1.1 3.4 4.5 0.5 0.3 0.5 0.3 0.4 0.4 1 0.9 0.8 1.7 0.8 0.9 1.1 Storage before eruption Recent geological, geophysical, and petrological studies proposed a series of magma storage zones below the ACC (Sosa-Ceballos et al. 2018; Avellán et al. 2020; PérezOrozco et al. 2021). These storage zones match well with the depths proposed by Townsend and Huber (2020) for magma reservoirs (150–250 MPa, approximately 6–10 km) and transitory-shallow magma accumulations (< 150 MPa), although they are slightly shallower than those proposed by Gualda and Ghiorso (2013) for rhyolites with the observed silica content. Petrologic evidence including exsolution textures and phenocryst sizes (assuming that phenocryst growth occurs at depth) found in the PEI deposits support a long, protracted interval of crystallization in a plutonic body. Furthermore, ambient noise seismic tomography of the anisotropic shear wave velocities reveals mafic intrusion below the E caldera border and NW outer zone (Perton et al. 2022). Geothermal exportation wells in Acoculco reach 1900 m below the surface and intersected an intrusive body at ~ 1600 m (Yáñez-García 1980; López-Hernández et al. 2009). Aeromagnetic studies found a series of magmatic intrusive bodies hosted by calcareous rocks at depths of 1000–2500 m below the ACC (López-Hernández et al. 2009; Avellán et al. 2020). Given their location, depth, and magnetic susceptibility, these intrusions were considered to be mafic. 13 We suggest a series of mafic magmas intruded the local basement, creating skarns and marbles, and underplated vast zones of the ACC surrounding areas; evidence for this is the magnetic anomaly interpreted for the NE border of the ACC (Fig. 13), below the PEI deposits, and the occurrence of monogenetic-mafic volcanism along the SE–NE vicinity of the ACC. We propose that these mafic magmas underplated the crystal-rich mush at depths of 100–250 MPa, progressively contributing to its partial melting (e.g., Sas et al. 2021; Eichelberger 2020). This eventually mobilized the mush sufficiently for the magma to erupt. Feldspar crystals in the PEI deposits contain rounded borders (Fig. 3) but are characterized by almost constant core–rim compositions (Figs. 7, 9). We propose that rim disequilibrium is produced by temperature changes during the ultimate processes that led to the eruption. Constant core–rim compositions could be produced by crystallization over a narrow range in P–T (e.g., such as might happen near the eutectic which is where we would expect to get such a high vol. fraction of crystals at constant composition), which drive melt-crystal evolution, buffered by volatile flux (e.g., Blundy et al. 2006). This process would maintain essentially constant crystal composition. Volatiles could be fluxed by magmas underplating the intrusive bodies (e.g., Huber et al. 2010; Sas et al. 2021). The interaction of the mush with new mafic pulses and its volatiles (Fig. 12), such as that occurs across other caldera systems (i.e., Sohn et al. 2009; Watts et al. 2016; Repstock et al. 2018), could initiate mush mobilization prior to eruption; Pérez-Orozco et al. (2021) and Peiffer et al. (2014) reported anomalous high fluxing of CO2 within the caldera perhaps as the result of mafic magmas degassing at depth. Thus, the mass transfer of volatiles from the new intrusion brings heat to the system. In consequence, this process can enhance “defrosting” and possibly could generate the subsequent remobilization of the mush to the surface (Huber et al. 2010). Mobilization to the surface Despite the processes of extraction and pre-eruptive storage are often treated in tandem (Bachmann and Bergantz 2004a, 2008b), it is important to separate the two processes (Gualda et al. 2019). We discuss the mobilization of PEI as a crystalrich magma in the context of ACC and its potential transportation to the surface. The transport of crystal-rich magma from the upper crust to the surface represents a rheological and structural problem as the high magma viscosity and yield strength hinders easy transport. To address this issue, we suggest that the PEI mush had a critical melt fraction prior to eruption and, most importantly, it was in an eruptible state. The melt fraction promoted a continuum transition from a partially molten solid to a mushy liquid at about Contributions to Mineralogy and Petrology (2023) 178:63 Page 17 of 24 63 Fig. 11  a Raman spectra of micro-crystals surrounding phenocrysts from different stratigraphic positions, as well as, a characteristic feldspar spectrum (Fritz et al. 2005). b FTIR spectra of tridymite and anorthoclase crystals from the bottom of the PEI. Note that the OH– H2O bands (3000–3800 ­cm−1) show no structural dislocation/defects of water or hydroxyl. This could suggest the diffusion of OH–H2O due to hydrothermal alteration 50–55 vol% crystals (e.g., (Marsh 2000). Furthermore, the bulk viscosity could change orders of magnitude between 40 and 60% crystallization (Lejeune and Richet 1995). Above the critical crystallinity of ~ 50 vol%, the crystals form a strong interlocking network that restricts flow, whereas viscosity and yield strength drop rapidly as crystallinity declines (Marsh 2000). Moreover, the transient state from Newtonian to non-Newtonian conduct occurs between 30 and 50% of crystals (Petford 2003), and the drop in the viscosity could be promoted by the recharge of hotter mafic magmas. For that reason, these critical rheological shifts arise close to the critical crystallinity for the very large volume of rhyolites described here. The lower viscosities produced through the hot magma recharges allows felsic magmas reach the surface in form of dikes, with width less than 20 m (Petford et al. 1994). These dikes are probably located both inside and outside the caldera (Fig. 1), following post-collapse faults generated during the lithospheric extensional episodes in the region (García-Palomo et al. 2002, 2018) let the magmas to infill the faults through the dike propagation (Petford et al. 1994; Weinberg 1996; Karlstrom et al. 2012; Cruden and Weinberg 2018). Dike propagation can experience a rheological transition when the eruption starts (Karlstrom et al. 2012) In the case of the PEI, we propose that extension could generate structural conditions that facilitate mafic recharge. Recharge events would 13 63 Page 18 of 24 in turn decrease the PEI viscosity and yield strength, at least temporarily, and enhance the likelihood of eruption. Consequently, the PEI eruption occurred in the eastern part of the 13 Contributions to Mineralogy and Petrology (2023) 178:63 ACC, where the displacement or reactivation of a regional fault, likely the ancient caldera rim, led to the eruption. This Contributions to Mineralogy and Petrology (2023) 178:63 Page 19 of 24 63 ◂Fig. 12  a) Magnetic field reduced to the pole (MFRP) map of the ACC region modified from Avellan et al. (2020); EAC1 and EAC2 represent the CFE boreholes. We propose that the anomalies characterized by -25 to 25 [nT] (center and NE of the caldera) represent mafic intrusions that contributed to form the crystal mush and expulsed it to the surface; the mafic magmas not only stalled at depth, but they also formed the monogenetic volcanos on the eastern side of the caldera. b Simplified cross-section profile (A–A’) showing the location of shallow intrusions and leftovers interpreted from the MFRP map. Given the crystal-rich nature of the PEI, we suggest that the crystal mush was very shallow and previously deformed by normal faults that facilitated the extraction to the surface fault is associated with negative magnetic anomalies and lineaments as described by Avellán et al. (2020). Geochemistry of the rhyolites and their relationship with the crystal mush The relationship between crystal mushes, melts that generate crystal-rich ignimbrites, and their impact on subsequent effusive volcanism has not been studied in detail. Because of the compositional diversity of its felsic eruptive products, the ACC (Fig. 1a) is an ideal location for studying how highly silicic, viscous magmas impact upper crustal evolution (Fig. 4). Rhyolite lavas are the most common rock in the ACC (Fig. 1b). Overall, these rocks can be classified into two groups: 1) I-type, classic hydrous calc-alkaline arc rocks, and 2) A-type ‘dry’ rocks (Fig. 6a) that have classical subduction-related trace element anomalies such as negative anomalies of Nb and Ta, also positive Pb anomalies (Fig. 5). This classification cannot determine if a rhyolite was formed by a mush. Hence, we classified the rhyolitic lavas into two groups as a function of their trace elements: the first group (G1) shows variable Eu/Eu*, low Ba contents, and variable #Mg, and the second group (G2) shows low Eu/Eu*, variable Ba contents, and low #Mg (Fig. 6). In addition, G1 lavas have variable concentrations of Sr, P, Zr, and REE, whereas G2 lavas are less variable, although the negative anomalies shown in Ba, P, Sr, and Ti are more pronounced (Fig. 5). G1 trace elements can be explained by events of fractional crystallization of plagioclase and ferromagnesian minerals. Whereas trace elements in G2 rocks are relatively constant and mimic the patterns and anomalies found in the PEI rocks (anomalies in Ba, Sr, P and Ti are characteristic in peralkaline rhyolites). However, we cannot neglect that some fractional crystallization occurred to form G2 rocks. Thus, rocks related to crystal mush extracts probably show features of partial melting caused by the distributed dissolution of crystals. Crystal dissolution is a process described in crystal mush systems (Bachmann et al. 2002; Deering et al. 2011; Pamukcu et al. 2013; Sliwinski et al. 2017; Szymanowski et al. 2017; Foley et al. 2020; Lubbers et al. 2020; Tavazzani et al. 2020) promoted, in general, by the injection of hot mafic magma that can lead to minimal chemical interaction and unevenly distributed Fig. 13  Geochemical models consider diverse degrees of partial melting (0.5–40% of partial melting) of a granitic xenolith (orange lines), lower (red line), middle (light red), and upper crust (pink line) showing the behavior of ACC rhyolitic samples in a La/Eu vs Rb, b La/Zr, c REE concentrations thermal flux (Foley et al. 2020). Consequently, the rocks produced by the partial melting of crystal mush will inherit certain geochemical patterns sensitive to modeling. 13 63 Page 20 of 24 Contributions to Mineralogy and Petrology (2023) 178:63 Fig. 14  Schematic model of rhyolite petrogenesis at the ACC. The generation of multiple melt compositions at the ACC is driven by fractional crystallization of subduction-related intermediate magmas and the mush extracts from A-type peralkaline intrusive bodies. The PEI represent mush (melt + crystals) mobilized to the surface, other ACC rhyolites seem to be pure melts extracted from the crystal mush and other I-type magmas. The underplating of mafic magmas drives some parts of the hot plutonic body to a crystal mush state. The thermal and volatile input also reduces density and viscosity through partial melting of mushed parts of the pluton increasing its buoyancy As observed in G1 and G2, some rocks fit the partial melt models; G2 rocks fit better the model performed with the xenolith as end member, > 40% of partial melting (Fig. 13), whereas G1, as seen in the model using PEI as end member, could represent relatively high degrees of PEI partial melting (Fig. 13). Although G1 and G2 overall fit the models, we cannot exclude the possibility that multiple crystal mushes serve as magma sources; idealizing such models suffer well-recognized limitations: lack of samples, uncertainty about the composition abundance of comprising phases and subsurface alteration. Hence, given the differences shown by G1 and G2 rocks, our best approximation is to suggest that rhyolites in the ACC were produced by at least two different sources, one of which could be related to the partial melting of the rocks remobilized to be deposited as the PEI. in this part of the TMVB. Thus, mafic magmas that erupted through monogenetic vents around the ACC, but these features are equally or even more common in western TMVB, served as the heating element melting the intrusive body, forming the crystal mush and the G2 rhyolites. Melting and extracting the crystal mush dispersed the pyroclastic density currents that emplaced the PEI deposits at the surface. Rhyolite formation model We propose a model for rhyolite formation in the ACC (Fig. 14). Primary melts are formed in a subduction zone for the TMVB (e.g., Gómez-Tuena et al. (2018); Parolari et al. (2021). The G1 rhyolites were formed by fractional crystallization of intermediate magmas, probably of tholeiitictransitional compositions that arrived at the upper crust and assimilated some country rock (e.g., Sosa-Ceballos et al. 2021). Given the distance from the trench ~ 450 km, the crustal thickness ~ 40 km, and the depth of the slab below the ACC ~ 300 km (Pérez-Campos et al. 2008; Ferrari et al. 2012; Castellanos et al. 2018), some melts formed A-type peralkaline intrusive bodies in a back-arc type environment (e.g., Karsli et al. 2012; Jiang et al. 2006). Perhaps, the influence of the mantle flow in this region, produced by the rollback and slab tear of the Cocos plate (Castellanos et al. 2018), generated the tectonic extension that promoted the ascent of the melts 13 Conclusion The ultra-crystalline rhyolitic PEI was erupted ~ 1.2 Ma ago through the eastern portion of the ACC. The PEI is composed primarily of compositionally homogeneous K-felspar and quartz polymorphs. The high crystallinity (close to ~ 100%) in the PEI is the result of two independent main processes: magmatic crystallization and hydrothermal recrystallization of the matrix. Protracted crystallization at shallow depths produced up to the 50% of crystals and postdepositional hydrothermal alteration causes devitrification of the remaining matrix to produce cristobalite + tridymite crystals. The presence of this ultra-crystal-rich rhyolitic ignimbrites in ACC suggests that some felsic eruptions are the product of the remobilization of a crystalline mush. Injection or recharge of hot magma resulted in the transport of the mush and its extracts from the crust to the surface. We suggest that the mafic magmas that erupted through monogenetic vents around the ACC served as the heating element that promoted the partial melting and mobilization of underplate intrusive bodies. Thus, the geochemical diversity of the ACC rhyolites reflects a complex combination of different magmatic processes. The ACC rhyolites can be divided into G1 (“I-type”) calkaline and G2 (“A-type”) peralkaline. We suggest that most of the G1 and G2 rhyolites Contributions to Mineralogy and Petrology (2023) 178:63 could represent relatively high degrees of partial melting and were produced by at least two different sources, one of which could be related to melt extracts of the crystal mush in which PEI magma were formed. Supplementary Information The online version contains supplementary material available at https://d​ oi.o​ rg/1​ 0.1​ 007/s​ 00410-0​ 23-0​ 2043-7. Acknowledgements Funding for this research was provided to JLM by Project#15 of CEMIE-Geo and WP4.4 of GEMEX. We thank Gabriela Reyes Agustín for sample preparation and his assistance with DRX, FTIR, and Raman analysis; Felipe García Tenorio for field mapping, sampling, thin section, and crystal preparation; Fabiola Mendiola for XRD sample preparation; Dario Torres for their support and guidance in the geochemical model; John Stix, Lydia Harmon, Gordon M. Moore, and unrevealed revisors for their comments that improve the manuscript. Data availability All the data was setted in the manuscript and supplementary information. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://​creat​iveco​mmons.​org/​licen​ses/​by/4.​0/. References Anderson AT, Davis AM, Lu F (2000) Evolution of Bishop Tuff Rhyolitic Magma Based on Melt and Magnetite Inclusions and Zoned Phenocrysts. J Petrol 41:449–473. https://d​ oi.o​ rg/1​ 0.1​ 093/​ PETRO​LOGY/​41.3.​449 Antao SM, Hassan I, Wang J et al (2008) State-of-the-art high-resolution powder x-ray diffraction (HRPXRD) illustrated with Rietveld structure refinement of quartz, sodalite, tremolite, and meionite. Can Mineral 46:1501–1509. https://​doi.​org/​10.​3749/​ CANMIN.​46.5.​1501 Arzilli F, Stabile P, Fabbrizio A et al (2020) Crystallization Kinetics of Alkali Feldspar in Peralkaline Rhyolitic Melts: Implications for Pantelleria Volcano. Front Earth Sci (lausanne). https://​doi.​ org/​10.​3389/​FEART.​2020.​00177/​BIBTEX Avellán DR, Macías JL, Layer PW et al (2019) Geology of the late Pliocene-Pleistocene Acoculco caldera complex, eastern TransMexican Volcanic Belt (México). J Maps 15:8–18 Avellán DR, Macías JL, Layer PW et al (2020) Eruptive chronology of the Acoculco caldera complex: A resurgent caldera in the eastern Trans-Mexican Volcanic Belt (México). J South Am Earth Sci. https://​doi.​org/​10.​1016/J.​JSAMES.​2019.​102412 Bachmann O, Bergantz GW (2003) Rejuvenation of the Fish Canyon magma body: A window into the evolution of large-volume silicic magma systems. Geology 31:789–792. https://​doi.​org/​ 10.​1130/​G19764.1 Page 21 of 24 63 Bachmann O, Bergantz GW (2004a) On the origin of crystal-poor rhyolites: extracted from batholithic crystal mushes. J Petrol 45:1565–1582 Bachmann O, Bergantz G (2008a) The Magma Reservoirs That Feed Supereruptions. Elements 4:17–21. https://​doi.​org/​10.​2113/​ GSELE​MENTS.4.​1.​17 Bachmann O, Bergantz GW (2008b) Rhyolites and their Source Mushes across Tectonic Settings. J Petrol 49:2277–2285. https://​doi.​org/​10.​1093/​PETRO​LOGY/​EGN068 Bachmann O, Dungan MA, Lipman PW (2002) The Fish Canyon Magma Body, San Juan Volcanic Field, Colorado: Rejuvenation and Eruption of an Upper-Crustal Batholith. J Petrol 43:1469–1503. https://​d oi.​o rg/​1 0.​1 093/​P ETRO​LOGY/​4 3.8.​ 1469 Bindeman I (2008) Oxygen Isotopes in Mantle and Crustal Magmas as Revealed by Single Crystal Analysis. Rev Mineral Geochem 69:445–478. https://​doi.​org/​10.​2138/​RMG.​2008.​69.​12 Bindeman IN, Valley JW (2003) Rapid generation of both high- and low-δ18O, large-volume silicic magmas at the Timber Mountain/Oasis Valley caldera complex, Nevada. GSA Bull 115:581– 595. https://​doi.​org/​10.​1130/​0016-​7606(2003)​115%​3c0581:​ RGOBHA%​3e2.0.​CO;2 Blundy J, Cashman K (2008) Petrologic Reconstruction of Magmatic System Variables and Processes. Rev Mineral Geochem 69:179– 239. https://​doi.​org/​10.​2138/​RMG.​2008.​69.6 Blundy J, Cashman K, Humphreys M (2006) Magma heating by decompression-driven crystallization beneath andesite volcanoes. Nature 443:76–80 Brophy JG (1991) Composition gaps, critical crystallinity, and fractional crystallization in orogenic (calc-alkaline) magmatic systems. Contrib Miner Petrol 109:173–182. https://​doi.​org/​10.​ 1007/​BF003​06477/​METRI​CS Browne PRL (1978) Hydrothermal alteration in active geothermal fields. Annu Rev Earth Planet Sci 6:229–248 Calcagno P, Trumpy E, Gutiérrez-Negrín LC, Norini G, Macías JA, Carrasco-Núñez G, Liotta D, Garduño-Monroy VH, Hersir GP, Vaessen L, Evanno G (2019) Updating the 3D Geomodels of Los Humeros and Acoculco Geothermal Systems (Mexico)-H2020 GEMex Project Canet C, Hernández-Cruz B, Jiménez-Franco A et al (2015) Combining ammonium mapping and short-wave infrared (SWIR) reflectance spectroscopy to constrain a model of hydrothermal alteration for the Acoculco geothermal zone, Eastern Mexico. Geothermics 53:154–165. https://​doi.​org/​10.​1016/J.​GEOTH​ERMICS.​2014.​ 05.​012 Capaccioni B, Coniglio S (1995) Varicolored and vesiculated tuffs from La Fossa volcano, Vulcano Island (Aeolian Archipelago, Italy): evidence of syndepositional alteration processes. Bull Volcanol 57:61–70. https://​doi.​org/​10.​1007/​BF002​98708/​METRI​CS Cashman KV, Giordano G (2014) Calderas and magma reservoirs. J Volcanol Geoth Res 288:28–45. https://​doi.​org/​10.​1016/J.​ JVOLG​EORES.​2014.​09.​007 Cashman KV, Sparks RSJ, Blundy JD (2017) Vertically extensive and unstable magmatic systems: A unified view of igneous processes. Science. https://​doi.​org/​10.​1126/​SCIEN​CE.​AAG30​55 Castellanos JC, Clayton RW, Pérez-Campos X (2018) Imaging the Eastern Trans-Mexican Volcanic Belt With Ambient Seismic Noise: Evidence for a Slab Tear. J Geophys Res Solid Earth 123:7741–7759. https://​doi.​org/​10.​1029/​2018J​B0157​83 Clemens JD, Bryan SE, Mayne MJ et al (2022) How are silicic volcanic and plutonic systems related? Part 1: A review of geological and geophysical observations, and insights from igneous rock chemistry. Earth Sci Rev. https://d​ oi.o​ rg/1​ 0.1​ 016/J.E ​ ARSCI​ REV.​ 2022.​104249 13 63 Page 22 of 24 Crovisier JL, Honnorez J, Fritz B, Petit JC (1992) Dissolution of subglacial volcanic glasses from Iceland: laboratory study and modelling. Appl Geochem 7:55–81. https://​doi.​org/​10.​1016/​S0883-​ 2927(09)​80064-4 Cruden AR, Weinberg RF (2018) Mechanisms of Magma Transport and Storage in the Lower and Middle Crust—Magma Segregation, Ascent and Emplacement. In: Steffi Burchardt (ed) Volcanic and igneous plumbing systems 13–53. https://​doi.​org/​10.​1016/​ B978-0-​12-​809749-​6.​00002-9 Deer WA, Howie RA, Wise WS, Zussman J (2004) Framework silicates : silica minerals, feldspathoids and the zeolites. In: Rockforming minerals, vol 4B, 2nd edn. The Geological Society, London, pp 493–523 Deering CD, Bachmann O, Vogel TA (2011) The Ammonia Tanks Tuff: Erupting a melt-rich rhyolite cap and its remobilized crystal cumulate. Earth Planet Sci Lett 310:518–525. https://​doi.​org/​10.​ 1016/J.​EPSL.​2011.​08.​032 Downs RT, Palmer DC (1994) The pressure behavior of c cristobalite. American Mineralogisl 79:9–14 Eichelberger J (2020) Distribution and Transport of Thermal Energy within Magma–Hydrothermal Systems. Geosciences 212(10):212. https://​doi.​org/​10.​3390/​GEOSC​IENCE​S1006​0212 Engvik AKP, Gerald AF, John D, Austrheim H (2008) Albitization of granitic rocks: the mechanism of replacement of oligoclase by albite. Can Mineral 46:1401–1415. https://​doi.​org/​10.​3749/​ canmin.​46.6.​1401 Ferrari L, Orozco-Esquivel T, Manea V, Manea M (2012) The dynamic history of the Trans-Mexican Volcanic Belt and the Mexico subduction zone. Tectonophysics 522:122–149 Foley ML, Miller CF, Gualda GAR (2020) Architecture of a Supersized Magma Chamber and Remobilization of its Basal Cumulate (Peach Spring Tuff, USA) P E T R O L O G Y. J Petrol 61:20. https://​doi.​org/​10.​1093/​petro​logy/​egaa0​20 Fritz J, Greshake A, Stöffler D (2005) Micro-Raman spectroscopy of plagioclase and maskelynite in Martian meteorites: Evidence of progressive shock metamorphism. Antarct Meteor Res 18:96 Furnes H (1978) Element mobility during palagonitization of a subglacial hyaloclastite in Iceland. Chem Geol 22:249–264. https://​doi.​org/​10.​1016/​0009-​2541(78)​90034-7 García-Palomo A, Macías JL, Tolson G et al (2002) Volcanic stratigraphy and geological evolution of the Apan region, eastcentral sector of the Trans-Mexican Volcanic Belt. Geofísica Internacional 41:133–150 García-Palomo A, Macías JL, Jiménez A et al (2018) NW-SE Pliocene-Quaternary extension in the Apan-Acoculco region, eastern Trans-Mexican Volcanic Belt. J Volcanol Geoth Res 349:240–255. https://​doi.​org/​10.​1016/J.​JVOLG​EORES.​2017.​ 11.​005 Giggenbach WF (1984) Mass transfer in hydrothermal alteration systems—A conceptual approach. Geochim Cosmochim Acta 48:2693–2711. https://​doi.​org/​10.​1016/​0016-​7037(84)​90317-X Giordano G, Caricchi L (2022) Determining the state of activity of transcrustal magmatic systems and their volcanoes. Ann Rev Earth Planet Sci 50:231-259. https://​doi.​org/​10.​1146/​ANNUR​ EV-​EARTH-​032320-​084733 Glazner AF, Johnson BR (2013) Late crystallization of K-feldspar and the paradox of megacrystic granites. Contrib Miner Petrol 166:777–799. https://​doi.​org/​10.​1007/​S00410-​013-​0914-1/​ FIGUR​ES/​16 Gómez-Alvarez F, Garduño-Monroy VH, Sosa-Ceballos G et al (2021) New constraints on tectonism and magmatism from the eastern sector of the Trans-Mexican Volcanic Belt (Chignahuapan Horst. J South Am Earth Sci, Puebla, México). https://​doi.​org/​ 10.​1016/J.​JSAMES.​2021.​103468 Gómez-Tuena A, Mori L, Straub SM (2018) Geochemical and petrological insights into the tectonic origin of the Transmexican 13 Contributions to Mineralogy and Petrology (2023) 178:63 Volcanic Belt. Earth Sci Rev 183:153–181. https://​doi.​org/​10.​ 1016/J.​EARSC​IREV.​2016.​12.​006 Gottsmann J, Lavallée Y, Martí J, Aguirre-Díaz G (2009) Magma– tectonic interaction and the eruption of silicic batholiths. Earth Planet Sci Lett 284:426–434. https://​doi.​org/​10.​1016/J.​EPSL.​ 2009.​05.​008 Gualda GAR, Ghiorso MS (2013) Geological note low-pressure origin of high-silica rhyolites and granites. J Geol 121:537–545. https://d​ oi.o​ rg/1​ 0.1​ 086/6​ 71395/A ​ SSET/I​ MAGES/L ​ ARGE/F ​ G4_​ ONLINE.​JPEG Gualda GAR, Gravley DM, Deering CD, Ghiorso MS (2019) Magma extraction pressures and the architecture of volcanic plumbing systems. Earth Planet Sci Lett 522:118–124. https://​doi.​org/​10.​ 1016/J.​EPSL.​2019.​06.​020 Harlow GE (1982) The anorthoclase structures: the effects of temperature and composition. Am Miner 67:975–996 Hay RL, Iijima A (1968) Nature and Origin of Palagonite Tuffs of the Honolulu Group on Oahu. Haw Mem Geol Soc Am 116:331– 376. https://​doi.​org/​10.​1130/​MEM116-​P331 Heaney PJ, Veblen DR, Post JE (1994) Structural disparities between chalcedony and macrocrystalline quartz. Am Miner 79:452–460 Hildreth W (1981) Gradients in silicic magma chambers: Implications for lithospheric magmatism. J Geophys Res Solid Earth 86:10153–10192. https://​doi.​org/​10.​1029/​JB086​IB11P​10153 Hildreth W, Wilson CJN (2007) Compositional Zoning of the Bishop Tuff. J Petrol 48:951–999. https://​d oi.​o rg/​1 0.​1 093/​P ETRO​ LOGY/​EGM007 Hirose T, Kihara K, Okuno M et al (2005) X-ray, DTA and Raman studies of monoclinic tridymite and its higher temperature orthorhombic modification with varying temperature. J Mineral Petrol Sci 100:55 Huber C, Bachmann O, Manga M (2010) Two Competing Effects of Volatiles on Heat Transfer in Crystal-rich Magmas: Thermal Insulation vs Defrosting. J Petrol 51:847–867. https://​doi.​org/​ 10.​1093/​PETRO​LOGY/​EGQ003 Ibrahim K, Hall A (1996) The authigenic zeolites of the Aritayn Volcaniclastic Formation, north-east Jordan. Miner Depos 31:514– 522. https://​doi.​org/​10.​1007/​BF001​96131/​METRI​CS Jiang YH, Jiang SY, Zhao KD, Ling HF (2006) Petrogenesis of Late Jurassic Qianlishan granites and mafic dykes, Southeast China: implications for a back-arc extension setting. Geol Mag 143:457– 474. https://​doi.​org/​10.​1017/​S0016​75680​50016​52 Jolles JSR, Lange RA (2021) Origin of Compositional Gradients with Temperature in the High-SiO2 Rhyolite Portion of the Bishop Tuff: Constraints on Mineral–Melt–Fluid Reactions in the Parental Mush. J Petrol 62:1–25. https://​doi.​org/​10.​1093/​PETRO​ LOGY/​EGAB0​87 Karlstrom L, Rudolph ML, Manga M (2012) Caldera size modulated by the yield stress within a crystal-rich magma reservoir. Nat Geosci. https://​doi.​org/​10.​1038/​ngeo1​453 Karsli O, Caran Ş, Dokuz A et al (2012) A-type granitoids from the Eastern Pontides, NE Turkey: Records for generation of hybrid A-type rocks in a subduction-related environment. Tectonophysics 530–531:208–224. https://​doi.​org/​10.​1016/J.​TECTO.​2011.​ 12.​030 Kihara K, Matsumoto T, Imamura M (1986) Structural change of orthorhombic-I tridymite with temperature: A study based on second-order thermal-vibrational parameters. Z Kristallogr Cryst Mater 177:27–38. https://​doi.​org/​10.​1524/​ZKRI.​1986.​177.​14.​27 Kimball KL (1990) Effects of hydrothermal alteration on the compositions of chromian spinels. Contrib Miner Petrol 105:337–346. https://​doi.​org/​10.​1007/​BF003​06543/​METRI​CS Kousehlar M, Weisenberger TB, Tutti F, Mirnejad H (2012) Fluid control on low-temperature mineral formation in volcanic rocks of Kahrizak. Iran Geofluids 12:295–311. https://​doi.​org/​10.​1111/​ GFL.​12001 Contributions to Mineralogy and Petrology (2023) 178:63 Lee S, Xu H (2019) Using powder XRD and pair distribution function to determine anisotropic atomic displacement parameters of orthorhombic tridymite and tetragonal cristobalite. Acta Crystallogr B Struct Sci Cryst Eng Mater 75:160–167. https://​doi.​org/​ 10.​1107/​S2052​52061​90009​33/​DK507​5SUP4.​PDF Lejeune AM, Richet P (1995) Rheology of crystal-bearing silicate melts: An experimental study at high viscosities. J Geophys Res Solid Earth 100:4215–4229. https://​doi.​org/​10.​1029/​94JB0​2985 Levien L, Prewitt CT, Weidner DJ (1980) Structure and elastic properties of quartz at pressure. Am Mineral 65(9-10):920–930 López-Hernández A, García-Estrada G, Aguirre-Díaz G et al (2009) Hydrothermal activity in the Tulancingo-Acoculco Caldera Complex, central Mexico: Exploratory studies. Geothermics 38:279– 293. https://​doi.​org/​10.​1016/J.​GEOTH​ERMICS.​2009.​05.​001 Lopez-Hernandez A, Castillo-Hernandez D (1997) Exploratory drilling at Acoculco, Puebla, Mexico: A hydrothermal system with only nonthermal manifestations (Conference) | OSTI.GOV. https://​ www.​osti.​gov/​biblio/​602731 Lubbers J, Deering C, Bachmann O (2020) Genesis of rhyolitic melts in the upper crust: Fractionation and remobilization of an intermediate cumulate at Lake City caldera. J Volcanol Geother Res, Colorado, USA. https://​doi.​org/​10.​1016/J.​JVOLG​EORES.​2019.​ 106750 Macías JL, Sheridan MF, Espíndola JM (1997) Reappraisal of the 1982 eruptions of El Chichón Volcano, Chiapas, Mexico: New data from proximal deposits. Bull Volcanol 58:459–471. https://​doi.​ org/​10.​1007/​S0044​50050​155/​METRI​CS Malfait WJ, Seifert R, Petitgirard S et al (2014) Supervolcano eruptions driven by melt buoyancy in large silicic magma chambers. Nat Geosci. https://​doi.​org/​10.​1038/​ngeo2​042 Marsh BD (1981) On the crystallinity, probability of occurrence, and rheology of lava and magma. Contrib Miner Petrol 78:85–98. https://​doi.​org/​10.​1007/​BF003​71146/​METRI​CS Marsh BD (2000) Reservoirs of magma and magma chambers. In: Sigurdsson H (ed) Encyclopedia of volcanoes. Academic, vol 1607. pp 191–206 Mathieu L (2018) Quantifying Hydrothermal Alteration: A Review of Methods. Geosciences 245(8):245. https://​doi.​org/​10.​3390/​ GEOSC​IENCE​S8070​245 Nadoll P, Mauk JL (2011) Wüstite in a hydrothermal silver-lead-zinc vein, Lucky Friday mine, Coeur d’Alene mining district, U.S.A. Am Miner 96:261–267. https://​doi.​org/​10.​2138/​AM.​2011.​3553/​ MACHI​NEREA​DABLE​CITAT​ION/​RIS Pamukcu AS, Carley TL, Gualda GAR et al (2013) The Evolution of the Peach Spring Giant Magma Body: Evidence from Accessory Mineral Textures and Compositions, Bulk Pumice and Glass Geochemistry, and Rhyolite-MELTS Modeling. J Petrol 54:1109–1148. https://​doi.​org/​10.​1093/​PETRO​LOGY/​EGT007 Pandarinath K, García-Soto AY, Santoyo E et al (2020) Mineralogical and geochemical changes due to hydrothermal alteration of the volcanic rocks at Acoculco geothermal system, Mexico. Geol J 55:6508–6526. https://​doi.​org/​10.​1002/​GJ.​3817 Pardo M, Suarez G (1995) Shape of the subducted Rivera and Cocos plates in southern Mexico: Seismic and tectonic implications. J Geophys Res Solid Earth 100:12357–12373. https://​doi.​org/​10.​ 1029/​95JB0​0919 Parolari M, Gómez-Tuena A, Errázuriz-Henao C, Cavazos-Tovar JG (2021) Orogenic andesites and their link to the continental rock cycle. Lithos. https://​doi.​org/​10.​1016/J.​LITHOS.​2020.​105958 Pauly BD, Schiffman P, Zierenberg RA et al (2011) Environmental and chemical controls on palagonitization. Geochem Geophys Geosyst 12:12017. https://​doi.​org/​10.​1029/​2011G​C0036​39 Peiffer L, Bernard-Romero R, Mazot A et al (2014) Fluid geochemistry and soil gas fluxes (CO2–CH4–H2S) at a promissory Hot Dry Rock Geothermal System: The Acoculco caldera, Mexico. Page 23 of 24 63 J Volcanol Geoth Res 284:122–137. https://​doi.​org/​10.​1016/J.​ JVOLG​EORES.​2014.​07.​019 Pérez-Campos X, Kim YH, Husker A et al (2008) Horizontal subduction and truncation of the Cocos Plate beneath central Mexico. Geophys Res Lett. https://​doi.​org/​10.​1029/​2008G​L0351​27 Pérez-Orozco JD, Sosa-Ceballos G, Macías JL (2021) Tectonic and magmatic controls on the evolution of post-collapse volcanism. Lithos, Insights from the Acoculco Caldera Complex, Puebla, México. https://​doi.​org/​10.​1016/J.​LITHOS.​2020.​105878 Perton M, Maldonado Hernández LT, Figueroa-Soto A et al (2022) The magmatic plumbing system of the Acoculco volcanic complex (Mexico) revealed by ambient noise tomography. J Volcanol Geother Res. https://​doi.​org/​10.​1016/J.​JVOLG​EORES.​2022.​ 107704 Petford N, Lister JR, Kerr RC (1994) The ascent of felsic magmas in dykes. Lithos 32:161–168. https://​doi.​org/​10.​1016/​0024-​ 4937(94)​90028-0 Petford N (2003) Rheology of Granitic Magmas During Ascent and Emplacement. Ann Rev Earth Planet Sci 31(1):399–427 Post JE, Bish DL (1989) Rietveld refinement of crystal structures using powder X-ray diffraction data. Mod Powder Diffr 20:277–30 Repstock A, Breitkreuz C, Lapp M, Schulz B (2018) Voluminous and crystal-rich igneous rocks of the Permian Wurzen volcanic system, northern Saxony, Germany: physical volcanology and geochemical characterization. Int J Earth Sci 107:1485–1513. https://​doi.​org/​10.​1007/​S00531-​017-​1554-X/​FIGUR​ES/​14 Rollinson HR (2014) Using Geochemical Data: Evaluation. Using Geochemical Data, Presentation, Interpretation. https://​doi.​org/​ 10.​4324/​97813​15845​548 Rowe MC, Ellis BS, Lindeberg A (2012) Quantifying crystallization and devitrification of rhyolites by means of X-ray diffraction and electron microprobe analysis. Am Miner 97:1685–1699. https://​ doi.​org/​10.​2138/​AM.​2012.​4006 Sánchez-Córdova MM, Canet C, Rodríguez-Díaz A et al (2020) Waterrock interactions in the Acoculco geothermal system, eastern Mexico: Insights from paragenesis and elemental mass-balance. Geochemistry. https://​doi.​org/​10.​1016/J.​CHEMER.​2019.​06.​003 Sas M, Shane P, Kuritani T et al (2021) Mush, Melts and Metasediments: a History of Rhyolites from the Okataina Volcanic Centre, New Zealand, as Captured in Plagioclase. J Petrol 62:1–26. https://​doi.​org/​10.​1093/​PETRO​LOGY/​EGAB0​38 Scaillet B, Holtz F, Pichavant M (2016) Experimental Constraints on the Formation of Silicic Magmas. Elements 12:109–114. https://​ doi.​org/​10.​2113/​GSELE​MENTS.​12.2.​109 Schwartz GM (1959) Hydrothermal alteration. Econ Geol 54:161–183. https://​doi.​org/​10.​2113/​GSECO​NGEO.​54.2.​161 Seifert W, Thomas R, Rhede D et al (2010) Origin of coexisting wüstite, Mg–Fe and REE phosphate minerals in graphite-bearing fluorapatite from the Rumburk granite. Eur J Mineral 22:495– 507. https://​doi.​org/​10.​1127/​0935-​1221/​2010/​0022-​2034 Sliwinski JT, Bachmann O, Dungan MA et al (2017) Rapid pre-eruptive thermal rejuvenation in a large silicic magma body: the case of the Masonic Park Tuff, Southern Rocky Mountain volcanic field, CO, USA. Contrib Miner Petrol 172:1–20. https://​doi.​org/​ 10.​1007/​S00410-​017-​1351-3/​FIGUR​ES/9 Sohn YK, Son M, Jeong JO, Jeon YM (2009) Eruption and emplacement of a laterally extensive, crystal-rich, and pumice-free ignimbrite (the Cretaceous Kusandong Tuff, Korea). Sediment Geol 220:190–203. https://​doi.​org/​10.​1016/J.​SEDGEO.​2009.​04.​020 Sosa-Ceballos G, Macías JL, Avellán DR et al (2018) The Acoculco Caldera Complex magmas: Genesis, evolution and relation with the Acoculco geothermal system. J Volcanol Geoth Res 358:288– 306. https://​doi.​org/​10.​1016/J.​JVOLG​EORES.​2018.​06.​002 Sosa-Ceballos G, Boijseauneau-Lopez ME, Perez-Orozco JD et al (2021) Silicic magmas in the Michoacan-Guanajuato volcanic 13 63 Page 24 of 24 field: An overview of plumbing systems, crustal storage, and genetic processes. Rev Mex Cienc Geol 38:210–225 Stroncik NA, Schmincke HU (2001) Evolution of palagonite: Crystallization, chemical changes, and element budget. GeochemGeophy Geosyst. https://​doi.​org/​10.​1029/​2000G​C0001​02 Sun SS, McDonough WF (1989) Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. Geol Soc Spec Publ 42:313–345. https://d​ oi.o​ rg/1​ 0.1​ 144/​ GSL.​SP.​1989.​042.​01.​19 Szymanowski D, Wotzlaw JF, Ellis BS et al (2017) Protracted nearsolidus storage and pre-eruptive rejuvenation of large magma reservoirs. Nat Geosci. https://​doi.​org/​10.​1038/​ngeo3​020 Tavazzani L, Peres S, Sinigoi S et al (2020) Timescales and Mechanisms of Crystal-mush Rejuvenation and Melt Extraction Recorded in Permian Plutonic and Volcanic Rocks of the Sesia Magmatic System (Southern Alps, Italy). J Petrol 61:49. https://​ doi.​org/​10.​1093/​PETRO​LOGY/​EGAA0​49 Toby BH, Von Dreele RB (2013) GSAS-II: the genesis of a modern open-source all purpose crystallography software package. J Appl Crystallogr 46(2):544–549. https://​doi.​org/​10.​1107/​S0021​ 88981​30035​31 Townsend M, Huber C (2020) A critical magma chamber size for volcanic eruptions. Geology 48:431–435. https://​doi.​org/​10.​1130/​ G47045.1 Troch J, Huber C, Bachmann O (2022) The physical and chemical evolution of magmatic fluids in near-solidus silicic magma reservoirs: Implications for the formation of pegmatites. Am Miner 107:190–205. https://​doi.​org/​10.​2138/​AM-​2021-​7915/​ASSET/​ GRAPH​IC/J_​AM-​2021-​7915_​FIG_​004.​JPG Vigneresse JL, Barbey P, Cuney M (1996) Rheological transitions during partial melting and crystallization with application to felsic magma segregation and transfer. J Petrol 37:1579–1600. https://​ doi.​org/​10.​1093/​PETRO​LOGY/​37.6.​1579 Walton AW, Schiffman P, Walton AW, Schiffman P (2003) Alteration of hyaloclastites in the HSDP 2 Phase 1 Drill Core 1. Description 13 Contributions to Mineralogy and Petrology (2023) 178:63 and paragenesis. Geochem Geophys Geosys. https://​doi.​org/​10.​ 1029/​2002G​C0003​68 Watts KE, John DA, Colgan JP et al (2016) Probing the VolcanicPlutonic Connection and the Genesis of Crystal-rich Rhyolite in a Deeply Dissected Supervolcano in the Nevada Great Basin: Source of the Late Eocene Caetano Tuff. J Petrol 57:1599–1644. https://​doi.​org/​10.​1093/​PETRO​LOGY/​EGW051 Weinberg RF (1996) Ascent mechanism of felsic magmas: news and views. Earth Environ Sci Trans R Soc Edinb 87:95–103. https://​ doi.​org/​10.​1017/​S0263​59330​00065​19 Whalen JB, Currie KL, Chappell BW (1987) A-type granites: geochemical characteristics, discrimination and petrogenesis. Contrib Miner Petrol 95:407–419 Whitney JA, Naldrett AJ, Robertson JM (1989) Origin and Evolution of Silicic Magmas. Ore Dep Ass Magmas. https://​doi.​org/​10.​ 5382/​REV.​04.​11.​k, Yáñez-García C (1980) Informe geológico del proyecto geotérmico Los Humeros-Las Derrumbadas, estados de Puebla y Veracruz. CFE-GPG Int Rep 01:80 Zhang JH, Yang JH, Chen JY et al (2018) Genesis of late Early Cretaceous high-silica rhyolites in eastern Zhejiang Province, southeast China: A crystal mush origin with mantle input. Lithos 296– 299:482–495. https://​doi.​org/​10.​1016/J.​LITHOS.​2017.​11.​026 Zhou Z, Fyfe WS (1989) Palagonitiz tion of basaltic glass from DSDP Site 335, Leg 37: Textures, chemical composition, and mechanism of formation. Am Miner 74:1045–1053 Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.