Vertical Fault Research Articles

Benchmarking numerical simulation of induced fault slip with semi-analytical solutions

Pore pressure fluctuation in subsurface reservoirs and its resulting mechanical response can cause fault reactivation. Numerical simulation of such induced seismicity is important to develop reliable seismic hazard and risk assessments. However, modeling of fault reactivation is quite challenging, especially in the case of displaced faults, i.e., faults with non-zero offset. In this paper, we perform a systematic benchmarking study to validate two recently developed numerical methods for fault slip simulation. Reference solutions are based on a semi-analytical approach that makes use of inclusion theory and Cauchy-type singular integral equations. The two numerical methods both use finite volume discretizations, but they employ different approaches to represent faults. One of them uses a conformal discrete fault model (DFM) while the other employs an embedded (non-conformal) fault model. The semi-analytical test cases cover a vertical frictionless fault, and inclined displaced faults with constant friction and slip-weakening friction. It was found that both numerical methods accurately represent pre-slip stress fields caused by pore pressure changes. Moreover, they also successfully cope with a vertical frictionless fault. However, for the case with an inclined displaced fault with a constant friction coefficient, the embedded method can not converge for the post-slip phase, whereas the DFM successfully coped with both constant and slip-weakening friction coefficients. In its current implementation, the DFM is therefore the model of choice when accurate simulation of local faulted systems is required.

The 2021 and 2022 North Coast California Earthquake Sequences and Fault Complexity in the Vicinity of the Mendocino Triple Junction

ABSTRACT The Mendocino Triple Junction (MTJ), one of the most tectonically active and complex regions of California, has damaging earthquakes on the San Andreas and Mendocino faults, within the oceanic and subducting regions of the Gorda section of the Juan de Fuca plate, and within the overriding North American plate (NAP). Two recent earthquake sequences in the MTJ region, starting on 20 December 2021 and 20 December 2022, highlight the complex interactions of regional faults. We explore these sequences to better define the deep faults in the MTJ region, and their rupture modes. Our finite-source analysis shows the 2021 sequence began with two M ∼6.0 earthquakes separated by ∼11 s in time and 30 km in distance. The first earthquake occurred offshore on the Mendocino fault at a depth of 16.5 km. Its S waves triggered an “onshore” intraplate Gorda event at a depth of 27 km, which ruptured a vertical fault toward the northeast. Finite-source analysis of the mainshock of the 2022 sequence, M 6.4, indicates the rupture started offshore north of Cape Mendocino within the Gorda plate and propagated east-northeast, toward populated areas. Damage to towns and infrastructure was exacerbated by directivity and the sediment-filled valleys, as well as by a large aftershock (M 5.4) centered 20 km south-southeast of the mainshock rupture plane. The depths and mechanisms of the onshore 2021 and the 2022 earthquakes and their aftershock sequences indicate that they occurred on different strike-slip faults within the subducted portion of the Gorda plate. The faults active in these earthquakes are unrelated to mapped surface faults in the overriding NAP and are oblique to the tectonic trends seen at the surface. The 2021 and 2022 earthquakes are close to the boundary between two distinct regions of the Gorda plate, where offshore north–south horizontal compression transitions to east–west downslab tension.

Effects of Dip Angle on Rupture Propagation Along Branch Fault Systems

ABSTRACT An important consideration in assessing seismic hazards is determining what is likely to happen when an earthquake rupture encounters a geometric complexity such as a branch fault. Previous studies showed parameters such as branch angle, stress orientation, and stress heterogeneity as key factors in the self-determined rupture path on branch faults. However, most of these studies were conducted in 2D or 3D with perfectly vertical faults. Therefore, in this study, we investigate the effects of dipping angle on rupture propagation along a branch fault system. We construct 3D finite-element meshes where we vary the dip angles (nine geometries in total) of the main and secondary faults, the stressing angle (Ψ=20°, 40°, and 65°), and the hypocenter location with nucleation on both the main and secondary segments. We find that for Ψ=40°, a rupture on the main fault is most likely to propagate across the branch intersection when the secondary fault is dipping. In addition, for Ψ=65°, a rupture on the secondary fault is most likely to propagate to the main fault when the secondary fault is shallowly dipping. This is caused by a fast rupture speed on the secondary fault and the dynamic stress effect that develops with the interaction of the free surface and the dipping secondary fault. These results indicate that dip angle is an important parameter in the determination of rupture path on branch fault systems, with potentially significant impact for seismic hazard, and should be considered in future dynamic rupture modeling studies.

Influence of velocity pulse directivity on seismic response of cross-fault bridges

Pulse-type ground motions have special destructive effects on structures, and directivity is one of the fundamental characteristics of velocity pulses. In this paper, the impact of the velocity pulse directivity of seismic motions on both sides of the fault on the seismic response of fault-crossing bridges was quantified using the original records from the Chi-Chi earthquake in Taiwan. First, an automatic baseline correction method was applied to restore the permanent ground displacement, and the direction corresponding to the strongest pulse energy on both sides of the fault was identified based on continuous wavelet transformation. On this basis, the variability of the peak intensity of ground motions on the hanging wall and footwall in different directions, as well as the differences in pulse characteristics between the strongest pulse component and the two initial horizontal components, were analysed. Finally, a typical 5 × 30 m continuous girder bridge was considered as the research object to reveal the influence of the velocity pulse directivity on the seismic response of key parts on both sides of the fault. The results showed that the variability in the peak ground-motion intensity on the hanging wall of the fault was higher than that on the footwall. The displacement variability was up to 863 %, velocity was 262 %, and acceleration was 124 %. The pulse characteristics of the strongest pulse component were stronger than those of the parallel and vertical fault components, and its acceleration, velocity, and displacement response spectra were larger. The amplification effects of the velocity pulse directivity on the transverse bending moment, longitudinal bending moment, and torque of the pier bottom reached 1.3, 1.1, and 1.6, respectively. However, the amplification effect of the velocity pulse directivity on the relative displacement of the pier top, displacement of the main beam, and displacement of the bearing was more than 1.6 times, which was more significant than the internal force of the pier bottom.

3D structural modelling of the Kopili fault zone in North-East India: Seismotectonic analysis utilising focal mechanism solutions of small-to-moderate earthquakes

This work presents new findings from recent seismic data which have never been exploited so far for the weak-to-moderately energetic (ML 2.3–6.4) earthquakes and to depict a 3D scheme for the Kopili Fault (KF) Zone. For the first time, we figure out different fault kinematics in the upper crust such as reverse and strike-slip mechanisms located along the KF zone. Moment Tensor (MT) solutions are successfully obtained through full-waveform inversions for 60 seismic events. The stress inversions are estimated based on the MT results as well as previous findings from 16 small to moderate earthquakes (MW 3.6–6.3). The MT solutions demonstrate that strike-slip kinematics is the most prevalent while a few thrust/normal earthquake events are observed in the area. The seismic cross-sections display that these earthquakes took place at a depth of 40 km and are mostly generated on the eastern side of the fault. According to the 3D spatial and temporal distribution of relocated events, the KF zone is not a completely vertical fault but rather one with a slight eastward dip. Our stress inversion results also indicate that this fault is equivalent to a fault plane with a dip angle of around 80˚, which is compatible with the previously estimated dip angle (75˚). We interpret that the research region is under a transpressional stress regime, which is likely to be the reason for the current activity of the KF zone. We postulate that our findings can be useful to enhance seismic hazard assessment in the densely populated region, which is also covered by vital infrastructure.

Constraints on the Geometry of Peripheral Faults above Mafic Sills in the Tarim Basin, China: Kinematic and Mechanical Approaches

Host rock deformation associated with sill emplacement is used to constrain magma transfer and storage within the upper crust. In contrast to classic models suggesting that the host rock above mafic sills is dominated by elastic bending, recent studies show that bounding faults that limit the uplift area can occur at the peripheries of a mafic sill. However, the accurate dip of this type of fault, named peripheral faults here, is still not well constrained. Their origin is also controversial in some cases. In this study, kinematic modeling and limit analysis are performed to better constrain the structure and mechanical properties of the peripheral faults based on seismic interpretation of a mafic sill from the Tarim Basin, China. The trishear kinematic model successfully reproduces peripheral faulting and associated folding of the host rock by performing a displacement of 58 m on a vertical fault plane with a fault propagation (P) to fault slip (S) ratio of 2.5. The limit analysis also predicts vertical damage at the sill tip by sill inflation. These results suggest that the dip angle of the fault in the case study is 90°, which is more accurate than that from the seismic interpretation with an 88° inward dip. This value may vary in other cases as it depends on the sill geometry (such as diameter and inclination), thickness, depth, and mechanical properties of the host rock. The study supports that peripheral faulting and associated folding can occur at the tips of the mafic sill due to the vertical uplift of the host rock caused by sill inflation. It is also suggested that trishear kinematic modeling and limit analysis are effective methods for studying the geometry of peripheral faults.

An approach for fault-related monitoring variables selection based on dual-layer correlation networks

PurposeFault-related monitoring variables selection is a process of obtaining a subset of variables from the original set, which is of great significance for reducing information redundancy and improving the performance of the fault diagnosis models. This paper aims to propose a novel variables selection approach based on complex networks.Design/methodology/approachFirstly, a dual-layer correlation networks (DlCN) which consists of mechanism-oriented correlation sub-network (MoCSN) and data-oriented correlation sub-network (DoCSN) is constructed. Secondly, an algorithm for identifying critical fault-related monitoring variables based on dual correlations is introduced. In the algorithm, the topological attributes of the MoCSN and correlation threshold of the DoCSN are used successively.FindingsIn the experiments of vertical elevator fault diagnosis, the critical fault-related monitoring variables selected by the DlCN-based approach is more effective than the traditional approaches. It indicates that fusion mechanism-oriented correlation can enhance the comprehensiveness of variable correlation analysis. Moreover, the approach has been proved to be adaptable to different fault diagnosis models.Originality/valueIn the DlCN-based variables selection approach, the mechanism-oriented correlation and data-oriented correlation are comprehensively considered. It improves the precision of variables selection. Meanwhile, it is an unsupervised and model-agnostic approach which addresses the shortcomings of some conventional approaches that require data labels and have insufficient adaptability for fault diagnosis models.

Recognition of crustal extension in the Basin and Range Province: A history

The Basin and Range Province of the western United States has inspired or augmented many important geological ideas but is most famous for the concept of crustal extension. Acceptance of the concept was a slow process. Prior to the recognition of plate tectonics, extension was unexpected and hard to understand; when mapping of the region began in the 1860s, mountain ranges were assumed to be the products of crustal compression due to contraction of the Earth. The first mapping suggested that ranges are anticlinal folds, but the idea of uplift along vertical range-bounding faults arose shortly thereafter. There was recognition starting in the 1880s that at least some faults were normal-sense rather than vertical, which implied extension, but there were decades of dissent by geologists who still promoted the anticlinal-folding idea, thought that inclined faults were reverse faults, or believed that the basins and ranges were created by water or wind or glacial erosion rather than tectonics. Normal faulting gradually became widely accepted, but comprehension of its meaning was hindered by apparent thrust faults in parts of the province that suggested contractional deformation. These faults would later be interpreted as low-angle normal faults, and crustal extension became widely, though not universally, accepted in the 1960s as magnetic anomalies confirmed seafloor spreading. Roughly a century elapsed between the first geologic mapping and general acceptance of crustal extension. A parallel evolution of thought occurred in Africa as a consensus arose that rift valleys were products of extension. "Even those who have more sympathy with man’s endeavor than with the affairs of Nature may take an interest in the Science of Tectonics…. as a rule, the knowledge of the structure of a mountain chain comes as the reward of glorious struggle, both physical and mental." —Edward Battersby Bailey, Tectonic Essays (1935, p. 1)

A Numerical Study of Fault Reactivation Mechanisms in CO2 Storage.

Due to extensive industrialization and ongoing fossil fuel consumption, CO2 emissions have significantly contributed to climate change and rising greenhouse gas levels. The collection and storage of CO2 in subsurface geological formations have been proposed as a feasible alternative. The well-documented In Salah storage site is the subject of the chosen case study. For the numerical analysis, a two-dimensional finite element simulation of fault reactivation processes was conducted in the context of the CO2 storage. The goal of this research is to conduct a mechanistic numerical analysis of a typical CO2 storage condition. The selected analysis domain is 4 km (width) × 2 km (depth) and includes all important domains and formations (overburden, main caprock, lower caprock, and underburden) of the In Salah site. The simulation results indicate that the influence of fault reactivation under 32 MPa of base injection pressure results in peak vertical deformation of 0.044 m in the caprock and a vertical displacement magnitude of 0.015-0.020 m at the surface level (Z = 0 m). The derived vertical deformation findings at the surface level are in agreement with the data obtained from the in situ InSAR monitoring system in 2009. The effects of the changes in the fault dip angle, key caprock mechanical parameters, and in situ stress ratio on the displacement profile are evaluated within the parametric study. In comparison to the benchmark numerical run, the scenario ratio of k = 0.5 led to a significant reduction in the displacement. The simulation in which the fault dip angle was 30° produced a more pessimistic result with a larger displacement field. This could be an indication of heightened fault reactivation risks associated with low-angle faults in storage sites with strong horizontal stress regimes due to the combined effect of increased shear stress and reduced inherent frictional resistance on the fault plane. Considering that a vertical fault dip angle (90°) and an additional three 20 m long vertical fractures above the reservoir produced similar vertical displacement observed with the fault dipping at a 60° angle, this indicated that the vertical faults in the vicinity of the storage site pose limited safety risks to the integrity of the sealing rock.

Three-dimensional fault model and features of chained hazards of the Luding MS 6.8 earthquake, Sichuan Province, China

The MS 6.8 Luding earthquake in 2022 is located on the NNW-trending Moxi segment of the Xianshuihe fault with left-lateral strike-slip behavior. This area is where the Xianshuihe, Anninghe, Daliangshan and Longmenshan faults intersect. China Earthquake Administration has identified that intersection area, among the Moxi segment of the Xianshuihe fault, the Anninghe fault, the Daliangshan fault and the southern part of the Longmenshan fault, as a high-magnitude earthquake hazard area. According to existing data on the Luding earthquake, including the focal parameters, the spatial distribution of re-located aftershocks, dominated azimuth of the earthquake intensities and earthquake-induced ground fissures, we built a 3D earthquake fault model. We found that two discontinuous NNW-trending vertical strike-slip faults with left stepping were the seismogenic faults of the Luding earthquake. Its coseismic left-lateral dislocation triggered transtensional slips and aftershocks on the NW-trending secondary faults at its northernmost tensile area. Meanwhile, local crustal coseismic shortening on the side of Mt. Gongga triggered the aftershocks on the NE- and NW-trending secondary conjugated strike-slip faults, which were confirmed by GNSS observations and InSAR deformation field around the epicenter. This earthquake rupturing pattern also controlled the spatial distribution of the earthquake intensity IX area and earthquake chain hazards. The Coulomb stress calculation shows that the Luding earthquake increases the risk of high-magnitude earthquake occurrence on the southernmost part of the Xianshuihe fault and the Anninghe fault. Finally, we suggested doing good monitoring of the Anninghe fault and the southernmost part of the Xianshuihe fault and avoiding active faults with seismogenic capacity and areas prone to earthquake-chained hazards during the site selection and planning of reconstruction.

Sealing faults and fluid-induced seismicity.

Sealing faults are nearly impermeable barriers that can form boundaries between subsurface pore-pressure domains. In hydrocarbon systems, sealing faults commonly form part of a structural trap; they are thus important elements for future storage of CO2 and other gases in depleted reservoirs. The Triassic Montney Formation in western Canada hosts low-permeability gas reservoirs containing sealing faults that have previously been assumed to compartmentalize pressure domains. In this study, we show that the distribution of induced seismicity associated with hydraulic fracturing (HF) exhibits a statistically significant spatial correlation with zones of high lateral gradient in pore pressure. These high-gradient zones are interpreted as sealing fault systems. The largest induced seismicity sequence, including a 4.5 ML mainshock on 30 November 2018, occurred during HF treatments in two horizontal wells, between which there is an exceptionally large contrast (~10 MPa) in measured pore pressure. Numerical simulation of a simplified model of a hydraulic fracture intersecting a nearby vertical fault, followed by fault rupture using rate-and-state friction rheology, generates results that are in good agreement with observed strike-slip faulting near one of the HF wells. Our study demonstrates that sealing faults exhibit previously unrecognized behaviour that may be important for understanding induced seismicity risk. This article is part of the theme issue 'Induced seismicity in coupled subsurface systems'.

Preliminary analysis of the indicators and accumulation models of gas hydrates in the central Zhongjiannan Basin, western South China sea

The Zhongjiannan Basin is located in the western continental margin of the South China Sea (SCS). It was recently reported that there are geophysical indicators of gas hydrates in the northern Zhongjiannan Basin. However, due to the lack of in-depth studies, it remains unknown whether geophysical indicators of gas hydrates are also available in the central and southern Zhongjiannan Basin. Based on the latest geophysical data, we have identified the seismic anomaly features of gas hydrates and observed an obvious bottom simulating reflector (BSR) in the seismic section in the central Zhongjiannan Basin. Based on a comprehensive analysis of seismic characteristics, geological structures (e.g., faults, diapirs and pockmarks) and regional sedimentary conditions, we propose three accumulation models of gas hydrates in the central Zhongjiannan Basin. The first model is that gas develops in deep depression areas, migrates upwards from deep source rocks along vertical faults and forms reservoirs in shallow formations. The second model is that gas also forms in deep depression areas and rises along mud diapirs or gas chimneys from source rocks to shallow formations, with reservoirs formed above diapirs and/or mud volcanoes. The third model is that gas is mainly generated in slope areas where strata are relatively thin and fault structures are well developed; gas migrates upwards from source rocks along vertical faults, and reservoirs are formed beneath pockmarks. The geophysical indicators of gas hydrates are found for the first time in the central Zhongjiannan Basin of the SCS, yielding an important reference for further exploration and development.

Multi-scale hydrostructural approach for karst environment. Application to the Arcier hydrosystem (eastern France)

Based on a multi-scale and hydrostructural approach, this study presents the most relevant methodology to be applied to a karst hydrosystem in order to get a full understanding of underground water flow. It implies a complete structural analysis, from the hydrosystem scale to the outcrop scale, including the intermediate scale of the major geological structures. We illustrate the method in the Arcier hydrosystem, in the northwestern border of the Jura fold-and-thrust belt (Eastern France).Field mapping and structural analysis allow to update the geological vision of the hydrosystem with two kink-type fault propagation folds, including a trishear kinematic model, on either side of a plateau presenting a hollow-and-dome configuration. Fracturing analysis reveals a fault-fracture network that we infer governs the entire hydrosystem. A Riedel pattern is highlighted, characterized by a N–S-striking (N355° ± 5), sinistral strike-slip, regional shear zone. Then, two 3D geological models, at different scales, constructed with MOVE and Visual Karsys softwares are combined with water levels and artificial tracer tests. It reveals a multilayer aquifer and a redefinition of groundwater circulations for the Arcier hydrosystem.The results demonstrate a strong geological control of karstic hydrosystems on groundwater circulations, proving that classical hydrogeological methods, such as natural and/or artificial tracers, must be combined with rigorous geological analysis. Moreover, the multi-scale approach provides an explanation of groundwater circulation based on the intersection between 3D geometry of impervious layers delimiting the aquifers and their base water level, instead of the 2D view (section or map) requiring systematic recourse to inferred vertical faults to cross permeability barriers vertically or laterally. This study also brings a new vision to the local protection of the water resource.

Developing meshing workflows in Gmsh v4.11 for the geologic uncertainty assessment of high-temperature aquifer thermal energy storage

Abstract. Evaluating uncertainties of geological features on fluid temperature and pressure changes in a reservoir plays a crucial role in the safe and sustainable operation of high-temperature aquifer thermal energy storage (HT-ATES). This study introduces a new automated surface fitting function in the Python API (application programming interface) of Gmsh (v4.11) to simulate the impacts of structural barriers and variations of the reservoir geometries on thermohydraulic behaviour in heat storage applications. These structural features cannot always be detected by geophysical exploration but can be present due to geological complexities. A Python workflow is developed to implement an automated mesh generation routine for various geological scenarios. This way, complex geological models and their inherent uncertainties are transferred into reservoir simulations. The developed meshing workflow is applied to two case studies: (1) Greater Geneva Basin with the Upper Jurassic (“Malm”) limestone reservoir and (2) the 5° eastward-tilted DeepStor sandstone reservoir in the Upper Rhine Graben with a uniform thickness of 10 m. In the Greater Geneva Basin example, the top and bottom surfaces of the reservoir are randomly varied by ± 10 and ± 15 m, generating a total variation of up to 25 % from the initially assumed 100 m reservoir thickness. The injected heat plume in this limestone reservoir is independent of the reservoir geometry variation, indicating the limited propagation of the induced thermal signal. In the DeepStor reservoir, a vertical sub-seismic fault juxtaposing the permeable sandstone layers against low permeable clay-marl units is added to the base case model. The fault is located in distances varying from 4 to 118 m to the well to quantify the possible thermohydraulic response within the model. The variation in the distance between the fault and the well resulted in an insignificant change in the thermal recovery (∼ 1.5 %) but up to a ∼ 10.0 % pressure increase for the (shortest) distance of 4 m from the injection well. Modelling the pressure and temperature distribution in the 5° tilted reservoir, with a well placed in the centre of the model, reveals that heat tends to accumulate in the updip direction, while pressure increases in the downdip direction.

Determining complex petroleum migration pathways along faults and sandbodies: a new method and application

Introduction: The reservoir-forming rule in the paleo-slope zone outside petroliferous basins’ source kitchens of the eastern basin in China is unclear, manifested in the distribution of oil and gas in complex multi-layers. This study focuses on petroleum migration along faults and sandbodies to clarify this question and improve the exploration efficiency of petroleum in areas with similar geological conditions.Method: We summarized two models of petroleum migration based on two factors: the connected sandbodies and the vertical unsealing faults. Model I is conducive to the migration of hydrocarbons within a single layer, while Model II is conducive to the vertical migration of hydrocarbons across layers. Combining the means for evaluating fault sealing and sandbody connectivity, integrating these models could form a set of characterization methods for the distribution of the complex fluid migration pathways.Results: This method is applied to characterize the complex multilayer migration pathways distribution in the first and second members of the Shahejie Formation in the central-southern area of the Wen’an Slope in the Baxian Sag, Jizhong Depression, Bohai Bay Basin, China. The results show that, in the first and second members of the Shahejie Formation in the central-southern area of the Wen’an Slope, the area for Model II is mainly distributed in the central-eastern regions. The second member has a Model I area distributed in its western region, and the first member has a Model I area distributed in its centraleastern parts. The latter distribution area is larger than the former. According to hydrocarbon migration analysis, the first member is more conducive to large regional migration and accumulation of hydrocarbons than the second member.Discussion: The predicted result of petroleum migration is consistent with the hydrocarbon distribution found in the central-southern areas of the Wen’an Slope, which is that oil and gas distribution in the first member is significantly greater than that in the second member. Thus, our method is feasible for characterizing the regional distribution of complex multi-layer migration pathways in simple tectonic setting and oil source conditions. However, further research is still needed to increase the effectiveness of present petroleum migration analysis.

The 12 October 2021, MW\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$M_W$$\\end{document}=6.4, Zakros, Crete earthquake

We study the strong 12 October 2021, MW\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$M_W$$\\end{document}=6.4, offshore Zakros, Crete earthquake, and its seismotectonic implications. We obtain a robust location (azimuthal gap equal to 17∘\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{\\circ }$$\\end{document}) for the mainshock by combining all freely available local, regional and teleseismic phase arrivals (direct and depth phase arrivals). Based on our location and the spatial distribution of the poor aftershock sequence we parameterise the fault area as a 30 km ×\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ imes $$\\end{document} 20 km planar surface, and using three-component strong motion data we calculate slip models for both earthquake nodal planes. Our preferred solution shows a simple, single slip episode on a NE-SW oriented, NW shallow-dipping fault plane, instead of a N-S oriented, almost vertical nodal plane. An anti-correlation of the aftershocks spatial distribution versus the maximum slip (∼\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sim $$\\end{document} 27 cm) of our model further supports this, although the accuracy of the aftershock hypocentral locations could be somewhat questionable. Coulomb stress changes calculated for both kinematic models do not show substantial differences, as the aftershock seismicity within the first 3 months after the mainshock is distributed along the stress shadow zone and over the stress enhanced areas developed at the southern fault edge, induced by the mainshock. The Kasos island tide gauge record analysis shows a small signal after the earthquake, but it can hardly demonstrate the existence of tsunami waves due to the low signal-to-noise ratio. Tsunami simulations computed for the two nodal planes do not yield conclusive evidence to highlight whether the causative fault plane is NE-SW oriented, NW shallow-dipping plane, or the N-S oriented plane, nevertheless, the power spectrum analysis of the NW shallow-dipping nodal plane matches the spectral peak at 8 s period and is overall closer to the spectrum of the tide gauge record. A USGS Shakemap was also produced with all available local strong motion data and EMSC testimonies. This was also investigated in an effort to document the responsible fault. The overall analysis in this study, slightly suggests a rather westward, shallow-dipping offshore fault zone, being antithetic to the main Zakros almost vertical normal fault which shapes the coast of eastern Crete and is perpendicular to the direction of Ptolemy Trench in this area. This result agrees with seismotectonic and bathymetric evidence which support the existence of approximately N-S trending grabens, east and northeast of Crete.

Optical flow-based vertical angular rate fault detection on UAVs

This paper discusses optical flow-based vertical angular rate fault detection based-on real fight IMU data and camera image sequences. The fight test platform utilizes a gimbal stabilized downward looking camera that is why only the angular rate around the inertial vertical axis can be estimated. Improvements relative to the previous work of the authors are the application of forward-backward optical flow to filter outliers, the direct detection of angular rate faults without the need for Euler angles or GPS velocities and tuning and verification on real images and fight data. After presenting the optical flow equations and the method for image-based vertical angular rate estimation the fight test scenarios and video processing steps are introduced. Four different Artificial errors are added to the IMU measurements and the fault detection applies residual thresholding and up-down counters. Tuning is based on false alarm and missed detection rates and the application of receiver operating characteristics. Finally, the estimation results and possible further improvements are presented.

Geomechanical modelling of ground surface deformation induced by CO2 injection at In Salah, Algeria: Three wells, three responses

This study investigates the ground surface deformation above three injection wells throughout CO2 injection between mid-2004 to early-2011 in the In Salah CO2 storage project at Krechba, Algeria. A coupled three-dimensional finite element model of three wells (KB501, KB502, and KB503) was developed. The model accounts for the CO2 diffusion in the reservoir and structural features including permeable fractures and faults, as well as the poroelastic deformation of the reservoir and overburden layers which result in the ground surface deformations during CO2 injection. The character of the ground surface deformation above each injection well was different, suggesting a different CO2 plume shape in the subsurface. An ellipse-shaped heave above KB501 indicates that the CO2 plume shape followed the anisotropic permeability in the reservoir, KB502’s double-lobe surface heave pattern is associated with the activation and pressurization of a vertical fault and finally, the elongated surface uplift pattern above KB503 is linked to the presence of permeable fractured zone in vicinity of the well, as indicated by seismic surveys. The two distinctive structural features in the vicinity of KB502 and KB503 were included in the numerical model, accounting for a hydraulically conductive fault in the caprock and a fractured zone in the reservoir, respectively. The observed surface deformation was utilised to tune the orthotropic mechanical properties of the fault and the fractured zones. The simulation results are in good agreement with the measured InSAR data. In addition, assigning anisotropic permeability to the reservoir plays a significant role in successfully simulating surface uplift patterns.

Study on Evolution Characteristics of Stress Field of Stepped Fault

With the development of mining to the deep, the ground stress increases correspondingly, and the geological conditions become more and more complicated. Mine dynamic disasters mainly occur near faults. Numerous studies have been conducted on single faults, but few studies have been done on stepped faults formed by double faults under complex geological conditions. Based on this, with the 3110-working face of Zhong Tai Mining in He bi as the engineering background, Single fault model a and different step fault models b, c, d, e and f were established by Flac3D numerical simulation software. The characteristics of stress field evolution of single fault and step fault under the influence of deep mining are studied. The results show that the stress distribution in the upper wall of the fault forms a stress concentration zone, and in the lower wall forms a stress weakening zone. Under the influence of excavation disturbance, the vertical fault distance and the horizontal distance between faults will affect the stress evolution characteristics of the stepped fault.

Lithospheric plane strain deformation caused by a vertical dip-slip fault

Lithospheric plane strain deformation caused by a vertical dip-slip fault