Plane Dipping Research Articles

Analysis of failure properties of red sandstone with structural plane subjected to true triaxial stress paths

In order to study the mechanical behavior and failure characteristics of rock mass with the structural plane in complex deep stress environment, the QKX-YB200 true triaxial rockburst test system, acoustic emission, and digital video recorder were adopted to conduct true triaxial loading and unloading indoor tests on the red sandstone cube specimens. The results show that: (1) Under true triaxial loading conditions, all specimens propagate in the form of anti-wing crack, and the dip angle of the structural plane is linearly related to the average value of the anti-wing crack propagation angle; under true triaxial unloading conditions, the specimens with structural plane dip angles of 30°, 45°, and 60° propagate in the form of wing crack, while those with a 90° angle propagate anti-wing cracks, and the closure of structural plane is related to the the crack type at the tips of structural plane. (2) The unstable failure stages in the stress–strain curve under true triaxial unloading are longer and exhibit more pronounced fluctuations; the peak strength of the specimens under the conditions of true triaxial loading and unloading decrease first and then increase with the increase of the dip angle, reaching the lowest at 45°. (3) During the true triaxial unloading tests, the rockburst degree of the intact rock specimens is more violent than that of the specimens with a structural plane. Additionally, the dip angle of the structural plane influences the severity of rockbursts under the same stress path. (4) The crack propagation mechanism and the closure degree of structural planes under different true triaxial stress paths and dip angles are discussed. It is observed that the crack initiation mode around the structural plane tips and the crack mode near free surface are closely related to the stress environment in different areas of the specimens. (5) Based on the indoor test results and previous experimental data, the differences in peak strength characteristics between open structural planes and closed structural planes are preliminarily analyzed and compared.

A rapid analysis of aftershock processes after a moderate magnitude earthquake with ML methods

SUMMARY Moderate magnitude earthquakes and seismic sequences frequently develop on fault systems, but whether they are linked to future major ruptures is always ambiguous. In this study, we investigated a seismic sequence that has developed within a portion of the stretching region of the Apennines in Italy where moderate to large earthquakes are likely to occur. We captured a total of 2039 aftershocks of the 2023 September 18, Mw 4.9 earthquake occurred during the first week, by using machine-learning (ML) based algorithms. Aftershocks align on two 5–7 km long parallel faults, from a length that exceeds what is expected from the main shock magnitude. The segments are ramping at about 6 km depth on closely spaced N100 striking 70 N dipping planes, at a distance of some kilometres from the main shock hypocentre. Our results indicate that even moderate magnitude events trigger seismicity on a spread set of fault segments around the main shock hypocentre, revealing processes of interaction within the crustal layer. The possibility that larger earthquakes develop during seismicity spread is favoured by pore pressure diffusion, in relation with the closeness to criticality of fault segments. Based on the very rapid activation of seismicity on the entire system and a back-front signal from the hypocentre of the main event, we infer that fluid pressure, initially high within the crustal layer, rapidly dropped after the main shock. Our study reinforces the importance of timely extracting information on fault geometry and seismicity distribution on faults. ML-based methods represent a viable tool for semi-real-time application, yielding constraints on short-time forecasts.

Contribution for the knowledge of the Alcochete-Pinhal Novo-Setúbal fault (Portugal) using ReMi and H/V Tests

The city of Setúbal is situated 30 km southeast of Lisbon, Portugal, and has experienced significant impacts from various earthquakes, including the one that occurred in 1858. This earthquake serves as an example of a nearby earthquake that caused significant social and structural damage, reaching a intensity of IX-VIII on the Modified Mercalli scale. Despite an average recurrence period of 3000–11000 years, as suggested by several authors, the Alcochete-Pinhal Novo-Setúbal fault (APSF) could potentially cause significant social, structural and economic impacts in future events. Several authors have reported on the APSF with slight variations in its outline, as it traverses multiple urbanized areas. The aim of this study was to locate the APSF in the Setúbal region and, to the best extent possible, determine its orientation and vertical offsets. 24 passive refraction microtremor tests (ReMi) and 145 horizontal-to-vertical (H/V) spectral ratios were employed to compute the variation of transverse wave velocity (VS) with depth and to compute the fundamental frequency (F0) and other natural frequency (F1) of the soils when present. Horizontal sections of VS and the spatial distribution of peak frequencies reveal alignments that can be related to the presence of the APSF fault zone. Considering all the alignments identified in this paper, two areas with distinct orientations for the APSF zone were proposed: one with a NNE-SSW strike and another with a NNW-SSE orientation. The lower values calculated for the fault plane dip (maximum of 5º) in the APSF area suggest a minimal vertical displacement, potentially corresponding to a strike-slip fault.

Contractional strains and maximum displacement-length ratios of lunar wrinkle ridges in four Maria of basalt

Compressional stresses from the basalt basins on the lunar are responsible for the formation of wrinkle ridges in lunar mare basalts. According to the wide angle multispectral camera (WAC) mosaic image, we selected 62, 75, 73, and 58 single wrinkle ridges in Mare Imbrium, Mare Serenitatis, Mare Fecunditatis, and Mare Tranquillitatis, respectively, for this paper. Several topographic profiles near the midpoint of each wrinkle ridge are generated to measure the maximum displacement (Dmax) and height of the wrinkle ridges using digital elevation model (LOLA) data. After that, we create 2D plots of displacement-length (L) for the ridge population in the four Maria and compare the results. A linear fit method derives the Dmax/L ratios (γ) from the D-L data. We calculated the contractional strains in each mare area based on the cap D sub m a. x over cap L data. Moreover, each mare's gravity pattern, mare thickness, and formation age are also presented. The ridges in Mare Imbrium and Mare Serenitatis have a higher γ value (1.83 × 10−2 and 1.98 × 10−2) than the ridges in Mare Fecunditatis and Mare Tranquillitatis, which have γ values of (1.67 × 10−2 and 1.75 × 10−2). Finally, the contractional strains (ε) in Mare Imbrium, Mare Serenitatis, Mare Fecunditatis, and Mare Tranquillitatis are estimated to be 0.23 %, 0.41 %, 0.39 %, and 0.18 % (considering 25° is the fault plane dip θ), respectively. The maximum values of the free-air gravity anomalies in Mare Fecunditatis range from −30 to 250 mGal, while minimum gravity anomalies in Mare Serenitatis range from −80 to 140 mGal. Mare Imbrium, Mare Serenitatis, Mare Fecunditatis, and Mare Tranquillitatis have an average thickness of 300 m, 910 m, 652 m, and 760 m, respectively. Furthermore, the Mare Imbrium ridge group is older than the Lunar Wrinkle Ridges in Mare Serenitatis. Mare Tranquillitatis ridge group formation takes longer than Mare Imbrium ridge group formation. Therefore, we believe that it has thicker basaltic units, a longer wrinkle ridge formation time, and higher gravity anomaly than the Mare Imbrium and Mare Serenitatis basins, even though the formation of the Mare Tranquillitatis and Mare Fecunditatis basins occurred earlier.

Brittle rock mass failure in deep tunnels: The role of infilled structural plane with varying dip angles

Filled rock discontinuities are commonly encountered in underground engineering excavations, and their mechanical behavior has a significant impact on the stability of surrounding rocks. However, their performance under three-dimensional stress conditions has rarely been explored in the past. In this study, we summarized the typical failure characteristics of filled structural planes in deep-buried tunnels and conducted laboratory true triaxial tests on granite specimens to investigate the effect of dip angle and infill thickness on the failure characteristics of high-stress brittle surrounding rocks. The study monitored the characteristics of vertical stress-displacement, acoustic emission, acceleration, and failure process of granite specimens during the loading process and revealed the failure mechanism of brittle rock masses under the influence of infill thickness and structural plane dip angle. The results showed that the dip angle of structural planes has a significant influence on the failure mode of rock masses in deep tunnels. There exists a critical range of dip angles, within which the rock mass undergoes shear slip failure along the pre-existing structural plane, manifested by rock spalling near the structural plane. Beyond this critical range, the rock mass experiences strain-type brittle failure, characterized by violent rock ejection. Furthermore, the filling of structural planes can change the failure mode of rock masses and lead to more severe damage when the friction coefficient of the structural plane falls below a specific threshold. Additionally, it was revealed that a perpendicular ejection force exists on the free surface of the rock mass under loading, which is approximately vertical to the normal direction and is the primary factor responsible for the formation of vertical parallel cracks within the rock mass. The findings of this work are of paramount importance for comprehending the failure characteristics of deep-buried tunnels containing infilling joints.

Quantitative seismic fracture characterization of a sandstone reservoir — Decatur, Illinois Basin

The Illinois Basin Decatur Project, a carbon capture and sequestration task, was undertaken to sequester 1 million tonnes of CO2 into a sandstone reservoir. A 3D seismic survey was conducted to characterize the reservoir. A geomodel was developed from seismic data, inversion results, and well data to geostatistically map the storage potential of the reservoir. However, no fracture model was created or utilized in this exercise. Fractures inherently influence the porosity and permeability of a reservoir. Ignoring them in reservoir characterization is not an optimal reservoir management practice. The image-log interpretation from a few vertical wells drilled in the area shows the bedding plane dips, but no fracture has been identified. However, the lack of fracture crossings in a few vertical wells does not imply that a formation is devoid of fractures altogether. Hence, seismic fracture characterization (leveraging the dense 3D seismic data) is necessary for a reservoir characterization exercise. We utilized the publicly available Decatur 3D seismic data set to run a seismic fracture characterization workflow to delineate potential fracture corridors present in the reservoir. We calculated three edge detection attributes (structural tensor, structure-oriented semblance, and structural dip) in combination to delineate the fracture lineaments. Our workflow extracts several quantitative measures of the seismic lineaments such as dip, azimuth, area, and length, which can be analyzed statistically. The principal focus of this work is to find a way forward to integrate the fractures from seismic data in a geologic model that can be utilized in simulations. Based on our interpretation of seismic fractures, we created a discrete fracture network that can be a building block for creating a finer-resolution fracture model. We also explored the fractal characteristics of seismic-derived fracture lineaments as a way forward for generating discrete fracture networks.

EARTHQUAKE of APRIL 24, 2018 with KR=10.5, Мw3.9 and I0=5 in the ANAPA-NOVOROSSIYSK ZONE

The article presents instrumental and macroseismic data on the earthquake of April 24, 2018 with Мw=3.9, which occurred at the Black Sea coast near the Abrau Peninsula. The earthquake was felt in 15 settlements at distances up to 52 km from the epicenter with an intensity of 3 to 5 points. The earthquake occurred in the Black Sea zone of occurrence of earthquakes. The scalar seismic moment М0 and the spectral moment magnitude Mw=3.9 were calculated based on the spectra of seven regional stations. The focal mechanism was constructed based on the signs of the first movement in the P-wave at 68 stations. It is established that a thrust has occurred in the source along a northwest-trending plane dipping to the northeast, which is consistent with the orientation and dip of the fault planes of tectonic structures in the area of the source. The macroseismic manifestations of the April 24, 2018 earthquake made it possible to choose the most acceptable depth value and calculate the intensity at the epicenter (h=10 km, I0=5)

Mechanical characteristics and energy evolution in a rock mass with a weak structural plane

The present aims to investigate the mechanical characteristics and energy evolution in rock masses containing weak structural planes under conventional triaxial loading conditions. Using a fluid–solid coupling test system of coal rock, numerous conventional triaxial compression tests were performed on rock masses at various dip angles of the structural plane. The obtained empirical outcomes revealed that the deviatoric stress–strain curve of the weak structural plane rock mass with an inclination angle greater than 20° rises step-by-step. On the macro level, slip-stability occurs on the upper and lower parts of the rock mass on the weak structural plane. Then mechanism of the slip-stability phenomenon is explored by analyzing the stress level in the rock mass with various inclination angles. It is found that the energy evolution during deformation and failure reflects the damaged state of the rock. Accordingly, the concept of ‘slip dissipation energy’ is proposed, and the values of each energy are calculated. The results have a good correspondence with the deviatoric stress–strain curve. Furthermore, it was found that the energy evolution of rock mass with a weak structural plane can be primarily classified into four stages, including storage of the initial energy, slip dissipation, abrupt increase in the pre-peak dissipation energy, and sudden drop in post-peak energy. Rock masses with various levels of dip angles exhibit similar elastic strain energy and dissipation energy at the peak point, demonstrating that energy evolution is dominated by energy storage and dissipation. At the same time, a negative correlation is observed between the structural plane dip angle and the occurrence of instantaneous impact instability failure in rock masses, indicating that a greater dip angle makes the rock mass less prone to experiencing instantaneous impact instability failure. This article provides a new idea for analyzing the geological disasters caused by external disturbances.

A new critical strain energy release rate failure criterion for anisotropic rock under high confining pressure

To research the influences of confining pressure on critical strain energy release rate (Gc) and energy transformation of anisotropic intact rocks, by combining a new modified Hoek–Brown failure criterion for anisotropic rock under high confining pressure (M-HB-HC) with the energy release rate failure criterion, the new critical strain energy release rate failure criterion is established. Triaxial compression tests for marble and sandstone samples with different bedding plane dip angles under different confining pressures have been performed. The study shows that Gc depends on both the bedding plane dip angle α and confining pressures, and the variation of Gc with the dip angle under low confining pressure is not evident, but is significant at the high confining pressure. Under the same confining pressure, the value of bedding marble’s Gc is greater than that of bedding sandstone. By comparing prediction results of several failure criteria with that of the proposed failure criterion and the experimental data, we found that the prediction of the new critical strain energy release rate failure criterion has better accuracy. The proposed failure criterion is applicable for different type of anisotropic intact rocks. This paper only researches the transverse isotropic case of the anisotropic rock. The other anisotropic will be studied in the future.

Structure-adapted Multichannel Matching Pursuit for Seismic Trace Decomposition

While the single-channel matching pursuit decomposes a seismic trace into a series of wavelets, the multichannel matching pursuit examines the lateral coherence of the seismic traces as a constraint to improve the lateral continuity of the decomposition. However, the presence of structures in the subsurface negatively affects the performance of the multichannel matching pursuit. We proposed a structure-adapted matching pursuit method that combines with a dynamic time-warping (DTW) algorithm to estimate the similarity between adjacent seismic traces and extract an optimal wavelet along the dip plane. This structure-adapted implementation would significantly speed up the convergence of the decomposition process. We modified the DTW algorithm by combining the Euclidean distance of the seismic trace and the first-order temporal derivative of the seismic trace. We also updated the amplitudes of all extracted wavelets simultaneously based on the least-squares principle. This DTW-based, structure-adapted, least-squares, multichannel matching pursuit method would improve the robustness and accuracy of seismic trace decomposition.

Structural Origin of Sinjar Anticline, NW Iraq

The Sinjar anticline is a double plunging, trending almost E-W in the northwestern part of Iraq. It extends in Syria for about 42 km, whereas in Iraq, its length is about 91 km, and the width is about 31 km. The northern limb (45° - 80°) is steeper than the southern limb (15° - 25°), with average plunges dip of 35° and axial plane dipping of 47.5° southwards. The exposed rocks in the anticline range in age from Upper Cretaceous, represented by the Shiranish Formation, to Upper Miocene, represented by the Injana Formation. Google Earth image was used to calculate structural data, which were used to indicate the structural origin of the Sinjar anticline. This was achieved by calculating the Aspect Ratio (AR), Fold Symmetry Index (IFS or IFS), and length of the mountain front (FS). Accordingly, it was found that the structural origin of the Sinjar anticline is a fault-bend fold.

Numerical modelling of the growth of polygonal fault systems

Despite over three-decades of active research and wide debate in the published literature, the mechanisms that govern the growth of polygonal faults are poorly understood. Here we investigate the growth of polygonal faults using a suite of geomechanical finite element forward models that couple dynamic fault propagation, sedimentation, and the mechanical compaction of unconsolidated granular sediment. We undertook a suite of numerical model simulations to explore the relationships between varying fault plane dip, residual friction of the fault, and the bulk material properties of the sedimentary sequence hosting the polygonal fault system. We find that the growth of polygonal faults within laterally-pinned sedimentary tiers can be explained by gravity-driven differential compaction and does not require additional causative elements to explain the gross pattern of strain accumulation. We also find that the magnitude of fault throw is influenced by the material properties and the original fault plane dip, but is most sensitive to the residual friction angle. Our models yield values for maximum throw versus height for the faults that fall within the range of global values compiled for polygonal faults, and throw rates are comparable to those recently measured in naturally occurring polygonal faults.

Precise aftershock distribution of the 2019 Yamagata-oki earthquake using newly developed simple anchored-buoy ocean bottom seismometers and land seismic stations

An earthquake with a magnitude of 6.7 occurred in the Japan Sea off Yamagata on June 18, 2019. The mainshock had a source mechanism of reverse-fault type with a compression axis of WNW–ESE direction. Since the source area is positioned in a marine area, seafloor seismic observation is indispensable for obtaining the precise distribution of the aftershocks. The source area has a water depth of less than 100 m, and fishing activity is high. It is difficult to perform aftershock observation using ordinary free-fall pop-up type ocean bottom seismometers (OBSs). We developed a simple anchored-buoy type OBS for shallow water depths and performed the seafloor observation using this. The seafloor seismic unit had three-component seismometers and a hydrophone. Two orthogonal tiltmeters and an azimuth meter monitored the attitude of the package. For seismic observation at shallow water depth, we concluded that an anchored-buoy system would have the advantage of avoiding accidents. Our anchored-buoy OBS was based on a system used in fisheries. We deployed three anchored-buoy OBSs in the source region where the water depth was approximately 80 m on July 5, 2019, and two of the OBSs were recovered on July 13, 2019. Temporary land seismic stations with a three-component seismometer were also installed. The arrival times of P- and S-waves were read from the records of the OBSs and land stations, and we located hypocenters with correction for travel time. A preliminary location was performed using absolute travel time and final hypocenters were obtained using the double-difference method. The aftershocks were distributed at a depth range of 2.5 km to 10 km and along a plane dipping to the southeast. The plane formed by the aftershocks is consistent with the focal mechanism of the mainshock. The activity region of the aftershocks was positioned in the upper part of the upper crust. Focal mechanisms were estimated using the polarity of the first arrivals. Although many aftershocks had a reverse-fault focal mechanism similar to the focal solution of the mainshock, normal-fault type and strike–slip fault type focal mechanisms were also estimated.Graphical

Experimental Investigation on the Initiation of Hydraulic Fractures from a Simulated Wellbore in Laminated Shale

Abstract Understanding the formation mechanisms of complex fracture networks is vitally important for hydraulic fracturing operations in shale formation. For this purpose, a hydraulic fracturing experiment under a core-plunger scale is conducted to investigate the impact of the bedding plane angle, borehole size, and injection rate on fracture initiation behaviors of laminated shale rock. The results on rock properties demonstrate that the anisotropic characteristics of shale rock are reflected not only in elastic modulus but also in tensile strength. The results of fracturing experiments show that the bedding plane dip angle and borehole size have significant effects on fracture initiation behaviors, in that fracture initiation pressure (FIP) decreases with the increase of those two factors. The impact of injection rate, by contrast, has no obvious variety regulation. The above data is further used to validate our previously proposed fully anisotropic FIP model, which shows better agreement with experimental results than those using other models under various parameter combinations. Finally, a postfracturing analysis is performed to identify the fracture growth patterns and the microstructures on the fracture surfaces. The results show that the hydraulic fractures (HFs) always grow along mechanically favorable directions, and the potential interaction between HFs and bedding planes mainly manifests as fracture arrest. Meanwhile, the roughness of fracture surfaces is physically different from each other, which in turn results in the difficulties of fluid flow and proppant migration. The findings of this study can help for a better understanding of the fracture initiation behavior of laminated shale rock and the corresponding fracture morphology.

Research on the anisotropic fracture behavior and the corresponding fracture surface roughness of shale

Understanding the anisotropic fracture mechanism of shale is of significance to assess the development of a stimulated reservoir volume. In this study, three-point bending tests were carried out using notched deep beam specimens of shale with 7 different bedding dip angles to investigate the anisotropic fracture behavior. The anisotropic characteristics of the fracture toughness, energy release rate, inelastic zone and crack propagation path were comprehensively analyzed and highlighted. The results showed that the fracture toughness decreases with increasing bedding dip angle. When the bedding dip angle is 0° and 90°, the fracture toughness reaches the maximum and minimum values, respectively, with a ratio of 1.60. The variation tendency of the energy release rate with the bedding inclination angle is consistent with that of the fracture toughness, but the ratio of the maximum value to the minimum value of the energy release rate is about 3.24. The size of the inelastic zone first decreases and then increases at different bedding inclination angles from 0° to 90° and reaches the maximum size of 1.51 mm when the bedding plane inclination angle is 45°. The crack propagates along the bedding plane rather than straight along the initial crack direction when the bedding dip angle is higher than 45°. As the bedding plane dip angle increases from 0° to 45°, the crack surface roughness of the tested shale first increases, then decreases and finally reaches the maximum. The discussion indicated that the presence of the Mode II fracture has a considerable influence on the crack surface roughness and crack propagation path, although it has little effect on the failure load. Compared with the maximum energy release rate (MERR) criterion, the maximum tangential stress (MTS) criterion seems more accurate in predicting the anisotropic Mode I fracture toughness. The applicability of these two classical crack initiation criteria of fracture mechanics has been analyzed based on our experimental results, which demonstrates the validity of the MTS criterion to assess the anisotropic fracture toughness of shales.

Identification of rock and fracture kinematics in high alpine rockwalls under the influence of elevation

Abstract. In alpine environments, tectonic processes, past glaciation and weathering processes fracture rock and prepare or trigger rockfalls, which are important processes of rock slope evolution and natural hazards. In this study, I quantify thermally and ice-induced rock and fracture kinematics and place these in the context of their role in producing rockfall and climate change. I conducted laboratory measurements on intact rock samples and installed temperature loggers and crackmeters at four rockwalls reaching from 2585 to 2935 m in elevation in the Hungerli Valley, Swiss Alps. My laboratory data show that thermal expansion followed three phases of rock kinematics, which resulted in a hysteresis effect. In the field, control crackmeters on intact rock reflected these temperature phases, and based on thermal expansion coefficients of these observed phases, I modelled thermal stress. Model results show that thermal stress magnitudes were predominantly below rock strengths. Crackmeters across fractures revealed fracture opening during cooling and reverse closing behaviour during warming on daily timescales. Elevation-dependent snow cover controlled the number of daily temperature changes and thermal stresses affecting both intact and fractured rock, while the magnitude is controlled by topographic factors influencing insolation. On a seasonal scale, slow ice-segregation-induced fracture opening can occur within lithology-dependent temperature regimes called frost cracking windows. Shear plane dipping controlled whether fractures opened or closed irreversibly with time due to thermally induced block crawling on an annual scale. Climate change will shorten snow duration and increase temperature extremes and will, therefore, affect the number and the magnitude of thermal changes and associated stresses. Earlier snowmelt in combination with temperature increase will shift the ice-induced kinematic processes to higher elevations. In conclusion, climate change will affect and change rock and fracture kinematics and, therefore, change rockfall patterns in alpine environments. Future work should quantify rockfall patterns and link these patterns to climatic drivers.

The Entire Crust can be Seismogenic: Evidence from Southern Malawi

AbstractThe Bilila‐Mtakataka Fault (BMF), at the southern end of the western branch of the East African Rift System (EARS), has been used in various scaling relation studies and arguments about the strength of the lithosphere. We present evidence for a similar, though more degraded, frontal scarp on the graben‐bounding synthetic Chirobwe‐Ntcheu Fault (CNF), showing that this fault is active simultaneously with the BMF. We deployed 17 geophones for ∼60 days around the southern end of Lake Malawi, across the footwall and hangingwall of the BMF. Continuous microseismicity can be seen from the surface to ∼35 km depth highlighting a plane dipping ∼42°E. Lower‐crustal earthquakes have previously been found in the EARS, and based on location and focal mechanism have been hypothesized to occur on planes that line up with the surface traces of large faults. However, no previous study of the EARS has revealed a fault plane throughout the crust that shows seismicity along its full length from the surface to the base of the crust. Rather, the lack of seismicity seen at mid‐lower crustal depths, has led some people to the “jelly sandwich” hypothesis. Our results show that the entire crust is seismogenic, so support the “crème brûlée” model. In our two month deployment we recorded 22 aftershocks ML ≥ 2 from the March 8th, 2018 earthquake 200 km south of our array, 7 months after the mainshock, confirming that aftershock sequences in regions of low strain have a long duration, and could be the main component of seismicity in slowly straining regions.

Earthquake Source Investigation of the Kanallaki, March 2020 Sequence (North-Western Greece) Based on Seismic and Geodetic Data

The active collision of the Apulian continental lithosphere with the Eurasian plate characterizes the tectonics of the Epirus region in northwestern Greece, invoking crustal shortening. Epirus has not experienced any strong earthquakes during the instrumental era and thus there is no detailed knowledge of the way the active deformation is being expressed. In March 2020, a moderate size (Mw 5.8) earthquake sequence occurred close to the Kanallaki village in Epirus. The mainshock and major aftershock focal mechanisms are compatible with reverse faulting, on NNW-ESE trending nodal planes. We measure the coseismic surface deformation using radar interferometry and investigate the possible fault geometries based on seismic waveforms and InSAR data. Slip distribution models provide good fits to both nodal planes and cannot resolve the fault plane ambiguity. The results indicate two slip episodes for a 337° N plane dipping 37° to the east and a single slip patch for a 137° N plane dipping 43° to 55° to the west. Even though the area of the sequence is very close to the triple junction of western Greece, the Kanallaki 2020 activity itself seems to be distinct from it, in terms of the acting stresses.

Investigation on the fracture and mechanical behaviors of simulated transversely isotropic rock made of two interbedded materials

Transversely isotropic rock composed of two interbedded layers has been rarely studied, although landslides frequently occur in bedded rock masses. With reference to the interbedded rock mass in the Miaowei Reservoir area of Yunnan Province, China, this study prepares transversely isotropic specimens made of two interbedded materials to investigate the fracture and mechanical behaviors of interbedded rock. The effect of the bedding plane dip angle and the application of the findings to an anti-dip interbedded rock slope are also studied. Specifically, two groups of isotropic specimens and five groups of transversely isotropic specimens are tested in a series of triaxial compressive tests; based on these tests, the crack types, failure modes and mechanical parameters are analyzed, and one damage model and its application to the failure probability analysis of an anti-dip interbedded rock slope are investigated. The experimental results show that twelve crack types, eight failure modes and five axial stress-strain curve types can be summarized and are significantly affected by the dip angle. The recognition of these crack types and failure modes can improve the qualitative instability analysis of bedded rock masses in the reservoir area, and five mechanical parameters, including peak stress, strain related to peak strength, residual strength, elastic modulus and Poisson's ratio, exhibit fluctuating variation trends versus the dip angle and confining stress due to the influence of the dip angle. Two groups of theoretical curves from the proposed damage model correspond with the pre-peak features of the test stress-strain curves. The proposed failure probability model considering rock damage shows good practicality when applied to an anti-dip interbedded rock slope; the slope displays discrete failure probabilities in various parts, e.g., the slope toe, slope shoulder and weaker layer show higher failure probabilities, and the relaxation behavior increases the failure probabilities and gradually eliminates the discreteness. Relevancy analysis shows that the strength of transversely isotropic rock is more sensitive to the elastic modulus, dip angle of bedding plane, confining stress and composition of the rock.

Numerical Simulation of the Interaction Between Normal Fault and Bedding Planes Using PFC

The effects of dip angles and thicknesses of bedding layers at the interaction line in between the planes of a normal fault and bedding layers are numerically studied by a discrete element modelling technics. The calibration of the numerical method is accomplished by using an inverse-modelling approach. The laboratory results of Brazilian tensile tests may be used to measure the micro-mechanical parameters of the intact rocks. The shear test results are also performed on the rock joints to determine the micro-mechanical parameters of the bedding interfaces in the present simulation work. The numerical simulation of the shear behaviour of normal fault is performed by making box models with dimensions of 100 mm × 72 mm to represent the bedding layers with different dips and thicknesses. The bedding plane dip angles change from 0° to 165° with an increment of 15° and the thicknesses of these layers change from 6 mm to 12 mm with an increment of 6 mm. The shear testing conditions are exposed to these models to form the normal faults for the numerical modelling of the physical problem under a lateral pressure of 2 MPa. The numerical modelling results show that the shear cracks initiate from all of the bedding interfaces. In this numerical modelling, three major oriented parallel tensile bands are developed and one vertical tensile fracture also propagates through the shear band so that the model is failed. As the bedding layers’ dips are increased, the angle in between the shear bands and the bedding layers is decreased, and as the bedding thickness is decreased, the lengths of shear bands are increased. However, when the bedding angles are 120° and 135° the minimum shear strength of the model is occurred.