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 near- and 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.
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.
The Norini and Groppelli (2020) criticism of the Urbani et al. (2020) 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 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).
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, 2014; 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.
Trace of the A–A
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. Van Loon and Wiggers, 1975, 1976; Branney and Kokelaar, 1994). Norini and Groppelli (2020) fail to discuss any possible alternative origin for their small-scale 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 first-order, 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. (2013, 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. Karaoğlu 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.
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 scenario is further supported by the recent work of
Jentsch et al. (2020; see also Deliverable 4.3 of GEMex, 2019a), where soil
temperature anomalies (
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) interpretation. Again, it is unclear why these shapes are reported as counterproof. Asymmetric bulges are common characteristics in many volcanic regions worldwide, in resurgent calderas (e.g. Ischia, Pantelleria, Sierra Negra and Alcedo; Galetto et al., 2017, and references therein) or associated with shallow intrusions, such as Usu (Goto and Tomiya, 2019), Chaine de Puys (van Wyk de Vries et al., 2014; Petronis et al., 2019), Bezymianny (Gorshkov, 1959) and Mt St. 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).
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 Urbani et al. (2020) 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.
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 large-scale 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.
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). One of the most famous examples of felsic hypabyssal intrusions in intracaldera ignimbrite deposits is in Long Valley Caldera (California). At Long Valley, the well logs revealed ca. 300 m cumulative thick succession of aphanitic to phyric rhyolitic intrusions emplaced during the post-caldera stage, into the older, ca. 1200 m thick, intracaldera Bishop Tuff (McConnell et al., 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.
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) completely misinterpreted and misquoted a recent work by
Cavazos-Álvarez et al. (2020), which only deals with the reinterpretation of
andesitic layers within the Xaltipan intracaldera ignimbrite (see blue
ellipses in wells: H10 – Fig. 3a; H20 – Fig. 3b; H42 – Fig. 3e) and does not
question the interpretation of the rhyolite bodies proposed by Urbani et al. (2020) as small intrusions located above and below the Xaltipan ignimbrite.
With regard to these rhyolite bodies, Cavazos-Álvarez et al. (2020) not only
confirm their existence in wells H20 and H26 (red ellipses in Fig. 3b and
d of this reply) but also identify previously unrecognized (i) ca. 400 m
cumulative thick rhyolite layers (between ca. 500–1000 m below the surface)
in well H25 (Fig. 3c) and (ii) ca. 50 m thick rhyolite layers (between
850–900 m below the surface) in well H42 (Fig. 3e). The depths of these
rhyolitic layers are compatible with the estimated intrusion depth of
The statement by Norini and Groppelli (2020) on the absence of an in-depth
sharp increase in the temperature and geothermal gradient in the H4 well
(considered to remain constant at ca. 20
Lithostratigraphic columns of the wells
Norini and Groppelli (2020) question the presence of domes younger than 7.3 ka based on stratigraphic and radiometric data presented in Carrasco-Núñez et al. (2018). In agreement with the geological map of Carrasco-Núñez et al. (2017a), Fig. 5 shows a perspective view of the Los Potreros caldera floor across the Maxtaloya and Los Humeros faults. The images show the presence of lava domes and flows of variable composition both covered by and emplaced above the 7.3 ka Cuicuiltic Member. Older lavas include those associated with the “resurgent phase” (50.7–44.8 ka; “Qr1” in Carrasco-Núñez et al., 2018, and references therein). Younger lavas show the absence of the 7.3 ka Cuicuiltic Member cover and a morphology poorly or unaffected by evidence of faulting. 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.
Perspective view from a satellite image of the Los Potreros caldera floor (Image Landsat from © Google Earth Pro, 2020, Inegi-Maxar Technologies; courtesy of Google). The dashed blue lines outline the lava domes and flows (Qta2, Qta3, Qb1, Qta4, Qt1) mapped by Carrasco-Núñez et al. (2017a), whereas the dashed white lines outline the mapping of the Cuicuiltic Member (Qtc) from Urbani et al. (2020).
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.
Urbani et al. (2020) simply summarize field evidence stating “the recent (post-caldera collapse) uplift in the Los Potreros caldera moved progressively northwards, from the south and north-eastern 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 decompression–degassing, 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 (
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
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
All the data presented in this paper are available upon request.
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.
The authors declare that they have no conflict of interest.
The original research (Urbani et al., 2020) 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.
This research has been supported by Horizon 2020 (grant no. GEMex (727550)) and the Mexican Energy Sustainability fund CONACYT-SENER (WP 4.5 of project no. 2015-04-268074).
This paper was edited by Joachim Gottsmann and reviewed by John Browning and one anonymous referee.