Articles | Volume 13, issue 2
Data description paper
23 Feb 2021
Data description paper |  | 23 Feb 2021

Petrophysical and mechanical rock property database of the Los Humeros and Acoculco geothermal fields (Mexico)

Leandra M. Weydt, Ángel Andrés Ramírez-Guzmán, Antonio Pola, Baptiste Lepillier, Juliane Kummerow, Giuseppe Mandrone, Cesare Comina, Paromita Deb, Gianluca Norini, Eduardo Gonzalez-Partida, Denis Ramón Avellán, José Luis Macías, Kristian Bär, and Ingo Sass

Petrophysical and mechanical rock properties are key parameters for the characterization of the deep subsurface in different disciplines such as geothermal heat extraction, petroleum reservoir engineering or mining. They are commonly used for the interpretation of geophysical data and the parameterization of numerical models and thus are the basis for economic reservoir assessment. However, detailed information regarding petrophysical and mechanical rock properties for each relevant target horizon is often scarce, inconsistent or distributed over multiple publications. Therefore, subsurface models are often populated with generalized or assumed values resulting in high uncertainties. Furthermore, diagenetic, metamorphic and hydrothermal processes significantly affect the physiochemical and mechanical properties often leading to high geological variability. A sound understanding of the controlling factors is needed to identify statistical and causal relationships between the properties as a basis for a profound reservoir assessment and modeling.

Within the scope of the GEMex project (EU H2020, grant agreement no. 727550), which aims to develop new transferable exploration and exploitation approaches for enhanced and super-hot unconventional geothermal systems, a new workflow was applied to overcome the gap of knowledge of the reservoir properties. Two caldera complexes located in the northeastern Trans-Mexican Volcanic Belt – the Acoculco and Los Humeros caldera – were selected as demonstration sites.

The workflow starts with outcrop analog and reservoir core sample studies in order to define and characterize the properties of all key units from the basement to the cap rock as well as their mineralogy and geochemistry. This allows the identification of geological heterogeneities on different scales (outcrop analysis, representative rock samples, thin sections and chemical analysis) enabling a profound reservoir property prediction.

More than 300 rock samples were taken from representative outcrops inside the Los Humeros and Acoculco calderas and the surrounding areas and from exhumed “fossil systems” in Las Minas and Zacatlán. Additionally, 66 core samples from 16 wells of the Los Humeros geothermal field and 8 core samples from well EAC1 of the Acoculco geothermal field were collected. Samples were analyzed for particle and bulk density, porosity, permeability, thermal conductivity, thermal diffusivity, and heat capacity, as well as ultrasonic wave velocities, magnetic susceptibility and electric resistivity. Afterwards, destructive rock mechanical tests (point load tests, uniaxial and triaxial tests) were conducted to determine tensile strength, uniaxial compressive strength, Young's modulus, Poisson's ratio, the bulk modulus, the shear modulus, fracture toughness, cohesion and the friction angle. In addition, X-ray diffraction (XRD) and X-ray fluorescence (XRF) analyses were performed on 137 samples to provide information about the mineral assemblage, bulk geochemistry and the intensity of hydrothermal alteration.

An extensive rock property database was created (Weydt et al., 2020;, comprising 34 parameters determined on more than 2160 plugs. More than 31 000 data entries were compiled covering volcanic, sedimentary, metamorphic and igneous rocks from different ages (Jurassic to Holocene), thus facilitating a wide field of applications regarding resource assessment, modeling and statistical analyses.

1 Introduction

The knowledge of petrophysical and mechanical rock properties of the deep subsurface is essential for reservoir exploration and assessment of the reservoir potential for a variety of industrial applications such as petroleum reservoir engineering, geothermal heat extraction, mining or nuclear waste disposal. The data are most commonly used for interpreting geophysical data, creating conceptual geological models or populating numerical models (Lévy et al., 2018; Scott et al., 2019; Deb et al., 2019a, b; Árnason, 2020). Depending on the scale of investigation (e.g., local, regional or continental scale), highly accurate spatial predictions of relevant rock properties are required to increase the success and accuracy of reservoir operations and to reduce economic risks.

Rock formations are usually characterized by a heterogeneous internal structure, mineral composition, and pore and fracture distribution resulting in great variability in petrophysical and mechanical properties (Schön, 2015). Thereby, tectonic events, diagenetic or metamorphic processes, and hydrothermal alteration significantly affect the rock properties (Pola et al., 2012; Aretz et al., 2015; Weydt et al., 2018a; Mordensky et al., 2019; Durán et al., 2019, Heap et al., 2020), leading to a high geological heterogeneity often observed within hundreds-of-meter to sub-meter scales (e.g., Canet et al., 2010). Although most exploration methods or geological models are aligned to the reservoir scale, the controlling factors within the reservoir need to be understood and quantified at different scales to estimate the heterogeneity of each relevant formation and to assess the uncertainty in the input parameters for different modeling approaches. However, on the one hand, detailed information about rock properties for the relevant target formations is often not available, is inconsistent or is distributed over the literature. On the other hand, important metadata such as petrographic descriptions, details on sample locations and applied methods for data acquisition are missing (Bär et al., 2020). Without sufficient information, it is often not possible to evaluate and profit from existing laboratory data from specific locations or reservoir formations for future modeling approaches or studies related to similar geological settings. Consequently, most reservoir models are based on assumed or generalized data sets and local geological heterogeneities are often not considered (Mielke et al., 2015). While most studies focus on a single parameter (Clauser and Huenges, 1995) or a small set of samples, extensive data sets are required, which contain data of numerous different analyses performed on each sample in order to constrain statistical and causal relationships between the parameters (Linsel et al., 2020).

Addressing these challenges, the GEMex project (Horizon 2020, grant agreement no. 727550) embedded the petrophysical and mechanical rock characterization of the target formations in a comprehensive workflow providing the basis for different modeling approaches, geophysical surveys, ongoing and future volcanological studies. The GEMex project is a European–Mexican collaboration which aims to develop new transferable exploration and exploitation approaches for enhanced geothermal systems (EGSs) and super-hot unconventional geothermal systems (SHGSs). For this purpose, the Acoculco and Los Humeros geothermal fields have been selected as demonstration sites. Both fields are linked to caldera complexes located in the northeastern part of the Trans-Mexican Volcanic Belt (TMVB). Extensive geological, geochemical, geophysical and technical investigations were performed to improve the reservoir understanding and to facilitate future drilling operations.

Up until the beginning of the project in 2016, information on rock properties of the different geological units in the study area was scarce or not available. Previous studies focused on the investigation of reservoir core samples of both geothermal fields (Contreras et al., 1990; García-Gutiérrez and Contreras, 2007; Canet et al., 2015). However, the existing data were not sufficient for the definition and parameterization of model units within the reservoir due to the limited core material available (six pieces for Acoculco; Canet et al., 2015) or the lack of petrographic descriptions and chemical data for individual samples (Contreras et al., 1990).

Therefore, outcrop analog studies and reservoir core studies were performed in order to characterize all relevant key units from the basement to the cap rock (Weydt et al., 2018b; Bär and Weydt, 2019). Geological heterogeneities were investigated on different scales: (1) macroscale (outcrops), (2) mesoscale (rock samples) and (3) microscale (thin-section and chemical analysis). Analog studies of the geological units exposed in outcrops around the investigated geothermal fields offer a cost-effective opportunity to investigate and correlate facies, diagenetic and metamorphic processes, and lithofacies-related rock properties from outcrops down to the subsurface (Howell et al., 2014). The definition of thermo-facies units (Sass and Götz, 2012) and the quantification of uncertainties for each parameter enable a reliable prediction of rock properties in the subsurface.

A comprehensive database was developed including petrophysical, thermophysical, magnetic, electric, dynamic and static mechanical properties combined with chemical and mineralogical data. In total 34 parameters were determined on more than 2160 plugs retrieved from 306 outcrop samples from both caldera complexes and 66 reservoir core samples of the Los Humeros geothermal field as well as 8 core samples of the Acoculco geothermal field covering volcanic, sedimentary, metamorphic and igneous rocks from Jurassic to Holocene age. Here, we present the workflow and current status of the GEMex rock property database (Weydt et al., 2020; These data not only provide the basis for ongoing research in the study area but also facilitate a wide field of applications in different disciplines, for example, a first assessment of the subsurface properties at early exploration stages (Bär et al., 2020), different modeling approaches, geostatistical and stochastic analyses, or the validation of different measurement methods.

GEMex project framework and sampling

The geothermal system in Los Humeros is steam dominated and has been under production since 1990, operated by the Comisión Federal de Electricidad (CFE). With a production of 94.8 MWe in 2018 it is the third-largest geothermal field in Mexico (Romo-Jones et al., 2019) with 65 wells drilled so far, of which 28 are productive and 5 are used as injection wells. With temperatures above 380 C encountered below 2 km depth in the northern part of the field, the Los Humeros caldera complex was characterized as a suitable target for the development of a SHGS within GEMex. In Acoculco two exploration wells have been drilled so far, which have encountered temperatures of approximately 300 C at a depth of about 2 km (Canet et al., 2015). Although a well-developed fracture network exists within the area, both wells were dry (López-Hernández et al., 2009). Thus, the GEMex project aims to develop a deep EGS in Acoculco in order to connect the existing wells to proximal fluid-bearing fracture zones.

The project comprises a multidisciplinary approach based on three milestones which are (1) resource assessment, (2) reservoir characterization and (3) concepts for site development (Jolie et al., 2018). The first milestone focused on a comprehensive understanding of structurally controlled permeability and the fluid flow in the reservoir including extensive fieldwork regarding stratigraphy and structural geology, fracture distribution, hydrological and geochemical studies of natural springs, comprehensive soil–gas studies (e.g., CO2 flux; Jentsch et al., 2020), and airborne thermal imaging. The second milestone includes several geophysical surveys (e.g., passive and active seismic, gravity and magnetotelluric surveys) to characterize active faults and to identify deep structures. In addition, extensive sampling campaigns were conducted for petrophysical, rock mechanical, chemical and mineralogical investigations of the key lithologies in the study area. Resulting data and models of all work groups are being combined in integrated reservoir models at a local, regional and superregional scale. The third milestone includes the investigation of transferable concepts for developing EGSs and the utilization of SGHSs; the identification of suitable materials and well designs, which can resist high temperatures and corrosive fluids in the reservoir; and the determination of possible drill pathways along with a comprehensive risk assessment and management.

Figure 1Schematic workflow of the GEMex project using the example of the El Dorado mine in Las Minas (d) with view on the footwall of the present fault (photo from Maximilian Bech). The quarry exposes exoskarn in many variations. Outcrop analysis included detailed investigation of kinematic indicators, mineralogy (a) and the main fracture pattern (e) to create numerical fluid flow models (f) as presented in Lepillier et al. (2019). Rock samples taken for lab investigation (b) and geochemical and thin-section analysis (c) (photo from Caterina Bianco). Cylindrical plugs drilled from the outcrop samples (g), which were distributed between the partners in order to determine rock properties, dating or high-T/P experiments (the experiments marked in blue are not included in this study).


The work presented in this study is part of milestone 2 (reservoir characterization) and focuses on the mineralogical, petrophysical and mechanical rock characterization of both geothermal systems. Several joint field campaigns with Mexican and European partners were conducted in order to cover and sample all relevant geological key units from the basement to the cap rock. In this context, work groups with different areas of expertise worked together in a joint approach (Fig. 1). Thus, structural geologists worked together with volcanologists, petrologists and petrophysicists on the same outcrops to, e.g., combine results of fracture pattern characterization and rock property analysis obtained from the same outcrops in a numerical fluid flow model (Lepillier et al., 2019). Likewise, samples for detailed mineralogical investigations were collected together with samples for petrophysical experiments. Over 300 representative samples were collected from more than 140 outcrops inside the caldera complexes and in the surrounding area (Fig. 2). In addition to outcrop analysis in the Acoculco and Los Humeros areas, particular attention was paid to the exhumed systems Zacatlán (east of Acoculco) and Las Minas (east of Los Humeros), where all units from the cap rock to the basement are exposed. These so-called “fossil systems” serve as proxies for the active geothermal fields and help in the understanding of the fluid flow and mineralization processes in the “active” geothermal reservoirs under discussion. Whenever possible, each geological unit was sampled several times at different outcrop locations to cover the unit's heterogeneity, and only samples with an overall fresh appearance unaffected by weathering were considered. Hydrothermal alteration of different intensities was observed in some outcrops in close proximity to fault zones and dikes. In these cases, hydrothermally altered samples were deliberately collected to analyze the effect of these processes on the rock properties. Besides analyzing outcrops and outcrop samples, the CFE granted extensive sampling of wellbore core material of both geothermal fields at the CFE camp in Los Humeros. In total 66 samples drilled from 37 core sections covering 16 wells drilled in Los Humeros and 8 core samples drilled from 6 core sections from well EAC1 of the Acoculco geothermal field were obtained. All samples were directly drilled within the field or sent as boulders to Europe or the Mexican institutes and subsequently distributed between the partners. This approach ensures that further work on the project, such as long-term flow experiments (Kummerow et al., 2020), high-T/P experiments, hydraulic fracture experiments (Deb et al., 2019c), detailed mineralogical analyses (thin-section and scattered electron microscope; Lacinska et al., 2020), isotope analyses or dating (Kozdrój et al., 2019), can be directly correlated with the results presented in this study. Furthermore, some parameters of the same sample set were analyzed by multiple institutes to compare and validate different analytical approaches.

Figure 2Geological map of the Acoculco and Los Humeros region including the sampling points of the outcrop samples (SGM, 2002a, b). The faults were recently mapped and characterized by Liotta et al. (2019) and Norini et al. (2019).

2 Geological setting

The Acoculco and Los Humeros caldera complexes are located in the northeastern part of the Trans-Mexican Volcanic Belt (TMVB), 125 and 180 km east of Mexico City, respectively. The E–W trending TMVB is a ∼1000 km long calc-alkaline arc which is directly linked to the subduction of the Rivera and Cocos plates beneath the North American Plate along the Middle America Trench (Ferrari et al., 2012; Macías et al., 2012; Avellán et al., 2018). The volcanic complexes are located over a ∼50 km thick continental crust (Pérez-Campos et al., 2008) and are situated ∼100 km north of the Popocatépetl and Pico de Orizaba volcanoes, which define the most active front of the TMVB in central eastern Mexico (Ferrari et al., 2012; Macías et al., 2012; Avellán et al., 2020).

Both volcanic complexes are emplaced on intensively folded Mesozoic sedimentary rocks (Mexican fold-and-thrust belt; Fitz-Díaz et al., 2017) belonging to the Sierra Madre Oriental comprising Jurassic sandstones, shales, hydrocarbon-rich limestones and dolomites overlain by Cretaceous limestones and shales (López-Hernández et al., 2009; Fitz-Díaz et al., 2017). The regional tectonic setting is characterized by Late Cretaceous–Eocene NW–SE striking thrusts and folds and subordinate NE-striking normal faults that are associated with an Eocene–Pliocene extensional deformation phase (Norini et al., 2019). Oligocene to Miocene granitic and syenitic plutons as well as andesitic and basaltic dikes intruded into the sedimentary sequences, leading to local metamorphism of marble, hornfels and skarn (Ferriz and Mahood, 1984; Fuentes-Guzmán et al., 2020). The sedimentary basement is exposed east and southeast of the Acoculco caldera close to Chignahuapan and Zacatlán as well as in the surroundings of the Los Humeros caldera. Furthermore, it was also cut at different depth levels in drill cores in both geothermal fields (López-Hernández, et al., 2009; Carrasco-Núñez et al., 2017a). The granitic plutons are spread over the study area, and new aeromagnetic data of the Acoculco caldera constrain the occurrence of at least four intrusive bodies hosted in the Cretaceous limestones at >1 km depth. Those were interpreted as a series of horizontal mafic intrusions providing the energy to maintain the geothermal field (Avellán et al., 2020).

The Acoculco caldera complex has an 18 km×16 km semi-circular shape (Avellán et al., 2018) and predominantly comprises Pliocene to Pleistocene basaltic to rhyolitic lavas, domes, cinder cones and ignimbrites. The caldera complex sits on an intersecting NE–SW and NW–SE fault system creating an orthogonal arrangement of grabens, half grabens and horsts (García-Palomo et al., 2002, 2018). Thereby the regional tectonic regime strongly affected the local tectonic behavior and structural deformation of the caldera (Sosa-Ceballos et al., 2018). The Acoculco caldera is located on the NE–SW Rosario-Acoculco horst and was built on top of Cretaceous limestones and the Zacatlán basaltic plateau (so far undated) as well as Miocene and Pliocene lavas and domes related to the regional volcanism of the TMVB (Avellán et al., 2018, 2020). Thereby the pre-caldera lavas and scoria cones exposed north and northeast of the Acoculco caldera complex were related to the Apan-Tezontepec Volcanic Field (Miocene and Pliocene), whereas Miocene andesitic and dacitic lavas are exposed west of the Acoculco caldera complex. Magmatic activity of the Acoculco caldera can be divided into five different eruptive phases, including recent deposits and hydrothermal altered areas inside the caldera (Avellán et al., 2018). It began with the emplacement of the Acoculco ignimbrite (∼2.7 Ma; 40Ar/39Ar), followed by several early (∼2.6–2.1 Ma) and late (∼2.0<0.016 Ma) post-caldera volcanic events producing basaltic to trachyandesitic and rhyolitic lava flows restricted within the caldera and rhyolitic lava domes, scoria cones and two ignimbrites that predominantly migrated to the caldera rim and periphery, respectively. The extra-caldera volcanism (2.4–0.19 Ma) comprises several basaltic trachyandesitic to basaltic andesitic lavas and scoria cones, related to the volcanism of the Apan-Tezontepec Volcanic Field. Products of the extra-caldera volcanism are interbedded with the lavas of the Acoculco caldera complex. It has to be emphasized that recent studies (Avellán et al., 2018, 2020) are not in line with previous volcanological studies performed by López-Hernández et al. (2009). In the study conducted by López-Hernández et al. (2009), the authors concluded that the Acoculco caldera (1.7–0.24 Ma) is nested within the older and larger Tulancingo caldera (∼3.0–2.7 Ma) forming the so-called Tulancingo–Acoculco caldera complex and that a third volcanic episode (1.8–0.2 Ma) occurred, which was related to monogenetic volcanism without a caldera collapse.

The younger Los Humeros caldera is the largest active caldera of the TMVB with a 21 km×15 km irregular shape and comprises predominantly Pleistocene to Holocene basaltic andesitic to rhyolitic volcanic rocks (Carrasco-Núñez et al., 2018; Norini et al., 2019). The oldest volcanic activity in this area is represented by a thick sequence of Miocene andesites, dacites and basaltic lava flows of the Cuyoaco and Alseseca andesite unit (∼10.5 Ma; Yáñez and García, 1982) and Pliocene to Pleistocene basaltic to andesitic lavas belonging to the Teziutlán andesite unit (dated between 1.44±0.312.65±0.43 Ma, 40Ar/39Ar; Carrasco-Núñez et al., 2017a). Miocene lavas have a cumulative thickness of up to 900 m and can be related to the Cerro Grande Volcanic Complex dated between 8.9–11 Ma (Carrasco-Núñez et al., 1997; Gómez-Tuena and Carrasco-Núñez, 2000), and Teziutlán andesite lavas have a reported thickness of up to 1500 m (López-Hernández, 1995). Both units are classified as “andesitic and basaltic volcanic basement” and form the currently exploited reservoir in the subsurface of the Los Humeros geothermal field (Carrasco-Núñez et al., 2018). The beginning of the magmatic activity of the Los Humeros volcanic complex is represented by rhyolitic lavas and abundant rhyolitic domes, mainly located at the western side of the volcanic complex (270±17 and 693±1.9 ka; Carrasco-Núñez et al., 2018). However, the caldera collapse itself is associated with the emplacement of the high-silica rhyolite Xaltipán ignimbrite at ∼160 ka with an estimated volume of 291 km3 and a thickness of up to 880 m (Carrasco-Núñez et al., 2018; Cavazos and Carrasco-Núñez, 2020). After the emplacement of the Xaltipán ignimbrite, which caused the characteristic trapdoor structure of the caldera, further explosive events led to the deposition of thick rhyodacitic Plinian deposits called Faby Tuff (Norini et al., 2015; Carrasco-Núñez et al., 2017a). Afterwards, a second caldera-forming eruption occurred at ∼69 ka and is related to the Zaragoza ignimbrite emplacement forming the Los Potreros caldera within the Los Humeros caldera. The post-caldera stage is represented by rhyolitic and dacitic domes within the center of the caldera (44.8±1.7 ka) and basaltic to trachyandesitic lava flows (8.9±0.03 ka), volcaniclastic breccias and fallout deposits (7.3±0.1 ka) with a highly variable lateral and vertical distribution (Carrasco-Núñez et al., 2017a, 2018).

3 Workflow

After the samples were distributed between the partners, cylindrical cores with diameters ranging from 25 to 65 mm were drilled and subsequently cut according to standards (ASTM D4543-19, 2019) for the required sample length, and the irregular and rough core ends were cut to be parallel. The laboratory tests were divided into three stages: (1) general petrophysical characterization including all non-destructive measurements, (2) mechanical rock characterization, and (3) chemical and mineralogical characterization. Non-destructive tests included particle density, bulk density, porosity, intrinsic matrix permeability, thermal conductivity in dry and saturated conditions, thermal diffusivity in dry and saturated conditions, P-wave velocity and S-wave velocity in dry and saturated conditions, specific heat capacity, magnetic susceptibility, and electric resistivity in dry and saturated conditions. Afterwards the destructive rock mechanical tests such as the Brazilian disc test, the chevron bend test, the point load test, and uniaxial and triaxial tests were performed to determine uniaxial compressive strength, Young's modulus, the Poisson ratio, tensile strength, fracture toughness, the friction angle and cohesion. Samples that were identified as suitable for destructive tests such as uniaxial or triaxial tests were ground plane-parallel prior to analysis. Quantitative and qualitative chemical analyses like X-ray fluorescence (XRF) and X-ray diffraction (XRD) as well as thin-section analyses were performed for the petrological and geochemical characterization. Figure 3 shows the schematic laboratory workflow of TU Darmstadt.

Figure 3Schematic workflow representing the measurement procedure at TU Darmstadt. The properties displayed in orange were determined on sample material and used to calculate those shown in red. Parameters marked with * were analyzed in dry and saturated conditions.


4 Structure of the database and sample classification

The database is publicly available under (Weydt et al., 2020) and contains petrophysical and rock mechanical properties as well as chemical data obtained by laboratory experiments within the scope of the GEMex project. This database is provided in a flat file Excel format and in CSV format to keep the handling as simple as possible. Its internal structure is based on the PetroPhysical Property Database – P3 – previously developed during the IMAGE project (Bär et al., 2020) with some project-specific modifications. The P3 database's internal design comprises multiple tables for petrography, stratigraphy, quality controls, chemical analyses and petrophysical properties and follows the concept of relational database management (Codd, 1970). As the database presented in this study is restricted to one study area, the P3 structure was simplified, and the sample's information has been compiled in two data sheets so far. The main objective was to provide the data in a user-friendly and well-structured form, allowing easy filtering and a transfer of data into other database formats like SQL (structural query language) to easily visualize it or to implement it for modeling approaches.

The first and main data sheet comprises all analyzed petrophysical parameters and sample information (metadata) compiled during this project. Each analyzed plug was provided with a sample ID, which acts as the primary key for all records. Sample information provided in the database is explained in the following sub-sections.

The second data sheet includes all chemical data, retrieved from composite sample material, and does not directly correspond to measurements on single plugs. The data are provided separately to increase handling and readability. Here, the sample name represents the primary key which links the data to the petrophysical measurements provided in the first table.

4.1 Metadata

The metadata include all additional sample information from sample ID to sample dimensions and can be used for rapid filtering and the precise categorizing of parameters.

Each analyzed plug or sample received a unique sample ID, which is derived from the sample name given in the field, the geothermal reservoir (LH or AC), the field trip (e.g., M17 for May 2017) and an abbreviation for the rock type (e.g., GD for granodiorite). This classification was developed within the project due to the high number of samples collected during different field trips. Furthermore, the sample ID provides information about the sample preparation. In hierarchical order the sample name, core name and plug name are provided. For each drilled core the sample name was complemented with C1 (i.e., core number 1), C2, C3 and so on. Whenever the core did not meet the requirements for destructive measurements (length-to-diameter ratio of 2:1 or too fragile), the core was cut into plugs. The core name was then complemented with capital letters A, B, C, etc. representing the way the core was cut (Fig. 4). The implementation of this hierarchical order allows for quick access of the parameters per plug, per core or per sample. Whenever a core was not cut into several plugs, the core and plug name are identical to avoid gaps in the database. For practical reasons the term “plug” was used for all cylindrical samples after sample preparation (cutting and grinding) ready to be analyzed. For the reservoir core samples, the existing core names were adopted. The ID begins with the well name (e.g., H23), followed by the core number (e.g., number 2), the core section (e.g., 14, or x for undefined) and the number of the drilled subcore (C1 or C2).

Figure 4Overview of the different preparation steps and sample labeling. Cores (b, various diameters; d, 40 mm in diameter) drilled from outcrop samples (a) and reservoir core samples (c) and subsequently cut into plugs (e) to meet the individual requirements of the measurement devices. The plugs were labeled with capital letters.


The samples were classified regarding their rock type and stratigraphic unit based on the recently published geological maps and volcanological studies conducted in Acoculco and Los Humeros (Avellán et al., 2018, 2020; Carrasco-Núñez et al., 2017a, b, 2018). Rock types were predominantly determined using macroscopic analyses complemented by thin-section analyses (whenever available). Additionally, bulk chemical analyses (XRF) were used to better characterize the volcanic rocks using the TAS classification (Le Maitre and Streckeisen, 2003). However, this classification is only applicable for unaltered sample material. The classification of the stratigraphic unit is based on the international chronostratigraphic chart of the IUGS (Cohen et al., 2013) according to international standardization. Whenever possible the local stratigraphic unit is given. The volcanological studies are still ongoing, and the ages of some units or areas are not yet well constrained.

Coordinates of the sampling locations are provided as latitude and longitude in decimal degrees (WGS84) and x and y coordinates (UTM WGS84). For the reservoir core samples, the coordinates of the well heads are included. All this information is given in meters above sea level (m a.s.l.) and represents the surface evaluation of the outcrops or the evaluation at reservoir depth for the reservoir core samples. The latter was provided in measured depth (MD) by the CFE, and the core sample material was obtained from vertically drilled wellbores.

Furthermore, the outcrop names and field trips are documented as project internal information and enable the placing of this work in relation to other work conducted within the study area. Samples from six field trips are provided in the database as shown in Table 1.

Table 1Overview of the field campaigns and related work.

Download Print Version | Download XLSX

The “location” was inserted in addition to the outcrop name and sample coordinates to classify the samples according to their sampling area, distinguishing between Acoculco, Los Humeros, and the exhumed systems Las Minas and Zacatlán–San Miguel Tenango (SMT). The column “institution” refers to the institution and authors that generated the data and indirectly links this to the applied methods described in Sect. 5.

Based on the rock type and stratigraphic classification, the samples were related to the model units of the regional and local geological models created within the GEMex project (Calcagno et al., 2018, 2020). The regional and local model units were defined to consider the most representative geological formations in the study area, the scale of the model and the objective of the project (Calcagno et al., 2018). For Los Humeros four regional and nine local model units were defined (Fig. 5). The classification is mostly based on recent work of Carrasco-Núñez et al. (2017a, b, 2018) and Norini et al. (2015, 2019) and information about formation depth, thickness and distribution provided by the CFE stratigraphic drilling profiles. Samples collected from basaltic and andesitic dikes as well as from intrusive bodies in Los Humeros and Las Minas were related to the basement (G4 and U9). The classification of the local units of the reservoir core samples represents the classification used for the latest update of the local model of Los Humeros (Calcagno et al., 2020).

Figure 5Regional and local model units of the 3D geological model of Los Humeros (slightly modified from Calcagno et al., 2018, 2020).

For the regional model of Acoculco, five units were defined (Fig. 6). All volcanic deposits were merged into one unit called AC5-Volcanites, whereas the basement rocks were split into four separate units: AC4-Limestones, AC3-Skarns, AC2-Granite and AC1-Basement. The description and stratigraphic classification is based on López-Hernández et al. (2009), Lorenzo-Púlido et al. (2010), Sosa-Ceballos et al. (2018) and Avellán et al. (2018).

Figure 6Regional model units of the 3D geological model of Acoculco (slightly modified from Calcagno et al., 2018).

As the last entities belonging to the metadata, sample descriptions and dimensions for each plug are provided. The sample description includes a brief macroscopic description and gives information about the occurrence of fractures, joints and fissures or other remarks (e.g., chert nodules or stylolites). Furthermore, the information is given as to whether thin sections were prepared or not. The section “sample dimensions” includes the length, diameter (exact and drilled diameter), weight (dry and saturated) and shape of the plug. Plug shapes were inserted for quality control and were distinguished between “ideal cylindrical plug”, “cylindrical plug with a broken edge”, “irregular shape” and “cuboid”. This information needs to be considered when the bulk density or volume is calculated by using the sample's dimensions. The exact sample dimensions provide the opportunity to analyze scale-dependent effects (Enge et al., 2007). Therefore, plugs with varying diameters and lengths were drilled and analyzed. Thus, small-scale samples (25 mm in diameter) for which the bulk volume reaches the minimal representative elementary volume (REV; e.g., Ringrose and Bentley, 2015) are included.

4.2 Rock properties

Rock properties provided are grouped into (1) classical petrophysical parameters such as density, porosity and permeability; (2) ultrasonic wave velocities; (3) thermal properties; (4) magnetic susceptibility; (5) electric resistivity; and (6) rock mechanical parameters. The results are provided as mean values with standard deviation (whenever possible) for each plug. For thermal conductivity and thermal diffusivity the maximum and minimum values were added. In total 34 different parameters were obtained following the recommendations of international standardization institutions and committees (e.g., ISRM, ASTM or DIN). Columns for specific remarks were included to provide further details whenever needed. Detailed information on methods and procedures is given in Sect. 5.

4.3 Chemical analyses

The results of chemical analyses (XRF and XRD) are provided in the second data sheet of the database. These data are retrieved from composite sample material, and a total of 131 samples (reservoir core samples and outcrop samples) were analyzed. The sample name acts as the primary key and allows for linking of chemical data with petrophysical data. Results of the XRF analyses are presented in weight percent for the major elements and in parts per million for the trace elements. For both analyses (XRF and XRD) the responsible institution is added to relate the data to the applied method.

5 Material and methods

The following sections briefly describe the applied methods conducted by the different partners. A more extensive description for the non-destructive measurements and the field trips can be found in project reports on the GEMex web page (Bär and Weydt, 2019;, last access: 21 October 2020). Sample material from TU Delft (field trip January 2017) and TU Darmstadt (field trip May 2017) were distributed to GFZ, RWTH Aachen and UniTO for non-destructive petrophysical measurements.

5.1 Sample preparation

Drill cores with diameters ranging from 25 to 65 mm were drilled from the outcrop samples and cut into plugs as described above. More than 2100 plugs and cores with an axial length ranging from ∼30 to 128 mm were prepared according to international standard ASTM D4543-19 (2019). The short plugs (diameter 25 to 40 mm, length 25 to ∼30 mm) were predominantly used for the non-destructive petrophysical measurements like bulk density, porosity and permeability due to the specific sample size requirements of the measurement devices. Remaining plugs were prepared to meet the requirements for the different destructive rock mechanical tests, which were conducted after the petrophysical characterization. For most of the rock mechanical tests a length-to-diameter ratio of 2:1 (uniaxial and triaxial tests) or 1:2 (Brazilian test) is required. Furthermore, the plane surfaces of the plugs had to be plane-parallel with a maximum angular misalignment of 0.05.

To ensure reproducibility of the results, the plugs were measured in oven-dry conditions (105 C for more than 24 h or 64 C for more than 48 h) and cooled down to room temperature in a desiccator (20 C). Microcracking or significant mass losses caused by mineralogical changes or the collapse of clay minerals during heating in the oven were not observed since a majority of the outcrop samples contain no clays and samples affected by hydrothermal or metamorphic processes contain mineral assemblages developed at higher temperatures.

In order to perform measurements in saturated conditions, the samples were evacuated in a desiccator and subsequently saturated with (de-ionized) water (TU Darmstadt and GFZ) or the samples were fully immersed in water for up to 4 weeks (RWTH Aachen and UniTO).

5.2 Non-destructive tests

At TU Darmstadt, density measurements were performed in a multi-step procedure using an AccuPyc helium pycnometer (ASTM D5550-14, 2014) and a GeoPyc powder pycnometer (Micromeritics, 1997, 1998, 2014), analyzing particle and bulk volume five times for each plug, respectively. Bulk density was then automatically calculated by dividing the dry weight of the plug by its measured volume. Afterwards porosities were calculated from the resulting differences in volume and represent the gas-effective porosity, also known as connected porosity. The accuracy of the method is 1.1 % (Micromeritics, 1998). Porosity measurements at TU Delft and UNAM were also performed using a helium gas pycnometer (Ultrapyc 1000 Version 2.12 and Ultrapyc 1200e gas pycnometers, respectively, both Quantachrome Corporation, USA) to determine the grain density (ASTM D5550-14, 2014), while bulk density was determined using caliper techniques according to ASTM D7263-16 (2016). Every plug was measured up to 20 times.

At GFZ and RWTH Aachen, particle density, bulk density and porosity were determined using the triple weighing method (ISRM, 1981). This method is based on the Archimedes principle, which uses the masses of the dry and fluid-saturated samples as well as that of the sample totally immersed in the fluid to calculate the pore volume and the porosity. The mass was determined with an accuracy of ±0.2 g. Usually, the accuracy is 1.5 % or better, but this especially depends on the surface condition for low-porosity samples. Thus, the measurements were performed up to three times per plug. A similar approach was used at UniTO by applying caliper techniques and the dry and saturated mass of each sample for the calculation of density and porosity (ISRM, 1979). Variations in particle and bulk density between the different methods applied on the same samples in this study range between 0.3 %–3 % (coefficient of variation) for limestones with porosities smaller than 3 % and 0.5 %–3.5 % for pyroclastic rocks with porosities between 11 % and 15 %, verifying the different methods and sample saturation procedures as sufficient to obtain data with the accuracy needed.

Matrix permeability was determined on cylindrical plugs (diameter and length ranging from 25 to 40 mm and ∼20 to 80 mm, respectively) with column gas permeameters constructed according to ASTM D4525-13e2 (2013) and ASTM D6539-13 (2013) standards at TU Darmstadt, GFZ and UNAM. The plugs were analyzed in a confined cell at constant differential pressure under a steady-state gas flow using at least five pore fluid pressure levels (Tanikawa and Shimamoto, 2008). Corresponding gas flow rates were measured with different flowmeters that allow for the detection of flow rates in the range between 10 and 10 000 cm min−1. This applied method is based on Darcy's law enhanced by factors for the compressibility and viscosity of gases in order to calculate the gas permeability (Scheidegger, 1974; Jaritz, 1999). The water equivalent permeability was derived from the gas permeability after the Klinkenberg correction (Klinkenberg, 1941). At TU Darmstadt the samples were analyzed with dried compressed air at five pressure levels ranging from 1 to 3 bar and 1 MPa confining pressure (Hornung and Aigner, 2004; Filomena et al., 2014). At GFZ a confining pressure of 8.5 MPa and five pressure levels ranging between 7.5 and 35 bar were applied (operated with argon), while at UNAM the permeability was determined using a confining pressure of 2.8 MPa and also five pressure levels up to 1 MPa (operated with nitrogen). Measurement accuracy of the TU Darmstadt permeameter varies from 5 % for highly permeable rocks (K>10-14 m2) to 400 % for impermeable rocks (K<10-16 m2) (Bär, 2012). The recorded flow rates were tested for turbulent fluid flow according to Kushnir et al. (2018) prior to the Klinkenberg correction to ensure a laminar fluid flow. A correction after Forchheimer (1901) was not required, since the corrected values were within the error range of the measurement device.

At TU Darmstadt, thermal conductivity and thermal diffusivity were measured simultaneously on oven-dried and saturated plugs using a thermal conductivity scanner (Lippmann and Rauen, Germany) after Popov et al. (1999, 2016). The device consists of a sample platform and an optical scanning system that moves along the sample surfaces, including a heat emitter and three infrared sensors facilitating a continuous profile. Samples are heated up by a defined heat flow, and the subsequent cooling rate is measured by the temperature sensors. Bulk thermal conductivity and thermal diffusivity were then calculated after Bär (2012) by using two reference standards. Both parameters were measured four to six times on each plug for saturated and dry conditions, respectively (two to three times on every planar surface including slight turning after every measurement to account for sample anisotropy). At RWTH Aachen, the same optical scanning method was used to determine thermal conductivity along the core axis of large cylindrical cores with diameters of 60 and 64 mm. To ensure uniform reflection conditions, the samples were painted with black acrylic paint on the planar surface (TU Darmstadt) and along the core axis (RWTH Aachen). According to Lippman and Rauen (2009), the measurement accuracy for thermal conductivity and thermal diffusivity is 3 % and 5 %, respectively.

Specific heat capacity was determined at TU Darmstadt using a heat-flux differential scanning calorimeter (C80; Setaram Instrumentation, 2009); crushed sample material was heated at a steady rate from 20 up to 200 C within a period of 24 h. Specific heat capacities were derived from the resulting temperature curves through heat flow differences. The accuracy is 1 % (Setaram Instrumentation, 2009). Volumetric heat capacity was calculated by multiplying the specific heat capacity with the associated bulk density of each sample. For direct comparison, specific heat capacity was calculated for each plug by dividing thermal conductivity by the product of bulk density and thermal diffusivity (Buntebarth, 1980).

Ultrasonic wave velocity was measured along the sample axis with pulse generators (TU Darmstadt – UKS-D including a USG-40 pulse generator and a digital PicoScope oscilloscope from Geotron-Elektronik, 2011; UniTO – Pundit Lab, Proceq, Switzerland according to ASTM D2845-08, 2008; GFZ – Panametrics HV pulser–receiver model 5058PR in combination with digital oscilloscope model DSO6012A from Agilent Technologies, USA) comprising point-source transmitter–receiver transducers. Thereby, the transducers were pressed against the parallel surfaces of the samples using a contact pressure of about 1 bar. Polarized pulses at high voltage in a frequency range from 20 kHz to 1 MHz for the USG-40 and Panametrics as well as from 54 to 250 kHz for the Pundit Lab were generated. The transmitted signals were recorded using digital oscilloscopes, and the arrival times of the P and S waves were picked manually and corrected for the dead time, which arises from the recording device (transducer, function generator, oscilloscope).

Bulk density and P- and S-wave velocities were used to determine dynamic elastic mechanical parameters, such as the dynamic shear modulus, Gdyn; dynamic Young's modulus, Edyn; and dynamic Poisson ratio, μdyn, after Zoback (2011):


where ρ is the bulk density [kg m−3], νp is the compressional wave velocity [m] and νs is the shear wave velocity [m].

Additional field measurements of P-wave velocities were performed by UniTO on irregularly shaped outcrop samples by using the same Pundit Lab Proceq device along different directions on the sample surfaces in order to identify anisotropy and the effect of fractures. Measurements were conducted following ASTM D2845-08 (2008) standard requirements. At TU Darmstadt both velocities were measured four to six times on each plug in both saturated and dry conditions. For analyzing the samples in saturated conditions, the samples were stored in degassed and de-ionized water to avoid desaturation. After preparing the device and measurement setup, the samples were immediately installed between the transducers and the transmitted signals were recorded until the sample started to desaturate. The data provided by GFZ represent average values from at least 4 to 10 individual measurements per plug (dry and saturated conditions), and at UniTO each sample was analyzed up to 20 times in order to depict the matrix heterogeneity of the larger cores and outcrop samples. The error in P-wave velocities is 3 % on average, whereas for S-wave velocities the average error is 8 % or higher, due to the higher attenuation and distortion of the S-wave signals.

Electric resistivity measurements were carried out on selected cylindric plugs at GFZ and UniTO and on outcrop samples in the field. At UniTO electric resistivity measurements were performed with a purpose-built square quadrupole (Syscal-Pro from Iris Instruments, France) after Clement et al. (2011). This consists of a rubber jacket with four steel electrodes (2 mm diameter and 40 mm length), arranged at the edges of two perpendicular diameters of the core sample at half of its longitudinal length. Electrical resistivity measurements were performed with a current injection between two subsequent electrodes and detection of the resulting electric potential between the remaining pair of electrodes. Current and potential electrodes were progressively reversed and rotated around the sample for a total of eight different potential measurements. The sequence was repeated three times, and each sample was tested in both dry and saturated (wet) conditions. Saturated conditions were reached by immersing the sample in a saline solution (with electrical conductivity equal to 1000 µS cm−1) for 24 h. A detailed description of the measurement procedure is also included in Vagnon et al. (2019).

Electric resistivity measurements at GFZ were executed with a four-electrode layout as well using an impedance spectrometer (Zahner Zennium electrochemical work station; Zahner Scientific Instruments, 2008), which supplied an AC voltage with an amplitude of 200 mV via disc-shaped current electrodes to the plane-parallel faces of the sample cylinders. The sample resistance was determined via detection of the impedance and the phase angle at distinct frequencies. Subsequently, the bulk resistivity was calculated from the sample resistance at 1 kHz, the cross-sectional area of the sample, and the distance between the potential electrodes that were pinned to the cylinder surface of the sample plugs. The measurements were performed on dry and on saturated samples. Oven-dry sample cores were saturated under vacuum with a NaCl solution with electrical conductivity equal to 1080 µS cm−1 and equilibrated for about 24 h. Prior to the measurements the samples were jacketed with a tight-fitting silicon sleeve to reduce the risk of desaturation. The accuracy of measurements in dry conditions is better than 3.5 %. In contrast, in saturated conditions for porous samples, the error increases to a maximum of 16 % if fluid evaporates or leaks from the pore space during the measurement interval.

The formation factor, F, of the samples was determined after Flovenz et al. (2005) from linear plots of bulk conductivities versus fluid conductivities at different brine concentrations, where F is the reciprocal of the linear fitting lines of the data points measured at fluid salinities varying between 0.56–10.42 S m−1.

Magnetic susceptibility was analyzed using the magnetic susceptibility meter SM30 (ZH Instruments, 2008), which consists of an oscillator with a pickup coil. An interpolating mode was applied including two air reference measurements and one measurement directly on the sample surface. The frequency change in the oscillator is proportional to the magnetic susceptibility of the rock sample. To ensure optimal contact of the sensor on the sample surface and to reduce the impact of air while measuring, only the plane surfaces of the plugs were analyzed.

Furthermore, a multi-sensor core logger (MSCL) from Geotek (2000) was used for measurements of gamma density, P-wave velocity, magnetic susceptibility and electrical resistivity at RWTH Aachen on whole cores with a diameter of 60–64 mm. Matrix density was calculated based on attenuation of gamma rays emitted from cesium-137, while porosity was calculated from the density measurements. P-wave velocity was measured using P-wave transducers (receiver and transmitter) mounted on opposite faces on the center sensor stand. A short pulse is produced at the transmitter, which propagates perpendicularly to the axis of the core and is detected by the receiver on the other side. The outer diameter of the core is measured with an accuracy of 0.1 mm. An absolute accuracy of ±3 m s−1 is achievable while computing the P-wave velocity. Magnetic susceptibility was determined using a Bartington loop sensor with a 5 % calibration accuracy. The sensor includes an oscillator circuit that generates a low-intensity alternating magnetic field at 0.565 kHz.

5.3 Destructive tests

Simple (non-cyclic) and cyclic uniaxial tests were performed to determine the rock's unconfined compressive strength and elastic rock mechanical properties, such as the static Young's modulus, Poisson's ratio, G modulus (also known as shear modulus) and bulk modulus. For the determination of the unconfined compressive strength (UCS) at TU Darmstadt, cylindrical plugs with a diameter of 40 mm and a length of 80 mm were introduced into a hydraulic uniaxial press (FORM+TEST Prüfsysteme, Germany) with a capacity of 1000 kN and a maximum loading rate of 0.5 kN s−1 until sample failure. The stress at this particular point represents the UCS, which was calculated according to ASTM D7012-14 (2014) and DIN 18141-1:2014-05:

(4) UCS = F A ,

where F is the load at failure [N] and A is the cross-sectional area of the sample [mm2]. Whenever the plugs were shorter than 80 mm and did not fulfill the required 2:1 length / diameter ratio, a correction function was applied as proposed by DIN 18141-1:2014-05:

(5) σ U ( 2 ) = 8 σ U 7 + 2 d l ,

where σU(2) is the corrected UCS [MPa] and σU the measured UCS [MPa], respectively, and d is the sample diameter [mm], while l denotes its length [mm]. At TU Darmstadt the destructive tests using the hydraulic uniaxial press were performed “force controlled” with a maximum loading rate of 0.5 kN s−1. The exceptions form very soft or fragile samples, such as ignimbrites, pumice or intensively fractured limestones. For these samples, the loading rate was individually reduced to 0.25 or 0.1 kN s−1 to meet the test requirements and to ensure the minimal test duration (e.g., 3 min for UCS and tensile strength).

For the determination of the static Young's modulus and Poisson's ratio, cyclic uniaxial tests were performed on three plugs (same dimension as described above) for each sample according to DIN 18141-1:2014-05 and Mutschler (2004). In order to record the axial displacement and lateral extension of the plug, three vertical and three lateral displacement transducers (LVDTs) were installed at an angle of 120 around the plug. The measurement was conducted in two cycles with the first cycle reaching 40 % and the second cycle reaching 60 % of the previously determined UCS from the same sample set. For intensively fractured limestones, the maximum load of the cycles was individually reduced to 30 % and 50 % of the previously determined UCS, respectively, to avoid an early rock failure and possible damage of the sensors. According to Mutschler (2004) a holding time of 5 min was set at the maximum value of each cycle. After the end of the holding time of the second cycle, the sensors were removed and the sample was loaded until failure to obtain the UCS. Using the results of the first unloading cycle, the static Young's modulus (average modulus) of each plug was calculated as the difference in stress divided by the difference in the vertical deformation according to ASTM D3148-02 (2002). Likewise, the static Poisson ratio was calculated as the ratio of lateral deformation and original diameter divided by the ratio of vertical deformation and original plug length. Subsequently, the G modulus, G, and bulk modulus, K, were calculated after ASTM D7012-14 (2014):


where E is the Young's modulus [N mm−2 or MPa] and μ is the Poisson ratio [–].

Furthermore, simple uniaxial tests were performed at TU Delft and UNAM to determine UCS, the static Young's modulus and the static Poisson ratio using a uniaxial stress–strain device with a capacity of 500 and 250 kN, respectively (GDSVIS load frame, GDS instruments, UK). Plugs with a dimension of 30 mm in diameter and a length of 75 mm drilled from marble, skarn, granodiorite and limestone samples from Las Minas were tested with a loading rate of 0.15 kN s−1 at TU Delft, while plugs with a dimension of 53 mm in diameter and a length of ∼110 mm drilled from volcanic rocks from Acoculco were analyzed at UNAM (displacement controlled with 0.05 mm min−1). Local axial and radial strains at UNAM were measured by the GDS LVDT local strain transducers, while at TU Delft axial displacement was recorded using two LVDTs and radial displacement was recorded using a radial chain with an LVDT sensor around the plugs. UCS, the static Poisson ratio and the static Young's modulus (TU Delft – tangent modulus; UNAM – secant modulus at 50 % of UCS) were calculated as described above following the ASTM guidelines (ASTM D3148-02; 2002).

Tensile strength of the sample material was determined at TU Darmstadt and TU Delft performing the indirect tensile test, also called the Brazilian test, according to ASTM D3967-16 (2016) and Lepique (2008). Cylindrical plugs with diameters of 55 and 40 mm (TU Darmstadt) and 30 mm (TU Delft) and a diameter / length ratio of 2:1 were loaded in a hydraulic uniaxial press by a linear distributed load until failure (diametrical compression). Afterwards the tensile strength of the plug was calculated using the following equation:

(8) σ t = 2 F π d l ,

where σt is the tensile strength [N mm−2 or MPa], F the load at failure [N], d the diameter [mm] and l the sample length [mm].

Fracture toughness was then calculated for granite, limestone, marble and skarn samples analyzed at TU Delft after Guo et al. (1993). In order to obtain more precise values, further chevron bend tests were performed on the same sample material at TU Delft. The tests were performed on cylindrical plugs with a length of 15 mm and a diameter of 30 mm using the uniaxial device following the methods proposed by ISRM (1988). Fracture toughness (KIc) of the sample material was determined first using a direct loading to failure (equal to KIc at Level I) and secondly using cyclic loading to calculate the correction of fracture toughness for non-linearity (equal to KcIc at Level II).

Additionally, point load tests were performed at UNAM in order to correlate the results to the tensile and uniaxial strength as proposed by ASTM D731-18 (2018). The tests were performed following the ISRM 325-89 (1984) and ASTM D5731-08 (2008) guidelines using a point load device from Controls (model 0550) with a maximum capacity of 100 kN. Therefore, cylindrical plugs with diameters of 25 mm and lengths ranging between 25 and 55 mm were jacked in a neoprene membrane during the test to confine the specimen and to avoid the fragmentation due to impacts with the ground.

Triaxial compression tests were performed on oven-dry samples at TU Darmstadt using a hydraulic triaxial press (Wille Geotechnik, Germany) with a capacity of 500 kN in order to determine the friction angle (φ), cohesion (c), shear (τ) and normal stress (σn) of the sample material. Depending on the availability, three plugs (diameter of 55 mm, length of 110 mm) for each sample were tested using different confining pressures (σ3) of 10, 20 and 30 MPa, respectively. According to ASTM D2664-04 (2004) the confining pressures and resulting vertical stresses (σ1) were transferred into a shear stress diagram to construct the Mohr–Coulomb criterion of failure to derive cohesion (intersection with the vertical axis) and the friction angle (the angle between the line and the horizontal axis). Whenever needed, the vertical stresses from UCS tests (with σ3=0) were considered to construct an additional circle in the shear stress diagram, thus enhancing the data evaluation.

5.4 Chemical analyses

In order to perform quantitative and qualitative chemical analyses, representative composite sample material from selected outcrop samples and the reservoir core samples was milled with a disc swing mill (Siebtechnik, Germany) for 2.5 min at 1000 rpm at TU Darmstadt and with a colloid mill (Mixer Mill MM301, Retsch GmbH, Germany) for about 1 min at TU Delft to obtain a grain size smaller than 63 µm.

XRD analyses at TU Delft and GFZ were performed using a Bruker D8 Advance diffractometer (Bruker, Karlsruhe, Germany) and the software DIFFRAC.EVA (TU Delft) and Match! (GFZ) for data evaluation. For XRF measurements at TU Delft, a Panalytical Axios Max WDXRF spectrometer was used and data evaluation was performed with SuperQ5.0i/Omnian software. In addition to the Omnian standards, many NIST SRM samples and pure compounds were used for calibration. At GFZ the XRF measurements were performed with a Panalytical Axios Advanced spectrometer in combination with the software SuperQ. For the analysis three reference standards (basalt ZGI-BM, granite ZGI-GM and shale ZGI-TB) were used. At TU Darmstadt, major and trace elements were analyzed with a Bruker S8 TIGER 4 WDXRF spectrometer using the Quant Express method. Accuracy is <5 % for the major elements and <10 % for the trace elements. The proposed limit of detection ranges between 400 ppm (Na) and 10 ppm (e.g., Rb, Sr, Nb). Further XRD analyses were performed at UniTO using a Siemens D5000 automatic X-ray diffractometer. The qualitative interpretation of the data has been realized with the software DIFFRACplus EVA Application (2001), by comparing the positions and intensity of the data with suitable databases (ICDD, previously JCPDS; ICSD; PCPDFWIN).

6 Status of the database

The database presented here comprises petrophysical and mechanical rock properties of outcrop samples and reservoir core samples of two caldera complexes located in the northeastern part of the TMVB. So far, the database comprises 31 982 data entries (Table 2) as a result of 34 properties determined for 2169 plugs and rock samples (2138 cylindrical plugs and 31 uncored samples). Destructive tests were conducted on more than 970 plugs. In addition, 133 XRF and 113 XRD analyses were performed.

In total 380 samples were analyzed covering volcanic rocks (950 plugs), sedimentary rocks (716 plugs), igneous rocks (147 plugs) and metamorphic rocks (356 plugs). Thereof, 80 outcrop samples were collected for Acoculco and 226 outcrop samples were collected for Los Humeros, resulting in 563 and 1606 analyzed plugs and samples including the reservoir core samples, respectively. The difference between the number of collected samples for Los Humeros and Acoculco is biased due to the purposes of the different field trips and the targets of the project. The main targets for the development of a deep EGS in Acoculco and SHGS in Los Humeros are marbles and skarns (AC3 and AC2) and the pre-caldera andesites and Cretaceous limestones and marbles (G3 and G4), respectively. As the basement rocks (AC1 to AC3) are not exposed in Acoculco, the exhumed systems were used as analogs. Therefore, the main attention was paid to Las Minas where 101 samples were collected (here associated with Los Humeros). In Las Minas it is possible to investigate the igneous bodies and their metamorphic products like skarn, hornfels or marble (Fuentes-Guzmán et al., 2020) as well as some outcrops belonging to the metamorphic basement below the Cretaceous and Jurassic units.

Figure 7Number of collected samples (outcrop and reservoir core samples) per model unit for the regional and local models of the Los Humeros (a) and Acoculco (b) geothermal systems.


The samples were classified regarding their model units as shown in Fig. 7. Following this approach almost all local model units for Los Humeros were covered. For some samples a classification is not possible at this stage of the project. Ongoing volcanological studies are underway and further dating is planned to overcome these knowledge gaps. The outcrop samples belonging to the pre-caldera group predominantly represent the Teziutlán andesite unit (U6) and the Cuyoaco andesite unit (U8). U5 comprises ignimbrites and pumice layers from the Xaltipán ignimbrite unit, while very recent basaltic lavas, ashfall deposits and ignimbrites collected within the Los Humeros caldera are associated with the post-caldera group (G1). The basement comprises a wide range of different rock types. G4 includes Jurassic sandstones and limestones; Cretaceous limestones, marls and shales; and Miocene granitic and granodioritic intrusive bodies and their metamorphic products marble and skarn. Regarding the regional model of Acoculco, outcrop samples from the two upper units AC5 and AC4 were collected. The uppermost unit comprises all volcanic deposits from the pre-caldera volcanics to the extra-caldera volcanism. Among others, samples from the Acoculco ignimbrite, Terrerillos andesite lava, Manzanito andesite and Perdernal rhyolitic lava were collected. The unit AC4 includes Jurassic limestones and sandstones and Cretaceous limestones. The reservoir core samples from well EAC1 cover ignimbrite (core 1), dacitic to rhyolitic lavas (core 2 and 3), skarn (core 4), marble (core 5) and granodiorite (core 5).

The number of measurements for each parameter resulted from the availability of measurement devices at the different institutes, required sample size, sample preparation and test duration as well as test setup. While most of the non-destructive parameters were analyzed on each plug, more time-intensive tests, such as specific heat capacity measurements or XRF and XRD analyses, were performed for each sample only (composite sample material). Likewise, rock mechanical tests are significantly more time-consuming as they require a specific sample size and sample preparation or in the case of triaxial tests a minimum number of samples to evaluate the test results. Although the total number of measurements significantly differs between some parameters, all parameters were analyzed on sample sets covering all relevant lithologies in the study area.

Table 2Number of measurements for each parameter.

Download Print Version | Download XLSX

7 Discussion

7.1 Data availability and data application

Rock properties are commonly used for reservoir exploration, assessment and modeling. While petrophysical, dynamic and static mechanical properties are the primarily used parameters for reservoir exploration, production and stimulation scenarios (Saller and Henderson, 1998; Rybacki et al., 2016; Gan and Elsworth, 2016; Ghassemi, 2017; Qu et al., 2019; Scott et al., 2019; Bohnsack et al., 2020), thermal properties are of great importance to assessing the subsurface temperature, the geothermal gradient, heat transport and heat storage (Weides et al., 2013; Weides and Majorowicz, 2014; Ebigbo et al., 2016; Franco and Donatini, 2017; Nurhandoko et al., 2019; Békési et al., 2020). Especially in active high-enthalpy hydrothermal systems, electric resistivity and magnetic susceptibility data are very useful to identify or map the cap rock and different lithologies or hydrothermally altered zones within the reservoir (Oliva-Urcia, 2011; Lévy et al., 2018, 2019), whereas high-T/P and detailed mineralogical studies help to estimate rock properties in reservoir conditions (Nono et al., 2020; Kummerow et al., 2020; Lacinska et al., 2020).

Within the scope of the GEMex project, petrophysical and rock mechanical data were used for various purposes. Deb et al. (2019a) used petrophysical and thermophysical properties to parameterize the structural model of Los Humeros and Acoculco (Calcagno et al., 2018) for simulating the initial state of the super-hot geothermal system. Several stimulation scenarios were investigated to evaluate the potential of the basement rocks in Acoculco for the development of an EGS (Deb et al., 2019b). Based on the fracture network characterization of outcrop analogs in Las Minas and petrophysical and rock mechanical data, Lepillier et al. (2019) created FEM models to calculate the fluid flow and heat exchange of fracture-controlled reservoirs in marble, skarn and limestone as an equivalent to the deep subsurface of Acoculco. Kruszewski et al. (2021) used rock mechanical parameters together with well parameters and geophysical logs to estimate the local stress field of the Acoculco geothermal field. Current studies focus on fracture propagation models and hydraulic fracture stimulation scenarios to estimate fracture geometries. The results of the petrophysical properties and volcanological studies are being used to interpret results of electric resistivity surveys (Benediktsdóttir et al., 2020), local earthquake tomography (Toledo et al., 2020), or gravity and magnetotelluric surveys (Cornejo, 2020).

Compared to siliciclastic or carbonate basins used for oil and gas exploitation, the number of petrophysical and mechanical rock property data for volcanic settings in the context of high-enthalpy geothermal systems is less documented.

So far, geothermal exploration studies in volcanic settings have provided rock properties analyzed on outcrop (e.g., Lenhardt and Götz, 2011; Pola et al., 2014; Mielke et al., 2016; Heap and Kennedy, 2016; Navelot et al., 2018; Mordensky et al., 2019; Eggertson et al., 2020) or reservoir core (Stimac et al., 2004; Siratovich et al., 2014; Ólavsdóttir et al., 2015; Mielke et al., 2015; Cant et al., 2018) samples. However, this study highlights the importance of the analysis of both outcrop and reservoir core samples. The comparison of reservoir samples, exhumed systems and outcrops in the surrounding area enables the identification of the processes that occurred within the reservoir and quantifying the impact on the properties correctly.

The need for valuable input data for reservoir modeling and assessment has recently led to an increased number of studies and publications (Bär et al., 2020). While several extensive national or global databases have already been developed and published for geothermal well data (National Geothermal data system NGDS, 2014; BritGeothermal, 2017; DOE Data Explorer, 2018); rock chemistry; geochronology; petrology; petrophysical data such as porosity, density or magnetic susceptibility derived from geophysical borehole data (Petlab, 2020; Sciencebase Minnesota, 2010; Georoc Mainz, 2020; Rock Properties Database British Columbia Canada, 2018; global whole-rock geochemical database compilation in Gard et al., 2019; National Geochemical Database USGS, 2014; the North American Volcanic and Intrusive Rock Database NAVDAT data base, 2020); lithology (the new global lithological map database GLiM in Hartmann and Moosdorf, 2012); mineralogy (BRITROCKS project, 2020); and petrography (RockViewer, 2020), a comprehensive and quality-proofed collection of laboratory rock properties has just recently been released by Bär et al. (2020; not considering fee-based and non-open-access databases that exist for oil and gas data like the AccuMap or geoSCOUT databases; IHS Markit, 2020; GeoScout, 2020). The PetroPhysical Property Database (P3) presented in Bär et al. (2020) collected rock property data from 316 research articles and comprises 75 573 data points of 28 different rock properties analyzed on a wide variety of lithologies worldwide. While the P3 database significantly increases the availability of standardized rock properties, it still contains a limited number of data points or parameters for each investigated area or formation. To increase the level of detail for the GEMex study area to the required spatial and stratigraphic coverage, the database presented in this paper contains more than 31 000 data points and 34 different parameters covering all important lithologies from the basement to the cap rock. The high number of analyzed plugs and samples enables detailed statistical and spatial geostatistical analyses on different scales (plug, sample, outcrop, formation or model unit), spatial evaluation of the results in 2D or 3D, and the validation of different analytical methods. Whenever possible, all parameters were analyzed on each plug. This approach allows the identification of statistical and causal relationships between the parameters and, thus, improves the accuracy of geostatistical predictions, which are crucial for upscaling or downscaling (see next section; Linsel et al., 2020). The usage of plugs with different dimensions (drilled diameter ranges from 25 to 65 mm with a length from ∼12 mm to 30 cm) enables the identification of scale effects, which need to be considered for the evaluation of dynamic mechanical properties (Bayuk and Tikhotsky, 2018). The level of detail presented in this study has not only significantly improved the geological understanding of both geothermal systems and super-hot geothermal systems in general but also helped in the better understanding of the relationship between different parameters and how they are affected by different processes (e.g., fracturing or hydrothermal alteration). The database not only provides the basis for ongoing research in the study area but also facilitates various applications in comparable geological settings within the TMVB or similar volcanic geothermal play types worldwide. Combined with other data sets (P3 in Bär et al., 2020, or Weinert et al., 2021), these data could be used to train machine learning algorithms to develop rock property prediction tools to improve and speed up parametrization of 3D geological models in the future.

7.2 Data processing and upscaling

The database presented in this study includes laboratory data analyzed on core and outcrop samples (centimeter to decimeter scale defined here as mesoscale), thus representing rock matrix properties only (with small-scale or single fractures in few samples). Oven-dried samples were analyzed under ambient laboratory conditions (room temperature of ∼21C and atmospheric pressure of 0.1 MPa) to standardize the test procedure and to ensure the comparability of the results for the different samples and rock types. Consequently, the data do not reflect in situ conditions such as high reservoir temperatures, overburden pressure, confining pressure and fluid properties at reservoir depth. Depending on the aim and scale of future applications, the data need to be corrected for reservoir conditions and transferred to the reservoir scale (macroscale). Hydraulic properties such as porosity and permeability tend to decrease with increasing stress and pressure at reservoir depth by closing fractures and compaction of the rock mass (rock compressibility; Zimmermann et al., 1986; Moosavi et al., 2014; Hatakeda et al., 2017; You et al., 2020), often also resulting in increased bulk density, heat conduction, electric resistivity and wave velocities (Horai and Susaki, 1989; Clauser and Huenges, 1995; Schön, 2015). However, the relationships between different properties related to temperature and pressure changes are complex. At higher temperatures thermal expansion of minerals can cause microfracturing, which again negatively affects thermal conductivity, ultrasonic wave velocities and rock strength (Heap et al., 2014; Vinciguerra et al., 2005) but increases hydraulic properties. Several analytical and empirical relationships and correction functions have been identified and developed in the past to transfer hydraulic (Zimmermann et al., 1986; Li et al., 2004; Zheng et al., 2015; Heap and Kennedy, 2016), thermal (Sass et al., 1971; Zoth and Hänel, 1988; Somerton, 1992; Vosteen and Schellschmidt, 2003; Hartmann et al., 2005; Whittington et al., 2009; Rühaak et al., 2015; Zhao et al., 2016; Merriman et al., 2018; Norden et al., 2020; Clauser, 2020), magnetic (Ohnaka, 1969; Ali and Potter, 2012; Zhang et al., 2020), electric (Shankland et al., 1997; Hatakeda et al., 2017; Kummerow and Raab, 2015; Kummerow et al., 2020; Nono et al., 2020) and mechanical (Mobarak and Somerton, 1971; Vinciguerra et al., 2005; Siratovich et al., 2011; Heap et al., 2014; Hassanzadegan et al., 2013; Vagnon et al., 2021) properties from laboratory to reservoir conditions. Transferring rock properties from core sample to the reservoir scale is challenging and has been the focus of numerous studies in the past (Christie, 1996; Farmer, 2002; Qi and Hesketh, 2005; Khajeh, 2013). Even though computer processing capacities have drastically increased over the past decades, the resolution (number of grids) and complexity of static geological models often tend to be too high to run numerical reservoir simulations, which solve complex, e.g., fluid or heat flow, equations. Thus, upgridding and upscaling techniques are required that retain as much of the original structure, geometry, petrophysical characteristics and facies heterogeneity as possible to deliver the vital information needed for reservoir assessment and operation (Walia and Leahy, 2014). Existing upscaling approaches can be grouped into direct or two-step and local or global upscaling methods (Wen and Gomez-Hernandez, 1996; Farmer, 2002). The most common upscaling techniques are simple cross correlations, (power-law) averaging (arithmetic, geometric or harmonic averaging often in combination with Monte Carlo techniques), renormalization, pressure-solver or tensor methods, and pseudofunctions (Qi and Hesketh, 2005). However, particularly the first-mentioned techniques tend to spatially smear out extremes within the reservoir, such as flow barriers or open fractures, and thus are not very useful for complex and heterogenous reservoirs (Ding et al., 1992; Qi and Hesketh, 2005). Geostatistical analyses and modeling using estimation algorithms (e.g., variogram analyses and kriging techniques) or sequential simulations (e.g., Gaussian simulation) have been applied to populate numerical models in geologically complex and/or fractured reservoirs (Hartanato, 2004; Bourbiaux et al., 2005; Ebong et al., 2019). However, integrating geological information regarding the geometry, distribution and connectivity of faults and fractures as well as linking fracture and matrix properties and fluid flow remains challenging (multiphase and dual-porosity modeling; Bourbiaux, 2010). Since hydrothermal alteration significantly influences the matrix properties (Heap et al., 2020), estimating the size and spatial distribution of hydrothermal aureoles along fractures in active volcanic settings becomes important to improving the accuracy of the reservoir model. While upscaling of hydraulic properties with application to oil and gas reservoirs has been intensively analyzed in the past (Wen and Gomez-Hernandez, 1996; Farmer, 2002; Sánchez-Vila et al., 2006), relatively little work has been done for thermal properties (Scheibe and Yabusaki, 1998; Hartmann et al., 2005; Rühaak et al., 2015). According to Rühaak et al. (2015) upscaling thermal conductivity can be fundamentally different from upscaling hydraulic or other transport parameters in porous media and rocks. The authors found that harmonic- and geometric-mean upscaled values most accurately reflect local values. Rühaak et al. (2014) and Gu et al. (2017) recommend kriging with external drift (KED) to interpolate subsurface temperature and thermal conductivity, respectively.

7.3 Limitations with respect to modeling the Los Humeros and Acoculco geothermal systems

Besides the many advantages described above, a number of limiting factors have to be considered prior to using this data set for modeling the Los Humeros and Acoculco geothermal systems. The fieldwork and the results of the petrophysical measurements revealed the complexity of both geothermal systems. Composition, lateral extension and distribution of the volcanic sequences are very variable within the study area. Furthermore, the basement rocks showed a high geological heterogeneity comprising several different rock types including shales, limestones, sandstones, intrusive bodies, marble and skarn. The definition of the preliminary model units is predominantly based on the local stratigraphy of the study area (Calcagno et al., 2018), and some model units comprise multiple different rock types. The results of the petro- and thermophysical properties however reveal high variability and a wide parameter range for individual units leading to high uncertainties during modeling. For this reason, the results for each lithostratigraphic unit were weighted with respect to their relative contribution in the study area for the population of the geological model of Los Humeros (Deb et al., 2019a), which was mainly based on lithostratigraphic well descriptions provided by the CFE. As this is not known in detail for every model unit, the relative contribution of each rock type was based on field observations.

The number of samples per unit strongly depended on the quality, availability and accessibility of representative outcrops in the field or reservoir core samples in the core storage. Thus, it was not possible to cover all local model units for Los Humeros.

Likewise, the number of measurements for each parameter was strongly affected by the availability of measurement devices, sample preparation and test duration. Although the data for each parameter cover all key lithologies in the study area, future work should focus on additional electric resistivity and rock mechanical tests (fracture toughness and triaxial tests) to better support the interpretation of MT, TEM and DC surveys or 3D geomechanical models. Furthermore, further research is needed on HT–HP experiments reaching supercritical conditions to better evaluate the processes within the reservoir and to transfer rock properties from laboratory to reservoir conditions of super-hot geothermal systems.

The core samples of the Los Humeros geothermal field were predominantly retrieved from the reservoir pre-caldera andesite units. They show high matrix variability due to hydrothermal alteration of different intensities, which caused significant differences regarding petrophysical and thermophysical properties compared to the equivalent outcrop samples. For about one-quarter of the samples, intensive hydrothermal alteration prevents a clear identification of the original rock type and correlation to equivalent units in the outcrops. This suggests that a comprehensive identification and characterization of the hydrothermal alteration aureoles in the geothermal fields is also required for the accurate assessment and modeling of these systems (e.g., by MT sounding or other direct or indirect analyses). Current studies on the reservoir core samples including detailed petrographic analyses and ICP-MS measurements aim to provide a better sample description and classification (Weydt et al., 2020, 2021). Only a few reservoir core samples were available representing the overlaying cap rock (Xaltipán ignimbrite) or the basement below. While the Xaltipán ignimbrite unit can be investigated in several outcrops around the Los Humeros caldera, the deeper part of the basement remains mostly unknown. The high number of collected samples in the exhumed systems and in the surrounding area of the caldera complexes greatly depicts the heterogeneity of the basement. However, the analyses of outcrops and the few reservoir core samples only cover the upper limited parts of the basement (approximately tens to hundreds of meters). Thus, in the field it is not possible to investigate the spatial extension of the intrusive bodies within the (meta)sedimentary basement. However, Urbani et al. (2020) concluded that the recent uplift within the Los Proteros caldera was caused by multiple intrusive bodies at a very shallow depth (425±170 to <1000 m). Likewise, in Acoculco several intrusive bodies had already been identified at 1000 m depth (below ground level; Avellán et al., 2020).

Regarding the regional model of Acoculco, only rocks of the two upper units are exposed in the field. For the parameterization of the remaining units, the project emphasizes using the exhumed system in Las Minas as an analog. Regarding the results of the petrophysical measurements, this concept can be applied for almost all units. However, the sedimentary sequences reveal the highest variability compared to other units comprising argillaceous mudstones to dolomitic marbles. The properties of the limestones and marbles resemble the different facies and diagenetic or metamorphic overprint. In Las Minas the limestones and marbles comprise dolomite, while the reservoir core samples from Los Humeros and most of the limestones collected from the outcrops in the surrounding area of both systems represent undolomitized, marine, fine-grained mudstones to wackestones. In addition, the reservoir core samples from the upper part of the carbonatic basement show intensive fracturing and recrystallization as a result of the complex tectonic activity caused by caldera collapses, uplift and ascending lavas. Furthermore, the term “skarn” has been widely used in the literature (related to the study area) without a precise description. The skarns in Las Minas commonly resemble Fe-rich ore deposits in close proximity to intrusive bodies. In contrast the units classified as skarn within the upper parts of the geothermal reservoirs (López-Hernández et al., 2009) formed instead due to intensive metasomatic processes caused by Ca-rich fluids migrating into the overlaying lavas. Once more, the physical properties reflect the different mineralogical composition of both skarn types.

8 Data availability

The data repository is available at (Weydt et al., 2020).

9 Conclusions

Within the scope of the GEMex project, an extensive rock property database was created comprising more than 31 000 data entries covering a great variety of different rock types and lithologies of Jurassic to Holocene age. The database includes petrophysical, thermophysical, magnetic, electric, and dynamic and static mechanical properties complemented by the results of XRF and XRD analyses. In total 34 properties were determined on 2169 plugs retrieved from more than 300 outcrop samples collected from the Acoculco and Los Humeros caldera complexes, 66 reservoir core samples drilled from 37 core sections from 16 wells of the Los Humeros geothermal field, and 8 core samples drilled from 6 core sections obtained from well EAC1 of the Acoculco geothermal field. The database was created in a simple and transparent format including comprehensive metainformation to facilitate application in various geoscientific disciplines worldwide.

The compiled data set allows for the

  • prediction of rock properties of target formations in the subsurface at early exploration stages or in the case of low data density,

  • assessment of the reservoir potential and estimation of economic risks and uncertainties,

  • population of 3D geological models (numeric thermo-hydraulic–mechanical–chemical (THMC) models),

  • statistical evaluation to identify relationships between the properties and trends required for upscaling approaches, and

  • validation of different analytical methods.

The data and workflow presented here will improve the planning and execution of future research projects. Outcrop analyses and the characterization of petrophysical and mechanical properties of outcrop and reservoir core samples are paramount for profound reservoir characterization and should in general be considered in future geoscientific studies to a greater extent to enable a more precise prediction of reservoir properties. Hereby, an integration of shallow geophysical and classical (e.g., scan-line) or state-of-the-art (lidar) fracture network characterization methods has great potential to further enhance 3D reservoir characterization.

The current structure of the database allows for easy modification and extension. It is planned to create an outcrop catalogue of all field campaigns conducted within GEMex and to improve it by adding the results of ongoing ICP-MS and detailed petrographic analyses.

Author contributions

This study was conducted by LMW, who analyzed samples, evaluated and compiled data, and wrote the manuscript. AARG, AP, BL, JK, GM, CC and PD contributed to the database by providing results of rock property measurements and chemical analyses. LMW, AARG, AP, BL, GM, CC, KB, GN, EGP, DRA and JLM were involved in the sampling campaigns. GN, EGP, DRA and JLM provided the main input to the geological interpretation and sample classification. The concept of the study as well as the manuscript outline and composition was designed by LMW, KB and IS. Rock characterization was coordinated by LMW, KB and IS and AP for the European and Mexican consortium, respectively. All authors contributed to this study and reviewed the manuscript.

Competing interests

The authors declare that they have no conflict of interest.


We thank Miguel Angel Ramírez Montes, Subgerencia de Estudios, Gerencia de Proyectos Geotermoeléctricos and the Comisión Federal de Electricidad (CFE) team for providing us with access to the core storage and for their help during sampling and drilling the CFE core samples. We highly appreciated working at the Los Humeros camp.

Furthermore, we thank Jana Perizonius, Thomas Kramer, Maximilian Bech and Roland Knauthe (master's and bachelor's graduates at TU Darmstadt) for their contribution to this project.

Special thanks go to Cord Peters, Patrick Höfler, Fatih Ekinci, Ruud Hendrikx, Gabriela Schubert, Dirk Scheuvens, Rainer Seehaus, Reimund Rossmann, Georg Wasmer, Angelo Agostino, Jessica Chicco, Chiara Colombero, Anna Ferrero, Sergio Vinciguerra and Federico Vagnon for their great support in the laboratory.

We also thank Christopher Rochelle, Domenico Liotta and his team, Víctor Hugo Garduño-Monroy and his students, Geovanny Hernández-Avilés, Luís and Daniel González-Ruiz, Irais Franco, and Gerardo Carrasco-Núñez for their support during the field campaigns. Many thanks go to Caterina Bianco for providing Fig. 1c.

Financial support

This project has received funding from the European Union's Horizon 2020 research and innovation programme (grant no. 727550) and the Mexican Energy Sustainability Fund CONACYT-SENER (project 2015-04-68074).

Review statement

This paper was edited by Kirsten Elger and reviewed by Léa Lévy and one anonymous referee.


Ali, A. and Potter, D. K.: Temperature dependence of the magnetic properties of reservoir rocks and minerals and implications for in situ borehole predictions of petrophysical parameters, Geophysics, 77, WA211–WA221,, 2012. 

Aretz, A., Bär, K., Götz, A. E., and Sass, I.: Outcrop analogue study of Permocarboniferous geothermal sandstone reservoir formations (northern Upper Rhine Graben, Germany): Impact of mineral content, depositional environment and diagenesis on petrophysical properties, Int. J. Earth Sci., 135, 1431–1452,, 2015. 

Árnason, K.: New Conceptual Model for the Magma-Hydrothermal-Tectonic System of Krafla, NE Iceland, Geosciences, 10, 27 pp.,, 2020. 

ASTM D2664-04: Standard Test Method for Triaxial Compressive Strength of Undrained Rock Core Specimens Without Pore Pressure Measurements (Withdrawn 2005), ASTM International, West Conshohocken, PA, USA, 4 pp.,, 2004. 

ASTM D2845-08: Standard Test Method for Laboratory Determination of Pulse Velocities and Ultrasonic Elastic Constants of Rock, ASTM International, West Conshohocken, PA, USA, 7 pp.,, 2008. 

ASTM D3148-02: Standard Test Method for Elastic Moduli of Intact Rock Core Specimens in Uniaxial Compression, ASTM International, West Conshohocken, PA, USA, 6 pp.,, 2002. 

ASTM D3967-16: Standard Test Method for Splitting Tensile Strength of Intact Rock Core Specimens, ASTM International, West Conshohocken, PA, USA, 5 pp.,, 2016. 

ASTM D4525-13e2: Standard Test Method for Permeability of Rocks by Flowing Air, ASTM International, West Conshohocken, PA, USA, 5 pp.,, 2013. 

ASTM D4543-19: Standard Practices for Preparing Rock Core Specimens and Determining Dimensional and Shape Tolerances, ASTM International, West Conshohocken, PA, USA, 13 pp.,, 2019. 

ASTM D5550-14: Standard Test Method for Specific Gravity of Soil Solids by Gas Pycnometer, ASTM International, West Conshohocken, PA, USA, 5 pp.,, 2014. 

ASTM D5731-08: Standard Test Method for Determination of the Point Load Strength Index of Rock and Application to Rock Strength Classifications, ASTM International, West Conshohocken, PA, USA, 12 pp.,, 2008. 

ASTM D6539-13: Standard Test Method for Measurement of the Permeability of Unsaturated Porous Materials by Flowing Air, ASTM International, West Conshohocken, PA, USA, 10 pp.,, 2013. 

ASTM D7012-14: Test methods for compressive strength and elastic moduli of intact rock corespecimens under varying states of stress and temperatures, ASTM International, West Conshohocken, PA, USA, 9 pp.,, 2014. 

ASTM D7263-16: Standard Test Methods for Laboratory Determination of Density (Unit Weight) of Soil Specimens, ASTM International, West Conshohocken, PA, USA, 7 pp.,, 2016. 

ASTM D731-18: Standard Test Method for Molding Index of Thermosetting Molding Powder, ASTM International, West Conshohocken, PA, USA, 4 pp.,, 2018. 

Avellán, D. R., Macías, J. L., Layer, P. W., Sosa-Ceballos, G., Cisneros, G., Sanchez, J. M., Martha Gómez-Vasconcelos, G., López-Loera, H., Reyes Agustín, G., Marti, J., Osorio, S., García-Sánchez, L., Pola-Villaseñor, A., García-Tenorio, F., and Benowitz, J.: Geology of the Pleistocene Acoculco Caldera Complex, eastern Trans-Mexican Volcanic Belt (México), J. Maps, 15, 8–18,, 2018. 

Avellán, D. R., Macías, J. L., Layer, P. W., Sosa-Ceballos, G., Gómez-Vasconcelos, M. G., Cisneros-Máximo, G., and Benowitz, J.: Eruptive chronology of the Acoculco caldera complex – A resurgent caldera in the eastern Trans-Mexican Volcanic Belt (México), J. S. Am. Earth Sci., 98, 102412,, 2020. 

Bär, K.: Untersuchung der tiefengeothermischen Potenziale von Hessen, PhD thesis, Technische Universität Darmstadt, Germany, 265 pp., 2012. 

Bär, K. and Weydt, L. M.: Comprehensive report on the rock and fluid samples and their physical properties in the Acoculco and Los Humeros regions, Deliverable D6.1, WP6, GEMex H2020 project, European Comission, 304 pp., available at: (last access: 25 May 2020), 2019. 

Bär, K., Reinsch, T., and Bott, J.: The PetroPhysical Property Database (P3) – a global compilation of lab-measured rock properties, Earth Syst. Sci. Data, 12, 2485–2515,, 2020. 

Bayuk, I. and Tikhotskiy, S.: Upscaling and downscaling of reservoir elastic properties – Rock Physics approach, in: SEG International Exposition and 88th Annual Meeting, 14–19 October 2018, Anaheim, CA, USA, 3653–3657,, 2018. 

Békési, E., Struijk, M., Bonté, D., Veldkamp, H., Limberger, J., Fokker, P. A., Vrijlandt, M., and van Wees, J.-D.: An updated geothermal model of the Dutch subsurface based on inversion of temperature data, Geothermics, 88, 101880,, 2020. 

Benediktsdóttir, A., Arango-Galván, C., Páll Hersir, G., Held, S., Romo-Jones, J. M., Luis-Salas, J., Avils, T., Ruiz-Aguilar, D., and Már Vilhjálmsson, A.: The Los Humeros superhot geothermal resource in Mexico: Results from an extensive resistivity survey, GEMex Final Conference, Potsdam, Germany, 18–19 February 2020, S6.1, 2020. 

Bohnsack, D., Potten, M., Pfrang, D., Wolpert, P., and Zosseder, K.: Porosity-permeability relationship derived from Upper Jurassic carbonate rock cores to assess the regional hydraulic matrix properties of the Malm reservoir in the South German Molasse Basin, Geotherm. Energy, 12, 1–47,, 2020. 

Bourbiaux, B.: Fractured Reservoir Simulation: a Challenging and Rewarding Issue, Oil Gas Sci. Technol., 65, 227–238,, 2010. 

Bourbiaux, B., Basquet, R., Daniel, J. M., Hu, L. Y., Jenni, S., Lange, G., and Rasolofosaon, P.: Fractured reservoirs modelling a review of the challenges and some recent solutions, First Break, 23, 33-40,, 2005. 

BritGeothermal: Deep geothermal energy research in the UK, available at: (last access: 21 October 2020), 2017. 

BRITROCKS project: BRITROCKS Rock collections: BGS Mineralogy and petrology collection database, available at: (last access: 21 October 2020), 2020. 

Buntebarth, G.: Geothermie, Springer, Berlin ans Heidelberg, Germany,, 1980. 

Calcagno, P., Evanno, G., Trumpy, E., Gutiérrez-Negrín, L. C., Macías, J. L., Carrasco-Núñez, G., and Liotta, D.: Preliminary 3-D geological models of Los Humeros and Acoculco geothermal fields (Mexico) – H2020 GEMex Project, Adv. Geosci., 45, 321–333,, 2018. 

Calcagno, P., Trumpy, E., Gutiérrez-Negrín, L. C., Liotta, D., Carrasco-Núñez, G., Norini, G., Brogi, A., Garduño-Monroy, V. H., Benediktsdóttir, A., Gaucher, E., Toldeo, T., Páll Hersir, G., Manzella, A., Santilano, A., Gola, G., Macías, J. L., Vaessen, L., Evanno, G., and Arango-Galván, C.: 3D Geomodels of Los Humeros and Acoculco geothermal systems (Mexico) – H2020 GEMex Project: Methodology, products and feedback, GEMex Final Conference, Potsdam, Germany, 18–19 February 2020, S1.2, 2020. 

Canet, C., Arana, L., González-Pertida, E., Pi, T., Prol-Ledesma, R. M., Franco, S. I., Villanueva-Estrada, R. E., Camprubí, A., Ramírez-Silva, G., and López-Hernández, A.: A statistics-based method for the short-wave infrared spectral analysis of altered rocks: An example from the Acoculco Caldera, Eastern Trans-Mexican Volcanic Belt, J. Geochem. Explor., 105, 1–10,, 2010. 

Canet, C., Trillaud, F., Prol-Ledesma, R., González-Hernández, G., Peláez, B., Hernández-Cruz, B., and Sánchez-Córdova, M. M.: Thermal history of the Acoculco geothermal system, eastern Mexico: Insights from numerical modeling and radiocarbon dating, J. Volcanol. Geoth. Res., 305, 56–62,, 2015. 

Cant, J. L., Siratovich, P. A., Cole, J. W., Villeneuve, M. C., and Kennedy, B. M.: Matrix permeability of reservoir rocks, Ngatamariki geothermal field, Taupo Volcanic Zone, New Zealand, Geotherm. Energy, 2, 1–28,, 2018. 

Carrasco-Núñez, G., Gómez-Tuena, A., and Lozano, V. L.: Geologic map of Cerro Grande volcano and surrounding area, Central Mexico, Geological Society of America, Map and chart series, MCH 081, available at: (last access: 9 December 2020), 1997. 

Carrasco-Núñez, G., Hernández, J., De Léon, L., Dávilla, P., Norini, G., Bernal, J. P., Jicha, B., Jicha, B., Navarro, M., and López-Quiroz, P.: Geologic Map of Los Humeros volcanic complex and geothermal field eastern Trans-Mexican Volcanic Belt, terra digitalis, 1, 1–11,, 2017a. 

Carrasco-Núñez, G., López-Martínez, M., Hernández, J., and Vargas, V.: Subsurface stratigraphy and its correlation with the surficial geology at Los Humeros geothermal field, eastern Trans-Mexican Volcanic Belt, Geothermics, 67, 1–17,, 2017b. 

Carrasco-Núñez, G., Bernal, J. P., Dávilla, P., Jicha, B., Giordano, G., and Hernández, J.: Reappraisal of Los Humeros Volcanic Complex by New U/Th Zircon and 40Ar/39Ar Dating: Implications for Greater Geothermal Potential, Geochem. Geophy. Geosy., 19, 132–149,, 2018. 

Cavazos, J. and Carrasco-Núñez, G.: Anatomy of the Xáltipan ignimbrite at Los Humeros Volcanic Complex; the largest eruption of the Trans-Mexican Volcanic Belt, J. Volcanol. Geoth. Res., 392, 106755,, 2020. 

Christie, M.: Upscaling for reservoir simulation, J. Pet. Technol., 48, 1004–1010,, 1996. 

Clauser, C.: Thermal Storage and Transport Properties of Rocks, II: Thermal Conductivity and Diffusivity, in: Encyclopedia of Solid Earth Geophysics. Encyclopedia of Earth Sciences Series, edited by: Gupta, H. K., Springer, Cham, Switzerland, 1–20,, 2020. 

Clauser, C. and Huenges, E.: Thermal conductivity of rocks and minerals. Rock Physics & Phase Relations, in: Rock Physics and Phase relations: A Handbook of Physical Constants, edited by: Ahrens, T. J., American Geophysical Union, Washington, USA, 3, 105–126,, 1995. 

Clement, R., Bergeron, M., and Moreau, S.: COMSOL Multiphysics modelling for measurement device of electrical resistivity in laboratory test cell, in: Proceedings of the 2011 COMSOL Conference, 26–28 October 2011, Stuttgart, 6 pp., 2011. 

Codd, E. F.: A Relational Model of Data for Large Shared Data Banks, Communications of the ACM, 13, 377–387,, 1970. 

Cohen, K. M., Finney, S. C., Gibbard, P. L., and Fan, J.-X.: The ICS International Chronostratigraphic Chart, Episodes 36, 199–204, available at: (last access: 09 Decmeber 2020), 2013. 

Contreras, L. E., Domínguez, A. B., and Rivera, M. O.: Mediciones petrofisicas en nucleos de perforacion del campo geotermico Los Humeros, Geotermia, 6, 9–42, 1990. 

Cornejo, N.: Towards visualization ofthe reservoir settings in the Los Humeros and Acoculco geothermal fields using gravity, GEMex Final Conference, , Potsdam, Germany, 18–19 February 2020, S1.6, 2020. 

Deb, P., Knapp, D., Clauser, C., and Montegrossi, G.: Modeling Natural Steady-State of Super-Hot Geothermal Reservoir at Los Humeros, Mexico, in: Proceedings of the European Geothermal Congress 2019, Den Haag, the Netherlands, 11–14 June 2019, 6 pp., 2019a. 

Deb, P., Knapp, D., Marquart, G., and Clauser, C.: Report on the numerical reservoir model used for the simulation of the Acoculco reservoir in Mexico, Deliverable 6.2, WP6, GEMex H2020 project, European Comission, 1.2, available at: (26 May 2020), 2019b. 

Deb, P., Salimzadeh, S., Dübner, S., and Clauser, C.: Laboratory experiments and numerical simulations of hydraulic fracturing for enhanced geothermal systems, in: Proceedings of the European Geothermal Congress 2019, Den Haag, the Netherlands, 11–14 June 2019, 4 pp., 2019c. 

DIN 18141-1: 2014-05: Baugrund-Untersuchung von Gesteinsproben – Teil 1: Bestimmung der einaxialen Druckfestigkeit, Beuth,, 2014. 

Ding, L. Y., Mehra, R. K., and Donnelly, J. K.: Stochastic Modeling in Reservoir Simulation, SPE Reservoir Engineering, 7, 98–106,, 1992. 

DOE Data Explorer: U.S. Department of Energy Office of Scientific and Technical Information, DOE Data Explorer Utah Forge, available at: (last access: 21 October 2020), 2018. 

Durán, E. L., Adam, L., Wallis, I. C., and Barnhoorn, A.: Mineral Alteration and Fracture Influence on the Elastic Properties of Volcaniclastic Rocks, J. Geophys. Res.-Sol. Ea., 124, 4576–4600,, 2019. 

Ebigbo, A., Niederau, J., Marquart, G., Dini, I., Thorwart, M., Rabbel, W., Peching, R., Bertani, R., and Clauser, C.: Influence of depth, temperature, and structure of a crustal heat source on the geothermal reservoirs of Tuscany: numerical modelling and sensitivity study, Geotherm. Energy, 5, 1–29,, 2016. 

Ebong, E. D., Akpan, A. E., and Ekwok, S. E.: Stochastic modelling of spatial variability of petrophysical properties in parts of the Niger Delta Basin, southern Nigeria, Journal of Petroleum Exploration and Production Technology, 10, 569–585,, 2019. 

Eggertsson, G. H., Lavallée, Y., Kendrick, J. E., and Markússon, S. H.: Improving fluid flow in geothermal reservoirs by thermal and mechanical stimulation: The case of Krafla volcano, Iceland, J. Volcanol. Geoth. Res., 391, 106351,, 2020. 

Enge, H. D., Buckley, S. J., Rotevatn, A., and Howell, J. A.: From outcrop to reservoir simulation model: workflow and procedures, Geosphere, 3, 469–490,, 2007. 

Farmer, C. L.: Upscaling: a review, Int. J. Numer. Meth. Fl., 40, 63–78,, 2002. 

Ferrari, L., Orozco-Esquivel, T., Manea, V., and Manea, M.: The dynamic history of the Trans-Mexican Volcanic Belt and the Mexico subduction zone, Tectonophysics, 522/523, 122–149,, 2012. 

Ferriz, H. and Mahood, G.: Eruptive rates and compositional trends at Los Humeros volcanic center, Puebla, Mexico, J. Geophys. Res., 89, 8511–8524,, 1984. 

Filomena, C. M., Hornung, J., and Stollhofen, H.: Assessing accuracy of gas-driven permeability measurements: a comparative study of diverse Hassler-cell and probe permeameter devices, Solid Earth, 5, 1–11,, 2014. 

Fitz-Díaz, E., Lawton, T. F., Juárez-Arriaga, E., and Chávez-Cabello, G.: The Cretaceous-Paleogene Mexican orogen: structure, basin development, magmatism and tectonics, Earth-Sci. Rev., 183, 56–84,, 2017. 

Flovenz, O. G., Spangenberg, E., Kuhlenkampff, J., Arnason, K., Karlsdottir, R., and Huenges, E.: The role of electrical interface conduction in geothermal exploration, in: Proceedings of the World Geothermal Congress 2005, Antalya, Turkey, 24–29 April 2005, 9 pp., 2005. 

Forchheimer, P.: Wasserbewegung durch Boden, Z. Ver. Dtsch. Ing., 45, 1782–1788, 1901. 

Franco, A. and Donatini, F.: Methods for the estimation of the energy stored in geothermal reservoirs, J. Phys., 796, 012025,, 2017. 

Fuentes-Guzmán, E., González-Partida, E., Camprubí, A., Hernández-Avilés, G., Gabites, J., Ruggieri, G., Iriondo, A., and López-Martínez, M.: The Miocene Tatatila-Las Minas IOCG skarn deposits (Veracruz) as a result of adakitic magmatism in the Trans-Mexican Volcanic Belt Short running title: Miocene IOCG deposits and adakitic magmas in the Trans-Mexican Volcanic Belt, Boletín de la Sociedad Geológica Mexicana, 73, A1105020,, 2020. 

Gan, Q. and Elsworth, D.: Production optimization in fractured geothermal reservoirs by coupled discrete fracture network modelling, Geothermics, 62, 131–142,, 2016. 

García-Gutiérrez, A. and Contreras, E.: Measurement of Thermal Conductivity and Diffusivity of Drill Core Samples from the Los Humeros Geothermal Field, Mexico, by a Line-Source Technique, GRC Transactions, 31, 555–559, 2007. 

García-Palomo, A., Macías, J. L., Tolson, G., Valdez, G., and Mora, J. C.: Volcanic stratigraphy and geological evolution of the Apan región, east-central sector of the Trans-Mexica Volcanic Belt, Geofís. Int., 41, 133–150, 2002. 

García-Palomo, A., Macías, J. L., Jiménez, A., Tolson, G., Mena, M., Sánchez-Núñez, J. M., Arce, J. L, Layer, P. W., Santoyo, M. A., and Lermo-Samaniego, J.: NW-SE Pliocene-Quaternary extension in the Apan-Acoculco region, eastern Trans-Mexican Volcanic Belt, J. Volcanol. Geoth. Res., 349, 240–255,, 2018. 

Gard, M., Hasterok, D., and Halpin, J. A.: Global whole-rock geochemical database compilation, Earth Syst. Sci. Data, 11, 1553–1566,, 2019. 

Georoc Mainz: GEOROC database, Geochemistry of Rocks of the Oceans and Continents, available at: (last access: 21 October 2020), 2020. 

GeoScout database: geoLogicSystems, geoScout database, available at: (last access: 11 November 2020), 2020. 

Geotek: Multi-Sensor Core Logger, Manual, Geotek, Nene House, Drayton Fields, Daventry, Northants, UK, 2000. 

Geotron-Elektronik: LightHouse UMPC V1.02, Installations- und Bedienungshandbuch, 1.6, Pirna: Geotron-Elektronik, 2011. 

Ghassemi, A: Application of rock failure simulation in design optimization of the hydraulic fracturing, in: Porous Rock Fracture Mechanics with Application to Hydraulic Fracturing, Drilling and Structural Engineering, edited by: Shojaei, A. K. and Shao, J., Elsevier, Duxford, UK, 3–23,, 2017. 

Gómez-Tuena, A. and Carrasco-Núñez, G.: Cerro Grande volcano: the evolution of a Miocene stratocone in the early Trans-Mexican Volcanic Belt. Tectonophysics, 318, 249–280,, 2000. 

Gu, Y., Rühaak, W., Bär, K., and Sass, I.: Using seismic data to estimate the spatial distribution of rock thermal conductivity at reservoir scale, Geothermics, 66, 61–72,, 2017. 

Guo, H., Aziz, N. I., and Schmidt, L. C.: Rock fracture-toughness determination by the Brazilian test, Eng. Geol., 33, 177–188,, 1993. 

Hartanato, L.: Different Scales and Integration of data in Reservoir Simulation, PhD thesis, Curtin University of Technology, Perth, Australia, 135 pp., 2004. 

Hartmann, A., Rath, V., and Clauser, C.: Thermal conductivity from core and well log data, Int. J. Rock Mech. Min. Sci., 42, 1042–1055,, 2005. 

Hartmann, J. and Moosdorf, N.: The new global lithological map database GLiM: A representation of rock properties at the Earth surface, Geochem. Geophys. Geosyst., 13, Q12004,, 2012. 

Hassanzadegan, A., Blöcher, G., Milsch, H., Urpi, L., and Zimmermann, G.: The Effects of Temperature and Pressure on the Porosity Evolution of Flechtinger Sandstone, Rock Mech. Rock Eng., 47, 421–434,, 2013. 

Hatakeda, K., Lin, W., Hirose, T., Tanikawa, W., Hamada, Y., and Tadai, O.: Electrical resistivity measurements of rocks under confining pressure condition, JAMSTEC Report of Research and Development, 24, 1–9,, 2017. 

Heap, M. J. and Kennedy, B. M.: Exploring the scale-dependent permeability of fractured andesite, Earth Planet. Sci. Lett., 447, 139–150,, 2016. 

Heap, M. J., Lavallée, Y., Petrakova, L., Baud, P., Reuschlé, T., Varley, N. R., and Dingwell, D. B.: Microstructural controls on the physical and mechanical properties of edifice-forming andesites at Volcán de Colima, Mexico, J. Geophys. Res.-Sol. Ea., 119, 2925–2963,, 2014. 

Heap, M. J., Gravley, D. M., Kennedy, B. M., Gilg, H. A., Bertolett, E., and Barker, S. L.: Quantifying the role of hydrothermal alteration in creating geothermal and epithermal mineral resources: The Ohakuri ignimbrite (Taupo Volcanic Zone, New Zealand), J. Volcanol. Geoth. Res., 390, 106703,, 2020. 

Horai, K.-I. and Susaki, J.-I.: The effect of pressure on the thermal conductivity of silicate rocks up to 12 kbar, Phys. Earth. Planet. Inter., 55, 292–305,, 1989. 

Hornung, J. and Aigner, T.: Sedimentäre Architektur und Poroperm-Analyse fluviatiler Sandsteine: Fallbeispiel Coburger Sandstein, Franken, Hallesches Jahrb. Geowiss., 18, 121–138, 2004. 

Howell, J. A., Allard, W. M., and Good, T. R.: The application of outcrop analogues in geological modeling: a review, present status and future outlook, Geol. Soc. Lond. Spec. Publ., 387, 1–25,, 2014. 

IHS Markit: The AccuMap database, available at: (last access: 11 November 2020), 2020. 

ISRM: Suggested methods for determining the uniaxial compressive strength and deformability of rock materials, in: The Complete ISRM Suggested Methods for Rock Characterization, Testing and Monitoring: 1974–2006, edited by: Ulusay, R. and Hudson, J. A., Ankara, Turkey, Pergamon Press, 137–138, 1979. 

ISRM: Suggested Methods – Rock characterization testing and monitoring, edited by: Brown, E. T., Pergamon Press, Oxford, UK, 211 pp., 1981. 

ISRM: Suggested methods for determining point load test, RTH 325-89, in: The Complete ISRM Suggested Methods for Rock Characterization, Testing and Monitoring: 1974–2006, edited by: Ulusay, R. and Hudson, J. A., Pergamon Press, Ankara, Turkey, 53–60, 1984. 

ISRM: Suggested methods for determining the fracture toughness of rock, Int. J. Rock Mech. Min. Sci. Geomech., 25, 71–96,, 1988. 

Jaritz, R.: Quantifizierung der Heterogenität einer Sandsteinmatrix am Beispiel des Stubensandstein (Mittlerer Keuper, Württemberg), Tübinger Geol. Abhandlungen, 48, 104 pp., 1999. 

Jentsch, A., Jolie, E., Jones, D. G., Taylor-Curran, H., Peiffer, L., Zimmer, M., and Lister, B.: Magmatic volatiles to assess permeable volcano-tectonic structures in the Los Humeros geothermal field, Mexico, J. Volcanol. Geoth. Res., 394, 106820,, 2020. 

Jolie, E., Bruhn, D., López Hernández, A., Liotta, D., Garduño-Monroy, V. H., Lelli, M., Páll Hersir, G., Arango-Galván, C., Bonté, D., Calcagno, P., Deb, P., Clauser, C., Peters, E., Hernández Ochoa, A. F., Huenges, E., González Acevedo, Z. I., Kieling, K., Trumpy, E., Vargas, J., Gutiérrez-Negrín, L. C., Aragón-Aguilar, A., Halldórsdóttir, S., González Partida, E., van Wees, J.-D., Ramírez Montes, M. A., Diez León, H. D., and the GEMex team: GEMex – A Mexican-European Research Cooperation on Development of Superhot and Engineered Geothermal Systems, in: Proceedings of the 43rd Workshop on Geothermal Reservoir Engineering, Stanford, CA, USA, 12–14 February 2018, 10 pp., 2018. 

Khajeh, M. M.: Heterogeneity Consideration and Upscaling of Elastic Properties in Coupled Geomechanical Flow Simulation of SAGD, PhD thesis, University of Alberta, Edmonton, Canada, 131 pp., 2013. 

Klinkenberg, L. J.: The permeability of porous media to liquids and gas, Drilling and Production Practice, 200–213, API-41-200, 1941. 

Kozdrój, W., Nawrocki, J., Pańczyk-Nawrocka, M., ZiółkowskaKozdrój, M., Wójcik, K., Kumek, J., and González-Partida, E.: Stratigraphic, petrological features and datings of Los Humeros rocks from outcrops and boreholes, in: Final report on active systems: Los Humeros and Acoculco, Deliverable 4.1, WP4, GEMex H2020 project, European Comission, available at: (26 May 2020), 2019. 

Kruszewski, M., Hofmann, H., Gomez-Alvarez, F., Bianco, C., Jimenez-Haro, A., Garduno, V. H., Liotta, D., Trumpy, E., Brogi, A., Wheeler, W., Bastesen, E., Parisio, F., and Saenger, E. H.: Integrated Stress Field Estimation and Implications for Enhanced Geothermal System Development in Acoculco, Mexico, Geothermics, 89, 101931,, 2021. 

Kummerow, J. and Raab, S.: Temperature dependence of electrical resistivity – Part II: A new experimental set-up to study fluid-saturated rocks, Energy Proced., 76, 247–255,, 2015. 

Kummerow, J., Raab, S., and Spangenberg, E.: The impact of reactive flow on electrical and hydraulic rock properties in supercritical geothermal settings, GEMex Final Conference, Potsdam, Germany, 18–19 February 2020, 32, 2020. 

Kushnir, A. R. L., Heap, M. J., and Baud, P.: Assessing the role of fractures on the permeability of the Permo-Triassic sandstones at the Soultz-sous-Forêts (France) geothermal site, Geothermics, 74, 181–189,, 2018. 

Lacinska, A. M., Rochelle, C., Kilpatrick, A., Rushton, J., Weydt, L. M., Bär, K., and Sass, I.: Evidence for fracture-hosted fluid-rock reactions within geothermal reservoirs of the eastern Trans-Mexican Volcanic Belt, GEMex Final Conference, Potsdam, Germany, 18–19 February 2020, 33, 2020. 

Le Maitre, R. W. and Streckeisen, A.: Igneous rocks: a classification and glossary of terms – recommendations of the International Union of Geological Sciences Subcommission on the Systematics of Igneous Rocks, Cambridge University Press, Cambridge, UK,, 2003. 

Lenhardt, N. and Götz, A. E.: Volcanic settings and their reservoir potential: an outcrop analog study on the Miocene Tepoztlán Formation, Central Mexico, J. Volcanol. Geoth. Res., 204, 66–75,, 2011. 

Lepillier, B., Daniilidis, A., Doonechaly Gholizadeh, N., Bruna, P.-O., Kummerow, J., and Bruhn, D.: A fracture flow permeability and stress dependency simulation applied to multi-reservoirs, multi-production scenarios analysis, Geotherm. Energy, 24, 1–16,, 2019. 

Lepique, M.: Empfehlung Nr. 10 des Arbeitskreises 3.3 “Versuchstechnik Fels” der Deutschen Gesellschaft für Geotechnik e. V.: Indirekter Zugversuch an Gesteinsproben – Spaltzugversuch, Bautechnik, 85, 623–627,, 2008. 

Lévy, L., Gibert, B., Sigmundsson, F., Flóvenz, Ó. G., Hersir, G. P., Briole, P., Pezard, P. A.: The role of smectites in the electrical conductivity of active hydrothermal systems: electrical properties of core samples from Krafla volcano, Iceland, Geophys. J. Int., 215, 1558–1582,, 2018. 

Lévy, L., Maurya, P. K., Byrdina, S., Vandemeulebrouck, J., Sigmundsson, F., Árnson, K., Ricci, T., Deldicque, D., Roger, M., Gibert, B., and Labazuy, P.: Electrical resistivity tomography and time-domain induced polarization field investigations of geothermal areas at Krafla, Iceland: comparison to borehole and laboratory frequency-domain electrical observations, Geophys. J. Int., 218, 1469–1489,, 2019. 

Li, C., Chen, X., and Du, Z.: A New Relationship of Rock Compressibility with Porosity, SPE 88464-MS, in: SPE Asia Pacific Oil and Gas Conference and Exhibition, Perth, Australia, 18–20 October 2004, 2004. 

Linsel, A., Wiesler, S., Hornung, J., and Hinderer, M.: High-resolution analysis of the physicochemical characteristics of sandstone media at the lithofacies scale, Solid Earth, 11, 1511–1526,, 2020. 

Liotta, D., Brogi, A., Garduño-Monroy, V. H., Gomez, F., Wheeler, W. H., Bastesen, E., Torabi, A., Bianco, C., Jimenez-Haro, A., Olvera-Garcia, E., and Zucchi, E.: Regional Geological Structures, in: Final report on active systems: Los Humeros and Acoculco, Deliverable 4.1, WP4, GEMex H2020 project, European Comission, available at: (last access: 26 May 2020), 2019. 

Lippmann, E. and Rauen, A.: Measurements of Thermal Conductivity (TC) and Thermal Diffusivity (TD) by the Optical Scanning Technology, Lippmann and Rauen GbR, Schaufling, Germany, 2009. 

López-Hernández, A.: Estudio Regional Volcánico y Estructural del Campo Geotérmico de Los Humeros, Puebla., México, Geotermia Revista Mexicana de Geoenergía, 11, 17–36, 1995. 

López-Hernández, A., García-Estrada, G., Aguirre-Díaz, G., González-Partida, E., Palma-Guzmán, H., and Quijano-Léon, J.: Hydrothermal activity in the Tulancingo-Acoculco Caldera Complex, central Mexico: Exploratory studies, Geothermics, 38, 279–293,, 2009. 

Lorenzo-Pulido, C., Armenta-Flores, M., and Ramírez-Silva, G.: Characterization of the Acoculco Geothermal Zone as a HDR System, GRC Transactions, 34, 369–372, 2010. 

Macías, J. L., Arce, J. L., García-Tenorio, F., Layer, P. W., Rueda, H., Reyes-Agustin, G., and Avellán, D.: Geology and geochronology of Tlaloc, Telapón, Iztaccíhuatl, and Popocatépetl volcanoes, Sierra Nevada, central México, Field Guides 25, 163–193,, 2012. 

Mandrone, G., Comina, C., and Vacha, D.: Faults characterization aimed at geothermal fluid path identification and quantification, GEMex Final Conference, Potsdam, Germany, 18–19 February 2020, 41, 2020. 

Merriman, D., Hofmeister, A. M., Roy, D. J., and Whittington, A. G.: Temperature-dependent thermal transport properties of carbonate minerals and rocks, Geosphere, 14, 1961–1987,, 2018. 

Micromeritics: AccuPyc 1330 Pycnometer, V2.02, Part No. 133-42808-01, Micromeritics GmbH, Munich, Germany, 1997. 

Micromeritics: GeoPyc 1360, V3., Part 136-42801-01, Micromeritics GmbH, Munich, Germany, 1998. 

Micromeritics: AccuPyc II 1340, Product Broschure, Micromeritics Instrument Corporation, Norcross, available at: (last access: 9 December 2020), 2014. 

Mielke, P., Nehler, M., Bignall, G., and Sass, I.: Thermo-physical rock properties and the impact of advancing hydrothermal alteration – a case study from the Tauhara geothermal field, New Zealand, J. Volcanol. Geoth. Res., 301, 14–28,, 2015. 

Mielke, P., Weinert, S., Bignall, G., and Sass, I.: Thermo-physical rock properties of greywacke basement rock and intrusive lavas from the Taupo Volcanic Zone, New Zealand, J. Volcanol. Geoth. Res., 324, 179–189,, 2016. 

Mobarak, S. A. and Somerton, W. H.: The Effect Of Temperature And Pressure On Wave Velocities In Porous Rocks, in: Proceedings of the Fall Meeting of the Society of Petroleum Engineers of AIME, New Orleans, Louisiana, 18 pp.,, 1971. 

Moosavi, S. A., Goshtasbi, K., Kazemzadeh, E., Bakhtiari, H. A., Esfahani, M. R., and Vali, J.: Relationship between porosity and permeability with stress using pore volume compressibility characteristic of reservoir rocks, Arab. J. Geosci., 7, 231–239,, 2014. 

Mordensky, S. P., Heap, M. J., Kennedy, B. M., Gilg, H. A., Villeneuve, M. C., Farquharson, J. I., and Gravley, D. M.: Influence of alteration on the mechanical behaviour and failure mode of andesite: implications for shallow seismicity and volcano monitoring, B. Volcanol., 44, 1–12,, 2019. 

Mutschler, T.: Neufassung der Empfehlung Nr. 1 des Arbeitskreises “Versuchstechnik Fels” der Deutschen Gesellschaft für Geotechnik e. V.: Einaxiale Druckversuche an zylindrischen Gesteinsprüfkörpern, Bautechnik, 81, 825–834,, 2004. 

National Geochemical Database USGS: Geochemical database of the US Geological Survey, USA, available at: (last access: 21 October 2020), 2014. 

National Geothermal data system NGDS: Database of various reports and data sets related to geothermal and geological projects in the USA, available at: (last access: 11 November 2020), 2014. 

NAVDAT data base: The North American Volcanic and Intrusive Rock Database (NAVDAT), available in the EarthChem Portal, available at: (last access: 21 October 2020), 2020. 

Navelot, V., Géraud, Y., Favier, A., Diraison, M., Corsini, M., Lardeaux, J.-M., Verati, C., de Lépinary, J. M., Legendre, L., and Beauchamps, G.: Petrophysical properties of volcanic rocks and impacts of hydrothermal alteration in the Guadeloupe Archipelago (West Indies), J. Volcanol. Geoth. Res., 360, 1–21,, 2018. 

Nono, F., Gibert, B., Parat, F., Loggia, D., Cichy, S. B., and Violay, M.: Electrical conductivity of Icelandic deep geothermal reservoirs up to supercritical conditions: Insight from laboratory experiments, JVGR, 391, 106364,, 2020. 

Norden, B., Förster, A., Förster, H.-J., and Fuchs, S.: Temperature and pressure corrections applied to rock thermal conductivity: impact on subsurface temperature prognosis and heat-flow determination in geothermal exploration, Geotherm. Energy, 1, 1–19,, 2020. 

Norini, G., Gropelli, G., Sulpizio, R., Carrasco-Núñez, G., Dávila-Harris, P., Pellicioli, C., Zucca, F., and De Franco, R.: Structural analysis and thermal remote sensing of the Los Humeros Volcanic Complex: Implications for volcano structure and geothermal exploration, J. Volcanol. Geoth. Res., 301, 221–237,, 2015. 

Norini, G., Carrasco-Núñez, G., Corbo, F., Lermo, J., Hernández, J., Castro, C., Bonini, M., Montanari, D., Corti, G., Moratti, G., Piccardi, L., Chavez, G., Zuluaga, M. C., Ramírez, M., and Cedillo, F.: The structural architecture of the Los Humeros volcanic complex and geothermal field, J. Volcanol. Geoth. Res., 381, 312–329,, 2019. 

Nurhandoko, B. E. B., Kurniadi, R., Susilowati, Tryoso, K., Widowati, S., Hadi, M. R. A., Abda, M. R., Martha, R. K., Fatiah, E., and Komara, R.: Integrated Subsurface Temperature Modeling beneath Mt. Lawu and Mt. Muriah in The Northeast Java Basin, Indonesia, Open Geosci., 11, 341–351,, 2019. 

Ohnaka, M.: Stability of Remanent Magnetization of Rocks under Compression – Its Relation to the Grain Size of Rock-forming Ferromagnetic Minerals, J. Geomagn. Geoelectr., 21, 495–505,, 1969. 

Ólavsdóttir, J., Andersen, M. S., and Boldreel, L. O.: Reservoir quality of intrabasalt volcaniclastic units onshore Faroe Islands, North Atlantic Igneous Province, northeast Atlantic, AAPG Bull., 99, 467–497,, 2015. 

Oliva-Urcia, B., Kontny, A., Vahle, C., and Schleicher, A. M.: Modification of the magnetic mineralogy in basalts due to fluid-rock interactions in a high-temperature geothermal system (Krafla, Iceland), Geophys. J. Int., 186, 155–174,, 2011. 

Pérez-Campos, X., Kim, Y., Husker, A., Davis, P. M., Clayton, R. W., Iglesias, A., Pacheco, J., Singh, S., Constantin, M. V., and Gurnis, M.: Horizontal subduction and truncation of the Cocos plate beneath central Mexico, Geophys. Res. Lett. 35, 1–6,, 2008. 

Petlab database: GNS Science, Petlab database, available at: (last access: 21 October 2020), 2020. 

Pola, A., Crosta, G., Fusi, N., Barberini, V., and Norini, G.: Influence of alteration on physical properties of volcanic rocks, Tectonophysics, 566/567, 67–86,, 2012. 

Pola, A., Crosta, G. B., Fusi, N., and Castellanza, R.: General characterization of the mechanical behaviour of different volcanic rocks with respect to alteration, Eng. Geol., 169, 1–13,, 2014. 

Popov, Y., Beardsmore, G., Clauser, C., and Roy, S.: ISRM Suggested Methods for Determining Thermal Properties of Rocks from Laboratory Tests at Atmospheric Pressure, Rock Mech. Rock Eng., 49, 4179–4207,, 2016. 

Popov, Y. A., Sass, P. D., Williams, C. F., and Burkhardt, H.: Characterization of rock thermal conductivity by high resolution optical scanning, Geothermics, 28, 253–276,, 1999. 

Qi, D. and Hesketh, T.: An Analysis of Upscaling Techniques for Reservoir Simulation, Pet. Sci. Technol., 23, 827–842,, 2005. 

Qu, H., Yang, B., Tian, X., Liu, X., Yang, H., Dong, W., and Chen, Y.: The primary controlling parameters of porosity, permeability, and seepage capability of tight gas reservoirs: a case study on Upper Paleozoic Formation in the eastern Ordos Basin, Northern China, Pet. Sci., 16, 1270–1284,, 2019. 

Ringrose, P. and Bentley, M.: Reservoir Model Design, Springer, Dordrecht, the Netherlands,, 2015. 

Rock Properties Database British Columbia Canada: Rock Properties Database, a collaborative Geoscience BC project among Mira Geoscience, the Geological Survey of Canada and CAMIRO, Open file resource, available at: (last access: 21 October 2020), 2018. 

RockViewer: Software for Organization of Petrographic Photomicrographs, available at: (last access: 21 October 2020), 2020. 

Romo-Jones, J. M., Gutiérrez-Negrín, L. C., and Canchaela-Félix, I.: 2018 México Country Report, IEA Geothermal, available at: (last access: 28 May 2020), 2019. 

Ruggeri, G., Morelli, G., Zucchi, M., Braschi, E., Agostini, S., Ventruti, G., Brogi, A., Liotta, D., Boschi, C., and Gonzalez-Partida, E.: Insight into the fluids occurring in the super-hot reservoir of the Los Humeros geothermal system from fluid inclusions and isotopic data of the Las Minas exhumed system (Mexico), GEMex Final Conference, Potsdam, Germany, 18–19 February 2020, 53, 2020. 

Rühaak, W., Bär, K., and Sass, I.: Combining numerical modeling with geostatistical analysis for an improved reservoir exploration, Energy Proced., 59, 315–322,, 2014. 

Rühaak, W., Guadagnini, A., Geiger, S., Bär, K., Gu, Y., Aretz, A., Hohmuth, S., and Sass, I.: Upscaling thermal conductivities of sedimentary formations for geothermal exploration, Geothermics, 58, 49–61,, 2015. 

Rybacki, E., Meier, T., and Dresen, G.: What controls the mechanical properties of shale rocks? – Part II: Brittleness, J. Petrol. Sci. Eng., 144, 39–58,, 2016. 

Saller, A. H. and Henderson, N.: Distribution of Porosity and Permeability in Platform Dolomites: Insight from the Permian of West Texas, AAPG Bull., 82, 1528–1550,, 1998. 

Sánchez-Vila, X., Guadagnini, A., and Carrera, J.: Representative hydraulic conductivities in saturated groundwater flow, Rev. Geophys., 44, RG3002,, 2006. 

Sass, I. and Götz, A. E.: Geothermal reservoir characterization: a thermofacies concept, Terra Nova, 24, 142–147,, 2012. 

Sass, J. H., Lachenbruch, A. H., Munroe, R. J., Greene, G. W., and Moses Jr., T. H.: Heat flow in the Western United States, J. Geophys. Res., 76, 6376–6413,, 1971. 

Scheibe, T. and Yabusaki, S.: Scaling of flow and transport behavior in heterogeneous groundwater systems, Adv. Water Resour., 22, 223–238,, 1998. 

Scheidegger, A. E.: The Physics of Flow Through Porous Media, 3rd ed., 372 pp., University of Toronto Press, Toronto, Canada, (last access: 25 May 2020), 1974. 

Schön, J. H.: Physical properties of rocks: Fundamentals and principles of petrophysics, Developments in petroleum science, Elsevier, Amsterdam, the Netherlands, 512 pp., 2015. 

Sciencebase Minnesota: US Geological Survey (USGS) – Collection of rock properties database from Minnesota USA, available at: (last access: 21 October 2020), 2010. 

Scott, S. W., Covell, C., Júlíusson, E., Valfells, Á., Newson, J., Hrafnkelsson, B., Pállson, H., and Gudjónsdóttir, M.: A probalistic geological model of the Krafla geothermal system constrained by gravimetric data, Geotherm. Energy, 29, 1–30,, 2019. 

Setaram Instrumentation: K/C80-1A C80 Commissioning, Setaram Instrumentation KEP Technologies, Caluire, France, 52 pp., 2009. 

SGM: CARTA GEOLÓGICO – MINERA, E14-2, Servicio Geológico Mexicano, Ciudad de Mexico, Mexico, first edition, 2002a. 

SGM: CARTA GEOLÓGICO – MINERA, E14-3, Servicio Geológico Mexicano, Veracruz, Mexico, first edition, 2002b. 

Shankland, T. J., Duba, A. G., Mathez, E. A., and Peach, C. L.: Increase of electrical conductivity with pressure as an indicator of conduction through a solid phase in midcrustal rocks, J. Geophys. Res., 102, 14741–14750,, 1997. 

Siratovich, P., Heap, M. J., Villeneuve, M., Cole, J., and Reuschlé, T.: Physical property relationships of the Rotokawa Andesite, a significant geothermal reservoir rock in the Taupo Volcanic Zone, New Zealand, Geotherm. Energy, 10, 1–31,, 2014. 

Siratovich, P. A., Sass, I., Homuth, S., and Bjornsson, A.: Thermal Stimultaion of Geothermal Reservoirs and Laboratory Investigation of Thermally Induced Fractures, GRC Transactions, 35, 1529–1536, 2011. 

Somerton, W. H.: Thermal properties and temperature-related behavior of rock-fluid systems, Dev. Pet. Sci., 37, Elsevier, Amsterdam, the Netherlands, 257 pp., 1992. 

Sosa-Ceballos, G., Macías, J. L., Avellán, D. R., Salazar-Hermenegildo, N., Boijseauneau-López, M. E., and Pérez-Orozco, J. D.: The Acoculco Caldera Complex magmas: genesis, evolution and relation with the Acoculco geothermal system, J. Volcanol. Geoth. Res., 358, 288–306,, 2018. 

Stimac, J. A., Powell, T. S., and Golla, G. U.: Porosity and permeability of the Tiwi geothermal field, Philippines, based on continuous and spot core measurements, Geothermics, 33, 87–107,, 2004. 

Tanikawa, W. and Shimamoto, T.: Comparison of Klinkenberg corrected gas permeability and water permeability in sedimentary rocks, Int. J. Rock Mech. Min. Sci., 46, 229–238,, 2008. 

Toledo, T., Gaucher, E., Jousset, P., Maurer, H., Krawczyk, C., Calò, M., and Figueroa, Á.: Local earthquake tomography at the Los Humeros geothermal field, GEMex Final Conference, Potsdam, Germany, 18–19 February 2020, 6.3, 2020. 

Urbani, S., Giordano, G., Lucci, F., Rossetti, F., Acocella, V., and Carrasco-Núñez, G.: Estimating the depth and evolution of intrusions at resurgent calderas: Los Humeros (Mexico), Solid Earth, 11, 527–545,, 2020. 

Vagnon, F., Colombero, C., Colombo, F., Comina, C., Ferrero, A. M., Mandrone, G., and Vinciguerra, S. C.: Effects of thermal treatment on physical and mechanical properties of Valdieri Marble, NW Italy, Int. J. Rock Mech. Min. Sci., 116, 75–86,, 2019. 

Vagnon, F., Colombero, C., Comina, C., Ferrero, A. M., Mandrone, G., Missagia, R., and Vinciguerra, S. C.: Relating physical properties to temperature induced damage in carbonate rocks, Geotechnique Letters, accepted, 2021. 

Vinciguerra, S., Trovato, C., Meredith, P. G., and Benson, P. M.: Relating seismic velocities, thermal cracking and permeability in Mt. Etna and Iceland basalts, Int. J. Rock Mech. Min. Sci., 42, 900–910,, 2005. 

Vosteen, H.-D. and Schellschmidt, R.: Influence of temperature on thermal conductivity, thermal capacity and thermal diffusivity for different types of rock, Phys. Chem. Earth, 28, 499–509,, 2003. 

Walia, S. and Leahy, G.: Addressing challenges in petrophysical modelling, HARTENERGY, available at: (last access: 6 November 2020), 2014. 

Weides, S. and Majorowicz, J.: Implications of Spatial Vari-ability in Heat Flow for Geothermal Resource Evalua-tion in Large Foreland Basins: The Case of the West-ern Canada Sedimentary Basin, Energies, 7, 2573–2594,, 2014. 

Weides, S., Moeck, I., Majorowicz, J., Palombi, D., and Grobe, M.: Geothermal exploration of Paleozoic forma-tions in Central Alberta, Can. J. Earth Sci., 50, 519–534,, 2013. 

Weinert, S., Bär, K., and Sass, I.: Database of Petrophysical Properties of the Mid-German Crystalline High, Earth Syst. Sci. Data Discuss. [preprint],, in review, 2021. 

Wen, X.-H. and Gomez-Hernandez, J.: Upscaling hydraulic conductivities in heterogeneous media: an overview, J. Hydrol., 183, ix–xxxii,, 1996. 

Weydt, L. M., Heldmann, C.-D. J., Machel, H. G., and Sass, I.: From oil field to geothermal reservoir: assessment for geothermal utilization of two regionally extensive Devonian carbonate aquifers in Alberta, Canada, Solid Earth, 9, 953–983,, 2018a. 

Weydt, L. M., Bär, K., Colombero, C., Comina, C., Deb, P., Lepillier, B., Mandrone, G., Milsch, H., Rochelle, C. A., Vagnon, F., and Sass, I.: Outcrop analogue study to determine reservoir properties of the Los Humeros and Acoculco geothermal fields, Mexico, Adv. Geosci., 45, 281–287,, 2018b. 

Weydt, L. M., Ramírez-Guzmán, Á. A., Pola, A., Lepillier, B., Kummerow, J., Mandrone, G., Comina, C., Deb, P., Norini, G., Gonzalez-Partida, E., Avellán, D. R., Macías, J. L., Bär, K., and Sass, I.: Petrophysical and mechanical rock property database of the Los Humeros and Acoculco geothermal fields (Mexico), TU Darmstadt datalib,, 2020. 

Weydt, L. M., Lucci, F., Lacinska, A., Scheuvens, D., Carrasco-Núñez, G., Giordano, G., Rochelle, Schmidt, S., Bär, K., and Sass, I.: Modeling of volcanic systems – part I: Petrophysical characterization and the affect of hydrothermal alteration - the case of the Los Humeros geothermal field (Mexico), in preparation, 2021. 

Whittington, A. G., Hofmeister, A. M., and Nabelek, P. I.: Temperature-dependent thermal diffusivity of the Earth's crust and implications for magmatism, Nature, 458, 319–321,, 2009. 

Yáñez, C. and García, S.: Exploración de la región geotérmica Los Humeros-Las Derrumbadas, Estados de Puebla y Veracruz (Internal report), Comisión Federal de Electricidad, Mexico City, Mexico, 96 pp., 1982.  

You, B., Xu, J., Shi, S., Liu, H., and Li, H.: Effect of Stress and Water Pressure on Permeability of Fractured Sandstone Based on Response Surface Method, Front. Earth. Sci., 11, 1–8,, 2020. 

Zahner Scientific Instruments: Zahner-Zennium Electrochemical Workstation, User's manual, 212 pp., Zahner Scientific Instruments, Kronach, Germany, 2008. 

Zhao, X. G., Wang, J., Chen, F., Li, P. F., Ma, L. K., Xie, J. L., and Liu, Y.: Experimental investigation on the thermal conductivity characteristics of Beishan granitic rocks for China's HLW disposal, Tectonophysics, 683, 124–137,, 2016. 

Zhang, H., Sun, Q., and Ge, Z.: Analysis of the characteristics of magnetic properties change in the rock failure process, Acta Geophys., 68 289–302,, 2020. 

Zheng, J., Zheng, L., Liu, H.-H., and Ju, Y.: Relationships between permeability, porosity and effective stress for low-permeability sedimentary rock, Int. J. Rock Mech. Min., 78, 304–318,, 2015. 

ZH Instruments: ZH Instruments Magnetic suceptibility meter SM30, User's manual, Brno, Czech Republic, 2008. 

Zimmermann, R. W., Somerton, W. H., and King, M. S.: Compressibility of porous rocks, J. Geophys. Res., 91, 765–777,, 1986. 

Zoback, M. D.: Reservoir geomechanics, Cambridge University Press, Cambridge, UK, 2011. 

Zoth, G. and Hänel, R.: Thermal conductivity, in: Handbook of terrestrial heat flow density determination, edited by: Hanel, R., Rybach, L., and Stegena, L., Kluwer, Dordrecht, the Netherlands, 449–453, 1988. 

Short summary
Petrophysical and mechanical rock properties are essential for reservoir characterization of the deep subsurface and are commonly used for the population of numerical models or the interpretation of geophysical data. The database presented here aims at providing easily accessible information on rock properties and chemical analyses complemented by extensive metadata (location, stratigraphy, petrography) covering volcanic, sedimentary, metamorphic and igneous rocks from Jurassic to Holocene age.