Fault Zone Permeability Research Articles

Fault‐Valve Instability: A Mechanism for Slow Slip Events

AbstractGeophysical and geological studies provide evidence for cyclic changes in fault‐zone pore fluid pressure that synchronize with or at least modulate slip events. A hypothesized explanation is fault valving arising from temporal changes in fault zone permeability. In our study, we investigate how the coupled dynamics of rate and state friction, along‐fault fluid flow, and permeability evolution can produce slow slip events. Permeability decreases with time, and increases with slip. Linear stability analysis shows that steady slip with constant fluid flow along the fault zone is unstable to perturbations, even for velocity‐strengthening friction with no state evolution, if the background flow is sufficiently high. We refer to this instability as the “fault valve instability.” The propagation speed of the fluid pressure and slip pulse, which scales with permeability enhancement, can be much higher than expected from linear pressure diffusion. Two‐dimensional simulations with spatially uniform properties show that the fault valve instability develops into slow slip events, in the form of aseismic slip pulses that propagate in the direction of fluid flow. We also perform earthquake sequence simulations on a megathrust fault, taking into account depth‐dependent frictional and hydrological properties. The simulations produce quasi‐periodic slow slip events from the fault valve instability below the seismogenic zone, in both velocity‐weakening and velocity‐strengthening regions, for a wide range of effective normal stresses. A separation of slow slip events from the seismogenic zone, which is observed in some subduction zones, is reproduced when assuming a fluid sink around the mantle wedge corner.

Transverse faulting in the Paradox Basin, Colorado Plateau, USA: Active, intermittent faulting over a 300 million year history of strain accommodation and fluid flow

The Paradox Basin (Colorado Plateau, USA) is dominated by major, first-order northwest-trending structures commonly 40 km or more in trace length, including (1) regional salt wall corridors related to passive diapirism during Pennsylvanian to Jurassic time, (2) gentle upright folds produced by Laramide shortening during the Late Cretaceous through early Cenozoic, (3) late Laramide normal faults, and (4) normal faults representing Neogene salt dissolution collapse. Less obvious at a regional scale are the fault zones aligned perpendicular (northeast-trending) to the dominant structural grain. There are 16 such faults zones, marked by short trace lengths (3−12 km), small offsets (10−100 m), and predominantly extensional kinematics. Based on published geological maps, field observations of fault zone properties (including fluid flow indicators), and kinematic analysis, we interpret these structures as transverse accommodation faults. Co-spatial structural associations reveal the transverse fault zones were active intermittently, likely as a means of minimizing along-strike strain incompatibility that accrued along the first-order structures as they evolved during late Paleozoic to Late Jurassic halokinesis, (mild) Late Cretaceous to Eocene folding, late Laramide normal faulting, and Neogene to recent collapse faulting. Locally, transverse faulting was influenced by reactivation of northeast-striking faults that offset sub-salt formations, including basement. Active, intermittent transverse faulting over the past ∼300 m.y. is consistent with the synthesis of published interpretations and age determinations focusing on the timing of diverse fluids that exploited the permeability of the transverse fault zones. The Paradox Basin, with its enormous subsurface salt volume and enduring fluid flow, has been an ideal dynamic environment for producing intermittently active transverse accommodation faulting.

Development of fluid pathways and associated diagenetic variations in a natural CO2 leaking fault zone.

Abstract Permeability within fault zones can vary through time due to repeated deformation events and rock–fluid interactions. Understanding the history of fault zone alteration is critical when building hydrogeologic models and evaluating the risk of mechanical rock failure during subsurface storage. Newly acquired drill core recovered within the fault damage zone and fault core of the Little Grand Wash fault (LGWF) is combined with observations from outcrop, optical petrography, computerized tomography image analysis, and ultrasonic velocity measurements to characterize the rock types and preserved structural deformation features within this fault zone. These data are used to understand the history of mechanical rock failure and mineralization associated with this fluid-charged fault system. We identify multiple structural features and use their cross-cutting relationships to understand the history of deformation and their association with changes in fault zone permeability and rock mechanical properties. At the LGWF zone, structural deformation features vary temporally and are used to recognize a decoupling between fault slip and fluid flow. The formation of most open-mode veins occurred after shear failure within the LGWF zone. Early-developed shear bands are cut by carbonate veins, that are in turn cut by shear fractures, followed by a second phase of vein formation and ultimately by folding in the fault core. This sequence of formation reflects changes in mechanical rock properties due to subsurface rock–fluid interactions and indicates that the alteration of rock matrix by secondary carbonate cement results in increased unconfined compressive strength, decreased permeability, and increased ultrasonic velocity. In this fault zone, and possibly other fault zones, the changes in rock properties associated with deformation can be detected through their geophysical response.

Physico-chemical constraints associated with hypogene hydrothermal alteration and supergene oxidation at the Farallón Negro intermediate-sulfidation epithermal Au-Ag deposit, NW Argentina

The Farallón Negro deposit is located in the Farallón Negro Volcanic Complex (FNVC; ∼9.40 to 4.88 Ma) in NW Argentina and includes intermediate-sulfidation, epithermal Au-Ag mineralization in the Alto de la Blenda deposit. Detailed mineralogical and geochemical investigations of the Esperanza-Esperanza Sureste, Cecilia veins, both hosted in the Alto de la Blenda monzonite, and the Laboreo vein, hosted in the volcaniclastic Morada Central andesitic breccia, reveal four hypogene vein stages. Each of these four stages are comprised of a texturally early quartz, followed by a sulfide-sulfosalt and a late Mn-carbonate sub-stage. Stage 1 is characterized by a quartz, base metal and Mn-carbonate sub-stage showing high δ18OVSMOW (δ18OQuartz = 13.7 ± 0.4‰ and δ18OMn-Calcite = 11.6 ± 0.2‰) and low δ13CVPDB values (δ13CMn-Calcite = –4.7 ± 0.2‰). Stage 2 includes most gold and silver minerals in a tennantite-tetrahedrite-sphalerite-pyrite-galena±electrum sub-stage. Fluid inclusions trapped in vein-anhydrite, only associated with stage 2, record the highest temperatures (TEntrapment = 327–335°C) and highest salinity (3.8NaCleq.wt.%) when compared to all other stages. Stage 2 quartz (δ18OQuartz = 9.9 to 11.9‰) and carbonate display lower δ18OVSMOW (δ18OMn-Calcite = 7.9 to 10.3‰), as well as lower δ13CVPDB values (δ13CMn-Calcite = -5.2 to 4.7‰) than all other stages. Stages 3 and 4 consist of barren quartz and carbonate sub-stages that show lower temperatures (TEntrapment = 187–262°C) and higher δ18OVSMOW (δ18OQuartz = 13.4 to 13.9‰ and δ18OMn-Calcite = 10.7 to 11.8‰) and δ13CVPDB values (δ13CMn-Calcite = –3.5 to −3.9‰) when compared to earlier stages.Petrographic, sulfide mineral chemistry, microthermometric, and stable isotope data combined with distinct alteration mineral assemblages allow the reconstruction of a model for the Alto de la Blenda deposit. During the hypogene hydrothermal fluid flow stage, around 7.16 Ma, at ∼2.2 km depth (PHYDROSTATIC = 220 bars), the high temperature alteration mineralogy, low fluid salinities and 18O and 13C isotope compositions indicate cooling of ascending magmatic-hydrothermal fluids. These fluids were derived from a contracted magmatic vapor phase that exchanged with unexposed sedimentary rocks of the Sierras Pampeanas. Hypogene quartz-R3 illite/smectite-chlorite-calcite-pyrite alteration in the andesitic-breccia and quartz-illite-R3 illite/smectite-chlorite-calcite-pyrite alteration in the monzonite are temporally and spatially to the hydrothermal vein minerals deposited between 187 and 335°C. During the post-hydrothermal supergene weathering event, around 2.70 Ma, sulfuric acids generated by meteoric fluids and supergene oxidation of hypogene sulfides, particularly pyrite, descended along normal fault zones and related vein and fracture networks to at least 450 m depth below the surface level, forming the supergene alteration assemblage of dioctahedral smectite-kaolinite-Mn-/Fe-oxides/-hydroxides-gypsum±interstratified chlorite/smectite±interstratified R1 illite/smectite in both wall rock types at temperatures <50°C.The dominant hypogene alteration assemblage, present in the monzonite and volcaniclastic andesitic breccia, may be indicative of deeper-seated, transitional magmatic-hydrothermal systems (>1 km depth), generally controlled by fluid cooling, as is the case at Farallón Negro. Pervasive supergene alteration can conceal these systems but the identified supergene alteration assemblage of dioctahedral smectite-kaolinite-interstratified chlorite/smectite-interstratified R1 illite/smectite that replaces the hypogene alteration assemblage along the permeable fault, vein and fracture zone network may present a useful vectoring tool towards concealed, deeper-seated, epithermal Au-Ag mineralization in arid climates.

Implications of Fault‐Valve Behavior From Immediate Aftershocks Following the 2023 Mj6.5 Earthquake Beneath the Noto Peninsula, Central Japan

AbstractThe Mj6.5 (Mw6.2) event that occurred on 5 May 2023 near the northern shoreline of the northeastern tip of the Noto Peninsula, central Japan, is the largest event to date in a long‐lasting, intense earthquake swarm. Here we have created a more precise aftershock catalog associated with the 2023 Mj6.5 and the second‐largest 2022 Mj5.4 sequence to understand the rupture process of this largest earthquake. Most of the aftershocks are aligned along a ∼45° SE‐dipping plane. The mainshock initially ruptured the same deep section of the fault zone that had been ruptured by the 2022 Mj5.4 event, before propagating rapidly to shallow depths and to offshore along the ruptured fault plane. The aftershock front migrated at a speed of ∼20 km/hr. This rapid upward migration of the immediate aftershocks might be driven by upwelling of crustal fluids along the intensely fractured and permeable fault zone via mainshock dynamic rupture.

Along‐Strike Segmentation of Seismic Tremor and Its Relationship With the Hydraulic Structure of the Subduction Fault Zone

AbstractAlong the strike of subduction zones, tectonic tremor episodicity is segmented on a geologic scale. Here, we study how this segmentation reflects large‐scale variations of the structure and conditions of the fault interface where tremor is generated. We try to understand which properties of the hydraulic system of the fault allow elementary tremor sources to synchronize, leading to the emergence of long‐period, large‐scale episodic activity. We model tremor sources as being associated with rapid openings of low‐permeability valves in the fault zone, which channels the upward flow of metamorphic fluids. Valve openings cause pressure transients that allow interaction between sources. In such a system, tremor activity is thus controlled by unsteady fluid circulation. Using numerical simulations of fluid flow, we explore the impact of valve spatial distribution and fluid flux on the emergence of large‐scale patterns of tremor activity. We show that when valves are densely distributed and submitted to near‐critical input flux, they synchronize and generate more episodic activity. Based on our model, the most periodic and spatially coherent tremor bursts should thus be emitted from segments densely populated with valves, and therefore of lower permeability than less synchronized segments. The collective activity of their valve population is responsible for fluid‐pressure cycling at the subduction scale. In the tremor zone of Shikoku, Japan, the most temporally clustered segment coincides with a downgoing seamount chain, suggesting that the segmentation of the fault zone permeability, and hence of tremor activity, could be inherited from the topography of the subducting oceanic plate.

3D topological analysis in deformation bands: Insights for structural characterization and impact on permeability

This study uses a 3D exposure of a fault zone in the Rio do Peixe Basin (NE Brazil) to investigate the impact of deformation band connection on permeability in porous sandstones. We analyze the topological deformation band distribution pattern and in situ permeability measurements. The results show that more Y-nodes dominate the fault core, while the damage zone registers more Y- and X-nodes. The I-nodes are associated with the highest permeability values, while X-nodes have lower values close to the fault core. In the damage zone, permeability acquisition on segment intersections (nodes) varies according to the type of node and directions (horizontal and vertical) of measurements. The vertical section (horizontal permeability) shows maximum values for I-nodes and minimum for X-nodes, whereas the horizontal section (vertical permeability) exhibits higher permeabilities for I-nodes and the lowest for Y-nodes. Therefore, this study contributes with new perspectives on the influence of deformation band connections on the permeability of fault zones and has implications for the modeling and flow simulation of siliciclastic reservoirs with deformation bands in the damage zones.

Geophysical characteristics of a fault system in the northern Central Range of Taiwan and its applications for geothermal energy exploration

The northern Central Range of Taiwan is a high-potential geothermal region. Since the formations are mainly tight metasandstone and slate, permeable structures associated with faults are commonly considered as conduits of geothermal fluids. This study determines the characteristics and orientations of the permeable fault zones by analyzing the geophysical logs and microresistivity formation image log (FMI) of the JT-4 well in Jentse, an important geothermal area in the northern Central Range. Between 720 and 1480 m measured depth (MD), the effective porosity of the intact host rock is mostly below 3% calculated by the geophysical log. Zones with porosity greater than 5% are only clustered within a few thin intervals. The FMI interpretations show these porous zones are in the interior of the fractured and faulted intervals. These porous fault zones comprise fault damage zones with a high density of open fracture planes and fault cores with porous fault breccias. There is a highly brecciated fault core in 1334–1339 m MD, which would be the most permeable interval of the well. Additionally, some healed fault zones with sealed fractures are observed. The picked drilling-induced tensile fractures signify that the direction of the present-day maximum horizontal principal stress is N40–50°E, and most of the open fractures also strike parallel to the NE–SW direction. The study results show that the open fractures are concentrated in the four fault zones belonging to one major normal fault system. After integrating the orientations and locations of the fault zones, we propose that the permeable normal fault system is about 200 m wide, trends N50–70°E, and dips 70–80° to the NW. The development of the open fractures and the permeable fault system in the northern Central Range may be controlled by the current rifting of the Okinawa Trough offshore northeastern Taiwan. The study exhibits the characteristics of fractured fluid conduits of the regional geothermal system, which will benefit future geothermal exploration in northeastern Taiwan.

3D coupled hydro-mechanical modeling of multi-decadal multi-zone saltwater disposal in layered and faulted poroelastic rocks and implications for seismicity: An example from the Midland Basin

Knowledge of subsurface coupled pore pressure and stress perturbations due to saltwater disposal (SWD) is crucial for understanding and managing induced seismicity. Current studies choose Biot poroelasticity to implement fluid-rock coupling and perform numerical modeling in relatively simple settings. In this study, we first introduce a general poromechanical framework to capture a broader range of hydro-mechanical coupling and then apply it to model real-world SWD problems with geological and operational complexities. We demonstrate a viable workflow through the case study of a 2500 km2, layered, faulted, and recently seismically awaken area that has undergone multi-decadal multi-zone injection by nearly 200 SWD wells in the Midland basin. The example provides a unique opportunity to study the respective impact of shallow and deep injections across depths and the role of inter-zone interactions, among others. Our modeling suggests over 500 and 250 psi increases in pore pressure and Coulomb stress, respectively, in the deep disposal zone in the last decade, which were 5 to 10 times greater than those accumulated over the past four decades in the two shallower disposal zones. Vertically, shallower injections impacted mostly shallower zones, whereas the deep injection produced changes across depths and dominated especially in the deep zone and the upper segment of the basement. Towards the mid-segment of the basement, changes due to shallow and deep injections were comparable, with the latter 2 to 5 times higher. Laterally, in-zone changes were dominated by pore pressure diffusion in the near field and poroelastic stressing in the far field in shallower zones but were dominated by pore pressure diffusion everywhere in the deep zone, owing to their different hydraulic diffusivities. Faults acted as pore pressure and Coulomb stress sinks for disposal zones and sources for non-disposal zones, and fault-zone permeability structure had a minor impact. Finally, seismological modeling suggests that seismicity in the area was most likely driven by the deep injection.

Characterizing along- and across-fault fluid-flow properties for assessing flow rates and overburden fluid migration along faults: a case study from the North Sea

Assessing fault zones as fluid-migration pathways requires the characterization of permeability both across and along faults, as well as the adjacent volume. The hydraulic properties of the Vette Fault Zone, North Sea, are described by modelling the mixing of host-rock lithologies into the fault zone, and the fault width is derived from empirical relationships as a function of throw and clay content. To better understand the sensitivity related to the uncertainties in overburden lithologies and fault-width correlations, a parametric study with 1125 model realizations were solved in a 2D steady-state, single-phase, subsurface flow model. The fault zone, included as a discrete permeable structure, significantly alters the flow field compared to a model that only considers lithological juxtaposition. The most prominent hydraulic communication in the Vette Fault Zone is downwards from the storage reservoir where sand is mixed into the fault zone. Increasing the host-rock permeability in the overburden also increases the fault permeability and shifts the inflection point for down-fault flow, causing the pressurized reservoir to drain towards the overburden and the top surface. For CO 2 storage application, the models highlighted the potential for downward communication along the fault for brine, and the CO 2 capillary sealing towards the overburden. Thematic collection: This article is part of the Fault and top seals 2022 collection available at: www.lyellcollection.org/topic/collections/fault-and-top-seals-2022

Episodic transient deformation revealed by the analysis of multiple GNSS networks in the Noto Peninsula, central Japan

Since November 30, 2020, an intense seismic swarm and transient deformation have been continuously observed in the Noto Peninsula, central Japan, which is a non-volcanic/geothermal area far from major plate boundaries. We modeled transient deformation based on a combined analysis of multiple Global Navigation Satellite System (GNSS) observation networks, including one operated by a private sector company (SoftBank Corp.), relocated earthquake hypocenters, and tectonic settings. Our analysis showed a total displacement pattern over 2 years shows horizontal inflation and uplift of up to ~ 70 mm around the source of the earthquake swarm. In the first 3 months, the opening of the shallow-dipping tensile crack had an estimated volumetric increase of ~ 1.4 × 107 m3 at a depth of ~ 16 km. Over the next 15 months, the observed deformation was well reproduced by shear-tensile sources, which represent an aseismic reverse-type slip and the opening of a southeast-dipping fault zone at a depth of 14–16 km. We suggest that the upwelling fluid spread at a depth of ~ 16 km through an existing shallow-dipping permeable fault zone and then diffused into the fault zone, triggering a long-lasting sub-meter aseismic slip below the seismogenic depth. The aseismic slip further triggered intense earthquake swarms at the updip.

Characteristics of Slow Slip Events Explained by Rate‐Strengthening Faults Subject to Periodic Pore Fluid Pressure Changes

AbstractGeophysical observations suggest that temporal changes in pore fluid pressure correlate with slow slip events (SSEs) at some subduction zones, including the Hikurangi and Cascadia subduction zones. These fluctuations in pore fluid pressure are attributed to fluid migration before and during SSEs, which may modulate SSE occurrence. To examine the effect of pore fluid pressure changes on SSE generation, we develop numerical models in which periodic pore‐pressure perturbations are applied to a stably sliding, rate‐strengthening fault. By varying the physical characteristics of the pore‐pressure perturbations (amplitude, characteristic length and period), we find models that reproduce shallow Hikurangi SSE properties (duration, magnitude, slip, recurrence) and SSE moments and durations from different subduction zones. The stress drops of modeled SSEs range from ∼20–120 kPa while the amplitudes of pore‐pressure perturbations are several MPa, broadly consistent with those inferred from observations. Our results indicate that large permeability values of ∼10−14 to 10−10 m2 are needed to reproduce the observed SSE properties. Such high values could be due to transient and localized increases in fault zone permeability in the shear zone where SSEs occur. Our results suggest that SSEs may arise on faults in rate‐strengthening frictional conditions subject to pore‐pressure perturbations.

Constraints on upper crustal fluid circulation and seismogenesis from in-situ outcrop quantification of complex fault zone permeability

The permeability of fault zones plays a significant role on the distribution of georesources and on seismogenesis in the brittle upper crust, where both natural and induced seismicity are often associated with fluid migration and overpressure. Detailed models of the permeability structure of fault zones are thus necessary to refine our understanding of natural fluid pathways and of the mechanisms leading to fluid compartmentalization and possible overpressure in the crust. Fault zones commonly contain complex internal architectures defined by the spatial juxtaposition of “brittle structural facies” (BSF), which progressively and continuously form and evolve during faulting and deformation. We present the first systematic in-situ outcrop permeability measurements from a range of BSFs from two architecturally complex fault zones in the Northern Apennines (Italy). A stark spatial heterogeneity of the present-day permeability (up to four orders of magnitude) even for tightly juxtaposed BSFs belonging to the same fault emerges as a key structural and hydraulic feature. Insights from this study allow us to better understand how complex fault architectures steer the 3D hydraulic structure of the brittle upper crust. Fault hydraulic properties, which may change through space but also in time during an orogenesis and/or individual seismic cycles, in turn steer the development of overpressured volumes, where fluid-induced seismogenesis may localize.

Tectonostratigraphic controls on pore fluid pressure distribution across the Taranaki Basin, New Zealand

AbstractSignificant variations in pore pressure across the Taranaki Basin, New Zealand, are attributed to changes in lithofacies and structure, usefully illustrated in terms of ten areas that we term geopressure provinces, each displaying individual pore pressure trends. Cretaceous to Early Miocene formations in different parts of the basin can be either normally pressured (near or at hydrostatic) or significantly overpressured (up to 28 MPa) at the same depth. Variations in Eocene–Oligocene facies types and thicknesses both within and between geopressure provinces provide first‐order control on the magnitude, distribution and maintenance of overpressure across the basin. Examples of hydraulic compartmentalisation due to sealing faults and stratigraphic architecture are identified within the basin. Deep pore pressure transitions are sealed by diagenetic, structural or stratigraphic mechanisms in different places and are associated with an increase in mudrock volume (reduced permeability) or gas generation. Thus, pore pressure distribution in the Taranaki Basin is controlled by a combination of sediment loading, lithofacies variations, fault zone permeability and structural architecture. This work represents an appraisal of the pore pressure distribution across the whole of a multiphase structurally complex basin, and the approach taken provides a framework for better understanding the distribution of pore fluid pressures and pore fluid migration in other sedimentary basins.

Hydrochemical and stable isotopes (δ2H and δ18O) changes of groundwater from a spring induced by local earthquakes, Northwest China

It has been well reported that earthquakes can cause changes in groundwater chemistry and isotopes, and much of those changes were occurred in wells or hot springs; however, changes in cold spring caused by local small earthquakes have received less attention. Here, we collected continuous monitoring of the X10 spring (September 2018 to December 2019), investigated the hydrogeochemical characteristics of the spring by using water chemistry analysis and isotope methods. We compared the changes in water chemical ion concentrations and hydrogen-oxygen isotope ratios with the surrounding seismic activity, and the results show that 1) major chemical ion concentrations in X10 springs have an annual dynamic pattern of being high in winter and low in summer, and this change may be related to the seasonal effect of snowmelt; 2) the spring water originates from long-distant meteoric water, snowmelt, and bedrock fissure water and is affected by rock weathering and evaporation; 3) the hydrogen and oxygen stable isotope ratios and HCO3− concentration in groundwater are sensitive to local small seismic activity. We considered that small earthquakes can change the permeability in fault zones or aquifers, leading to mixing of groundwater with different chemical composition and isotopes. Our study demonstrates that the hydrogen and oxygen stable isotopes are more sensitive to seismic activity than the commonly used chemical constitutents, and that the sensitive constitutents vary in different observation wells or springs; therefore, combined monitoring of isotopes and water chemistry should be considered in the future to capture hydrogeochemical precursor signals caused by earthquakes.

Hindcasting injection-induced aseismic slip and microseismicity at the Cooper Basin Enhanced Geothermal Systems Project

There is a growing recognition that subsurface fluid injection can produce not only earthquakes, but also aseismic slip on faults. A major challenge in understanding interactions between injection-related aseismic and seismic slip on faults is identifying aseismic slip on the field scale, given that most monitored fields are only equipped with seismic arrays. We present a modeling workflow for evaluating the possibility of aseismic slip, given observational constraints on the spatial-temporal distribution of microseismicity, injection rate, and wellhead pressure. Our numerical model simultaneously simulates discrete off-fault microseismic events and aseismic slip on a main fault during fluid injection. We apply the workflow to the 2012 Enhanced Geothermal System injection episode at Cooper Basin, Australia, which aimed to stimulate a water-saturated granitic reservoir containing a highly permeable (k = 10^{-13} - 10^{-12}hbox {m}{^2}) fault zone. We find that aseismic slip likely contributed to half of the total moment release. In addition, fault weakening from pore pressure changes, not elastic stress transfer from aseismic slip, induces the majority of observed microseismic events, given the inferred stress state. We derive a theoretical model to better estimate the time-dependent spatial extent of seismicity triggered by increases in pore pressure. To our knowledge, this is the first time injection-induced aseismic slip in a granitic reservoir has been inferred, suggesting that aseismic slip could be widespread across a range of lithologies.

Contribution of Remote Sensing and GIS to the Inventory and Analysis of Factors influencing on the Development of Karst Features in the Middle Atlas, Morocco

The development of karst features is a complex and interactive process of many factors with varying influence over time. The following main factors seem to have influenced the occurrence and development of karst features in the Middle Atlas in Morocco: Precipitations after extreme weather events depending on the precipitation pattern modified by orographic effects and exposition of the terrain towards the main rain-bearing clouds are determining surface water input and infiltration as well as the outline of morphologic watersheds and drainage basins. The lithologic conditions such as the outcrops of lower Jurassic limestones, the geomorphologic disposition, the volcanic activity within limestones combined with hydrothermal dissolution processes, and the structural conditions influencing water permeability in fault and fracture zones and, thus, the dissolution intensity, are further factors. Ongoing aseismic vertical and horizontal neotectonic movements due to plate tectonic activities and shallow earthquake activities are causing the development and reactivation of fault zones and as consequence modifying groundwater flow conditions in the karst hydrologic systems. Seismic activity can lead as well to the development of sinkholes and the collapse of caves. The development of a systematic Karst-Geo Information System (GIS) is recommended that should be updated regularly.

Transport of Geothermal Fluids along Dikes and Fault Zones

Field observations of active and fossil natural geothermal fields indicate that geothermal fluids are primarily transported along dikes and fault zones. Fluid transport along dikes (commonly through fractures at their margins) is controlled by the cubic law where the volumetric flow rate depends on the aperture of the fracture in the 3rd power. Dikes (and inclined sheets) also act as heat sources for geothermal fields. In high-temperature fields in volcanoes in Iceland dikes and inclined sheets constitute 80–100% of the rock at crustal depths of 1.5–2 km. Holocene feeder-dikes are known to have increased the activity of associated geothermal fields. Fault zones transport geothermal fluids along their two main hydromechanical units, the core and the damage zone. The core is comparatively thin and primarily composed of breccia, gouge, and clay and related low-permeability porous materials. By contrast, the fault damage zone is characterised by fractures whose frequency is normally highest at the contact between the core and the damage zone. Fluid transport in the damage zone, and in the core following fault slip, is controlled by the cubic law. During non-slip periods fluid transport in the core is primarily controlled by Darcy’s law. Secondary mineralisation (forming mineral veins and amygdales) tends to reduce the fault-zone permeability. Repeated earthquake activity is thus needed to maintain the permeability of fault zones in active natural geothermal fields.

Borehole observation-based in situ stress state estimation of the Los Humeros geothermal field (Mexico)

In this study, we present results of an in situ stress state estimation for the Los Humeros geothermal field (Mexico) located in the eastern section of the Trans-Mexican Volcanic Belt. During more than 40 years of geothermal energy production, issues related to induced seismicity, reservoir depletion, subsidence, and wellbore stability have occurred in Los Humeros. The in situ stress tensor, being one of the major controls on the above-mentioned phenomena, remains highly uncertain. The high temperatures of the geothermal reservoir at Los Humeros, often exceeding the critical point of water, promote difficulties in performing in situ stress measurements. In this study, based on data commonly acquired from drilling operations from eleven high-temperature wells within the Los Humeros caldera, we constrain the in situ stress state of the geothermal field. For this purpose, we employ borehole observations including borehole breakouts, instances of fluid circulation loss, and injection tests. Results from this study indicate predominantly strike-slip regime with a strong reverse component and NE-SW acting maximum horizontal stress. Utilizing slip and dilation tendency analysis, permeability and seismic activity of major fault zones within the geothermal field could be aptly explained by the acting stress regime. The outcome of this study provides valuable insights into the in situ stress state of the Los Humeros caldera and supports the exploration of conventional and unconventional (i.e., supercritical) geothermal resources.

Origin and migration of fluoride in the area of the Aluto Volcanic Complex (Main Ethiopian Rift)

Fluoride-enriched ground and surface waters represent a major health risk for the local population in many areas along the East African Rift. The present study investigates the origin of fluoride and the reason for its accumulation in the rift waters, following two hypotheses: (i) fluid-rock-interactions release fluoride from minerals into the water and (ii) magmatic-derived fluoride-containing liquids and gases migrate along permeable fault zones until they dissolve in ground- and surface water or be released to the atmosphere.Rock-, gas, water-, soil-, and plant samples were collected from the area within and close by the Aluto Volcanic Complex, which is part of the Main Ethiopian Rift. Most analyzed waters showed fluoride concentrations above the drinking water limit (>1.5 mg/L) with the highest values in hot springs (up to 70 mg/L) and the geothermal well (76 mg/L). In the solid phase, a high fluoride content was found in acid volcanic rocks (ignimbrite: 4391 ppm; rhyolite: 3248 ppm) as well as in pumice (up to 1955 ppm). The fluoride content of soil samples collected within the volcanic complex varied between 82 and 1036 ppm, whereas former lake sediments from outside the Aluto Volcanic Complex showed higher fluoride contents ranging from 674 to 8171 ppm.Identified fluoride-rich minerals are various amphiboles (about 3 wt.-% F−), fluor-apatite (about 3 wt % F−), minerals of the fluor-caphite group (about 5 wt.-% F−), parisite (up to 9 wt.-% F−), and fluorite (CaF2). Elevated fluoride concentrations were also measured in some gas samples from fumaroles (up to 50 ppm) and in plant samples collected next to the fumaroles (up to 65 ppm). Leaching experiments of solid samples with deionized water demonstrated that fluoride can easily be mobilized from pumice and sediments but hardly from consolidated rocks. This fluoride release increased with temperature (up to 150 °C) and correlated roughly with dissolved silica indicating the binding of some fluoride to the amorphous or weakly crystalline silica fraction.Based on these results it was concluded that fluoride migrates via different pathways through the environment: At the depth of the magma chamber during magmatic differentiation fluoride enriches initially in the magmatic melt and accumulates in some late-crystallizing minerals of igneous rocks such as fluorite or mica. Upon volcanic eruption fluoride is predominantly incorporated in the glass (ignimbrite) and ash phase (pumice). On the surface, these rocks are exposed to weathering and fluoride leaches partly out into the aqueous phase. Soft and porous rocks such as pumice release fluoride first whereas hard extrusive/igneous rocks are less prone to weathering and retain the fluoride. Pumice and (lake) sediments might act both, as source and as sink for fluoride in the area. Although some fluoride might drain from the surface into the deeper subsoil, we conclude that magmatic fluids (liquid and gas) contribute more substantially to the overall fluoride ground water content because (i) of the much higher fluoride content in deep geothermal waters and hot springs as compared to surface-near waters; (ii) active geothermal surface manifestations located along fault zones indicate that fluids migrate from deep magmatic intrusions (as gas and liquid) towards the surface, where fluoride dissolves in groundwater; (iii) and the good correlation between bicarbonate (deriving from dissolution of magmatic CO2) and fluoride content in all analyzed water samples.