Applying Reactivation Tendency Analysis and Mohr-space to Evaluate Shear Strength Decrease and Anisotropy with Pre-existing Weakness(es) under Uniform Stress State
1State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing 102249, China
2Tectonic Laboratory, Department of Earth Sciences, Uppsala University, Uppsala, Sweden
*Corresponding author: Hengmao Tong, State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing 102249, China, Tel: +86 10 89734607, Fax: +86 10 89734158, E-mail: firstname.lastname@example.org ,email@example.com
Received: February 22, 2017 Accepted: April 7, 2017 Published: April 13, 2017
Citation: Tong H, Koyi H, Ramberg H. Applying Reactivation Tendency Analysis and Mohrspace to Evaluate Shear Strength Decrease and Anisotropy with Pre-existing Weakness(es) under Uniform Stress State. Int J Petrochem Res. 2017; 1(1): 31-39. doi: 10.18689/ijpr-1000107
Copyright: © 2017 The Author(s). This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Understanding the mechanical controls on shear strength decrease due to preexisting weakness is a fundamental problem in tectonic studies. In this study, by applying reactivation tendency analysis theory, a theoretical framework and defined Shear-strength Coefficient (fd) are developed for evaluating the shear-strength decrease and anisotropies due to the presence of preexisting weakness(es). The proposed study managed to overcome the restrictions of previous studies assumption that a pre-existing weakness plane contains the intermediate stress (σ2 ) and vertical or horizontal orientations of principal stresses (Andersonian stress state). A new graphical technique (Mohr-space) was utilized to predict the shear-strength decrease and anisotropies caused by preexisting weakness(es). The Mohr-space technique made easier to visualize the state of stress and results of shear strength changes and able to build the quantitative and intuitive relationship between Shear-strength Coefficient (fd) and weakness relativeorientation(θ’,φ›), weakness mechanical properties (Cw and μw) and relative in any uniform tri-axial stress state.. In this study, Shear-strength decrease and anisotropies of a rock sample are evaluated theoretically, and shear strength properties and deformation characteristics of a geological body with multiple pre-existing weaknesses are analyzed and predicted.
Key words: Shear strength; preexisting weakness; reactivation tendency; Mohr-space; sandbox experiment.
It is well documented that pre-existing weaknesses (fracture planes or faults, layering, fabrics etc.) can lead to decrease of shear strength and strength anisotropies [1-3], and the potential strength anisotropy created by pre-existing weakness is considerable . Ranalli (1990)  proposed a unified quantitative model to evaluate strength anisotropies by the pre-existing weakness in terms of the three tectonic faulting regimes when the weakness plane contains the intermediate stress (σ2 ). However, to our knowledge strength anisotropies evaluation with pre-existing weakness is limited to two-dimensional cases (weakness plane containing σ2 ) before slip-tendency theory was proposed . However, strength anisotropies evaluation with slip-tendency  is confined to the qualitative analysis only. In addition, there is also some limitation in slip-tendency theory, i.e. the principal stresses are oriented either vertically or horizontally, and the cohesive strength of all preexisting weakness is neglected . However, principal stress direction may depart significantly from vertical and horizontal with depth in the upper crust [8-9]. Furthermore, some weakness zones may possess cohesive strength, particularly in cemented faults, and properties of weaknesses may vary . For a example, Sibson (1974)  showed that cementation of a fault zone can create 1.0-MPa cohesive strength or more. In this study, Shear-strength Coefficient (fd) is defined, and Reactivation Tendency analysis theory , and Mohr-space  are applied to evaluate strength decrease and anisotropies due to the presence of pre-existing weakness(es) with arbitrary azimuth in any uniform tri-axial stress state. Strength decrease and anisotropies caused by weakness may be intuitively simplified and quantitatively analyzed and evaluated with Mohr space.
Mohr-space, pole (σn, τn) of any oriented plane in tri-axial stress state
Mohr diagrams, which was introduced by Otto Mohr (1882), is one of the most used and useful tools in structural geology , and has been used extensively in mechanical problems, such as stress analysis, failure envelopes, fractures opening and reactivation    . Although real three-dimensional Mohr diagrams do exist for any tri-axial stress state , Mohr-diagram is usually considered to be two-dimensional , which is well known as Mohr-circles. Mohr-cyclides, which can be used to represent any second rank tensor (including stress tensor) was introduced by Coelho and Passchier (2008) . However, the stress components (σn , τn ) of any given plane are the most important, and the general diagrams of Mohr-cyclides are not so convenient to be used and prepared. In contrast, the Mohr space, which was proposed by Tong and Yin (2011)  can be used to express the normal stress and shear stress of a plane with an arbitrary azimuth in an arbitrary three-dimensional stress state.
In any given stress state, the pole (σn , τn ) of any plane (i.e. pre-existing weakness plane, defined by dip direction θ and dip angle φ) is either located on the three Mohr-circles (i.e. P1 on large Mohr circle σ1 -σ3 , P2 and P3 on small Mohr circles σ1 - σ2 and σ2 -σ3 respectively, Figure 1) or in the area (grey area in Figure 1) between large Mohr-circle and two small Mohrcircles(i.e. P4 in Figure 2) in σn -τn coordinate system [2,18]. There is one to one correspondence relationship between any plane and its pole (σn , τn ) in σn -τn coordinate system. This space (i.e. the three Mohr-circles and the area between them) is called Mohr-space . With the contour lines (pink and green dotted lines in Figure 2) of plane angles (pseudo-dip direction θ› and pseudo-dip angle φ›, θ>= α> +90 °, where α> is the angle between the plane and the intersection line of the plane and σ2 -σ3 plane; φ> the angle between the plane and σ2 -σ3 plane), the pole of the plane can be found in the Mohrspace. The relationship between θ>, φ>, θ, and φ can be determined with transformation of coordinates , and the contour lines of angles can be compiled with equation (1) .
Coordination definition and transformation can be seen in Tong and Yin (2011) .
Thus, by applying Mohr-space, normal stress (σn ) and shear stress (τn ) of any plane can be conveniently and intuitively determined, and the changes of normal and shear stresses with plane relative-orientation (θ´,φ>) can be easily analyzed.
Reactivation Tendency Factor and its expressions in Mohr-space
(where τn , σn are shear and normal stresses on the weakness plane, respectively; τnw is the corresponding rupture shear strength; Cw, and μw are cohesiveness and frictional coefficient of weakness plane, respectively), which is extended from Sliptendency  and is used to evaluate the reactivation likelihood of pre-existing weakness, is proposed by Tong and Yin (2011). Reactivation Tendency Factor (fa ) is determined by its relative-orientation (θ´,φ´), mechanical properties (Cw, μw ) of the weak plane, and the stress tensor. fa =1.0 shows that the pre-existing weakness is in critical state of reactivation; when fa <1, it has reactivated and when fa >1, it is in a stable state. The weaknesses, which fa ≥1.0 in the critical stress state of Coulomb rupture, will reactivate one by one in progressive deformation according to their fa value order.
Applying Mohr-space, fa can be intuitively expressed accordingto the following steps. (1) The lines τnw=± (Cw + μwσn ) are two symmetrical straight lines (when μw is a constant)or curved lines(b and b´ in Figure 2, when μw is variable) in σn -τn coordinate system, which were called “weakness reactivation lines”  or “shear strength line of weakness ” , where (0, Cw ) is the starting point, μw is the slope of the line. If Cw = 0 (weakness without cohesion), the line cross the point of origin. (2) The pole (σn , τn ) of a weakness plane (P in Figure 2) can be easily plotted in Mohr-space according to its relative-orientation (θ´and φ´) and its σn , τn w can be intuitively determined (Figure 2). As such, fa (τn /τn w) of the weakness plane can be easily and intuitively determined (PP0 /PPR in Figure 2). (3) When the pole of the weakness plane is on its “weakness reactivation line” (P1 in Figure 2), fa =1.0 and the pre-existing weakness is in the critical state of reactivation. However, when the pole is located outside the two “weakness reactivation lines” (striped yellow area in Figure 2, i.e. P2 ), fa >1.0 and the weakness plane has reactivated. When the pole is located inside the two weakness reactivation lines (yellow area without stripes in Figure 2, i.e. P3 ), fa<1 indicating that the plane is in a stable state.
The definition of Shear-strength coefficient and its relationship to Reactivation Tendency Factor of pre-existing weakness
As a weakness plane may be reactivated when stress is below Coulomb rupturing stress state, a pre-existing weakness will lead to decrease in shear strength   . However, the weakness with fa <1.0 will not reactivate at critical stress of Coulomb rupture  according to Reactivation Tendency analysis theory. This means that not all weakness planes will lead to decrease in shear strength. －Therefore, in the following shear strength analysis, we only consider weaknesses with fa≥1.0.
Consider a weakness plane (Pw) with a normal n and its reactivation tendency factor fa ≥1.0 at critical uniform stress state (σ1 , σ2 and σ3 ) of Coulomb rupture (Figure 3). Pw has been reactivated at critical stress state of Coulomb rupture , and its pole in Mohr-space is (σn , τn ) (Point A in Figure 3). With the same normal stress σn , Pw will be at critical state of reactivation when its shear stress is τnw (Point B in Figure 3, σ´1 , σ´2 and σ´3 are the three principal stresses) and τnw<τn . The large Mohr circle (σ´1 -σ´3 circle) of critical stress state of weakness reactivation is smaller than the σ1 -σ3 Mohr circle of critical stress state of Coulomb rupture (Figure 3), that means σ´1 - σ´3 < σ1 - σ3 . Thus, in Mohr-space (Figure 3), it is very easy to understand that a weaknesswith fa >1.0will lead to decrease in shear strength.
Because shear strength is related to normal stress and increases with increasing of normal stress, relative instead of absolute shear strength is more useful to consider.In order to quantitatively evaluate the shear strength decrease due to the existence of a weakness, we define the parameter fd called Shear-strength coefficient (fd = τnw /τn , where τnw is the critical shear stress of weakness reactivation, τn is the shear stress on the same plane with the same normal stress σn at critical stress state of Coulomb rupture (σ1 , σ2 and σ3 stress state in Figure 3)). It is easy to understand that the smaller Shear-strength coefficient (fd ) is the larger decrease will be in shear strength.
When fa ≥1.0, τnw can never be greater than τn (as stress will drop down when weakness reactivates according to Reactivation Tendency analysis theory  . Furthermore, when fa ≤ 1.0, which means that weakness will not be reactivated at critical stress state of Coulomb rupture and will not lead to shear strength decrease(i.e., fd=1.0). Therefore, shear-strength coefficient can never be greater than 1.0 (fd ≤ 1.0).
It is easy to find the following relationship between fd and Reactivation Tendency Factor (fa ) at critical stress state of Coulomb rupturing according to the definition of shearstrength coefficient and Reactivation Tendency analysis theory 
The relationship between Shear-strength Coefficient (fd ) and the ratio of Differential stress can be derived at Coulomb rupture ((σ1 - σ3 ) C =σ1 - σ3 in Figure 3) and at critical state of weakness reactivation ((σ1 - σ3 ) W =σ1 ´- σ3 ´ in Figure 3) according to the definition of Shear-strength Coefficient and using Mohr-space, and it means that fd = (σ1 - σ3 )W / (σ1 - σ3 )C. As the relative position of the weakness pole in (σ1 , σ2 , σ3 )-Mohr-space (point A) and in (σ1 ´, σ2 ´, σ3 ´)-Mohr-space (point B) is the same, OA/OB is equal to the ratio of radius of the two big Mohr circle ( 2 σ1 −σ 3 and 2 σ1 ´−σ 3´), which means that OA/OB = (σ1 - σ3 )/(σ1 ´- σ3 ´) = (σ1 - σ3 ) C/ (σ1 - σ3 ) W. While two triangles AOC and BO´C are similar, as a result, BC /AC =OB /OA and fd =τn w /τn = BC /AC = OB / OA = (σ1 - σ3 ) C/ (σ1 - σ3 ) W, which means that fd = (σ1 - σ3 ) C/ (σ1 - σ3 ) W when fa ≥1.0. So, fd (1/ f a ) is related to weakness orientation (θ´and φ´), mechanical properties (Cw, μw ) and the stress tensor according to Reactivation Tendency analysis theory , and it can be quantitatively calculated using equations 1 and 2. Since fd is related to weakness relative-orientation, weakness can lead to shear-strength anisotropies.
Analysis of shear strength affection factors due to pre-existing weakness(es)
Shear-strength coefficient (fd) is used to discuss how the related factors (mechanical properties and orientation of weakness, and ) affect shear strength due to preexisting weakness in triaxial stress state. As fd = 1/fa , and fa can be intuitively and conveniently expressed in Mohr-space, Mohrspace is applied to do these analysis.
The relationship between fd and weakness relativeorientation (θ´,φ´)
It is easy to see that when θ´=90°, and φ´ =±φ´0 = ± (45°+0.5arctg (μw) (φ´0 is usually about 60°) (points A and A´ in Figure 4), fa reaches the highest value, whereas fd reaches the smallest value (Figure 4). Therefore, (90°, ± φ´0 ) (two points) are the two weakest relative-orientation of the weakness plane(s). As θ´ decreases from 90°, and φ´ deviates from φ´0 , fd will increase. In Mohr-space (Figure 4), the points of intersection between “weakness reactivation lines” and Mohr circles are demarcation points of weakness reactivation (B1 , B2 and B´1 , B´2 on σ1 -σ3 large Mohr circle, C1 , C2 and C´1 , C´2 on σ2 -σ3 small Mohr circle, Figure 4). The corresponding φ´ of B1 , B2 , B´1 and B´2 is φ´1 , φ´2 , φ´1 , φ´2 , and the corresponding θ´ of C1 , C2 and C´1 , C´2 is θ´1 , θ´2 , (180° - θ´1 ) and (180° -θ´2 ) respectively. It is easy to find that when φ´<φ´1 or φ´> φ´2 , or θ´<θ´1 or θ´2 <θ´<180°-θ´2 or θ´>180°-θ´1 , the weakness will all locate inside the two weakness reactivation lines (grey area in Figure 4) and fa will always <1.0, which means there exist critical angle φ´1 , φ´2 for φ´ and θ´1 ,θ´2 and (180° - θ´1 ),(180° -θ´2 ) for θ´, when φ´<φ´1 or φ´> φ´2 , or θ´<θ´1 or θ´2 <θ´<180°-θ´2 or θ´>180°-θ´1 , fa will always <1.0 and the weakness cannot reactivate and will not lead shear strength decrease (fd = 1.0).
Point A and A´ (Figure 4) are two Coulomb rupture points, and fd of the same oriented weakness plane is the lowest. Weakness plane with low fd value concentrates around point A (or A´) in the yellow area of Figure 4(θ´=90°~75° and φ´=φ´0 ± 15°). Change of fd shows that weakness relative-orientation is the predominant factor which leads to shear strength anisotropies. It is noteworthy tounderline that (θ´, φ´) of weakness is determined jointly by orientations of weakness and three axes of principal stresses.
The relationship between fd and weakness mechanical properties (Cw and μw)
In Mohr-space (Figure 5), it is easy to find that the change of mechanical property of a weakness plane (Cw and μw) will cause the change of the position (determined by Cw) and the slope (determined by μw) of weakness reactivation lines. It is easy to understand that fd will decrease with decreasing of Cw and μw.
Weakness mechanical properties are the predominant factors which lead to shear strength decrease. For example, low friction coefficient (μw< 0.2 ) and no cohesion along the weakness may lead to more than 90% shear strength decrease (fd< 0.1).
The relationship between fd and relative value of σ2
Shear strength is also related to relative value of . However, our analysis in Mohr-space provided intuitive results which are easy to follow. It is easy to see that relative determines the shape of Mohr-space (Fig. 6). The relationship between and fd is complicated and depends on θ´ and φ´, particularly φ´(Fig. 6). In general, φ´ can be divided into 3 intervals: φ´≥70°, 40°<φ´<70°, and φ´≤40°. When φ´≥70°fd will increase (Pw1 in Figure 6). When 40°<φ´<70°, f d decreases a little (Pw3, Pw4 and Pw5 in Figure 6), and when φ´≤40°,fd changes little (Pw6 and Pw7 in Figure 6). These three cases are valid whenσ2 decrease(or increases, Fig. 6). As such, the change of σ2 should not be ignored in shear strength analysis in the presence of pre-existing weakness. However, when φ´≥70° or φ´≤40° (particularly when φ´≤40°), most of the weaknesses usually cannot be reactivated at critical stress state of Coulomb rupturing and thus doesnot lead shear-stress decrease (fd = 1.0). As a result, the effect of relative σ2 is most prominant only in the interval 40°<φ´<70°; under normal circumstances, fd will decrease a little with decreasing σ2 (Fig. 6).
Mechanical properties and relative-orientation, which are the governing factors of shear strength decrease and strength anisotropies, respectively, are internal factors, while is an external factor and not so important according to the above analysis,. When the relative-orientation (θ´, φ´) and mechanical properties (Cw, μw ) of pre-existing weakness, and are given, the Shear-strength Coefficient can be quantitatively evaluated in Mohr-space.
Shear strength evaluation for rock samples and geological bodies: theoretical analysis and verification
There are many kinds of pre-existing weakness, such as faults, geologic contacts, bedding and foliation that affect shear strength   . Based on their structure, Morley (2002)  divided weaknesses into two types: “discrete” and “pervasive”. Based on the value of cohesive strength, Tong and Yin (2011)  divide weaknesses into “strong weakness” (with relatively large cohesive strength, such as bedding, foliation etc.) and “weak weakness” (with relatively small or zero cohesive strength (i.e. ignorable), such as faults, fracture planes etc.). “Discrete” weakness is usually a “weak weakness”, and “pervasive” weakness probably is a “strong weakness”. We will evaluate the effect of “strong weakness” and “weak weakness” on shear strength in rock samples, and analyze deformation sequences with multiple pre-existing weaknesses in geological bodies.
“Pervasive” weaknesses (or “strong weakness”, bedding in sedimentary rocks, foliation in metamorphic rocks) do exist in rocks. There may also exist“pervasive” weaknesses in magmatic rocks with ductile deformation or flow foliation. In homogeneous-looking rock samples there may exist “pervasive” (particularly in sedimentary or metamorphic rocks) and /or “discrete” weaknesses (“weak weakness”, internal small or micro-fracture) leading to shear strength anisotropies. In order to quantitatively evaluate shear strength decrease and anisotropies of rock samples with “pervasive” or “discrete” weakness, the situations of ①Cw =0.5C, 0.33C and 0.2C for “pervasive”weakness, and Cw =0 for “discrete” weakness, and ② = 0.5, 1.0 and 2.0 are considered (let μw = µ = 0.6, C = 20MPa). Applying Mohr-space, the calculated values of Shear-strength Coefficient can be seen in Table 1. For other C values, the result is completely the same, while the result will change a little bit as µ ( μw ) changes (if μw = µ ). However, if μw changes while µ isconstant, fd will change proportionally with μw .The results show thata “discrete” weakness may lead to more than 80% maximum drop of shear strength, while “pervasive” weakness usually leads to 20-60% maximum shear strength decrease for rock samples.
Table 1. Shear-strength Coefficient (fd ) value in different value of CW,relative σ2 , θ´ and φ´
Haimson and Rudnicki (2010)  conducted true tri-axial compression tests on siltstone samples, and showed that shear strength is related to σ2 . Although the orientation of siltstone bedding is not mentioned in the paper, we speculate that the shear strength change is caused by “pervasive” weakness (siltstone bedding). If bedding does not lead to mechanical anisotropy (i.e. when it is mechanically homogeneous), shear strength will not change with σ2 . The tri-axial compression tests of Haimson and Rudnicki (2010)  may be the verification for the above theoretical analysis of dependece of shearstrength on σ2 . Naturally, availability of more experiment data is needed.
In the presence of multiple pre-existing weaknesses in a rock sample, each weakness has its own Shear-strength Coefficient. The overall Shear-strength Coefficient is determined by and is equal to that of the weakest weakness, where the Shear-strength Coefficient is the smallest, if interactions of weaknesses are ignored.
Geological bodies and deformation sequence
There are the two main differences between geological bodies and rock samples when the size and locationof weakness are concerned. It is possible to prepare relatively homogeneous small rock samples (it is necessary for regular rock mechanics test). However, it is inevitable that there are pre-existing weaknesses more or less in the geological bodies, which will lead to shear strength decrease (shear strength is much smaller than that of rock samples   .
On a larger scale, there may or probably exist multiple weaknesses  in geological bodies (i.e. rift basins, orogenic belts, or a part of them). Different relative-orientation of the weaknesses and/or their mechanical properties will lead to different fd of weakness. The relative-orientation and/or mechanical properties may also vary greatly along large scale pre-existing weakness (i.e. big pre-existing faults) and will lead to different fd in different segments along the same weakness.
Unlike rock samples, in the presence of multiple pre-existing weaknesses in geological bodies, one of the weaknesses reactivation will lead to stress drop and form local stress field only in the area along and near the weakness   instead of the whole geological body. Differential stress within other areas of the body away from the weakness can increase with progressive deformation until another weakness reactivates or Coulomb rupturing occurs. In otherwords, a weakness can lead to shear strength decrease only within a part of geological bodies (along and near a “prefered” weakness). As a result, in the condition of homogeneous regional stress field, it is easy to understand that the weaknesses, which are away from each other and fd < 1.0, can reactivate according to their fd values (from small to large) in the progressive deformation: the weakest weakness with the smallest fd value reactivates first, then the second weakest, and so on. Coulomb rupturing will occur at last in the area away from the weaknesses.
In order to verify the above statements, a sandbox model is run, where multiple pre-existing weaknesses are built away from each other. The fd values of these weaknesses can be quantitatively determined, and the regional stress field can be regarded as homogeneous.
In the models, three pre-existing weakness planes (Pw1, Pw2 and Pw3), which are oblique to the extension direction (σ3 ), with (φ´, θ´) = (80°, 55°), (60°, 70°) and (30°, 49°) respectively, are set. The weakness planes are represented by a 80mm wide slice of thin sheets of paper inserted into the base of a homogeneous layer (8.0 cm thick) of dry quartz sand (Fig. 7). Using these sheets of paper allowed localizing fault initiation, and it was easy to determine friction angle between the paper and loose sand. The friction angle of dry quartz sand is 31°, while the friction angle of sand with paper is 20° . So the μ and μw is 0.60 (= tg31°) and 0.36 (= tg20°). The cohesion of dry sand is very small, so C and Cw are both assumed to be equal to 0. As a result, the relative Reactivation Tendency Factor of the three weakness planes is fa1 =1.59>fa2 =1.42>1.0 >fa3 =0.85 in the critical stress state of Coulomb rupture (Fig. 8), and Shear-strength Coefficient fd1 =0.63<fd2 = 0.70<1.0= fa3.
With the above model (Figure 7), the weakness, with Shear-strength Coefficient fd1<1.0, reactivated to form several faults, one after another, in the progressive extension as predicted, and Coulomb ruptures faults form at last (Figure 9). The deformation in the progressive extension is summarized in the following:
Figure 7. Experiment design of sandbox modeling. (a) Skeleton figure of experimental model. It is 40cm wide, 50cm long and 8cm thick. Pw1, Pw2, and Pw3 are weakness planes which is oriented oblique to the extension direction, with (θ´, φ´) are (80°, 55°), (60°, 70°) and (30°, 49°) respectively, and represented by paper being inserted into the homogeneous dry sand. Elastic rubber is under the sand layer and connected to driving end. (b) Section of AA´ in figure (a).
The weakness plane with the smallest Shear-strength Coefficient (Pw1) reactivated to form a fault at first (fault number 1, Fig. 9b-1). Then, the second weakness plane (Pw2) reactivated to form a fault (fault number 2, Fig. 9b-1). In the middle stage (after the weakness reactivated and before the fault near perpendicular to extension direction formed), weakness related faults which are parallel (or sub-parallel) to and controlled by Pw1 or Pw2 (antithetic and synthetic faults, fault numbers 3-8in Figures 22b-2 and b-3) began to develop near the weakness plane area. In the final stage (Fig. 9b-3), faults near perpendicular to the extension direction (small faults in Figure 9b-3) began to develop. While the third weakness plane, which Shear-strength Coefficientfd3=1.0 (fa3 =0.85<1.0, Fig. 8), does not reactivate as predicted. Similar experiments carried out by Tong and Yin (2011)  and Tong et al. (2014)  showed the same results.
In general, a phase of deformation indicated by the initiation of a new fault trend is attributed to a specific stress regime. However, our analysis suggests that in a region where the magnitude of differential stresses progressively increase while their directions are kept constant, multiple phases of fault initiation with different trends can occur. This could be due to reactivation of pre-existing weaknesses. The sequence of reactivation of any preexisting weakness is predictable according to the relative-orientation and mechanical properties of the weak planes.
It is not easy to estimate the shear strength of pre-existing weaknesses in the field. Even though, the relative-orientation of a weakness may be determined accurately in some circumstances (when high quality 3D seismic data is available, and fault systems in shallow level can be mapped accurately),mechanical properties of the weakness (Cw, μw ) are not easily determined.However, if the sequence of faulting in the geologic record is known , this information may be reversely used for assessing the mechanical properties of the preexisting weaknesses at the time of their reactivation.
Although this study is expanding the work of Ranalli (1990)  and Morris et al , our newly defined Shearstrength Coefficient, newly developed graphical techniqueMohr-space and the underlying theoretically framework, which is based on Reactivation Tendency analysis theory (Tong and Yin, 2011), provide a much general and intuitive treatment of the shear strength decrease and anisotropies caused by pre-existing weakness(es). Specifically, the assumption that the weakness plane containing the intermediate stress (σ2 ) in Ranalli´s (1990)  analysis can now be neglected. Finally, we predicted that weaknesses will reactivate sequentially according to the Shear- strength coefficient values order (from small to large) and new fractures (Coulomb rupture) form at last in the progressive deformation. It was verified by a simple sandbox experiment.
Both in this study and that of Ranalli (1990) , it is assumed that the preexisting weakness must be planar. This assumption, however, does not prevent the application of the analysis developed in this study to non-planar weaknesses, as curved surfaces can always be divided into approximately small planar segments. Another important assumption in our study is that fault formation, propagation and activity do not affect the regional stress distribution. This assumption is clearly an oversimplification, as both detailed analysis of faultzone evolution and regional modeling show that frictional sliding on faults is capable of creating local stress fields near the faults that are different from the regional stress field  . However, the results of our sandbox experiment imply that the effect of activating preexisting weakness does affect local stress, but is minimal in changing regional stress. More research is clearly needed to address this problem.
Possible complex interactions between faults that lie near one another were not considered in our model. As stress concentrates at crack tips, it is expected that complex stress fields can be induced by the presence of cracks or weaknesses under uniform regional stress  (e.g., Gudmundsson et al., 2009). Thus, our model should be considered as an idealized conceptual guide, which can only be used in realistic situations when the above factors are considered.
Despite the complexities discussed above, our analysis does provide a new insight into the temporal evolution of multiple pre-existing weaknesses with different Shear– strength coefficient. Our analysis suggests that in a region with a progressive increase in the magnitude of differential stresses while the directions of the principal stresses are not changed, multiple phases of fault initiation with different trends can be generated. On the other hand, the results of theoretical analysis will provide some information and clues to understand actual shear strength decrease and anisotropies due to the pre-existing weaknesses.
Shear-strength coefficient (fd), which is defined to evaluate the shear strength decrease due to the presence of preexistingweaknesses, is determined by the orientation and mechanical properties of weakness(the intrinsic factors) and stress tensor (the external factors), and can be calculated quantitatively. It also can be expressed intuitively/graphically in Mohr-space. The results of theoretical analysis in Mohrspace show that:
(1)Weakness relative-orientation (θ´, φ>), which is determined jointly by orientation of the weakness and three principal stress axis, is the predominant factor leading to shear strength anisotropies. (0°, φ´0 ) and (180°, φ´0 ) are two relativeorientations with the lowest fd. As θ´ increases from 0° or decreases from 180°, and φ´ deviates from φ´0 , fd will increase. There are critical angles φ´1 , φ´2 and θ´1 , 180°-θ´1 , when φ´1 <φ´< φ´2 or θ´1 <θ´<180°-θ´1 , fd = 1.0. Low fd value of weakness plane concentrates in the area around θ´=0°~15° or 165°~180° and φ´=φ´0 ± 15°.
(2)Weakness mechanical properties (Cw and μw) are the predominant factorsthat lead to shear strength decrease. “Discrete” weaknesses may suffer more than 80% decrease in maximum shear strength, while “pervasive” weaknesses usually suffer 20-60% decrease in maximum shear strength.
(3)The effect of relative to shear strength is a little complicated, and is related to θ´ and φ´, particularly φ´. The effect of relative σ2 is most prominant only in the interval 40°<φ´<70°. Under normal circumstances, fd will decrease a little with increasing of .
In a region with a progressive increase in the magnitude of differential stresses acting on multiple pre-existing weaknesses, while the directions of the principal stresses maintain the same, multiple phases of fault initiation with different trends can be generated with predicted sequence according to their shear–strength coefficient (from small to large). This theoretical prediction wasverified by the results of a sandbox model.
We would like to thank graduate students Mingyang WANG and Hua wu HAO for drawing the figures of the Mohr space and their help to complete the sandbox model.This study was supported by China National Major Project of Oil and Gas (2016ZX05024-005-004 & 2011ZX05023-004- 0122011zx05006-006-02-01) and China Natural Science Foundation (Grant No. 41272160&40772086). HAK is supported by the Swedish Research Council.
- Jaeger JC. Shear fracture of anisotropic rocks. Geol. Mag. 1960; 97: 65-72.
- Jaeger J C, Cook NGW. Fundamentals of rock mechanics (third edition): London, United Kingdom, Chapman and Hall, 1979; 593.
- Ranalli G. Rheology of the Earth: Deformation and Flow Processes in Geophysics and Geodynamics. Allen & Unwin, London. 1987.
- Morley C K, Haranya C, Phoosongsee WS. J Structural Geology. 2004; 26(10): 1803-29. doi: 10.1016/j.jsg.2004.02.014
- Ranalli G, Yin ZM. Critical stress difference and orientation of faults in rocks with strength anisotropies: the two-dimensional case. J Structural Geology. 1990; 12(8): 1067-71. doi: 10.1016/0191-8141(90)90102-5
- Morris AP, David AF, Brent H. Slip-tendency analysis and fault reactivation. Geology. 1996; 24(3): 275-78. doi: 10.1130/0091-7613
- Morris AP, Ferrill DA. The importance of the effective intermediate principal stress (σ2 ) to fault slip patterns. J Structural Geology. 2009; 31(9): 950-59. doi: 10.1016/j.jsg.2008.03.013
- Spencer JE, Chase CG. Role of crustal flexure in initiation of low-angle normal faults and implications for structural evolution of the Basin and Range province. J. Geophys. Res. 1989; 94: 1765-75. doi: 10.1029/ JB094iB02p01765
- Yin A. Origin of regional, rooted low-angle normal faults: a mechanical model and its tectonic implications. Tectonics.1989; 8(3): 469-82. doi: 10.1029/TC008i003p00469
- Sibson RH. Frictional constraints on thrust, wrench, and normal faults. Nature. 1974; 249: 542-44. doi: 10.1038/249542a0
- Tong H, Yin A. Reactivation tendency analysis: A theory for predicting the temporal evolution of preexisting weakness under uniform stress state. Tectonophysics. 2011; 503(3-4): 195-200. doi: 10.1016/j. tecto.2011.02.012
- Tong HM, Wang JJ, Zhao HT, et al. Mohr space and its application to the activation prediction of pre-existing weakness. Science China: Earth Sciences. 2014; 57(7): 1-10. doi: 10.1007/s11430-014-4860-1
- Coelho S, Passchier C. Mohr-cyclides, a 3D representation of geological tensors: The examples of stress and flow. J Structural Geology. 2008; 30(5): 580-601. doi: 10.1016/j.jsg.2008.01.003
- Brace WF. Mohr construction in the analysis of large geological strain. Geolological Society of America Bulletin. 1961; 72(7): 1059-79. doi: 10.1130/0016-7606
- Ramsay JG. Folding and Fracturing of Rocks. McGraw-Hill, New York 1967.
- Delaney PT, Pollard DD, Zioney JI, MacKee EH. Field relations between dikes and joints: Emplacement processes and paleostress analysis. Journal of Geophysical Research. 1986; 91: 4920-38. doi: 10.1029/ JB091iB05p04920
- Jolly RJH, Sanderson DJ. A Mohr circle construction for the opening of a pre-existing fracture., J Structural Geology. 1997; 19(6): 887-892. doi: 10.1016/S0191-8141(97)00014-X
- Twiss RJ, Moores. Structural Geology. New York, W. H. Freeman and Company 1992.
- Tong H, Cai D, Wu Y, Li X, Li X, Meng L. Activity criterion of pre-existing fabrics in non-homogeneous deformation domain. Science China (Earth Science). 2010; 53(8): 1115-25. doi: 10.1007/s11430-010-3080-6
- Collettini C, Trippetta F. A slip tendency analysis to test mechanical and structural control on aftershock rupture planes. Earth and Planetary Science Letters. 2017; 255(3-4): 402-413. doi: 10.1016/j.epsl.2007.01.001
- Bailey W, Ben-Zion Y. Statistics of Earthquake Stress Drops on a Heterogeneous Fault in an Elastic Half-Space. Bulletin of the Seismological Society of America. 2009; 99(3): 1786–1800.
- Yin ZM, Ranalli G. Estimation of the frictional strength of faults from inversion of fault-slip data: a new method. J Structural Geology. 1995; 17(9): 1327-35. doi: 10.1016/0191-8141(95)00028-C
- Wallace RE. Geometry of shearing stress and relation to faulting. J. Geol. 1951; 59: 118-130.
- Bott MHP. The mechanics of oblique slip faulting. Geol.Mag. 1959; 96(2): 109-117.
- McKenzie DP. The relation between fault plane solutions for earthquakes and the directions of the principal stresses. Bull.Seismol. Soc. Am. 1969; 59: 591-601.
- Morley CK. A tectonic model for the Tertiary evolution of strike–slip faults and rift basins in SE Asia. Tectonophysics. 2002; 347(4): 189-215. doi: 10.1016/S0040-1951(02)00061-6
- Haimson B, Rudnicki JW. The effect of the intermediate principal stress on fault formation and fault angle in siltstone. J Structural Geology. 2010; 32(11): 1701-11. doi: 10.1016/j.jsg.2009.08.017
- Tong H, Koyi H, Huang S, Zhao H. The effect of multiple pre-existing weaknesses on formation and evolution of faults in extended sandbox models. Tectonophysics. 2014; 626: 197-212. doi: 10.1016/j. tecto.2014.04.046
- Gudmundsson A, Simmenes TH, Larsen B, Philipp SL. Effects of internal structure and local stresses on fracture propagation, deflection, and arrest in fault zones. J. Struct. Geol. 2010; 32(11): 1643-55. doi: 10.1016/j. jsg.2009.08.013
- Yin A. Mechanics of monoclonal systems in the Colorado plateau during the Laramide orogen. J. Geophys. Res. 1994; 99: 22043-58.