The importance of the effective intermediate principal stress (σ′2) to fault slip patterns

In normal faulting regimes, the magnitudes and orientations of the maximum and minimum principal compressive stresses may be known with some confidence. However, the magnitude of the intermediate principal compressive stress is generally much more difficult to constrain and is often not considered to be an important factor. In this paper, we show that the slip characteristics of faults and fractures with complex or nonoptimal geometry are highly sensitive to variation or uncertainty in the ambient effective intermediate principal stress (σ′₂). Optimally oriented faults and fractures may be less sensitive to such variations or uncertainties. Slip tendency (T_s) analysis provides a basis for quantifying the effects of uncertainty in the magnitudes and orientations of all principal stresses and in any stress regime, thereby focusing efforts on the most important components of the system. We also show, for a normal faulting stress regime, that the proportion of potential surfaces experiencing high slip tendency (e.g., T_s ≥ 0.6) decreases from a maximum of about 38% where σ′₂ = σ′₃, to a minimum of approximately 14% where σ′₂ is halfway between σ′₃ and σ′₁, and increases to another high of approximately 29% where σ′₂ = σ′₁. This analysis illustrates the influence of the magnitude of σ′₂ on rock mass strength, an observation previously documented by experimental rock deformation studies. Because of the link between fault and fracture slip characteristics and transmissivity in critically stressed rock, this analysis can provide new insights into stress-controlled fault transmissivity.     

An inverse method to determine the optimal stress from imperfect fault data

An inverse method is developed to determine the optimal stress inversion from imperfect fault data. It is specifically designed to process fault data at variable division—either with observed slip directions or with observed fault/slip senses or with both. The method has the flexibility of processing data with a few fault/slip senses incorrectly determined—a case we occasionally meet in the field. A vast number of artificial 10-datum examples are processed to illustrate the reliability of the method and to examine the effects of the data set at variable division on stress estimation. The results show that the precision of solution decreases with the increase in either of two factors, the percentage of the data limited to slip sense and the range of measurement error, both resulting in a broader dispersion of estimated stress vectors relative to the assigned stress vector in the sigma space. Thus, in order to keep the dispersion within a range of 5° requires the percentage of data for which only slip sense is known to be less than about 20%. This in turn allows for a measurement error range in the data of 10° in real space, much larger than real measurement errors. We believe that the proposed method would be of practical value for stress inversion of fault/slip data.     

Constraining slip rates and spacings for active normal faults

Numerous observations of extensional provinces indicate that neighbouring faults commonly slip at different rates and, moreover, may be active over different time intervals. These published observations include variations in slip rate measured along-strike of a fault array or fault zone, as well as significant across-strike differences in the timing and rates of movement on faults that have a similar orientation with respect to the regional stress field. Here we review published examples from the western USA, the North Sea, and central Greece, and present new data from the Italian Apennines that support the idea that such variations are systematic and thus to some extent predictable. The basis for the prediction is that: (1) the way in which a fault grows is fundamentally controlled by the ratio of maximum displacement to length, and (2) the regional strain rate must remain approximately constant through time. We show how data on fault lengths and displacements can be used to model the observed patterns of long-term slip rate where measured values are sparse. Specifically, we estimate the magnitude of spatial variation in slip rate along-strike and relate it to the across-strike spacing between active faults.     

Comminution and fluidization of granular fault materials: implications for fault slip behavior

Whilst faulting in the shallow crust is inevitably associated with comminution of rocks, the mechanical properties of the comminuted granular materials themselves affect the slip behavior of faults. Therefore, the mechanical behavior of any fault progresses along an evolutionary path. We analyzed granular fault rocks from four faults, and deduced an evolutionary trend of fractal size frequency. Comminution of fault rocks starts at a fractal dimension close to 1.5 (2-D measurement), at which a given grain is supported by the maximum number of grains attainable and hence is at its strongest. As comminution proceeds, the fractal dimension increases, and hence comminution itself is a slip weakening mechanism. Under the appropriate conditions, comminuted granular materials may be fluidized during seismic slip events. In this paper, we develop a new method to identify the granular fault rocks that have experienced fluidization, where the detection probability of fragmented counterparts is a key parameter. This method was applied to four fault rock samples and a successful result was obtained. Knowledge from powder technology teaches us that the volume fraction of grains normalized by maximum volume fraction attainable is the most important parameter for dynamic properties of granular materials, and once granular fault materials are fluidized, the fault plane becomes nearly frictionless. A small decrease in the normalized volume fraction of grains from 1 is a necessary condition for the phase transition to fluidization from the deformation mechanism governed by grain friction and crushing by contact stresses. This condition can be realized only when shearing proceeds under unconstrained conditions, and this demands that the gap between fault walls is widened. Normal interface vibration proposed by Brune et al. [Tectonophysics 218 (1993) 59] appears to be the most appropriate cause of this, and we presented two lines of field evidence that support this mechanism to work in nature.     

Permeability upscaling of fault zones in the Aztec Sandstone,Valley of Fire State Park,Nevada, with a focus on slip surfaces and slip bands

The effect of open and filled slip surfaces on the upscaled permeability of two fault zones with 6 and 14 m strike-slip in an eolian Aztec Sandstone, Nevada, USA is evaluated. Each fault zone is composed of several fault components: a fault core, bounded by filled through-going slip surfaces referred to as slip bands, and a surrounding damage zone that contains joints and deformation bands. Slip band geometry, composition, and petrophysical properties are characterized. Measurements and modeling show that slip band permeabilities can vary over 12 orders of magnitude depending on the degree of fill within the slip bands. The slip bands along with other fault zone components are represented in finite volume numerical calculations and the impact of various slip-band representations on upscaled fault zone permeability is tested. The results show 2 orders of magnitude variation in upscaled fault zone permeability in the fault-normal direction and a factor of 2 variation in the fault-parallel direction. The numerical results presented here are compared to the earlier numerical results in which structured Cartesian grids were used for the numerical simulations, and are in qualitative agreement with earlier calculations but use about a factor of 250–400 fewer numerical cells.     

The Gubbio normal fault (Central Italy): geometry, displacement distribution and tectonic evolution

Normal faults within orogenic belts can be pre-, syn- or post-orogenic features. We studied the Gubbio normal fault (central Italy), which is an example of a pre-orogenic fault reactivated in a post-orogenic stage. The Gubbio Fault is a 22-km-long fault bordering a Quaternary basin and part of an active faults system in the Umbria–Marche region (Central Italy). The interpretation of a set of seismic profiles enables us to reconstruct the fault geometry in detail and to measure displacement and throw distributions along the fault strike. Seismic data indicate that the Gubbio Fault represents an example of multiple reactivation: at least a portion of the fault was active in the Miocene and only a part of the total displacement was achieved in the Quaternary. The reconstruction of the fault geometry at depth shows that the fault is characterised by listric geometry. The fault is also characterised by a bend along strike and structure contours show that this geometry is maintained at depth. As the fault is commonly addressed as presently active, the maximum fault dimensions are correlated to the maximum expected earthquake, and the presence of the fault bend is discussed as a possible barrier to seismic ruptures propagation.     

Finite-element 3D modeling of stress patterns around a dipping fault

Stresses in a block around a dipping fracture simulating a damage zone of a fault are reconstructed by finite-element modeling. A fracture corresponding to a fault of different lengths, with its plane dipping at different angles, is assumed to follow a lithological interface and to experience either compression or shear. The stress associated with the destruction shows an asymmetrical pattern with different distances from the highest stress sites to the fault plane in the hanging and foot walls. As the dip angle decreases,the high-stress zone becomes wider in the hanging wall but its width changes negligibly in the foot wall.The length of the simulated fault and the deformation type affect only the magnitude of maximum stress,which remains asymmetrical relative to the fault plane. The Lh/Lfratio, where Lhand Lfare the widths of high-stress zones in the hanging and foot walls of the fault, respectively, is inversely proportional to the fault plane dip. The arithmetic mean of this ratio over different fault lengths in fractures subject to compression changes from 0.29 at a dip of 80°to 1.67 at 30°. In the case of shift displacement, ratios are increasing to 1.2 and 2.94, respectively.Usually they consider vertical fault planes and symmetry in a damage zone of faults. Following that assumption may cause errors in reconstructions of stress and fault patterns in areas of complex structural setting. According geological data, we know the structures are different and asymmetric in hanging and foot walls of fault. Thus, it is important to quantify zones of that asymmetry. The modeling results have to be taken into account in studies of natural faults, especially for practical applications in seismic risk mapping, engineering geology, hydrogeology, and tectonics.     

四川汶川5.12大地震同震滑动断层泥的发现及意义

2008年汶川8.0级地震沿龙门山断裂带内的映秀—北川断裂和灌县—安县断裂产生了近300 km的同震地表破裂带。震后地质科学考察发现地表变形以逆冲为主,并伴有右旋走滑。地震地表破裂带大多沿古生代碳质泥岩、页岩和三叠系煤系地层内的滑动面出露地表,这些软弱地层为地震破裂带冲到地表提供了超低摩擦滑动带。我们发现在同震垂直和水平位错达6m左右的地表破裂带,地震的同震滑动发生在厚度约0.5~2cm 的狭窄滑动带内,以发育新鲜的灰色断层泥为特征,这些断层泥是地震断层快速滑动过程中岩石—流体相互作用的结果。     

四川汶川5.12大地震同震滑动断层泥的发现及构造意义

2008年汶川8．0级地震沿龙门山断裂带内的映秀-北川断裂和灌县-安县断裂产生了近300km的同震地表破裂带。震后地质科学考察发现地表变形以逆冲为主,并伴有右旋走滑。地震地表破裂带大多沿古生代碳质泥岩、页岩和三叠系煤系地层内的滑动面出露地表,这些软弱地层为地震破裂带冲到地表提供了超低摩擦滑动带。我们发现在同震垂直和水平位错达6m左右的地袁破裂带,地震的同震滑动发生在厚度约O．5—2cm的狭窄滑动带内,以发育新鲜的灰色断层泥为特征,这些断层泥是地震断层快速滑动过程中岩石-流体相互作用的结果。     

Weakness and mechanical anisotropy of phyllosilicate-rich cataclasites developed after mylonites of a low-angle normal fault (Simplon Line,Western Alps)

The Simplon Fault Zone is a late-collisional low-angle normal fault (LANF) of the Western Alps. The hanging wall shows evidence of brittle deformation only, while the footwall is characterized by a c. 1 km-thick shear zone (the Simplon Fault Zone), which continuously evolved, during exhumation and cooling, from amphibolite facies conditions to brittle-cataclastic deformations. Due to progressive localization of the active section of the shear zone, the thermal-rheological evolution of the footwall resulted in a layered structure, with higher temperature mylonites preserved at the periphery of the shear zone, and cataclasites occurring at the core (indicated as the Simplon Line). In order to investigate the weakness of the Simplon Line, we studied the evolution of brittle/cataclastic fault rocks, from nucleation to the most mature ones. Cataclasites are superposed on greenschist facies mylonites, and their nucleation can be studied at the periphery of the brittle fault zone. This is characterized by fractures, micro-faults and foliated ultracataclasite seams that develop along the mylonitic SCC′ fabric, exploiting the weak phases mainly represented by muscovite and chlorite. Approaching the fault core, both the thickness and frequency of cataclasite horizons increase, and, as their thickness increases, they become less and less foliated. The fault core itself is represented by a thicker non-foliated cataclasite horizon. No Andersonian faults or fractures can be found in the footwall damage zone and core zone, whilst they are present in the hanging wall and in the footwall further from the fault. Applying a stress model based on slip tendency, we have been able to calculate that the friction coefficient of the Simplon Line cataclasites was <0.25, hence this fault zone is absolutely weak. In contrast with other fault zones, the weakening effect of fluids was of secondary importance, since they accessed the fault zone only after an interconnected fracture network developed exploiting the cataclasite network.     

A grid-search inversion method looking for the best classification of polyphase fault/slip data

A novel, grid search-based stress inversion method is developed in this paper to find the global minima in the solution region for the classification of fault/slip data into many single-phase subsets. Exhaustively repetitious grid searches are taken to deal with possible local minima, in a departure from existing grid search-based inversion methods. Two stopping rules, to stop at the abrupt change of the objective function or at the least change of the classification, are adopted in the method to look for the best classification. Much calculation time is saved by using a modified version of conventional grid search. The feasibility of this proposed method is demonstrated by applying it to two artificial examples and two real examples. However, enormous time in calculation is still needed in the case of a data set either with a large number of data or for a large number of assigned subsets.     

Effect of sampling and linkage on fault length and length–displacement relationship

The power-law exponent (n) in the equation: D=cL ⁿ , with D = maximum displacement and L = fault length, would be affected by deviations of fault trace length. (1) Assuming n=1, numerical simulations on the effect of sampling and linkage on fault length and length–displacement relationship are done in this paper. The results show that: (a) uniform relative deviations, which means all faults within a dataset have the same relative deviation, do not affect the value of n; (b) deviations of the fault length due to unresolved fault tip decrease the values of n and the deviations of n increase with the increasing length deviations; (c) fault linkage and observed dimensions either increase or decrease the value of n depending on the distribution of deviations within a dataset; (d) mixed deviations of the fault lengths are either negative or positive and cause the values of n to either decrease or increase; (e) a dataset combined from two or more datasets with different values of c and orders of magnitude also cause the values of n to deviate. (2) Data including 19 datasets and spanning more than eight orders of fault length magnitudes (10⁻²–10⁵ m) collected from the published literature indicate that the values of n range from 0.55 to 1.5, the average value being 1.0813, and the peak value of n _d (double regression) is 1.0–1.1. Based on above results from the simulations and published data, we propose that the relationship between the maximum displacement and fault length in a single tectonic environment with uniform mechanical properties is linear, and the value of n deviated from 1 is mainly caused by the sampling and linkage effects.     

Modal and vectorial analysis for determination of stress axes associated with fault slip data

Based on Anderson's faulting theory, for a given set of fault slip data including the sense of slip on faults and fault planes, this paper provides two possible methods to reconstruct the principal stress axes using vectorial and modal analysis procedure. The vectorial analysis consists of computing eigenvectors of the orientation matrices defined by axes P, B, and T (P, B, and T being geometrical axes parallel to axes ₁, ₂, ₃ associated with a single striated fault) which are determined geometrically knowing the slip vector S, the normal to fault plane N, and the dihedral angle 2. A parameter R, related to the maximum eigenvalue of the orientation matrix and the size of data sample, is shown to be a good test value for the homogeneity of the data. A process of refinement of this parameter enables bad data (representing faults generated under different tectonic events) to be ignored in the final computation of the principal stress axes; results are thus significantly improved and the vectorial analysis procedure enables the orientation of stress to be numerically determined. On defining probability density regions around each P, B, and T axes that incorporate the variation of fault geometry, the modal analysis is carried out in two density functions established by a convolution process to locate modes which represent optimal locations of the principal stress axes in spherical space. Splitting of heterogeneous data samples into homogeneous subsets is achieved by a dynamical cluster procedure which enables the principal stress axes associated with each subset to be determined separately.     

Structural control of geothermal reservoirs in extensional tectonic settings: An example from the Upper Rhine Graben

In extensional tectonic settings major structural elements such as graben boundary faults are typically oriented subparallel to the maximum horizontal stress component S_Hmax. They are often structurally accompanied by transfer zones that trend subparallel to the extension direction. In the Upper Rhine Graben, such transfer faults are typically characterized by strike-slip or oblique-slip kinematics. A major re-orientation of the regional stress field by up to 90° of the Upper Rhine Graben in the Early Miocene led to the present-day normal and strike-slip faulting regimes in the North and South of the Upper Rhine Graben, respectively, and a transition zone in-between. Consequently, conditions for fault frictional failure changed significantly. Moreover, it has been observed during tracer and stimulation experiments that such transfer faults may be of major importance for the hydraulic field of geothermal reservoirs under the present stress condition, especially, when located between production and injection well.In this context we have investigated slip and dilation tendencies (T_S and T_D) of major structural elements at reservoir scale for two representative geothermal sites, Bruchsal (Germany) and Riehen (Switzerland), located close to the Eastern Main Boundary Fault of the Upper Rhine Graben. We have evaluated the quality and uncertainty range of both tendencies with respect to potential variation in S_Hmax orientation. Despite significant differences in orientation of the structures and the stress regimes, the resulting variation of T_S and T_D reveal major similarities concerning the reactivation potential of both, the graben-parallel structures and the transfer faults. The conditions of criticality for tensile failure and non-criticality for shear failure suggest that transfer faults are most likely naturally permeable structures with low stimulation potential. This is in agreement with the absence of both immediate tracer recovery and seismicity in the studied geothermal sites.     

Time dependent structural architecture of subsidiary fracturing and stress pattern in the tip region of an extensional growth fault system, Tarquinia basin, Italy

Fault tip regions, relay ramps and accommodation zones in between major segments of extensional fault systems provide zones of additional structural and stratigraphic complexities and also significantly affect their hydraulic behaviour. The great interest for both academic and industrial purposes encouraged specific studies of fault tip regions that, in some cases, produced controversial results. We approached the study of fault tip regions by integrating structural, AMS and stratigraphic analyses of the tip of an extensional growth fault system in the Tarquinia basin, on the Tyrrhenian side of the Northern Apennines. Detailed structural mapping indicates the occurrence of systematic relationships between the orientation of the main subsidiary fault zones, the orientation and position of the two main joint sets developed in the fault damage zones, and the overprinting relationships between the two main joint sets themselves. Microstructural analysis of fault core rocks indicates a progression of deformation from soft-sediment to brittle conditions. The AMS study supports the evolution of deformation under a constantly oriented stress field. By combining this multidisciplinary information we propose an evolutionary model for the tip of the extensional growth fault system that accounts for the progressively changing sediment rheological properties, and for the time dependent subsidiary deformation pattern by invoking the interplay between the regional stress field and the local, kinematically-derived one by fault activity. We also speculate on the overall implications for fluid flow of the proposed evolutionary model.     

Experimental frictional heating of coal gouge at seismic slip rates: Evidence for devolatilization and thermal pressurization of gouge fluids

High velocity (1 m/s) friction experiments on bituminous coal gouge display several earthquake-related phenomena, including devolatilization by frictional heating, gas pressurization, and slip weakening. Stage I is characterized by sample shortening and reduction in the coefficient of friction (μ) from  1 to 0.6. Stage II is characterized by high frequency ( 5 Hz) oscillations in stress and strain records and by gas emissions. Stage III is marked by rapid weakening (μ  0.1 to 0.35) and sample shortening, together with continued gas emissions. Stage IV produces stable stress records and continued weakness (μ  0.2), but without gas emission. Stage I shortening is due to compaction of the gouge and the weakening is attributed to mechanical or thermal effects. Stage II behavior is interpreted as due to coal gasification and fluctuations in fluid pressure, resulting in high frequency stick-slip type behavior. Dramatic reduction in shear stress in stage III is attributed to gas pressurization by pore collapse and corresponds to a frictional instability, analogous to nucleation of an earthquake. Microstructural observations indicate the deformation was brittle during stages I and II but ductile during stages III and IV. Time dependent finite element frictional heat models indicate the center of the samples became hot ( 900 °C) during stage II, whereas the edge of samples remained relatively cold (< 300 °C). Vitrinite reflectance of coal samples shows an increase in reflectance from  0.5 to  0.8% over the displacement interval 20–40 m (20–40 s), indicating that the reflectance responds to frictional heating on a short time scale. The energy expended per unit area in these low stress, large displacement experiments is similar to that of higher stress ( 50 MPa), short displacement ( 1 m) earthquakes ( 10⁷ J/m²).     

Neotectonics of the Marikina Valley fault system (MVFS) and tectonic framework of structures in northern and central Luzon, Philippines

Recognition of neotectonic features along the Marikina Valley fault system (MVFS) in central Luzon, Philippines indicates a dominantly dextral strike-slip motion during its most recent activity believed to be Late Pleistocene to Holocene in age. Variations in the ratios of vertical to horizontal displacements for the segments imply a dominantly dextral motion of the West Marikina Valley fault (WMVF) and oblique dextral motion for the East Marikina Valley fault (EMVF). The displacement data further suggest that rupturing along the EMVF involved multiple segments and occurred separately from the events along the WMVF segments. Estimated earthquake magnitudes for the WMVF and EMVF based on single-event offsets fall within the range M 7.3–7.7. The vertical slip component in the northern part of the Marikina Valley is associated with the development of a basin between the EMVF and WMVF while the large vertical component in the southernmost segment of the EMVF (Talim) is attributed to volcanism-related extension. Lateral advection of the block bounded by the MVFS and the Philippine fault zone (PFZ), rather than pure shear resulting from an assumed east–west compression, best explains the observed kinematics of the MVFS. This is the result of compression during the westward drift of the Philippine Sea Plate and northern Luzon and occurs through slip along the WMVF and EMVF at rates of 5–7 mm/yr.     

Fault growth and interactions in a multiphase rift fault network: Horda Platform,Norwegian North Sea

Physical models predict that multiphase rifts that experience a change in extension direction between stretching phases will typically develop non-colinear normal fault sets. Furthermore, multiphase rifts will display a greater frequency and range of styles of fault interactions than single-phase rifts. Although these physical models have yielded useful information on the evolution of fault networks in map view, the true 3D geometry of the faults and associated interactions are poorly understood. Here, we use an integrated 3D seismic reflection and borehole dataset to examine a range of fault interactions that occur in a natural multiphase fault network in the northern Horda Platform, northern North Sea. In particular we aim to: i) determine the range of styles of fault interaction that occur between non-colinear faults; ii) examine the typical geometries and throw patterns associated with each of these different styles; and iii) highlight the differences between single-phase and multiphase rift fault networks. Our study focuses on a ca. 350 km² region around the >60 km long, N–S-striking Tusse Fault, a normal fault system that was active in the Permian–Triassic and again in the Late Jurassic-to-Early Cretaceous. The Tusse Fault is one of a series of large (>1500 m throw) N–S-striking faults forming part of the northern Horda Platform fault network, which includes numerous smaller (2–10 km long), lower throw (<100 m), predominantly NW–SE-striking faults that were only active during the Late Jurassic to Early Cretaceous. We examine how the 2^nd-stage NW–SE-striking faults grew, interacted and linked with the N–S-striking Tusse Fault, documenting a range of interaction styles including mechanical and kinematic isolation, abutment, retardation and reactivated relays. Our results demonstrate that: i) isolated, and abutting interactions are the most common fault interaction styles in the northern Horda Platform; ii) pre-existing faults can act as sites of nucleation for 2^nd-stage faults or may form mechanical barriers to propagation; iii) the throw distribution on reactivated 1^st-stage faults will be modified in a predictable manner if they are intersected or influenced by 2^nd-stage faults; iv) sites of fault linkage and relay-breaching associated with the first phase of extension can act as preferential nucleation sites for 2^nd-stage faults; and v) the development of fault intersections is a dynamic process, involving the gradual transition from one style to another.     

Geometry and architecture of faults in a syn-rift normal fault array: The Nukhul half-graben, Suez rift, Egypt

The geometry and architecture of a well exposed syn-rift normal fault array in the Suez rift is examined. At pre-rift level, the Nukhul fault consists of a single zone of intense deformation up to 10 m wide, with a significant monocline in the hanging wall and much more limited folding in the footwall. At syn-rift level, the fault zone is characterised by a single discrete fault zone less than 2 m wide, with damage zone faults up to approximately 200 m into the hanging wall, and with no significant monocline developed. The evolution of the fault from a buried structure with associated fault-propagation folding, to a surface-breaking structure with associated surface faulting, has led to enhanced bedding-parallel slip at lower levels that is absent at higher levels. Strain is enhanced at breached relay ramps and bends inherited from pre-existing structures that were reactivated during rifting. Damage zone faults observed within the pre-rift show ramp-flat geometries associated with contrast in competency of the layers cut and commonly contain zones of scaly shale or clay smear. Damage zone faults within the syn-rift are commonly very straight, and may be discrete fault planes with no visible fault rock at the scale of observation, or contain relatively thin and simple zones of scaly shale or gouge. The geometric and architectural evolution of the fault array is interpreted to be the result of (i) the evolution from distributed trishear deformation during upward propagation of buried fault tips to surface faulting after faults breach the surface; (ii) differences in deformation response between lithified pre-rift units that display high competence contrasts during deformation, and unlithified syn-rift units that display low competence contrasts during deformation, and; (iii) the history of segmentation, growth and linkage of the faults that make up the fault array. This has important implications for fluid flow in fault zones.     

Analogue and numerical modeling of normal fault patterns produced due to slip along a detachment zone

In this paper, we investigate normal fault patterns produced by the sliding motion along a gently dipping normal fault by using analogue model tests and numerical modeling. The motivation for this study was offered by microseismic test data that indicate the existence of an active low-angle shear zone at a depth of 9–11 km in the extensional region of high seismic activity of the Gulf of Corinth (Greece). Both modeling techniques seem to support the hypothesis that the system of high-angle normal faults that are responsible for the final asymmetrical graben formation initiate at the tip of the active basal detachment nearest to the free surface. The normal faults propagate upwards with progressive sliding of the inclined basal plane, resulting in a first phase of symmetrical graben configuration that is delimited by a main synthetic fault and an antithetic fault forming a Rankine zone. Subsequent sliding on the inclined base induces a family of secondary antithetic normal faults, which are responsible for the asymmetry of the failure pattern and the diffusive character of deformation in that area. Shear deformation is more intense and localized along the synthetic normal fault than along the antithetic faults. Elaboration on the analogue test results has led to the phenomenological relations among four main parameters that describe the geometry of grabens, namely, (i) the width and (ii) the maximum subsidence of the graben, (iii) the dip angles of the conjugate normal faults, and (iv) the amount of sliding along the low-angle normal fault. However, analogue models do not produce the system of synthetic faults that is observed in the Gulf of Corinth. The effects of both friction angle variation along the detachment base and of the constitutive behavior of the model material on the configuration of the final structural pattern were also studied with a series of numerical continuum models. It was found that (a) the fault pattern of the Gulf of Corinth may be reproduced with either a strain-softening material with low elastic modulus or a constant strength material, and (b) two consecutive grabens, such as those of Gulfs of Corinth and Evia, may also be reproduced by an appropriate combination of variation of dip and frictional properties along the hypothesized detachment zone.