Slip Displacement Research Articles

Impact of reservoir pressure fluctuations on fracture shear slip behavior in shale gas hydraulic fracturing

To elucidate the effect of fluid pressurization rate on the shear slip behavior of shale fractures during reservoir hydraulic pressure fluctuations, shale slip experiments with constant axial stress were conducted at 0.5, 2, 8, and 30 MPa/min pressurization rates based on the pressure fluctuation characteristics of hydraulically fractured reservoirs of deep shale gas in the Sichuan Basin. The results indicated that quasi-static slip exhibited higher average velocities than creep slip by 2–3 orders of magnitude at equivalent pressurization rates. Additionally, at 30 MPa/min pressurization rates, the slip type transitioned directly from creep to dynamic slip. The slip velocity increases with increasing pressurization rate, and the average velocities of quasi-static slip are 7.10 × 10−4, 2.20 × 10−3, and 5.40 × 10−3 mm/s, respectively. Dynamic slip exhibited the highest critical slip pressure of 7.14 MPa and the largest friction coefficient increased by 55% at a pressurization rate of 2 MPa/min. Friction strength primarily increased by 42% during creep slip at a pressurization rate of 0.5 MPa/min, while the friction coefficient minimally decreased by 0.03% during dynamic slip. Notably, dynamic slip exhibited a significant increasing trend in the percentage of friction coefficient increment. The critical dynamic slip pressure and energy released during slip exhibited an initially increasing and then decreasing pattern over the range of pressurization rates in this experiment, which is a result of the degree of inhomogeneity of the fluid increasing and then decreasing with the pressurization rate. Consequently, this led to maximum accumulation and energy release at the fracture surface, resulting in minimal quasi-static slip displacement, and is not conducive to improved reservoir permeability characteristics and safe shale gas recovery.

Comminution‐Induced Transient Frictional Behavior in Sheared Granular Halite

AbstractGrain comminution is commonly observed in numerous geological settings. To elucidate the role of grain comminution in dry granular friction, we sheared breakable halite (NaCl) grains using a ring‐shear configuration at a constant slip rate under various normal stresses. We observed transient frictional behaviors: a constant regime exhibiting a high friction coefficient at small slip displacements, and a weakening regime showing a substantial decay in friction at large slip displacements. The characteristic slip lengths for both regimes decreased with normal stress and were characterized by similar exponents. Micro‐X‐ray tomography revealed the evolution of microstructure from distributed grain comminution to progressive shear localization for these two regimes. We propose that the filling processes of comminuted fine particles, during which fine particles saturate and then overflow the shear zone, define transient frictional behaviors. This study may hold significant implications for natural shear systems, given the ubiquity of comminution and localization phenomena.

Uncovering the Fracturing Mechanism of Granite Under Compressive–Shear Loads for Sustainable Hot Dry Rock Geothermal Exploitation

Shear-dominated hazards, such as induced earthquakes, pose an escalating threat to the sustainability and safety of the geothermal exploitation. Variations in fault orientations and compression–shear stress ratios exert a profound influence on the failure processes underlying these disasters. To better understand these effects on the shear failure mechanisms of hot dry rocks, mode-II fracturing tests on granites were conducted at varying loading angles (specifically, 55°, 60°, 65°, and 70°). These tests were accompanied by a comprehensive analysis of the mechanical properties, energy dissipation behavior, acoustic emission (AE) responses, and digital image correlation (DIC)-extracted displacement fields. The tensile–shear properties of stress-induced microcracks were discerned via AE characteristic parameter analysis and DIC displacement decomposition, and the mode-II fracture energy release rate was quantitatively characterized. The results reveal that with increasing compression–shear loading angles, the mechanical properties of granites are weakened, and the elastic strain energy at peak stress gradually decreases, while the slip-related dissipated energy increases. Throughout the fracturing process, the AE count progressively climbs and reaches a peak near catastrophic failure, with an upsurge in low-frequency and high-amplitude AE events. Microcrack distribution concentrates aggregation along the shear plane, reflecting the emergent displacement discontinuities evident in DIC contours. Both the AE characteristic parameter analysis and DIC displacement decomposition demonstrate that shear-sliding constitutes the paramount mechanism, and the fraction of shear-oriented microcracks and the ratio of tangential versus normal displacement escalate with increases in shear stress. This analysis is supported by the heightened propensity for transgranular microcracking events observed through scanning electron microscopy. As the shear-to-compression stress increases, the energy concentration along the shear band intensifies, with the gradient of the fitting line between cumulative AE energy and slip displacement steepening, indicative of a heightened mode-II energy release rate. These results contribute to a deeper understanding of the mode-II fracture mechanism of rocks, thereby providing a foundational basis for early warnings of shear-dominant geomechanical disasters, and improving the safety and sustainability of subsurface rock engineering.

The study on bond‐slip constitutive model of recycled concrete under different loading rates

AbstractIf recycled concrete is to be widely, reasonably, and safely applied in practical engineering, the study of bond‐slip performance between steel bars and concrete is crucial. In this paper, the central pull‐out test of recycled concrete under different loading rates was conducted. The characteristics of failure mode and crack directions of the specimens were observed and analyzed. The ultimate bond strength and slip displacement between rebar and concrete were analyzed. The effects of three parameters, namely concrete strength grades, replacement percentage of recycled coarse aggregate, and loading rates on the failure modes and bond‐slip mechanical properties of recycled concrete, were analyzed and summarized. The experimental results showed that there were significant differences in the influence trend and amplitude of the above three parameters on the ultimate bonding strength, slip amount, and slip curve. Based on the above experimental and theoretical analysis, a reasonable bond‐slip constitutive model was established in this paper and was in good agreement with the experimental results.

Transition From Reverse to Left‐Slip on the Eastern Haiyuan Fault, NE Tibetan Plateau, From the Structure and Age of the Ganyanchi Pull‐Apart Basin

AbstractBecause the active, left‐slip Haiyuan fault is a first‐order structure along the NE margin of the Tibetan Plateau, its tectonic evolution provides insight into deformation processes along this margin over time. Late Cenozoic deformation in the area of the eastern Haiyuan fault initiated as thrust faulting, followed by left‐lateral strike slip displacement. However, the time of this kinematic change has been long debated. Here we use the structural evolution and the age of the Ganyanchi basin to date the onset of this kinematic transition. Seismic reflection data, integrated with published structural mapping, indicate that the basin is a pull‐apart controlled by strands of the eastern Haiyuan fault. Previously reported magnetostratigraphy of a core from the basin interior indicates deposition started at ∼2.8 Ma. However, this age likely postdates initiation of the bounding strike‐slip faults. We estimate that an additional 0.2–1.0 Myr of lateral slip occurred before the dated section formed. Thus, we find that the age of the eastern Haiyuan fault is at least ∼3.5 Ma, which is significantly older than early estimates of less than 2 Ma. Integrating our new data with prior work reveals that the kinematic change from NE‐SW shortening to left‐slip along the Haiyuan fault youngs to the east, which we interpret to result from growth of the Haiyuan strike‐slip fault via progressive lateral propagation from west to east. Miocene‐Pliocene onset of the Haiyuan strike‐slip fault coincides with initiation ages of other major strike‐slip faults in NE Tibet, likely implying that the NE margin of the Tibetan Plateau had switched to a strike‐slip dominated mode of deformation by this time.

Interparticle friction behaviors of kaolinite: Insights into macroscale friction from nanoscale

The friction behavior of widespread clay minerals is a major concern in many geo-engineering problems, such as the stability of soft soil foundations and the induction of seismic fault zones. The present work aimed to study the friction behaviors of kaolinite at the particle level using the molecular dynamics method. The effects of normal force (Fn), shear velocity (v), and interfacial water film on nanofriction were discussed. The “stick-slip” phenomenon and periodic evolution of friction forces (Ff) were observed in dry friction and became less pronounced with water lubrication. The dry Ff of kaolinite was found to be insensitive to Fn. However, wet Ff exhibited a linear increase with Fn and then transitioned to a non-linear relationship as slip displacement increases due to the continuous loss of water molecules from the interface during friction. Notably, at high loads (Fn ≥ 30 nN), the peak friction of wet kaolinite showed characteristics similar to dry friction. A velocity-strengthening behavior of kaolinite at high velocities was observed in both dry and wet conditions. The macroscale friction coefficients of kaolinite were predicted from nanofriction data and results showed good agreement with experimental values. This study lays the foundation for bridging micro- and macro-mechanical behaviors, suggesting a new pathway for acquiring precise macroscopic friction of minerals through cross-scale studies.

Tangential contact modeling to seal considering elastic-plastic for lubrication transition mechanism

Researches on the reciprocating seals lubrication under elastic-hydrodynamic lubrication (EHL) are mature, and leakage of failure criterion has been described with Couette flow for more than 30 years, which could not reflect to leakage mechanism in multiscale. Thus, it is necessary to explore seal leakage failure to focus on micro convex bodies of asperities model under contact pressure and cycle tangential friction stress in mixed lubrication. This is because convex bodies behaviors determine sealing effect and lubrication behaviors in microscale with subjecting to the dominant pressure. Interaction between convex bodies fatigue hysteresis curve and elastic-plastic strain are also calculated to explore sealing effect with leakage rate mechanism. We then design water seal at high pressure of 100 MPa with the speed of 0.1–1 m/s and medium temperature from 5 ℃ to 40 ℃. IWAN model of the convex bodies are calculated. Probability density function and Jenkin’s elements critical slip displacement are obtained. Convex bodies of radial and tangential deformation are also calculated through loading of contact pressure and friction stress with cycle strength coefficient to obtain seal fatigue lifespan. Ratio to seal stick and slip friction is calculated for evaluating the loading boundary. Seal leakage failure mechanism with fatigue characteristics is illustrated with macro lubrication parameters and micro convex bodies mechanics performance.

Dynamic reliability and seismic fragility analysis for high concrete-faced rockfill dam slopes subjected to stochastic earthquake and parameter excitation via PDEM

The randomness of material parameters and earthquake excitation greatly impacts the seismic stability of high concrete-faced rockfill dam (CFRD) slopes. However, no research has comprehensively analyzed the impacts of stochastic earthquake excitations, random parameters and their coupled effects on seismic stability of CFRD slopes and selected a multi-indicator evaluation criterion to carry out performance-based safety assessment. This paper develops a novel and comprehensive framework for evaluating the seismic stability of CFRD slopes based on the generalized probability density evolution method (GPDEM) with high efficiency and accuracy. The effects of three kinds of randomness on seismic stability of CFRD slopes were comprehensively compared and analyzed on the basis of multi-indices (safety factor (FS), cumulative slip displacement and cumulative time with FS < 1.0). Firstly, the nonlinear shear strength parameters of 40 high CFRDs were statistically collected to obtain the statistical characteristics of rockfill material strength parameters of actual CFRD projects. Secondly, three kinds of random samples were generated using generalized F-discrepancy (GF-discrepancy) method and Spectral expression-Random function (SERF) method. Then, the GPDEM was introduced to combine with three indices to perform the stochastic analysis and reliability evaluation. Finally, the fragility analysis of the CFRD slope was conducted. The results reveal that FS is more sensitive to the ground motions randomness, while cumulative slip displacement and cumulative time with FS < 1.0 are more sensitive to the material parameters randomness. The seismic safety evaluation of CFRD slopes based on multiple indices with full consideration of coupled randomness effects is necessary.

Slip‐Dependence of Fault Frictional Stability Under Hydrothermal Conditions

AbstractIn the rate‐state friction law framework, the transition from velocity weakening (V‐W) to velocity strengthening (V‐S) behavior marks the base of the seismogenic crust. Here we investigate the role of fault slip displacement under hydrothermal conditions in controlling the V‐W to V‐S transition. We shear simulated gabbro gouges at slip velocities ranging from 16 nm/s (∼50 cm/year) to 10 μm/s (∼8 cm/day) under hydrothermal conditions (300–400°C temperature; 30 MPa pore fluid pressure). We observe that cumulative fault slip increases the critical velocity for the V‐W to V‐S transition. The transition is accompanied by localized to distributed deformation mode, the formation of smectite‐type clays and occurrence of intergranular mass transfer. Our results provide insights into understanding the deepening and shallowing of V‐W/V‐S boundary (lower limit of the seismogenic zone) following a mainshock. Besides strain rate effects, slip‐enhanced chemical alteration and grain size‐sensitive deformation may temporarily contribute to the shallowing process.

Slope displacement reliability analysis considering rock parameters spatial variability subjected to stochastic mainshock-aftershock earthquake

A new and efficient framework for analyzing the stochastic dynamic response and reliability of the rock slope was constructed based on the generalized probability density evolution method (GPDEM). The left bank slope of Jinping dam was taken as an illustrative example and the effects of mainshock-aftershock sequence randomness, spatial variability of geotechnical materials and coupled uncertainties on the slope dynamic response and reliability were compared and analyzed for the first time. Firstly, a series of mainshock-aftershock sequences and parameter random fields were generated as stochastic factors inputs. Then, the effectiveness of GPDEM was validated by comparing with Monte Carlo Simulation (MCS). The Newmark slip displacement method was combined with GPDEM to analyze the slope dynamic response and reliability subjected to single and coupled uncertainties from a probabilistic viewpoint. Finally, the fragility curves of three risk levels were established based on the reliabilities with different seismic intensities. The result suggests that parameter spatial variability affects the dynamic response and reliability of the slope to some extent, but this effect will be overshadowed by the uncertainty of mainshock-aftershock sequences when it is considered simultaneously. Besides, aftershocks will cause the decrease of slope reliability and the impact of aftershock cannot be overlooked.

Study of the UHPC–NC Interfacial Bonding Properties of Steel Tubes with Internally Welded Reinforcement Rings

Ultra-high-performance concrete (UHPC) has the advantages of high strength, excellent durability performance, etc., and its compressive strength is several times that of normal concrete (NC). Due to the role of materials such as steel fibers, the tensile strength of UHPC is higher than that of NC. For steel-tube concrete columns in corrosive seawater environments, UHPC–NC columns with welded ring reinforcement inside the steel tube are proposed to strengthen the interfacial bonding performance, and the effects of seawater corrosion of steel-tube concrete are studied. Eight steel-tube UHPC–NC specimens were designed for push-out tests. The steel tubes were internally constructed with glossy unconstructed and reinforcing rings, with the core concrete with UHPC used below and the C40 plain concrete used above. By examining push-out load, slip displacements, and steel-tube wall strains, this study analyzed the influence of different factors on the bond behavior and failure mechanism of bond-slip in shear-resistant reinforcing ring connectors. The push-out simulation of the steel-tube concrete was carried out using ABAQUS 2021 software, and the simulation results were compared with the experimental results, which showed good agreement. The results show that the bond strength of the steel tube–concrete column interface can be significantly improved by using the construction measures of internally welded reinforcement rings; for specimens with the same percentage of core concrete UHPC and C40 thickness, the bond strength of the two rings was significantly improved by approximately 33% over that of the one-ring reinforcement ring; corrosive environments will degrade the bond strength of the steel tube–concrete column interface.

Numerical investigation of the behavior of combination connections with pretensioned high-strength bolts and longitudinal fillet welds

Numerical analysis is conducted to investigate the behavior of double shear combination connections made with pretensioned high-strength bolts and longitudinal fillet welds. The finite element (FE) models were validated using the results of an experimental investigation and utilized to develop a better understanding of the influence of critical input parameters on the connection behavior. Slip-dependent surface frictional and weld ductile fracture models were incorporated into the numerical models to simulate the nonlinear behavior of the connections. The numerical simulations confirm that at low slip displacements, the force carried by the combination connection can be estimated as the summation of the friction resistance introduced by the pretensioned bolts and the shear force carried by the fillet welds. The results also indicate that the behavior and capacity of the combination connection are not sensitive with respect to the weld location compared to the center of gravity of the bolts as long as no eccentricity is introduced within the force transfer mechanism. The paper also investigates the stress distribution within the connection plates and studies the effect of fillet weld size on the behavior of combination connections.

Kinematic Analysis of the Nain and Darrech Deh Faults and Their Structural Position in the Western Side of the Central Iranian Range

The western margin of Central Iran – the Nain area – is composed of the Mesozoic ophiolitic mélanges deposited during the Late Cretaceous. The Nain and Darrech Deh faults whose trace indicators can be seen in the study area are among the largest faults therein. These faults have a north-northwest strike and are the right-lateral transpressional active lineaments. The main aim of this research is to perform structural mapping and geometric analysis of structural features along the Nain and Darrech Deh fault system. The research was conducted using satellite imagery and field data. The structural analysis showed that the faults comprise a complex network which formed during two stages: 1) thrust faulting in the Late Cretaceous and 2) reactivation of thrust faults in the form of the right-lateral reverse-slip structures as a result of transpressional movements in the Upper Miocene. The Nain and Darrech Deh faults are presently active and exhibit the right-lateral strike slip displacement with a reverse component.

Effect of mineral dissolution on fault slip behavior during geological carbon storage

Geological carbon dioxide storage is one of the most effective technologies to reduce greenhouse gas emissions to the atmosphere. The mineral dissolution is inevitable during the long-term storage process, while the chemical effect on fault slip behavior remains poorly understood. In this paper, we present a hydro-chemical–mechanical coupled fault reactivation study using the unified pipe-interface element method (UP-IEM). The impact of chemical dissolution on stress redistribution is considered. The accuracy of the numerical method to simulate reactive solute transport and fault frictional slip are confirmed by analytical solution and numerical verification. Two different fault tectonic scenarios, of which fault locating in the caprock and fault crosscutting the reservoir and caprock, are designed to investigate the fault slip behavior induced by mineral dissolution. The effect of dissolved CO2 solute concentration, fault aperture, reservoir permeability and chemical reaction rate are discussed. Results show that fault locating in the caprock is easier to be activated. The chemical reaction has a more significant effect on fault slip behavior than fluid pressure during the long-term low-velocity CO2 leakage. The increase of CO2 solute concentration and reservoir reaction rate cause an obvious growth on fault slip displacement. The decrease of fault initial aperture results in the fault slip velocity slower and the slipping area smaller. The high-permeability reservoir keeps the fault more stable than low-permeability reservoir under the influence of mineral dissolution.

Landslides triggered by the 2022 Ms. 6.8 Luding strike-slip earthquake: An update

In this study, we established an updated coseismic landslide inventory of the 2022 Ms. 6.8 Luding earthquake. This version contains over 16,000 hand-digitized landslides with a total landslide area of 43.2 km2. We compared the spatial distribution of landslides and topographic, seismic, and geological factors on both sides of the seismogenic fault. The results show that the coseismic landslides of the Luding earthquake were mainly distributed along both sides (SWW and NEE blocks) of the seismogenic fault, with more landslides concentrated on the west side (SWW block). Compared to the NEE block, landslides on the SWW block are more likely to occur in areas characterized by higher elevations and reliefs. Otherwise, we observe an exponential decrease in both landslide area density and number density as the distance from the seismogenic fault increases, and the proximity to the fault is positively correlated with landslide size. Finally, we explored the relationship between landslide distribution and slip patterns by examining the detailed rupture characteristics of the Luding earthquake. Our analysis indicates that density patterns of coseismic landslides closely track the distribution of fault slip. The landslide abundance area corresponds to the area with the largest slip displacement of the surface rupture. The more the rupture extends into the shallow portion of the fault plane, the more landslides occur.

Temporal evolution of shear-induced dilatancy of rock fractures: controls from surface roughness and normal stress

SUMMARY Understanding the shear-induced dilatancy of rock fractures is important for assessing the permeability evolution and seismic hazard in shale and geothermal reservoirs. The displacement dependence of fracture dilation has been well studied, while the influence of slip velocity is poorly constrained. In this study, we combined displacement- and velocity-dependent aperture models to reproduce the transient shear-induced dilatancy of fractures in sandstone in 16 normal stress unloading tests. Our results show that the combined aperture model can describe the transient fracture aperture evolution during accelerating slip induced by normal stress unloading better than the model dependent only on slip displacement. Slip velocity could enhance the aperture increase on smoother fractures at lower normal stresses and higher slip velocities. Both the dilation factor and characteristic slip distance decrease with increasing normal stress and surface roughness, signifying reduced contribution of slip velocity to transient shear-induced dilatancy at higher normal stresses and surface roughness. The dilation angle increases with the increase of surface roughness, and this increase diminishes at higher normal stresses primarily attributable to more severe asperity wear. These findings highlight the importance of slip velocity in controlling the transient evolution of aperture and permeability of a rock fracture. Our study also provides constraints on the constitutive parameters in the combined aperture model for describing transient shear-induced fracture dilatancy. We suggest that it is crucial to incorporate the velocity-dependent aperture model to simulate the nonlinear evolution of fracture aperture in future analytical and numerical models involving coupled hydromechanical processes in geoenergy systems.

Steadystate: A MATLAB-based routine for determining steady-state friction conditions in the framework of rate-and-state friction analysis

Reliable determination of rate-and-state friction (RSF) parameters depends on achieving steady-state (SS) friction conditions before and after experimental velocity-stepping friction tests. This operation, through nonlinear least squares fitting, is commonly preceded by the removal of any overall slip weakening/hardening after friction velocity steps (VSs) through a sufficiently large window of slip displacement at SS ( = linear detrend). However, to date, the identification of SS and thus the correct linear detrend is dependent on the user, which potentially results in differing RSF outputs from the same data set. Here, we demonstrate that large errors in the determination of the fitted RSF parameters can result if SS conditions are not reached before and after VSs. Such errors can be particularly relevant for materials characterized by long evolution of frictional resistance with slip, such as clay-rich gouge layers, in which identifying SS after VSs is not always obvious. To this end, we propose a methodology to accurately and consistently identify where SS is achieved after VSs. This methodology is coded into a new MATLAB-based routine, steadystate. We show the key features of the methodology, as well as how to use steadystate and read its output. We also illustrate the broad applicability of the approach to friction data with different noise levels and sampling frequencies referenced to slip velocity, by reviewing observations from synthetic data sets and specific examples of experiments from different laboratories involving various sheared materials.

Numerical Investigation of Stratified Slope Failure Containing Rough Non-Persistent Joints Based on Distinct Element Method

To address the critical issue of slope stability in jointed rock masses with complex and rough structural planes, a rough joint network model using the Fourier transform was proposed and applied to the Shilu open pit mine. The on-site structural plane survey results were combined with MATLAB and PFC2D to establish numerical models for slope stability analysis considering both linear-jointed and rough-jointed rock slopes. By comparing the slip body and fracture distribution, it was found that the rough-jointed slope was stabler than the linear-jointed slope. In addition, the fracture patterns and slip displacement were significantly influenced by the inclination and undulation of the bedding planes. Slip failure patterns occurred when the angle of inclination was set at 60°. The joints played a crucial role in inducing the shear strength of slope rock masses, and the slide area was mainly observed in the shallow slope surface for inclination angles of 0° and 45°, and in the middle deep area for 60° and 90°. These results demonstrated the importance of considering rough non-persistent fractures when developing a new numerical model for slope failure modes.

Assessment of critical parameters affecting fretting fatigue life of bridge stay cables at saddle supports

Saddle systems are a popular solution for supporting cables in cable-stayed bridges. Fretting fatigue failure of cables at saddle supports is a primary design consideration for these systems. Critical parameters that affect the fretting fatigue life of the cables are the contact forces and the slip displacements at the contact points between the cables and the saddle. Determining these critical parameters is the first step in evaluating the fretting fatigue life of the cables. Wear in the contact points between the cables and the saddle can affect the load distribution in saddle systems and consequently affect these critical parameters. However, the effect of wear has not been studied in the previous works. This paper first discusses the methods proposed in the literature for evaluating contact forces and slip displacements in the contact points between a cable and a saddle. Then, a finite element model of the problem and a framework for modelling the wear are presented. Finally, fretting fatigue life is determined based on the different studied approaches. The main highlights of the current study are considering the effect of wear in the simulation and employing an enhanced FE model for slip displacement calculations. The results of the basic model without wear effects showed that the first contact point between the cable and the saddle is critical for fatigue failure. However, by incorporating wear in the model, the contact force at the first point dropped and the second contact point became critical; this is in line with the observations in large-scale fatigue tests of saddle systems.

The relationship between kinematics and fault geometry for surface coseismic ruptures on across-strike faults: New observations of slip vectors and displacements along the Pisia and Skinos faults from the 1981 Eastern Gulf of Corinth, Greece earthquakes

The relationships between kinematics and fault geometry for the coseismic ruptures from the 24th and February 25, 1981 earthquake sequence in the eastern Gulf of Corinth (Ms 6.7 and 6.4) are analysed. The two earthquakes ruptured faults located across strike rather than along strike as typifies other earthquake sequences. In detail, surface ruptures formed on the sub-parallel Pisia and Skinos Faults, with an 8 km along-strike overlap zone, separated across strike by < 2 km. The largest coseismic offsets occurred in the overlap zone. The 41-year-old ruptures are still well preserved as bedrock fault plane lichen-free stripes and colluvial ruptures, allowing detailed structural mapping at 213 rupture localities. A comparison between our measurements and Jackson et al. (1982) showed no overall consistent signal of post-seismic slip as some of our measurements were greater and some smaller than those recorded in 1981. The ruptures produced a single maximum asymmetric profile (Pisia: maximum throw of 223 cm) and a double maxima profile (Skinos: maximum throw of 109 cm and 130 cm). The shapes of the profiles differed in previous earthquakes on these faults, as evidenced by an older lichen-free stripe, implying non-characteristic earthquakes. Summing the two overlapping throw profiles across-strike reveals a single maximum symmetric bell-like profile. Using the above observations on coseismic offsets, kinematic information, and the geometry of faults, a rupture scenario has been proposed in terms of fault bends and corrugation orientations which suggests that parts of each fault may have ruptured in each earthquake.