Fault Length Research Articles

Evolution of rift faulting in incipient, magma-poor divergent plate boundaries: New insights from the Okavango-Makgadikgadi Rift Zone, Botswana

Continental rifts are thought to transition from a power-law fault length distribution during the juvenile stages of extension, to an exponential distribution during break-up and oceanic spreading. However, fault scaling relationships have rarely been quantified in natural incipient rifts, particularly in the absence of magmatic influence. We address this knowledge-gap in the incipient Okavango-Makgadikgadi Rift Zone, Northern Botswana, consisting of the amagmatic Okavango Rift Basin (ORB) and Makgadikgadi Rift Basin (MRB). We utilize high-resolution satellite topography data (12.5 m TanDEM-X and 30 m SRTM) and aeromagnetic data to image and map faults across the rift zones and generate a new robust fault map for the region. Further, we analyze the length-frequency distribution of faults using a two-sample Kolmogorov-Smirnov test. The results show that the Makgadikgadi Rift Basin exhibits an exponential distribution associated with a nucleating rift as it exhibits diffuse faulting and lacks a border-fault, whereas the Okavango Rift Basin exhibits a power-law distribution that is consistent with its relatively more evolved rift structure with established border faults. Thus, we propose that continental divergent plate boundaries commonly nucleate with an initial exponential distribution of fault lengths which subsequently transition into a power-law distribution as the rift evolves into its stretching phase, and may again transition into an exponential distribution as break-up initializes.

The influence of the strength of pre-existing weak zones on rift geometry and strain localization

Continental rifts normally initiate within previously deformed lithosphere and thus their evolution and architecture can be largely controlled by inherited weak zones in the pre-rift crust. Here, we quantify the role of the strength and obliquity of pre-existing crustal-scale weak zones in the evolution of continental rift systems. We use a 3D numerical geodynamic model to assess strain localization, associated fault development, and rift segmentation during the early stages of tectonic extension. We find that both the strength and obliquity of the weak zones significantly influence the patterns of strain localization. A pre-existing very weak zone with low obliquity can promote the development of continuous and long-lived border faults parallel to the rift axis. Conversely, a comparatively strong weak zone with high obliquity leads to a staggered en-echelon rift geometry that lacks rectilinear laterally persistent strain localization. Furthermore, we find that rift obliquity and weak zone strength may modulate rift fault length, throw, and azimuth. These results provide new and compelling insights into the structure and evolution of natural active rifts that develop within orogenic basement terranes.

Tectonic geomorphology and paleoseismology of the Angelochori fault segment of the Anthemountas extensional detachment fault, Central Macedonia, Greece: Paleoseismic evidence from the 1677 CE earthquake

This study presents the results of a geomorphological and paleoseismological investigation on the Angelochori normal fault. This fault controls the southern boundary of the Anthemountas basin and constitutes the western segment of the Anthemountas extensional detachment fault (AEDF), a ∼48 km-long, E-W-trending, north-dipping low-angle normal fault in Central Macedonia, Greece. To assess the relative tectonic activity of the Angelochori fault and to determine the number, timing, and displacement of surface-faulting earthquakes, we conducted detailed field mapping and paleoseismological investigations. Our approach involved utilizing geomorphic indices, analyzing stratigraphic and structural relations, considering luminescence ages, and employing OxCal modeling techniques. The acquired data reveal, for the first time, the Middle-Late Holocene activity of the Anghelochori fault, providing constraints on the age of three strong (M > 6) surface-rupturing earthquakes over the last ∼7.7 ka. The youngest of these events aligns with the historical earthquake of the Anthemountas basin in 1677 CE. The observed minimum displacement per event ranges from 0.48 m to 0.65 m, consistent with a fault length estimated at about 24–28 km, suggesting an earthquake magnitude on the order of 6.6. The cumulative slip rate over the recorded period averages 0.21 ± 0.02 mm/a, with an average open recurrence interval of 2.6 ± 0.2 ka. Nonetheless, the displacement pattern displays alternating intervals of rapid slip at 0.60 ± 0.16 mm/a and slower slip at 0.08 ± 0.01 mm/a. Classified as Class 1, the Angelochori fault signifies an active tectonic structure posing a significant seismic hazard for the Thessaloniki metropolitan area.

Source parameters and aftershock pattern of the October 7, 2021, M5.9 Harnai earthquake, Pakistan

On October 7, 2021, a magnitude 5.9 earthquake struck the Harnai (Baluchistan) region of Pakistan, causing several fatalities and injuries within the epicentral area. First-order tectonic deformation in this region is caused by the convergence of the Indian Plate with respect to the Eurasian Plate. The Katwaz Block hinders the motion of the Indian Plate, resulting in the formation of strike-slip faults. In this study, the P-wave first-motion polarity technique was used to determine the mainshock faulting style. Cyclic scanning of the polarity solutions was applied to determine the most suitable focal mechanism solution among the available solutions generated by the FOCMEC (focal mechanism) software. The nodal planes correspond to different faulting styles (i.e., thrust and strike-slip faulting). A nodal plane oriented in the NW-SE direction corresponded to a strike-slip mechanism, which was considered to be the fault plane. Tectonically, this earthquake was associated with the Harnai-Karahi strike-slip fault zone owing to the fault strike and direction of slip. The apparent stress drop, fault length, and moment magnitude of the Harnai earthquake were 35.4 bar, 6.1 km, and 5.9, respectively. A lower b-value for the Gutenberg-Richter law was observed prior to the earthquake. Higher α- than b-values (α >b) indicate that this earthquake was governed by large events as opposed to small-magnitude events. The Harnai sequence had a decay exponent close to unity, lasted for 145 days, and produced few aftershocks. The study will help the future hazard mitigation in the region.

Authentic fault models and dispersive tsunami simulations for outer-rise normal earthquakes in the southern Kuril Trench

The southern Kuril Trench is one of the most seismically active regions in the world. In this study, marine surveys and observations were performed to construct fault models for possible outer-rise earthquakes. Seismic and seafloor bathymetric surveys indicated that the dip angle of the outer-rise fault was approximately 50°–80°, with a strike that was slightly oblique to the axis of the Kuril Trench. The maximum fault length was estimated to be ~ 260 km. Based on these findings, we proposed 17 fault models, with moment magnitudes ranging from 7.2 to 8.4. To numerically simulate tsunami, we solved two-dimensional dispersive wave and three-dimensional Euler equations using the outer-rise fault models. The results of both simulations yielded identical predictions for tsunami with short-wavelength components, resulting in significant dispersive deformations in the open ocean. We also found that tsunami generated by outer-rise earthquakes were affected by refraction and diffraction because of the source location beyond the trench axis. These findings can improve future predictions of tsunami hazards.Graphical

On the Interplay Between Distributed Bulk Plasticity and Local Fault Slip in Evolving Fault Zone Complexity

AbstractWe numerically investigate the role of plastic strain accumulation on the mechanical response of a planar strike‐slip fault. Our models show that fault‐zone strength significantly impacts the ensuing sequence of earthquakes. Weaker fault zones accumulating more plastic strain promote more complexity in the seismicity pattern through aperiodic earthquake occurrences and intermittent episodes of rupture and arrest. However, if the fault zone strength is high enough, the overall earthquake sequence is characterized by periodic fault‐spanning events. We find that both the fault normal stress and the fault geometric profile evolve throughout the earthquake sequence, suggesting a self‐roughening mechanism. Despite the significant impact of plasticity on the fault response, the width of the plastically deforming region in the fault zone is small compared to the fault length. Our results suggest a rich behavior in dynamically evolving fault zones and support the need for further high‐resolution studies of the highly non‐linear near‐fault region.

Tectono‐sedimentary evolution of high‐displacement crustal‐scale normal faults and basement highs on rifted margins: Klakk Fault Complex and Frøya High, Mid‐Norwegian Margin

AbstractCrustal‐scale high‐displacement (>10 km) normal faults are not captured in existing tectono‐sedimentary models of rift basins. We used 2D and 3D seismic reflection and well data to perform a structural and source‐to‐sink analysis of the southern part of the Klakk Fault Complex and the western part of the Vingleia Fault Complex, Mid‐Norwegian rifted margin. The north–south trending Klakk Fault Complex has a zig‐zag to sinuous plan‐view geometry, forming a series of structural recesses and salients along strike. In cross‐section, the fault complex has a listric to convex‐up or low‐angle planar geometry with displacements above 20 km. This fault complex exhumed basement highs, the Frøya High and Sklinna Ridge, in its footwall and created a series of supradetachment basins, for example, the Rås Basin, in its hanging wall. In contrast, the northeast‐southwest trending Vingleia Fault Complex has a zig‐zag geometry in plan view and planar to listric geometry in cross‐section and displacement of up to 5 km. This fault has the Frøya High in its footwall and the southern Halten Terrace in its hanging wall. Restoration of selected structural cross‐sections shows a prominent fault‐parallel ridge, up to 15 km east of the Klakk Fault Complex interpreted as a palaeodrainage divide. This divide separates steep drainages developed along the west‐dipping footwall scarp to the Klakk Fault Complex, from broader, gentler east‐dipping drainages up to ca. 10 km long developed on a back‐tilted dip slopes along the eastern side of the Frøya High and Sklinna Ridge. Progressive headward erosion of active flank catchments was enhanced around topographically elevated structural salients to the point of capturing previous dip‐slope‐directed drainages during the earliest Cretaceous. A network of submarine canyons develop down‐dip of the drainage catchments along the Klakk Fault Complex scarp, whose geometries and length are controlled by their location with respect to the structural salients or recesses, and the presence of fault terraces. The middle Jurassic‐earliest Cretaceous synrift deposits form two seismic sequences that are filled with five distinctive seismic facies that record the evolution from a linked normal fault during rift climax to a high‐displacement stage. During the high displacement stage, exhumed local continental core complexes formed structural salients, separated along strike by structural recesses at the heads of supradetachment basins. Key elements of the high‐displacement fault stage include (i) the development of structural salients at sites of rift climax displacement maxima, (ii) development of supradetachment basins in rift climax displacement minima and (iii) migration of major depocentres away from the centre of rift climax fault segments. We synthesise these observations into a generic tectono‐sedimentary model for high‐displacement faults.

Geomorphological traits of landscapes in continental rifts—From fault‐elastic rebound to sedimentary sinks

AbstractWe analyse 498 faults identified in satellite imagery and interpret the height and width of associated footwall ranges with respect to co‐seismic elastic rebound from tectonic and erosional unloading. The dynamics of footwall uplift link uplands to catchment patterns and interrelated hanging wall sedimentary fans. Height–length relations of some catchments and associated alluvial fans scale linearly whereas others, such as fault‐slope catchments and related down‐fault fans (building out from faults) show a significant scatter without an obvious trend. Perched basins abandoned in the footwalls of younger faults offer catchment‐fan height–length relations like watergap and dipslope‐related fans and, besides, hint at reduction of dip angle due to rollback of larger faults before abandonment. Analysis of the width‐to‐height ratio (W/h) of footwall ranges offer a robust linear statistical trend, h = 0.06 W and is identical between datasets. This trend is valid for both arid and tropical rifts, the latter offering smaller rebounds. Contributions of elastic rebound on fault throw in our data are simplistically considered through comparison to global trends on fault length versus throw. This allows consideration around maximum throw (Tmax) linked to the maximum height of footwall ranges (h) and to their width (W) above the reference level. Basic calculations indicate that co‐seismic rebound contributes from <1% to 17% of extensional fault throw. Width‐to‐height ratios for large faults (L > c. 50 km) show less spread than smaller faults. Such large faults expectedly dissect the brittle crust, indicating that these large faults which root in the ductile–brittle transition approach a balanced, steady‐state kinematic pattern. We speculate that significant crustal thinning associated with these large faults triggers the onset of isostatic adjustments that drive fault rotation, instigating fault abandonment and disconnected perched basins.

Friction models of one-dimensional earthquakes

A number of fundamental characteristics of earthquake physics are reproduced by simple stick-slip friction models. These include creep motion, stick/slip periodic earthquakes, reduced frequency of occurrence for more energetic quakes, ground velocity proportional to slip duration, and ground displacement proportional to fault length. By confining attention to one dimension, these results are obtained analytically. Though many earthquake characteristics cannot be reproduced in a one-dimensional model, our simple approach nevertheless has pedagogical value in providing insight into some of the complicated dynamics of earthquakes.

The N-S direction strike-slip activities in the Pamir hinterland under oblique convergence: the 2015 and 2023 earthquakes

SUMMARY The prominent Pamir plateau holds considerable significance in comprehending the processes of Asian continental collisional orogeny. However, due to harsh natural conditions and low seismic activity within the Pamir hinterland, our understanding of this region remains deficient. Recent major events and the accumulation of geodetic observations present a rare opportunity for us to get insights into the tectonic activities and orogenic processes occurring in this region. First, employing Sentinel-1 and Advanced Land Observation Satellite (ALOS)-2 Synthetic Aperture Radar (SAR) images, we acquire coseismic displacements associated with the most recent earthquakes in 2015 and 2023. Subsequently, we conduct the source models inversion with the constraints of surface displacements based on a finite-fault model. Our results reveal displacements ranging from −0.8 to 0.8 m for the 2015 Mw 7.2 Tajik earthquake and −0.25 to 0.25 m for the 2023 Mw 6.9 Murghob event, respectively. The optimal three-segment model for the 2015 event ruptured a fault length of 89 km with a surface rupture extending 59 km along the Sarez–Karakul fault (SKF), characterized predominantly by left-lateral strike-slip motion, with a maximum slip of 3.5 m. Meanwhile, our preferred uniform slip model suggests that the 2023 event ruptured an unmapped fault in the southern Pamir region with a strike angle of 31° and a dip angle of 76.8°. The distributed slip model indicates that the 2023 event ruptured a fault length of 32 km, resulting in an 8 km surface rupture. This event is characterized by left-lateral strike slip, with a peak slip of 2.2 m. Secondly, the Coulomb stress calculations demonstrate that the 2023 event was impeded by the 2015 event. Finally, interseismic Global Positioning System data revel a relative motion of 3.4–5.7 mm yr−1 in the N-S component and 3.2–3.8 mm yr−1 in the E-W component along the SKF in the Pamir hinterland, respectively. These N-S direction strike-slip activities and slip behaviours support an ongoing strong shear and extension in the Pamir regime, which is a response to the oblique convergence between the Indian and Eurasian plates.

Quantifying fault interpretation uncertainties and their impact on fault seal and seismic hazard analysis

Fault-horizon cut-off data extracted from seismic reflection datasets are used to study normal fault geometry, displacement distribution, and growth history. We assess the influence of three seismic interpretation factors (repeatability, measurement obliquity, and fault cut-off type) on fault parameter uncertainty. Two repeat interpretations resulted in mean differences of 5–15% for throw, 11–42% for heave, 9–31% for displacement, and 7–27% for dip across faults. Measurement obliquity, where faults are interpreted using non-perpendicular transects to fault strike, show increasing uncertainty with increasing obliquity. Uncertainty in throw is 14–24% at obliquities >20° and 6–13% where obliquities <20°. Continuous cut-offs, including non-discrete deformation, generally exhibit greater uncertainties compared to discontinuous (discrete) cut-offs. We consider the effect of interpretation factors on fault parameters used in seismic hazard assessment (SHA) and fault seal, using the established Shale Gouge Ratio (SGR). Even modest measurement obliquities and repeatability errors can affect inputs for SHA, causing large differences in throw- or slip-rate and inferred fault length. Measurement obliquity and repeatability have a variable impact on SGR calculations, highlighting the additional importance of sedimentary layer thickness and distribution. Our findings raise questions about the optimum workflow used to interpret faults and how uncertainties in fault interpretation are constrained and reported.

New Empirical Source-Scaling Laws for Crustal Earthquakes Incorporating Fault Dip and Seismogenic-Thickness Effects

Abstract New global source-scaling relations for the aspect ratio and rupture area for crustal earthquakes that include the width-limited effect and a possible free-surface effect are derived using a global dataset of finite-fault rupture models. In contrast to the commonly used scaling relations between moment magnitude (M), fault length (L), width (W), and area, we built self-consistent scaling relations by relating M to the aspect ratio (L/W) and to the fault area to model the change in the aspect ratio once the rupture width reaches the down-dip width limit of the fault. The width-limited effect of large-magnitude earthquakes depends on the fault dip and a regional term for the seismogenic thickness. The magnitude scaling of the aspect ratio includes a break in the magnitude scaling that is dip angle dependent. This dip angle-dependent magnitude scaling in the magnitude–area relation is modeled by a trilinear relation incorporating a dip-related transition range. The effect of the free surface was observed using a normalized depth term and parameterizing the source by the depth of the top of the fault rupture; it is more apparent in the area scaling relation. The scaling differences are related to the fault geometry, not to the rake angle, as commonly assumed. Finally, the corresponding L and W scaling relations obtained by converting the area and aspect ratio models to L and W models not only show good agreement with the previous regional scaling laws on average but also provide better fault-specific application due to the inclusion of a fault-specific dip angle and seismogenic thickness.

The Main Lineament Trends of Central Iraq, Depending on the THD of Residual Gravity and Magnetic Maps Obtained from Four Upward Continuation Levels

The upward continuation was used to obtain central Iraq's residual gravity and RTP magnetic maps at elevations 2, 12, 16, and 22 km. These maps are processed by the Total Horizontal Derivative (THD) technique to detect the main lineaments that represent the boundary of subsurface sources or faults. The obtained lineaments trends represented on the rose diagram and the main trends in the study area were found. The lineaments (faults) orders from the main trend descending, from gravity data, are N55W, N45W, N35W, N-S, and N65W, respectively, for the four considered upward elevations. The length of the faults obtained from gravity data ranges from 40-250km, while the length of faults obtained from magnetic data ranges from 20-100km. The lineament from the main trend descending, from magnetic data, are NS, N55W, N35W, N45W, N45E, N35E, N65W and N55E, respectively. The lineaments obtained from magnetic data indicate that the main trend in shallow sources is NW-SE, while for deep sources; the main trend is N-S. The NW-SE fault trends obtained from gravity data are related to the Najd faults system, while the N-S fault trends are related to the Nabitah system. Generally, it is indicated from the analysis of lineaments of gravity and magnetic data that the NW-SE is the main trend of shallow structures in most areas of central Iraq.

Analysis of Influencing Factors of Slippage and the Dynamic Process of Fault Slip Caused by Multi-Stage Fracturing

Casing deformation is evident during the development of shale oil and gas wells in the Sichuan and Junggar Basins in China. Their casing deformation characteristics, distribution law of deformation points, and main controlling factors were analyzed. According to the analysis results, shear is the main cause of casing deformation of shale oil and gas wells in the Sichuan and Junggar Basins in China and has the characteristics of “a dense heel end and a sparse toe end”. Faults account for 75% of casing deformation points, and fault slip caused by multi-stage fracturing is the primary factor responsible. The calculation model for fault slip that takes into account fracturing fluid invasion was established, and the dynamic variation law of fault slip was clarified: the fracturing fluid intruded into the fault, the relative dislocation of the damaged fault was caused by gravity, and the fault slippage was caused by the increase in fault activation length. This resulted in a linear increase in fault slippage, and the slippage reached its maximum when the fracturing fluid completely penetrated the fault and reached the fault boundary. The slip amount has a positive correlation with the fault length and the in situ stress difference; it increases first and then decreases with the increase in the fault dip angle. The slip amount reaches its maximum when the fault dip angle reaches 45°.

Stability analysis of the rock zone between the tunnel face and the fault fracture zone

Abstract Water and mud inrush is one of the main safety accidents that occur during tunnel construction in water-rich karst regions. Often, faulting occurs in front of the tunnel face, creating a conduit for water and inrush disasters can easily occur. Accurately predicting the safety distance between the tunnel face and the fault fracture zone allows for effectively avoiding water and mud inrush disasters during construction. First, an analytical model of the safety distance of water and mud inrush prevention is proposed, in which the rock zone between the tunnel face and the fault fracture zone is considered a thick rectangular plate with simple support on four sides. Subsequently, the proposed model is successfully verified through comparison with two existing models and engineering cases. Finally, the influence of main model parameters on the safety distance is further determined. This study shows that: (i) The safety distance increases with the increase in the cross-sectional height and width, and the burial depth of the tunnel; (ii) The safety distance increases with the increase in the effective gravity of the rock inside the fault fracture zone and the height of the groundwater table, and decreases in the dip angle of the fault; (iii) The safety distance increases with the increase in fault width, and the fault length has little influence on the safety distance.

Structural expression of the frontal thrust of an active fold-and-thrust belt: The Holocene 123-km-long Kur fault, Greater Caucasus, Azerbaijan

Here we present the main features of the frontal structure, known as Kur Fault, of the Plio-Quaternary Kura fold-and-thrust belt in the Greater Caucasus (Azerbaijan). The Kur Fault has been analysed thanks to geological-structural and geomorphological surveys of its whole length, integrated by a relocation of instrumental seismicity, data on historical seismicity, new focal mechanism solutions, and ambient vibration measurements across the fault trace. The in-depth study of the frontal structure can: i) provide insights into the shallow propagation of a regional reverse fault, ii) contribute to a better understanding of the earlier stage of development of a young continent-continent collision, and iii) have implications for seismic hazard assessment because the area is seismically active and hosts the most important infrastructure for energy production in the country. The results show that the fault deforms the surface for a total length of 123 km. The shallow expression is given by four main scarp segments, with a right-stepping arrangement, which have different structural significance; they are represented by an alternation of fault-propagation folds, folds with offset frontal limbs, and shallow faulting. Analyses of the age of deformed deposits and landforms suggest activity from Mid-Late Miocene times to the Holocene. The fault attitude and pure reverse kinematics are coherent with the Holocene and present-day state of stress, characterized by a N–S to NNE-SSW horizontal σ1, suggesting the capability for seismic reactivation. Calculation of potential Mw indicates values in the range 7.5–7.9 if we consider its entire fault length, 6.1–7.2 if we consider the single segments.

The NE-SW Sibari fault zone: A seismic hazard source in Ionian Northern Calabria (Italy)

A multidisciplinary approach including archaeological, geophysical, and geological/geomorphological surveys provided pieces of evidence that allowed us to identify the Sibari fault zone (SFZ) in Northern Calabria (Italy). The SFZ runs in a ∼ NE-SW direction for a length of ∼18 km from the Ionian coastline to Terranova da Sibari and has an oblique normal-dextral kinematics. The envelope of the SFZ is derived from several direct and indirect evidence resulting in subparallel and locally en-echelon fault traces over a maximum 500 m-wide band, running at different elevations across hills and flat lands. The SFZ was active since at least the Middle-Upper Pleistocene, producing faulting of alluvial deposits, marine terraces, drainage incisions, and the archaeological structures of Sybaris. Given the fault length and assuming a seismogenic behavior, the SFZ is a primary earthquake source possibly producing moderate to large earthquakes (M ≥ 6). We calculated the average slip rates along the SFZ based on the ages and on the accumulated displacements of offset streams and marine terraces. The estimates are of 0.05–0.18 mm/yr and 0.41–0.70 mm/yr for vertical and dextral slip, respectively. Based on both the measured (min. 30 cm) and the expected value (av. 40 cm) of lateral slip per event, we infer an average recurrence for surface faulting events on the SFZ of about 700–1000 yrs. The most recent surface faulting earthquake occurred on the fault is dated 1300–1100 yrs. ago, highlighting that the elapsed time approaches the estimated average recurrence. Considering these findings, the newly recognized SFZ should be included among the faults that contain a potential seismic hazard in this poorly known portion of the Ionian sector of northern Calabria.

Paleoseismology of the Sawtooth Fault and Implications for Fault Behavior in the Epicentral Region of the 2020 Mw 6.5 Stanley, Idaho, Earthquake

Abstract The 2020 moment magnitude (Mw) 6.5 Stanley, Idaho, earthquake raised questions about the history and extent of complex faulting in the northwestern Centennial Tectonic Belt (CTB) and its relation to the Sawtooth normal fault and Eocene Trans-Challis fault system (TCFS). To explore faulting in this area, we excavated a paleoseismic trench across the Sawtooth fault along the western margin of the CTB, and compared an early Holocene (9.1 ± 2.1 ka, 1σ) rupture at the site with lacustrine paleoseismic data and fault mapping in the 2020 epicentral region. We find: (1) a history of partial to full rupture of the Sawtooth fault (Mw 6.8–7.4), (2) that shorter ruptures (Mw≤6.9) are likely along distributed and discontinuous faults in the epicentral region, (3) that this complex system that hosted the 2020 earthquake is not directly linked to the Sawtooth fault, (4) that the northeast-trending TCFS likely plays a role in controlling fault length and rupture continuity for adjacent faults, and (5) that parts of the TCFS may facilitate displacement transfer between normal faults that accommodate crustal extension and rotation. Our results help unravel complex faulting in the CTB and imply that relict structures can help inform regional seismic hazard assessments.

Design displacement for lifelines at fault crossings: the code-based approach for Europe

The earthquake-resistant design of lifelines, such as pipelines, tunnels and bridges, is based on the reliable representation and estimation of the seismic loading. In the case of lifeline–fault crossings, the design fault displacement is typically derived from estimates based on fault dimensions via empirical fault scaling relations for a given “design” scenario event. This approach comes with an unknown level of safety because the fault productivity and the actual distribution of earthquake events are essentially disregarded. To overcome this challenge, a simplified approach is proposed by statistically analyzing the outcome of probabilistic fault displacement hazard analyses (PFDHAs). A selection of faults from the 2020 European Fault-Source Model is used to build the logic tree and to set the range of parameters considered in the PFDHAs. The methodology allows the (mostly conservative) approximation of the fault displacement corresponding to any given return period based on readily available data, namely fault productivity, fault mechanism, fault length, and lifeline crossing location on the fault. The proposed methodology has been proposed and adopted as an informative Annex in prEN 1998-4:2022.

Review of discrete fracture network characterization for geothermal energy extraction

Geothermal reservoirs are highly anisotropic and heterogeneous, and thus require a variety of structural geology, geomechanical, remote sensing, geophysical and hydraulic techniques to inform Discrete Fracture Network flow models. Following the Paris Agreement on reduction of carbon emissions, such reservoirs have received more attention and new techniques that support Discrete Fracture Network models were developed. A comprehensive review is therefore needed to merge innovative and traditional technical approaches into a coherent framework to enhance the extraction of geothermal energy from the deep subsurface. Traditionally, statistics extracted from structural scanlines and unmanned aerial vehicle surveys on analogues represent optimum ways to constrain the length of joints, bedding planes, and faults, thereby generating a model of the network of fractures. Combining borehole images with seismic attributes has also proven to be an excellent approach that supports the stochastic generation of Discrete Fracture Network models by detecting the orientation, density, and dominant trends of the fractures in the reservoirs. However, to move forward to flow modelling, computation of transmissivities from pumping tests, and the determination of hydraulically active fractures allow the computation of the hydraulic aperture in permeable sedimentary rocks. The latter parameter is fundamental to simulating flow in a network of discrete fractures. The mechanical aperture can also be estimated based on the characterization of geomechanical parameters (Poisson’s ratio, and Young’s modulus) in Hot Dry Rocks of igneous-metamorphic origin. Compared with previous review studies, this paper will be the first to describe all the geological and hydro-geophysical techniques that inform Discrete Fracture Network development in geothermal frameworks. We therefore envisage that this paper represents a useful and holistic guide for future projects on preparing DFN models.