Reply to Norini and Groppelli’s comment on “Estimating the depth and evolution of intrusions at resurgent calderas: Los Humeros (Mexico)” by Urbani et al. (2020)

Structural studies in active caldera systems are widely used in geothermal exploration to reconstruct volcanological conceptual models. Active calderas are difficult settings to perform such studies mostly because of the highly dynamic environment, dominated by fast accumulation of primary and secondary volcanic deposits, the variable and transient rheology of the shallow volcanic pile, and the continuous feedbacks between faulting, secondary porosity creation, and geothermal fluid circulation, alteration and cementation that tend to obliterate the tectonic deformation structures. In addition, deformation structures can be also caused by nearand far-field stress regimes, which include magmatic intrusions at various depths, the evolving topography and regional tectonics. A lack of consideration of all these factors may severely underpin the reliability of structural studies. By rebutting and providing a detailed discussion of all the points raised by the comment of Norini and Groppelli (2020) to the Urbani et al. (2020) paper, we take the opportunity to specify the scientific rationale of our structural fieldwork and strengthen its relevance for geothermal exploration and exploitation in active caldera geothermal systems in general and, particularly, for the Holocene history of deformation and geothermal circulation in the Los Humeros caldera. At the same time, we identify several major flaws in the approach and results presented in Norini and Groppelli (2020), such as (1) the lack of an appropriate ranking of the deformation structures considering an inventory method for structural analysis; (2) the misinterpretation and misquoting of Urbani et al. (2020) and other relevant scientific literature; and (3) irrelevant and contradictory statements within their comment.


Introduction
Structural studies in active calderas provide key elements for the exploration of geothermal systems and greatly contribute to the development of conceptual models for their exploitation. We herein reply to the comment by Norini and Groppelli (2020) on our paper Urbani et al. (2020) entitled "Estimating the depth and evolution of intrusions at resurgent calderas: Los Humeros (Mexico)", giving us the opportunity to further discuss our approach, results and the proposed reconstruction of the Holocene volcano-tectonic evolution of the Los Humeros volcanic complex (LHVC; Mexico) and its relevance for our understanding of the active geothermal system.

Reply to the criticism raised in the comment
The Norini and Groppelli (2020) criticism of the  paper concentrates on three main aspects: (1) supposed lack of structural field data and supposed geometric and structural inconsistency between the Holocene deformation and the proposed subsurface model; (2) supposed lack of validation of the obtained results with those available from Published by Copernicus Publications on behalf of the European Geosciences Union.
well log data; and (3) supposed contradictions with the available stratigraphic reconstruction and radiometric ages.
Here follows a point-by-point discussion of the critical points raised in Norini and Groppelli (2020).

Las Papas and Las Viboras structures: inventory
vs. selective method of structural analysis in active volcanic areas Norini and Groppelli (2020) question the reinterpretation made by Urbani et al. (2020) of Las Papas and Las Viboras structures as presently inactive morphological scarps, showing small-scale faults in the Cuicuiltic Member (see Fig. 2 of Norini and Groppelli, 2020). They are also critical of the supposedly few data presented. In contrast to the inventory method followed by Norini and Groppelli (2020), where all faults are mixed together without any hierarchy and discussed as unweighted data, in Urbani et al. (2020) we followed a selective method of structural analysis. By using this method, we ranked faults by adopting the following criteria in the field: (i) the topographic expression of the fault, (ii) the amount of displacement of individual fault strands and/or fault systems; (iii) the along-strike persistence of the fault trace; (iv) the presence of clear kinematic indicators; (v) the presence or absence of associated hydrothermal alteration, and (vi) the relative age with respect to the Holocene intracaldera Cuicuiltic Member fall deposit. Regarding this last criterion, the Cuicuiltic Member is regarded as a reference space-time marker to discriminate between faults older or younger than 7.3 ka due to its well-known age and distribution within the caldera (Dávila-Harris and Carrasco-Núñez et al., 2017a). Accordingly, in Urbani et al. (2020) we described only selected faults that exhibited clear metre-scale offsets, consistent lateral extent and/or evidence of hydrothermal alteration. We therefore strongly reject that our dataset is poor, because it deals with structures that, based on the above-listed criteria, allowed us to discriminate between and rank volcano-tectonic and hydrothermal processes, which are instead missed by the inventory method of an unweighted fault dataset adopted by Norini and Groppelli (2020). In terms of geothermal exploration, the faults presented in Norini and Groppelli (2020) are indeed questionable in terms of relevance. For example, the LH17106 and LH62 outcrops shown in Fig. 2c of Norini and Groppelli (2020) are in the same location of outcrop LH-08 shown in Fig. 5c of Urbani et al. (2020), where an erosional surface at the top of undeformed and unaltered pyroclastic deposits of the Xoxoctic Tuff, blanketed by the Cuicuiltic Member, is clearly visible. Figure 1a of this reply shows the same outcrop, where the erosional unconformity at the top of the underlying subhorizontal pyroclastics is sutured by the Cuicuiltic Member fall deposits. The large-scale blanketing geometry of the unaltered Cuicuiltic Member fall deposits across the Las Papas scarp is well visible in Fig. 1b of this reply. This indicates that Las Papas is currently an inactive morphological scarp without evidence of hydrothermal alteration. Whether or not this scarp was associated in origin (prior to 7.3 ka) with a fault trace is not evident in the field nor relevant for our study, focused on present-day relationships between faulting and geothermal circulation. It is noteworthy that even Norini et al. (2019;see Sects. 4 and 6.2) raise doubts regarding the relevance of the Las Papas structure within the Los Humeros geothermal field, suggesting a weak or no connection with the geothermal reservoir. The same holds for the Las Viboras structure. In our opinion, the small-scale faults shown by Norini and Groppelli (2020) in their Fig. 2d-e are not at all compelling and may be alternatively interpreted as small-scale normal faults generated by near-field (local) stresses affecting unlithified material (e.g. Wernicke and Birchfiel, 1982;Bridgwater et al., 1985;Branney and Kokelaar, 1994;Gao et al., 2020;Yang and Van Loon, 2016). In particular, Fig. 2d of Norini and Groppelli (2020) is unclear, whereas Fig. 2e does not show any evidence of displacement of the lower white and black pumice beds as well as on the upper brown beds, suggesting an intra-formational readjustment (Fig. 1c) rather than a fault. (e.g. Wiggers, 1975, 1976;Branney and Kokelaar, 1994). Norini and Groppelli (2020) fail to discuss any possible alternative origin for their smallscale faults, which, considering their location is in an active caldera floor, severely impinges on the reliability of the inventory dataset presented in Norini and Groppelli (2020) and its relevance for geothermal studies. Reinterpreting the small-scale offsets shown in Fig. 2 on Norini and Groppelli (2020) as minor gravitational structures (i) would justify why they have no connection with the geothermal circulation nor with any thermal anomaly and (ii) clarifies to the reader why the Urbani et al. (2020) paper instead focused only on faults that were ranked as being of first-order importance in terms of displacement, persistence in the field and age of the structurally controlled hydrothermal fluid circulation.
In summary, we (i) question the use of the inventory method for the structural analysis adopted by Norini and Groppelli (2020) when applied to active calderas, which may prove to be inappropriate and unable to discriminate firstorder, deep and geothermally relevant fault systems from small-scale, soft-state deformation structures that are also common in intracaldera domains (e.g. Branney and Kokelaar, 1994), and (ii) consider the selective method of structural analysis used by Urbani et al. (2020) as fully appropriate to rank deformation structures ( Fig. 1a-i) when the aim of structural fieldwork is to characterize the surface deformation related to the recent activity of a caldera, in order to constrain the morphotectonic fingerprints of the resurgence, and to evaluate its source and areal extent and, even more importantly, its relevance for the active geothermal system.  Norini and Groppelli (2020) state that "active or fossil alteration does not always allow for identifying faults or the age of faulting, because this also depends on their depth, life span of the hydrothermal system, spatial relationships, and fluid paths along primary permeability and fracture zones (e.g. Bonali et al., 2015;Giordano et al., 2016)". Although geothermal fields may also develop associated with primary porosity (Hurter and Shellschmidt, 2003), the two studies cited by Norini and Groppelli (2020) are in no way at odds with Urbani et al. (2020). The work of Bonali et al. (2016), on the active tectonics at Copahue (Argentina) points out that active fault systems in volcanic settings are responsible for driving hot fluids to the surface. Similarly, the work of Giordano et al. (2013Giordano et al. ( , 2016 on the Tocomar geothermal field (Puna Plateau, Argentina), investigated the evidence of a geothermal field based on the overlapping distribution of hot springs and active fault systems. We thank Norini and Groppelli (2020) for bringing to our attention these two very interesting papers because, along with mainstream literature (e.g. Karaoglu et al., 2019), they clearly indicate that hydrothermal fluids and associated alteration in volcanic settings are driven and controlled by active fault systems. The relationship between faulting and fluid circulation is well established also in exhumed systems, where it is clear how fault permeability is affected by the interplay between the far-field regional stress field and the near-field stress regime (e.g. Rossetti et al., 2011;Olvera Garcia et al., 2020). Therefore, the cited papers support the proposal of Urbani et al. (2020) to use the distribution and intensity of the hydrothermal alteration within the 7.3 ka Cuicuiltic Member marker beds that ubiquitously blanket the caldera floor and all the fault scarps as a valid space-time marker in the field to discriminate between active vs. inactive fault segments controlling the upwelling of geothermal fluids ( Fig. 1d-i).
In contrast to Norini and Groppelli (2020), we conclude that, in agreement with authoritative literature, hydrothermal alteration follows the space-time distribution of structurally controlled (fault-induced) secondary permeability pathways and that its distribution should be used, along with measured fault displacements, persistence and (relative) age, as an indication of fault activity and ranking for geothermal purposes. At Los Potreros, the presence of the 7.3 ka Cuicuiltic Member marker bed allowed us to track the type and intensity of deformation and its association with fluid circulation and alteration in space and time.

Surface thermal anomalies
Norini and Groppelli (2020) state that "The Maxtaloya fault trace is coincident with a sharp thermal anomaly identified by Norini et al. (2015) (Fig. 3). U2020 did not consider this positive (warm) anomaly when they discussed the thermal remote sensing results published by Norini et al. (2015) (Sect. 5.3 in U2020).". This statement is not correct, because we did discuss the published remote-sensing results as clearly written in Sect. 5.3 of Urbani et al. (2020). Moreover, the sharp and narrow temperature peaks, described by Norini et al. (2015), are spatially coincident with the Los Humeros and Loma Blanca faults described by Urbani et al. (2020). This

Identification and geometry of uplifted areas: topographic data and structural mapping
Norini and Groppelli (2020) criticize the location and geometry of the three uplifted areas of Los Humeros, Loma Blanca and Arroyo Grande identified by Urbani et al. (2020). However, in the topographic profiles across the bulges shown by Norini and Groppelli (2020) in their Fig. 4a-b, the uplifted areas at Loma Blanca, Arroyo Grande and Los Humeros are well visible and the existence of such bulges is unquestionable. Therefore, it is unclear on what basis Norini and Groppelli (2020) question the existence of such uplifted areas.
The asymmetry (Arroyo Grande) and tilt of the uplifted areas (Loma Blanca) detailed by Norini and Groppelli (2020) are in no way contradictory to the Urbani et al. (2020)  Helens (Lipman et al., 1981). Despite being stimulating for future work, the investigation of the exact origin of the bulge shapes was far beyond the scope of Urbani et al. (2020), who, for this reason, maintained the same initial and simplified geometric configuration for the analogue models. Therefore, the comment made by Norini and Groppelli (2020) is not relevant for the discussion presented in Urbani et al. (2020).

Apical depression of bulges
The model proposed by Urbani et al. (2020) predicts the formation of an apical depression at the top of a bulge induced by a shallow intrusion. Norini and Groppelli (2020) state that the topography of natural bulges identified by  does not show well-defined apical depressions in the asymmetric Arroyo Grande and Los Humeros uplifted areas, contradicting the model results. Analogue modelling in Urbani et al. (2020) injects symmetric intrusions, a condition appropriate for the morphology of the Loma Blanca bulge, where the apical depression is very well evident (Fig. 2) and measured in the field (Fig. 6f in Urbani et al., 2020). The Arroyo Grande and Los Humeros bulges are instead asymmetrical and likely developed as trapdoor uplifts (thus without apical depression) associated with asymmetric intrusions and with a deformation amount much larger than that at Loma Blanca and that considered in the analogue models. Therefore, the comment made by Norini and Groppelli (2020) is incorrect regarding the Loma Blanca bulge and not relevant in the other two cases, and, consequently, the predictive value of the model proposed by Urbani et al. (2020) is in no way compromised. Indeed, the experiments were designed to ensure the formation of an apical depression not considering the trapdoor uplift and faulting, and the asymmetric bulges are not directly modelled.

Reverse faults bounding uplifted areas
Norini and Groppelli (2020) state that Urbani et al. (2020) do not provide independent validation of the proposed multiple magmatic intrusion model, such as field evidence of reverse faults predicted by the analogue modelling results. Exposure of faults in active caldera floors depends on many factors: (i) elastic vs. anelastic response to deformation source, its location, intensity and duration; (ii) nucleation depth and surface propagation; and (iii) burial vs. exhumation rates. Therefore, while reverse faults accompanying both largescale resurgence and local uplifts are expected by analogue models, the scarcity of visible and measurable reverse faulting in no way disproves the hypothesis related to the intrusion of cryptodomes and resurgence (Bonanza, Lipman et al., 2015;Long Valley, Hildreth et al., 2017;Kutcharo, Goto and McPhie, 2018). Therefore, the statement by Norini and Groppelli (2020) claiming that the locations of such reverse faults "are a fundamental feature of their model" is incorrect.
In addition, Norini and Groppelli (2020) show the traces of inferred reverse faults at the periphery of the Loma Blanca bulge, just where the Urbani et al. (2020) model predicts (see Fig. 2), making Norini and Groppelli's (2020) statements contradictory.
2.5 Validation of the proposed model: geothermal wells log data

Lithology of intrusions
Norini and Groppelli (2020) invoke the thermal profile and the stratigraphy of just one well log (the H4 well, drilled at the top of the Loma Blanca bulge) to claim the lack of validation of the models proposed in Urbani et al. (2020). First, we would like to emphasize that the proposed reinterpretation of the subsurface stratigraphy presented in Urbani et al. (2020) is not just based on the H4 well. A great part of Sect. 2 ("Geological-structural setting") and Fig. 2a-b presented in Urbani et al. (2020) discuss in detail the published data from 12 well logs (including the H4 well log) as presented in Arellano et al. (2003) and in Carrasco-Núñez et al. (2017a, b). The model evaluation of the intrusion depths, as derived from the equation of Brothelande and Merle (2015), are valid within the modelling assumptions and are within the depth range of some rhyolitic-dacitic bodies drilled in geothermal wells, wherein they are simply described texturally as lavas (Carrasco-Núñez et al., 2017b, and references therein). The lithologic definition of "lava" is associated with aphanitic to phaneritic textures that are not only restricted to subaerial environments and may be impossible to distinguish from textures of sub-volcanic and hypabyssal bodies. Hypabyssal rocks are characterized by a rapid cooling and their textures are fine grained or glassy and mostly resemble those of volcanic rocks (Philpotts and Ague, 2009 , 1995). We therefore reject the criticism by Norini and Groppelli (2020) as the interpretation is based not only on the H4 well but also on stratigraphic reconstructions derived from 12 published well logs.

Geometry of caldera fill
The reinterpretation proposed by Urbani et al. (2020) of some of the rhyolitic-dacitic bodies of the Los Potreros subsurface as hypabyssal intrusives is not simply based on their lithology but also on their geometry, stratigraphic position and the whole geometry of the caldera fill; all are elements neither considered nor discussed in Norini and Groppelli (2020). When correlating the stratigraphic well logs, Urbani et al. (2020) documented (in Sect. 2 at p. 530 and Fig. 2) the irregular geometry of both the top of the Xaltipan intracaldera ignimbrite and the post-caldera units, as well as the lack of a clear topography-filling geometry: a stratigraphic setting that can hardly be reconciled with an intracaldera setting unless the emplacement of intrusive bodies has occurred in the shallow crust. It is noteworthy that the main geometric anomalies of the caldera fill appear directly in correspondence with the possible location of a felsic intrusion. For example, a 600 m thick rhyolitic-dacitic body showing all the petrographic features of a hypabyssal intrusion is reported to the west of Arroyo Grande in the H20 well at 470-1060 m depth from the surface (see also Carrasco-Núñez et al., 2017b). It is located at the top of the pre-caldera andesites, intrudes both the intracaldera and the post-caldera units, and shows no lateral continuity. Similar felsic bodies were also drilled in the H5, H26, H19 and H25 wells. Furthermore, Norini and Groppelli (2020) 2020) for what concerns the subsurface stratigraphy of the Los Potreros intracaldera domain, and therefore we reject the criticism of Norini and Groppelli (2020).

Thermal gradient
The statement by Norini and Groppelli (2020) Carrasco-Núñez et al., 2017b;Urbani et al., 2020). Measured downhole temperature profiles for well H4 (Torres-Rodriguez, 1995;, 1998Martínez-Serrano, 2002)  It should be noted in Fig. 5 that the obsidian dome is also the site of a sample whose magmatic zircon crystallization age was U-Th dated to 44.8 ka in Carrasco-Nuñez et al. (2017). At the same time, the dome is largely not covered by the Cuicuiltic Member. As the dated sample was taken at the base of the dome, this can be interpreted in various ways: one is that the dome is poly-phased and its upper part is younger than 7.3 ka, resting above the Cuicuiltic Member; the other is that the dome was exhumed after 7.3 ka. The emplacement and exhumation of the obsidian dome and the nearby faulting of the Cuicuiltic Member by tens of metres of displacement at the site of the Los Humeros Fault indicates that this section of the fault was active later than 7.3 ka. By contrast, the fault displacement drastically reduces southward along the Maxtaloya fault. This in our opinion supports our interpretation of the Maxtaloya-Los Humeros faults as segmented and diachronous during the Holocene, in agreement with the Urbani et al. (2020) interpretation.

Recent history of caldera floor uplift
Norini and Groppelli (2020) misquote Urbani et al. (2020), attributing to them the interpretation of a northward shift in volcanic activity within the Los Potreros caldera, which was neither declared nor intended in the paper.  simply summarize field evidence stating "the recent (post-caldera collapse) uplift in the Los Potreros caldera moved progressively northwards, from the south and northeastern sector of the caldera towards the north along the Los Humeros and Loma Blanca scarps". Urbani et al. (2020) did not discuss the causes of such northward shift and even less attributed it to a shift in "the volcanic feeding system" as erroneously and unjustifiably reported by Norini and Groppelli (2020). The fate of a magma intrusion, i.e. whether it will erupt or stop in the crust, depends on many factors, such as its buoyancy (density contrast with host rocks), the initial gas content, the rise speed and style of decompressiondegassing, and the rheology of the magma and of the intruded crust, including its layering, structure and so forth. The evolution over time and space of intrusions in a caldera may see different phases and have many different causes, partly depending on feedbacks existing between the evolving configuration of the magmatic plumbing system and the evolving rheology and structure of the caldera roof rocks. The plumbing system of the last 10 kyr at Los Humeros has been reconstructed in detail by Lucci et al. (2020). This study documents a multi-storey magmatic complex, which allows the eruption along the Los Potreros caldera floor of both deeply sourced (>30 km) olivine basalts and shallow differentiated (<3 km) felsic trachytes and rhyolites. The results of Lucci et al. (2020), curiously neither cited nor discussed by Norini and Groppelli (2020), highlight the absence of the classic large-volume, single magma chamber and suggest that the activation of magma sources at different depths appears not to have followed any specific pattern during the Holocene. The present absence below Los Humeros of a single large magma chamber or crystal mush able to form a rheological barrier to the rise of basalts directly from lower crustal depths severely impinges upon the model of classic resurgence supported by Norini and Groppelli (2020), which requires the existence of a voluminous viscous layer accommodating magma recharge and acting as a pressure source for resurgence (Galetto et al., 2017).

Summary and implications for the Los Humeros geothermal system
Understanding the anatomy of magma plumbing systems of active volcanic systems, from deeper reservoirs to subsurface ephemeral batches, is crucial to define temperature, depth and geometry of the heat sources for geothermal exploration. The Pleistocene-Holocene Los Humeros volcanic complex (LHVC), located in the eastern Trans-Mexican Volcanic Belt (central Mexico), represents one of the most important exploited geothermal fields in Mexico, with ca. 95 MW of produced electricity. Geological investigations at LHVC started at the end of the 1970s and culminated in the production of (i) the first comprehensive geological map ( Fig. 6a; after Ferriz and Mahood, 1984), (ii) a structural map of the intracaldera domain ( Fig. 6b; after Alcantara et al., 1988), (iii) the proposal of a petrological conceptual model of the plumbing system made of a single voluminous (ca. 1200 km 3 ) melt-dominated and zoned magma chamber at shallow depths (ca. 5 km, Fig. 6c; after Verma, 1985), and (iv) the proposal of an inflation-deflation caldera episodic and cyclic model ( Fig. 6d; after Campos-Enriquez and Arredondo-Fragoso, 1992) connected to the activity of the single voluminous conventional magma chamber of Verma (1985). Since these main studies, and up to the most recent published work, the understanding of the Los Humeros volcanic complex has been incremental, never questioning the consolidated model of the single zoned magma chamber where all petrologic, volcanologic and deformation processes originate (i.e., Ferriz and Mahood, 1984;Alcantara et al., 1988;Verma, 1985;Campos-Enriquez and Arredondo-Fragoso, 1992). Structural work by Norini et al. (2015Norini et al. ( , 2019, produced updates and refined versions ( Fig. 6e and f) of the original structural map by Alcantara et al. (1988). Based on the assumption of the existence of an active single voluminous magma chamber as proposed in the early 1980s (Verma, 1985), post-caldera deformation has been interpreted uniquely as due to a classic mechanism of resurgence (e.g. Fig. 6g; after Norini et al., 2019) that much (or completely) resembles the first proposal of Campos-Enriquez and Arredondo-Fragoso (1992). However, such a conceptual model is now under stress as the geothermal anomalies appear very localized, mainly confined along the NNW-SSEtrending Maxtaloya-Los Humeros-Loma Blanca-Los Conejos corridor and corresponding to the almost unique, narrow, rapidly declining, thermal anomaly recognized within the Los Potreros caldera (Norini et al., 2015;Peiffer et al., 2018;Jentsch et al., 2020). This geothermal configuration is  Verma (1985). (d) Schematic representation of the evolution of Los Humeros volcanic complex by Campos-Enriquez and Arredondo-Fragoso (1992) where magmatism, eruptive styles, inflation and deflation phenomena are all correlated to the activity of the single voluminous and shallow-seated magma chamber of Verma (1985).
(e) Morphostructural map of the Los Potreros caldera with interpretation of the sectorial resurgence as proposed by Norini et al. (2015). reflected in the relatively low number of productive geothermal wells (ca. 25 out of ca. 60; Gutiérrez-Negrín et al., 2019, 2020 but is difficult to reconcile with the existence of a single, deep-seated, large-volume magmatic source that should instead generate widespread and sustained thermal anomalies in the caldera floor, such as in active resurgent calderas like Ischia (Carlino et al., 2014). A step change in paradigm in the reconstruction of the Holocene magmatic plumbing system at Los Humeros has been proposed in Lucci et al. (2020) and GEMex (2019b), which despite their important implications for the understanding of the present-day geothermal system were not even cited by Norini and Groppelli (2020). Lucci et al. (2020) carried out a thermobarometric study of all exposed Holocene lavas, demonstrating that the scattered intracaldera monogenetic activity reflects the ascent of magmas from basaltic to trachytic in composition from sources located at depths comprised between >30 km (basalts) to <3 km (trachytes) and for variably evolved compositions with complex histories of ascent and stalling at various depths, depicting a multi-storey  Lucci et al. (2020), of the magmatic plumbing system feeding the Los Humeros post-caldera stage activity, beneath the Los Humeros Caldera, as derived by pressure-temperature estimates obtained from mineral-liquid thermobarometry models. The model is integrated with the crustal structure (see Lucci et al., 2020, for further explanations). (b) Schematic model, by Urbani et al. (2020), of the evolution and of the subsurface structure of the Los Potreros caldera floor. Multiple magmatic intrusions located at a relatively shallow depth (<1 km) are responsible for the localized bulging of the caldera floor (Arroyo Grande, Los Humeros and Loma Blanca uplifted areas). The Cuicuiltic Member eruption is assumed as a time marker in the evolution of the intracaldera domain.
plumbing system (e.g. Cashman and Giordano, 2014;Cashman et al., 2017;Sparks et al., 2019). This innovative reconstruction of the plumbing system suggests that the largevolume magma chamber at 5 km depth that produced the caldera collapses at the time of the eruption of the Xaltipan ignimbrite (164 ka) and Zaragoza ignimbrite (69 ka) does not exist anymore as a single melt-dominated volume, allowing the rise to the surface of mantle-derived magmas as well as differentiation processes at various depths of small batches of magma through the entire crust. Urbani et al. (2020) performed structural fieldwork based on a selective method approach combined with analogue models, showing that, at least during the Holocene, the classic resurgence model (e.g. Norini et al., 2019) does not explain the fault ranks and the spatio-temporal evolution of the deformation and alteration. This change in paradigm at Los Humeros implies (i) the inadequacy of the hypothesis of a single, large and voluminous shallow magmatic chamber homogeneously distributed beneath the caldera; (ii) the proposal of an innovative scenario, characterized by a complex magmatic plumbing system vertically distributed across the entire crust, from a deeper resi-dence zone for basalts to a shallower magmatic plexus made of small single-charge ephemeral pockets of heterogeneous magmas localized beneath the Los Humeros nested caldera (Fig. 7a;after Lucci et al., 2020); and (iii) the interpretation of the recent deformation at the Los Humeros volcanic complex not as a classical resurgence associated with the bulk inflation of a deep magma reservoir but as the response to the ascent and emplacement of multiple, small-volume magma batches at shallow crustal conditions (<1 km depth) ( Fig. 7b; after Urbani et al., 2020). These results have important consequences for the geothermal exploration and exploitation and the siting of future geothermal wells, where shallow magma bodies can act as scattered and localized short-lived heat sources complicating the pattern of isotherms related to deeper reservoirs. At the same time, the evidence of absence during the Holocene of an actively recharged, large and meltdominated magma chamber located at 5 km depth (i.e. the Xaltipan and Zaragoza magma chamber) may help understand the localized nature of the thermal anomalies at Los Humeros.
Data availability. All the data presented in this paper are available upon request.
Author contributions. SU, GG, FL, FR and GCN were equally responsible for the conceptualization of the study (research planning) and the investigation (field analysis, review of well logs and stratigraphic data). FL and SU wrote the original draft of the paper. All authors equally contributed reviewing and editing the paper.
Competing interests. The authors declare that they have no conflict of interest.
Acknowledgements. The original research  was funded by the European Union Horizon 2020 GEMex project (grant agreement no. 727550) and by the Mexican Energy Sustainability fund CONACYT-SENER, WP 4.5 of the project 2015-04-268074.
Review statement. This paper was edited by Joachim Gottsmann and reviewed by John Browning and one anonymous referee.