Reactivity of dolomitizing fluids and Mg source evaluation of fault-controlled dolomitization at the Benicàssim outcrop analogue (Maestrat Basin, E Spain)

Peer reviewedPostprin

The evolution of a gypsum vein network in compressive tectonic settings: the Montsant anticlinorium case (NE Spain)

The Montsant anticlinorium is part of the Pàndols-Cavalls-Montsant tectonic line situated in the southwestern area of the Catalan Coastal Ranges, adjacent to the contact with the Tertiary Ebro Basin. We have interpreted this Alpine structure as a triangular type I zone with two opposite faults. The centre of the anticlinorium is formed by middle Muschelkalk facies with intensively deformed gypsum layers and an intensively deformed zone with an associated framework of satin spar gypsum veins. A field structural analysis reveals that there are two sets of veins: one associated with a pre-folding stage (before the Alpine orogeny), and another one related to the development of the Montsant anticlinorium, and therefore syn-folding (and Alpine in age)

Fluid pressure drops during stimulation of segmented faults in deep geothermal reservoirs

Hydraulic stimulation treatments required to produce deep geothermal reservoirs present the risk of generating induced seismicity. Understanding the processes that operate during the stimulation phase is critical for minimising and preventing the uncertainties associated with the exploitation of these reservoirs. It is especially important to understand how the phenomenon of induced seismicity is related to the pressurisation of networks of discrete fractures. In this study, we use the numerical simulator CFRAC to analyse pressure drops commonly observed during stimulation of deep geothermal wells. We develop a conceptual model of a fractured geothermal reservoir to analyse the conditions required to produce pressure drops and their consequences on the evolution of seismicity, fluid pressure, and fracture permeability throughout the system. For this, we combine two fracture sets, one able to be stimulated by shear-mode fracturing and another one able to be stimulated by opening-mode fracturing. With this combination, the pressure drop can be triggered by a seismic event in the shear-stimulated fracture that is hydraulically connected with an opening-mode fracture. Our results indicate that pressure drops are not produced by the new volume created by shear dilatancy, but by the opening of the conjugated tensile fractures. Finally, our results reveal that natural fracture/splay fracture interaction can potentially explain the observed pressure drops at the Rittershoffen geothermal site

3DHIP-Calculator A New Tool to Stochastically Assess Deep Geothermal Potential Using the Heat-In-Place Method from Voxel-Based 3D Geological Models

The assessment of the deep geothermal potential is an essential task during the early phases of any geothermal project. The well-known 'Heat-In-Place' volumetric method is the most widely used technique to estimate the available stored heat and the recoverable heat fraction of deep geothermal reservoirs at the regional scale. Different commercial and open-source software packages have been used to date to estimate these parameters. However, these tools are either not freely available, can only consider the entire reservoir volume or a specific part as a single-voxel model, or are restricted to certain geographical areas. The 3DHIP-Calculator tool presented in this contribution is an open-source software designed for the assessment of the deep geothermal potential at the regional scale using the volumetric method based on a stochastic approach. The tool estimates the Heat-In-Place and recoverable thermal energy using 3D geological and 3D thermal voxel models as input data. The 3DHIP-Calculator includes an easy-to-use graphical user interface (GUI) for visualizing and exporting the results to files for further postprocessing, including GIS-based map generation. The use and functionalities of the 3DHIP-Calculator are demonstrated through a case study of the Reus-Valls sedimentary basin (NE, Spain)

Induced seismicity in pressurised single fractures: a numerical approach

The exploration and exploitation of deep geothermal reservoirs has significantly increased during the last years. These reservoirs use heat exchange to produce heat or electricity. The so-called Enhanced Geothermal Systems (EGS) are characterized by a stimulation phase that aims to increase fluid flow and heat transfer between wells by increasing the permeability and transitivity of the reservoir. This is achieved by injecting high-pressure fluids (normally water) in order to increase the apertures of existing fractures, enhancing their sliding and/or generating new ones. However, this technique induces low-magnitude seismicity that occasionally results in damage at the Earth's surface. Numerical simulations able to reproduce the hydro-thermo-mechanical behaviour of geological reservoirs are an essential tool for the evaluation and forecasting of induced seismicity in such systems. In this study, the numerical code CFRAC is used to systematically evaluate how the orientation of faults with respect to the stress field influences seismicity, the injection rate and the fracture sliding behaviour

Subgrain rotation recrystallization during shearing: insights from full-field numerical simulations of halite polycrystals

We present, for the first time, results of full-field numerical simulations of subgrain rotation recrystallization of halite polycrystals during simple shear deformation. The series of simulations show how microstructures are controlled by the competition between (i) grain size reduction by creep by dislocation glide and (ii) intracrystalline recovery encompassing subgrain coarsening by coalescence through rotation and alignment of the lattices of neighboring subgrains. A strong grain size reduction develops in models without intracrystalline recovery, as a result of the formation of high-angle grain boundaries when local misorientations exceed 15°. The activation of subgrain coarsening associated with recovery decreases the stored strain energy and results in grains with low intracrystalline heterogeneities. However, this type of recrystallization does not significantly modify crystal preferred orientations. Lattice orientation and grain boundary maps reveal that this full-field modeling approach is able to successfully reproduce the evolution of dry halite microstructures from laboratory deformation experiments, thus opening new opportunities in this field of research. We demonstrate how the mean subgrain boundary misorientations can be used to estimate the strain accommodated by dislocation glide using a universal scaling exponent of about 2/3, as predicted by theoretical models. In addition, this strain gauge can be potentially applied to estimate the intensity of intracrystalline recovery, associated with temperature, using quantitative crystallographic analyses in areas with strain gradients

Sulfur and lead isotope systematics: Implications for the genesis of the Riópar Zn-(Fe-Pb) carbonate-hosted deposit (Prebetic Zone, SE Spain)

The Zn-(Fe-Pb) deposits of the Riópar area (Prebetic Zone, SE Spain) are hosted by dolostones that replace Berriasian to Valanginian (Upper Jurassic-Lower Cretaceous) limestones. Mineralization consists of hypogene sphalerite, marcasite and galena, and supergene calamine zones. The hypogene ores are associated with a saddle dolomite gangue. The ore bodies occur as discordant and stratiform lenses, ore-cemented breccias, cm- to mm-wide veins and veinlets, disseminations and stylolite porosity filling within the host dolomites. The main ore controls include stratigraphy and/or lithology, tectonics (faults, fractures and breccias) and availability of metals and sulfur. The morphologies and epigenetic character of the hypogene ore bodies are consistent with the classification of this mineralization as a Mississippi Valley-type (MVT) deposit. The Ga/Ge geothermometer in sphalerite yielded a temperature range of 194-252ºC, which represents the temperature of the source region of the ore solution. This value is comparable to the temperature obtained in the ore deposition site, 159±15ºC from the Δ34S geothermometer in sphalerite galena pairs. This similitude points to a hydrothermal fluid that did not cool down significantly during flow from the fluid reservoir area to the precipitation site. δ34S values of base-metal sulfides (-7.5 to +3.5 ¿) are consistent with thermochemical reduction of Triassic sulfate (seawater and/or derived from dissolution of evaporites) by interaction with organic compounds (e.g., hydrocarbons, methane), which reduced sulfate to sulfide in the deposition site. The lead isotope ratios (206Pb/204Pb = 18.736-18.762; 207Pb/204Pb = 15.629-15.660; 208Pb/204Pb = 38.496-38.595) of galena suggest that Pb, and probably other metals as Zn, is derived from continental crustal rocks. On the other hand, these relations points to an unique metal source probably derived from the Paleozoic basement rocks. The relationship between bedding-parallel stylolites, dolomitization, sulfide precipitation and Alpine tectonic affecting the MVT ore, suggests a relative timing range for the mineralization in the Riópar area of 95 to 20 Ma (Upper Cretaceous-Tertiary). The sulfide mineralization and the associated dolomitization are thus explained by the contribution of two fluids that mixed in different proportions during dolomitization and mineralization: i) a fluid probably derived from Cretaceous seawater saturating Mesozoic sediments (Fluid A), characterized by being dilute and initially low temperature, which should have contained organic rich compounds in the ore deposition site (e.g., hydrocarbons and CH4 dissolved gas); and ii) a high salinity hydrothermal brine (Fluid B) rich in both metals and sulfate, circulated through the Paleozoic basement. During the pre-ore dolomitizing stage the fluid phase was dominated by the diluted fluid (Fluid A > Fluid B), whereas in a later fluid pulse, the proportion of the high salinity fluid increased (Fluid A < Fluid B) which allowed sulfide precipitation. MVT exploration in the Prebetic Zone should focus towards the SW of the Riópar mines, in the vicinity of the Alto Guadalquivir-San Jorge fault

Seismic Anisotropy of Temperate Ice in Polar Ice Sheets

We present a series of simple shear numerical simulations of dynamic recrystallization of two‐phase nonlinear viscous materials that represent temperate ice. First, we investigate the effect of the presence of water on the resulting microstructures and, second, how water influences on P wave (Vp) and fast S wave (Vs) velocities. Regardless the water percentage, all simulations evolve from a random fabric to a vertical single maximum. For a purely solid aggregate, the highest Vp quickly aligns with the maximum c‐axis orientation. At the same time, the maximum c‐axis development reduces Vs in this orientation. When water is present, the developed maximum c‐axis orientation is less intense, which results in lower Vp and Vs. At high percentage of water, Vp does not align with the maximum c‐axis orientation. If the bulk modulus of ice is assumed for the water phase (i.e., implying that water is at high pressure), we find a remarkable decrease of Vs while Vp remains close to the value for purely solid ice. These results suggest that the decrease in Vs observed at the base of the ice sheets could be explained by the presence of water at elevated pressure, which would reside in isolated pockets at grain triple junctions. Under these conditions water would not favor sliding between ice grains. However, if we consider that deformation dominates over recrystallization, water pockets get continuously stretched, allowing water films to be located at grain boundaries. This configuration would modify and even overprint the maximum c‐axis‐dependent orientation and the magnitude of seismic anisotropy

Petrography and geochemistry of fault-controlled hydrothermal dolomites in the Riópar area (Prebetic Zone, SE Spain)

The present paper reports the first detailed petrographical and geochemical studies of hydrothermal dolomites related to MVT Zn-(Fe-Pb) deposits in the Riópar area (Mesozoic Prebetic Basin, SE Spain), constraining the nature, origin and evolution of dolomitizing and ore-forming fluids. Mapping and stratigraphic studies revealed two stratabound dolostone geobodies connected by other patchy bodies, which replace carbonate units of Upper Jurassic to Lower Cretaceous ages. These dolostones are associated to the W-E trending San Jorge fault, indicating a main tectonic control for fluid flow. Seven different dolomite types were identified: i) matrix-replacive planar-s (ReD-I); ii) matrix-replacive planar-e (ReD-II); iii) planar-e sucrosic cement (SuD); iv) non-planar grey saddle dolomite cement (SaD-I) pre-dating Zn-(Fe-Pb) sulfides; v) non-planar milky to pinkish saddle dolomite cement (SaD-II) post-dating Zn-(Fe-Pb) ores; vi) ore-replacive planar-e porphyrotopic (PoD); and vii) planar-s cloudy cement (CeD). Meteoric calcite types were also recognized. The different dolomite types are isotopically characterized by: i) depleted δ18O (from +25.1 to +27.6¿ V-SMOW) and δ13C (from -2.3 to +0.9¿ V-PDB) values compared to Upper Jurassic to Lower Cretaceous limestone signature (δ18O: +27.6 to +30.9¿ V-SMOW; δ13C: +0.5 to +3.2¿ V-PDB); and ii) 87Sr/86Sr ratios for the main dolomitization phases (ReD and SuD: 0.70736-0.70773) close to the Jurassic and Cretaceous carbonate values (0.70723-0.70731) whereas more radiogenic values (0.70741-0.70830) for saddle dolomites (SaD) related to the Zn-(Fe-Pb) sulfide mineralization prevailed after fluid interaction with Rb-bearing minerals. Microthermometrical studies on two-phase liquid and vapor fluid inclusion populations in planar and non-planar dolomites and sphalerite show homogenization temperatures between 150 and 250ºC. These data indicate that both planar and non-planar dolomite textures formed at high-temperatures under hydrothermal conditions in deep-burial diagenetic environments. The main dolomitizing phase (ReD-I/ReD-II and SaD-I) shows low to moderate fluid inclusions salinity (5 to 14 wt.% eq. NaCl), whereas the dolomitization related to ore precipitation (sphalerite and SaD-II) spreads to higher salinity values (5 to 25 wt.% eq. NaCl). These data may respond to a mixing between a low salinity fluid (fluid A, less than 5 wt.% eq. NaCl) and a more saline brine (fluid B, more than 25 wt.% eq. NaCl) at different fluid proportions

Geology of the Cerro Quema Au-Cu deposit (Azuero Peninsula, Panama)

The Cerro Quema district, located on the Azuero Peninsula, Panama, is part of a large regional hydrothermal system controlled by regional faults striking broadly E-W, developed within the Río Quema Formation. This formation is composed of volcanic, sedimentary and volcano-sedimentary rocks indicating a submarine depositional environment, corresponding to the fore-arc basin of a Cretaceous–Paleogene volcanic arc. The structures observed in the area and their tectono-stratigraphic relationship with the surrounding formations suggest a compressive and/or transpressive tectonic regime, at least during Late Cretaceous–Oligocene times. The igneous rocks of the Río Quema Formation plot within the calc-alkaline field with trace and rare earth element (REE) patterns of volcanic arc affinity. This volcanic arc developed on the Caribbean large igneous province during subduction of the Farallon Plate. Mineralization consists of disseminations of pyrite and enargite as well as a stockwork of pyrite and barite with minor sphalerite, galena and chalcopyrite, hosted by a subaqueous dacitic lava dome of the Río Quema Formation. Gold is present as submicroscopic grains and associated with pyrite as invisible gold. A hydrothermal alteration pattern with a core of advanced argillic alteration (vuggy silica with alunite, dickite, pyrite and enargite) and an outer zone of argillic alteration (kaolinite, smectite and illite) has been observed. Supergene oxidation overprinted the hydrothermal alteration resulting in a thick cap of residual silica and iron oxides. The ore minerals, the alteration pattern and the tectono-volcanic environment of Cerro Quema are consistent with a high sulfidation epithermal system developed in the Azuero peninsula during pre-Oligocene time