Fracture and fault networks are characterised by complex spatial and dimensional properties that might affect the flow and accumulation of subsurface fluids. In geological applications, DFN models are commonly employed to compute the multiscale properties of fractured rock volumes in terms of porosity and equivalent permeability. In the present contribution, we focus on an outcrop-to-reservoir investigation of Mesozoic shallow-water carbonates exposed along the axial zone of the southern Apennines fold-and-thrust belt, FTB, Italy. The carbonates were originally deposited in lagoon-to-proximal ramp settings of the Paleo Apenninic Platform during Lower Jurassic–Upper Cretaceous times and include well-layered and massive associations of bed packages several m-thick. By integrating the results of both field and digital structural analyses, we build multiple DFNs to assess the hydraulic behaviour of geocellular volumes representative of the different scales of observation, and stochastically populated with high-angle fractures, small, and medium faults. Fractures are either strata-bound, SB, or non-strata-bound, NSB, and results compartmentalized within single-bed packages. Small faults crosscut multiple bed-packages and show a few cm-to-m throws, whereas medium faults displace multiple bed-package associations and have throws > 1m. Both small and medium faults exhibit high peaks of fracture density, P₂₀, and intensity, P₂₁, in correspondence with the releasing jogs disrupting mechanical interfaces such as bed package boundaries and pre-existing low-angle thrust faults. As data input for DFN modeling, the aperture values of the stochastic fractures are set as proportional to either fracture length (most favourable flow conditions) or to the square root of fracture length (least favourable conditions). At the outcrop scale, 5m-side DFN models show the highest values of fracture porosity among those considered in this work. These results are therefore consistent with both SB and NSB fractures forming the main repository for fluid accumulation. At larger scales, the 50m-side and 500m-side DFN models including small and medium faults, are characterized by higher values of equivalent permeability, which range between 10^-2 and 10^-1 mD. Considering the computed Kxx and Kyy values for the single geocellular volumes, near-isotropic horizontal conditions are assessed across all scales of investigation. Accordingly, altogether, high-angle SB and NSB fractures, small and medium faults form a system of well-connected network through the shallow-water carbonates. Interpreting these data in light of published values for platform carbonates in Italy, we interpret this multiscale horizontal permeability isotropy due to the severe exhumation (~ 4 to 5km) the studied carbonates went through during the Quaternary downfaulting of the southern Apennines FTB.

How to cite: Abdallah, I. B., Panza, E., Dastoli, S., Manniello, C., Prosser, G., and Agosta, F.: Fracture porosity and equivalent horizontal permeability computed for shallow-water carbonates of the southern Apennines fold-and-thrust belt, Italy, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-277, https://doi.org/10.5194/egusphere-egu24-277, 2024.

Diagenetic and tectonic processes taking place in platform carbonates produce significant textural and mineralogical modifications through time, controlling the pore space in terms of dimension, geometry, shape, and connectivity of single pores, and influencing both total and effective porosity. Focusing on the different types of solution surfaces, this study is conducted on Lower Jurassic, Cretaceous, and Eocene, carbonates exposed at the Viggiano Mt. and Raparo Mt., southern Apennines, Italy. Microscale analyses show that the primary porosity of these rocks was occluded by pervasive blocky cements, which precipitated during burial diagenesis of the carbonates. We aim to assess the role exerted by the roughness of bed-parallel and low-angle to bedding solution surfaces, on the pore properties and permeability values of a variety of carbonate lithofacies such as mudstones, packstones, grainstones and rudstones. Specifically, we show the results of Nuclear Magnetic Resonance (NMR), gas-porosimetry and water-permeability tests conducted on plugs cored either orthogonal or parallel to bedding interfaces. All the study plugs show an amount of effective porosity lower than 5%, with mean values of ca. 3%. Excluding larger microfractures, and sporadic intrafossil and intercrystal molds, among the various types of solution surfaces we document that the rough, seismogram-type stylolites localize secondary porosity, while smooth, wave-type stylolites do not. The seimogram type stylolites, due to the non-selective carbonate’s dissolution, form a poorly connected vuggy porosity, and the NMR results the pores are subspherical to tubular (r<3 µm), with low aspect ratios (stiff pores), differently from the pores associated to open fractures (soft pores). Connectivity in the seismogram-type stylolites-related pores is due mainly to small microfractures forming capillary porosity (pore throat ca. r=1 µm). The results of permeability measurements at room pressure indicate that the amount of bed-perpendicular permeability is generally low (10^-1 and 10^-3 mD). The results of permeability measured at increasing confining pressure show that the in stylolite-dominated the bed-parallel samples have slightly higher values with respect to the bed-orthogonal ones and, at increasing confining pressure conditions, the permeability decreases of less than one order of magnitude. Differently, the fracture-dominated plugs show a permeability decrease of two orders of magnitude. These results are therefore consistent with a pore connectivity affected by open fractures only at shallow depths (<25 MPa) and influenced by both stylolites and primary pores at depths. Results of ongoing X-ray tomography analyses will better clarify the pore distribution along stylolites and at the fracture-stylolite intersections.

How to cite: Manniello, C., La Bruna, V., Bezerra, H. F. R., Araùjo, R. E. B., Morais, X. M., Michie, E., Faulkner, D., Allen, M. J., Prosser, G., and Agosta, F.: Pore space properties of solution surfaces in shallow-water carbonates, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-536, https://doi.org/10.5194/egusphere-egu24-536, 2024.

On a global scale, the South Caspian is a basin where mud volcanoes are
most densely distributed. Numerous active volcanoes are registered here, both on
land and in the adjacent sea basin (Caspian Sea), which, as a result of their daily
activity, discharge water fluids to the Earth's surface.
In this work, using samples taken from coastal and island mud volcanoes, a
comparative study of the chemical compositions of the fluids expelled from the
salse and gryphon-type emissions, as well as their origin, including the
contribution of various stratigraphic units to the water formation processes.
According to our results, the contribution of the condensed waters of the
Caspian Sea was greater in the formation of Na-Cl type waters of mud volcanoes
in the South Caspian Basin. In addition, the origin of volcanic waters also formed
as a result of the contribution of shallow ‘low-mineralized’ pore fluids and deep-
seated ‘high-mineralized’ brines correlates well with the depths (1-5.5 km) of the
Productive Series strata.

How to cite: Bayramova, A.: A Comparative study of water fluids of coastal and island mud volcanoes: Acase study of the South Caspian Basin, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-569, https://doi.org/10.5194/egusphere-egu24-569, 2024.

The complexity of fault zone structure and faulting mechanisms significantly impact the permeability of fault zone rocks. Various factors, including the type and porosity of the host rock, the fault’s geometry, differential strain, the history of deformation, and the tectonic setting, can cause permeability to exhibit wide fluctuations. However, research into the permeability of complex fault zones is constrained by the scarcity of in situ measurements. Utilizing analytical relationships, laboratory tests, outcrop measurements, or numerical modeling often yields biased results, as they may not accurately represent real-world conditions. Moreover, the effects of fault branching and interconnections have not been thoroughly explored. This study compares the permeability of faulted carbonate rocks using a comprehensive database of in situ permeability tests. We investigate how a network of seven faults influenced permeability variations at different locations and depths by examining surface and subsurface data, including information from excavated tunnels. Our findings reveal that factors such as fault dips, length, fault structure, and rock characteristics can create diverse impacts on permeability. We observe permeability values in fault damage zones one to five orders of magnitude higher than those in the host rocks. The thickness and condition of damage zones shed light on the range of fault zone permeability. Furthermore, we find that faulted rocks with higher porosity and lower mechanical strength exhibit more substantial alterations in permeability. This study provides valuable insights into the behavior of faulted carbonate rocks.

How to cite: Akbariforouz, M., Zhao, Q., and Zheng, C.: Understanding the Heterogeneity and Anisotropy of Permeability in Carbonate Rocks Within a Fault Network, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-1240, https://doi.org/10.5194/egusphere-egu24-1240, 2024.

Fluid extraction from geological formations for purposes of subsurface utilization leads to pore pressure drop in reservoirs and subsequent compaction and seismicity, especially in porous sandstones. Petrography controls the geomechanical properties of the reservoir, crucial for predicting a reservoir's response to fluid extraction and understanding its lateral variability. This study focuses on the Groningen gas field in the Netherlands, addressing its compaction-induced surface subsidence and seismic events resulting from gas depletion. Core samples were examined to delineate spatial petrographic trends and develop a microscale model of the Groningen gas field. Optical microscopy (OM) and scanning electron microscopy (SEM) were used to determine mineralogical compositions, textural relationships and diagenetic processes. Predominantly constituted of sublitharenites, the sandstones exhibit dolomite and quartz cement as primary authigenic cements. Variations in clay types—kaolinite, illite, and chlorite—were observed, influencing localized pore-filling cementation processes.  Across the field, mineral relations revealed notable trends: depth-related feldspar decrease, correlation between kaolinite and feldspar abundance, and elevated chlorite content towards the northern sector together with the presence of an early quartz cementation phase, which is also observed within aquifer cores. The dissolution of feldspar potentially impacts the structural integrity of the sandstones, while authigenic mineralization appears intricately linked to depositional facies and localized fault-related fluid movements. The timing and extent of these diagenetic processes emerged as pivotal factors dictating sandstone stability within the reservoir. This comprehensive analysis enhances our understanding of Groningen's reservoir heterogeneity, offering critical insights to predict and manage subsurface responses to extraction-induced pressure changes. By providing predictive models, this study facilitates the evaluation of reservoir behavior and aids in mitigating risks associated with compaction-induced subsidence and seismicity.

How to cite: Mulder, S., Bublik, D., and Miocic, J.: Petrographic Heterogeneity and Sandstone Evolution in the Groningen Gas Field: Implications for Reservoir Geomechanics, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-2471, https://doi.org/10.5194/egusphere-egu24-2471, 2024.

Multiple episodes of brittle deformation tend to increase the structural complexity of fault zones. This commonly results in the development of juxtaposed and non-coeval distinct Brittle Structural Facies (BSF) formed at different times, depths, and temperature. Indeed, these BSFs are characterized by an irregular distribution of inherited, syn-kinematic, and post-kinematic minerals, whose study provides useful information about the temperature conditions of (de)formation, and the origin of fluids circulating within faults. We combined X-ray diffraction (XRD) analyses and polytype determinations of whole-rock and several grain-size fractions (6-10 μm, 2-6 μm, 0.4-2 μm, 0.1-0.4 μm, and <0.1 μm), with H-isotope data of 76 fault rocks and 6 protoliths from two structurally and well-characterized fault zones with different kinematics: the Carboneras strike-slip fault zone (Betic Cordilleras, SE Spain) and the Kornos-Aghios Ioannis extensional fault zone on the Lemnos Island (North Aegean Trough, Greece). Mineralogical and geochemical data allowed us to (1) reconstruct the distributions of syn/post-kinematic minerals in distinct BSFs, (2) constrain their formation temperature, and (3) unravel the origin of fluids involved during faulting. In the Carboneras fault rocks, we distinguished a protolithic mineralogical assemblage consisting of quartz, carbonates, K/Na-micas (2M₁ polytype), chlorite, kaolinite, and Fe/Ti-oxides, a syn-kinematic assemblage composed of mixed layers chlorite-smectite and illite-smectite (1M_d polytype), and a post-kinematic assemblage made up of smectite, chlorides, and sulphates. The coexistence of randomly (R0), short-range (R1), and long-range (R3) ordered illite-smectite in different BSFs indicates contrasting formation temperatures. Bulk samples generally display δ²H values (V-SMOW) between -15‰ and -60‰ (with a few exceptions at -90‰), while their respective <2μm fractions show δ²H values between -10‰ and -60‰. The combination of mineralogical and geochemical data from the Carboneras fault zone depicts a complex history of multiple brittle events occurring at different temperature conditions, wherein parental fluids of mostly meteoric origin infiltrated into the fault zone and interacted with the host rocks at various degrees and depths. In the fault rocks of Lemnos Island, we identified a host-rock mineralogical assemblage composed of quartz, feldspars, carbonates, K-mica (2M₁ polytype), chlorite, R0 illite-smectite, anatase, and Fe-oxides and hydroxides, a syn-kinematic assemblage made up of R3 illite-smectite (1M_d polytype), and kaolinite, and a post-kinematic assemblage characterized by halite and gypsum. Bulk samples display δ²H values (V-SMOW) between -70‰ and -140‰, while their respective <2μm fractions show δ²H values between -60‰ and -85‰. Such results indicate that Kornos-Aghios Ioannis fault rocks formed during a deformation event with predominantly hydrothermal fluids circulating into the fault zone (>160°C). This multidisciplinary approach represents an innovative point of view for studying fluid circulation, mineral crystallization, and temperature evolution in complex fault zones, and it can be applied to both orogen scale faults and smaller fault systems.

How to cite: Moretto, V., Del Sole, L., Curzi, M., Dallai, L., Vignaroli, G., Viola, G., Balsamo, F., Berio, L. R., Grathoff, G., Warr, L. N., and Aldega, L.: Tracing the origin of parental fluids and paleo-temperature distribution in fault zones: case studies from the Carboneras Fault (Spain) and Lemnos Island (Greece), EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-4381, https://doi.org/10.5194/egusphere-egu24-4381, 2024.

Deformation may exert a positive or negative impact on rock strength, stiffness and porosity depending on the initial properties of the rock. In porous rocks, deformation bands (DBs) reduce both porosity and permeability, and potentially increase rock strength and stiffness. In low porosity (tight) rocks, deformation occurs via shear and opening-mode fractures, thus enhancing permeability and deteriorating strength and stiffness.

In this study we describe field observations and laboratory analyses to characterize sub-vertical deformation bands in the porous carbonates of the Calcarenite di Gravina Fm. exposed in Matera and Gravina in Puglia, Southern Italy. Matera and Gravina in Puglia are located at the boundary between the Apulian foreland and the foredeep of the Southern Apennines thrust belt. Both study areas consist of an asymmetrical horst structure involving the Cretaceous (Senonian) tight carbonates of the Apulian platform (Calcare di Altamura Fm.) unconformably overlain by Plio-Pleistocene shallow-marine coarse-grained lithic sandstone and grainstone (Calcarenite di Gravina Fm.). The Calcarenite di Gravina Fm. is dominated by pervasive DBs organized into 2 main sets dipping at high angle and striking N-S and NNE-SSW, except in some limited areas where the DBs form a complex network with the presence of a secondary set, striking NW-SE, showing mutual crosscutting relationships. The DBs have a positive relief, due to their relatively higher resistance to erosion, and appear whitish, tabular with some slight undulations, 10’s of meters long in map view, continuous and steeply inclined (from 75° to 90°). These structures have thicknesses from few mm up to a few cm depending on the facies of the calcarenite. At the outcrop scale the DBs don’t seem to accommodate significant shear offsets. A total of 53 oriented samples were collected for thin sectioning, petrophysical, and microstructural analysis. We acquired porosity measurements using a Hg-intrusion porosimeter, which showed that the DBs have an average porosity of 11%, while for the host rock is 38%. 29 Blue-impregnated thin sections of host rocks and DBs were analysed by standard microscopy, cathodoluminescence (CL), and a scanning electron microscope (SEM). The analysis of the CL images shows that the clast size distributions are similar in the DBs and host rock, and they are not crushed or fractured. Furthermore, while in the host rock the clasts are randomly oriented, in the DBs they are iso-oriented with the long axis describing low angles to the band.

Our work shows that these DBs are characterized by the absence of grain size reduction with a strong preferred orientation of long axes without significant fragmentation. These data may indicate that DBs formed in high porosity conditions, at very shallow depths. Moreover, given the strong difference between the petrophysical properties of the host rock and DBs and the apparent abundance and continuity of the DBs, they would provide a significant anisotropy in the flow pathways of the calcarenite.

How to cite: Freda, G., Mittempergher, S., Balsamo, F., Pizzati, M., Di Cuia, R., and Ricciato, A.: Microstructural characterization of deformation bands in shallow porous carbonates of Apulian Platform, Southern Italy, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-5693, https://doi.org/10.5194/egusphere-egu24-5693, 2024.

Here we present the results of surface and subsurface (caves) geological mapping in the Lefka Ori Massif. The structural field data collection and the stratigraphy analysis were supported by the IGCP-715 project in collaboration of the Speleological Association of Crete (SPOK) and the Gourgouthakas 2023 exploration team. The overall structure of the Lefka Ori Massif consists of a roughly E-W trending open anticline. Regional uplift is related to compressional and extensional events during the Late Oligocene to present. The primary lithologies of the Lefka Ori Massif are comprised of the Plattenkalk Unit and the Trypali Unit. The Plattenkalk Unit contains (1) black brecciated karstified marbles/dolostones with interbedded stromatolites (Triassic), (2) white to gray marbles/dolostones (Triassic-Lias), and (3) thin bedded marbles (Dogger-Oligocene). Black dolostones and brecciated marbles (with stromatolite fragments) are building up the overlying Trypali Unit (Late Triassic-Lias).

The results from the structural analysis in the Lefka Ori Massif highlight 3 sets of faults, which are of (A) E-W striking thrust faults, (B) dextral NNW-SSE strike-slip faults and (C) E-W striking steep dip-slip faults. A regional thrust (hereafter "Pachnes Thrust") has a vertical displacement of 1000 m and duplicates the Plattenkalk Unit stratigraphy (sections 1 and 2) in the Lefka Ori Massif. The hanging wall consists of Plattenkalk units 1 and 2, while the footwall comprises the Plattenkalk unit 3. The Pachnes Thrust has an apparent horizontal displacement of approximately 12 km. The cross-cutting relationship of the thrust with the strike-slip faults indicate that the faults are coeval. The steep dip-slip faults are still of an unknown relative age but certainly most of them appear as normal, younger faults. The Pachnes Thrust is folded around a roughly E-W trending open fold axis. The Late Oligocene green and red, calcschist sediments in the footwall of the Pachnes Thrust indicates the maximum age of thrusting.

The third deepest cave in Lefka Ori Massif (Sternes Cave, depth -616 m) contains a system of galleries reaching a total length of 5.6 km. This karst system of subhorizontal galleries is well explained by the Pachnes Thrust as it parallels the thrust. Likewise, the Pachnes Thrust has been delineated in the deepest cave (Gourgouthakas) at a depth of -700 m. In addition, the accepted hydrology in the Lefka Ori Massif is defining an upper fast flowing reservoir and a lower slow flowing reservoir. The physical model for these two reservoirs is explained by the hanging wall 1 and 2 lithologies (fast flowing reservoir) and the footwall lithology 3 (slow flowing reservoir) of the Pachnes Thrust.

The Pachnes Thrust in Lefka Ori Massif probably correlates with the recorded southward thrusting and nappe stacking due to convergence between the African and Eurasian plates (Oligocene-Early Miocene). The deformed Pachnes Thrust plane is coeval with the southward thrusting, due to the uplift and exhumation phases during extension described in the literature. The present findings are offering new evidence of the tectonic evolution and the exhumation of the high pressure metamorphosed rocks in the Hellenic Subduction Channel.

How to cite: Moraetis, D., Leontaritis, A., Scharf, A., Fassoulas, C., Zacharias, S., Pavlopoulos, K., Mattern, F., Alnaqbi, A., Yu, X., Pennos, C., Adamopoulos, K., Hamdan, H., and Nikolaidis, N. P.: A new geological model for the Lefka Ori Massif in the accretionary wedge of the Hellenic Subduction Zone. , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9013, https://doi.org/10.5194/egusphere-egu24-9013, 2024.

Understanding the origin and distribution of damage within carbonate-hosted fault zones is crucial, yet it remains a complex challenge, which hampers the overall assessment of their mechanical and hydraulic structure. In carbonate-hosted fault zones, shattered to intensely brecciated non-cohesive rocks have been reported. Although their origin has been related to the propagation of multiple seismic ruptures, great uncertainties persist regarding their interpretation and distribution.

The NW-SE striking, approximately 15 km long Roccapreturo Fault, in the central Apennines of Italy, is an intriguing case study where non-cohesive fault rock domains occur within its footwall damage zone. These domains elongate in a NE-SW direction for ~200 meters from the main slip surface.

We employed a multiscale approach to better understand the distribution and origin of the non-cohesive fault rocks. The fault geometry and throw distribution along the main fault segments were characterized through fault-perpendicular geological cross-sections. Virtual outcrop models of key exposures, located in and around an abandoned quarry, were constructed using Structure from Motion-Multiview Stereo photogrammetry. These models utilized photos taken with a Mavic Mini 2 drone. The interpretation of virtual outcrop models, combined with classical fieldwork, allowed us to map the damage and minor fault strands.

The Roccapreturo Fault displaces Cretaceous rocks originally deposited in various depositional environments. Along its strike, from NW to SE, the fault intersects rocks from internal or restricted carbonate platform, margin, and proximal slope to basin depositional environments. Notably, non-cohesive fault rocks are exposed between the margin and proximal slope rocks. This area coincides with the maximum throw of the fault, which is ca. 600 meters, and with the intersection with a system of pre-existing NE-SW-striking steeply dipping faults.

At the outcrop scale, faults exhibit two preferred orientations, parallel and perpendicular to the main slip surfaces of the Roccapreturo Fault, respectively. The former ones show predominant dip-slip kinematics, while the latter ones show both dip-slip and strike-slip kinematics.

We interpret the distribution of non-cohesive fault rocks along the Roccapreturo Fault as influenced by its intersection with the NE-SW fault system, where most of the slip accumulated. Accordingly, the pre-existing NE-SW faults accommodated transtensional slip during latest extensional deformation and coeval rock exhumation from depth. The transition of the Cretaceous depositional environments, which was accommodated by the NE-SW-striking faults, therefore highlights the pivotal role of pre-existing anisotropies in dictating the distribution of damage, particularly of non-cohesive fault rocks, in carbonate hosted faults.

How to cite: Mercuri, M., Agosta, F., Fondriest, M., Smeraglia, L., Billi, A., Tavani, S., and Carminati, E.: Influence of pre-existing faults on damage distribution in carbonate fault zones: the case study of the Roccapreturo Fault, central Apennines, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-10510, https://doi.org/10.5194/egusphere-egu24-10510, 2024.

Networks of deformation bands in porous granular rocks represents potential low-permeable baffles to fluid flow in subsurface reservoirs. However, little work addresses the network properties of such networks, like spatial intensity, network connectivity and network geometry. Motivated by this, we here present an investigation of two-dimensional, horizontal exposures of deformation band networks within the Jurassic Entrada Sandstone in the San Rafael Desert (Utah). We analyse the geometry and topology (i.e. a network represented as nodes and branches) of the studied networks to: 1) characterise deformation band orientation, connectivity and areal intensity; 2) assess spatial topological variability; 3) elucidate large scale variation across the study area; 4) evaluate effective network permeabilities. Effective deformation band network permeability is calculated by incorporating a topological measure of network connectivity into the permeability calculations. Deformation band networks show distinct topological signatures, typically being dominated by Y-nodes, and IC- and CC-branches. Depending on the orientation of deformation bands and numbers of different sets of deformation bands within each studied network, both topology and areal intensity vary. Low proportion of isolated II-branches reflects the evolution of deformation bands through bifurcation and abutment, creating Y-nodes, to form interconnected networks. We document great spatial variability I network connectivity and topology within individual networks. Similarly, the effective permeability within well-connected (parts of) the studied deformation band networks (>1.5 connections per branch) significantly reduce effective permeabilities, whereas areas within the networks with low connectivity offer higher-permeable pathways for tortuous fluid flow.

How to cite: Heggernes, H., Rotevatn, A., Demurtas, M., Nixon, C. W., and Fossen, H.: Spatial variability in topology, connectivity and permeability within deformation band networks, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17174, https://doi.org/10.5194/egusphere-egu24-17174, 2024.