Geology Reports⌕ Search

Geology topics

J.D. Byerlee

Publications and source records attributed to J.D. Byerlee.

At least 19 recordsLinked to original sources
USGS2012

Fault stability inferred from granite sliding experiments at hydrothermal conditions

Seismicity on crustal faults is concentrated in the depth interval 1–3 to 12–15km. Tse and Rice (1986) suggested that the lower bound on seismicity is due to a switch with increasing temperature from velocity weakening (destabilizing) to velocity strengthening (stabilizing) friction. They inferred this transition from friction data for dry granite; however, pore fluids exist at elevated temperatures throughout the crust, and may strongly influence strength and sliding behavior. We present new data from sliding experiments on granite at elevated T (23° to 600°) plus elevated PH 2 O (100 MPa),

Geophysical Research Letters↗
USGS2012

Coulomb plasticity within the fault zone

We represent a well‐developed fault by a layer of granular material (fault gouge) confined between two competent fault blocks. Slip on such a fault involves plastic shearing of the fault gouge. That is, the fault gouge behaves as a Coulomb material, and the plastic flow is accomplished by slip on the two sets of Coulomb shears appropriate to the stress state and the frictional properties of the gouge. Neither set of Coulomb shears is coplanar with the fault zone. Generally, the initial plastic flow is self arresting, and increased shear stress is required to drive further shear. But as the shear stress is increased a stress state is reached within the gouge in which the plane of maximum shear stress is parallel to the plane of the fault zone. Then the fault gouge shears at whatever rate is required to keep the shear stress across the fault from increasing further as was demonstrated experimentally by Mandl et al. [1977]. In a mature fault previous slip cycles have presumably left the fault gouge in a stress state such that the plane of maximum shear is parallel to the fault plane. Then, with additional shear loading, unlimited shear deformation is possible at the shear stress corresponding to the onset of Coulomb failure, and the fault appears to slip as if friction on the plane of the fault were less than on planes of other orientations. Slip on the fault does not imply that the fault plane coincides with a Coulomb shear as is sometimes assumed.

Geophysical Research Letters↗
USGS2012

Frictional slip of granite at hydrothermal conditions

Sliding on faults in much of the continental crust likely occurs at hydrothermal conditions, i.e., at elevated temperature and elevated pressure of aqueous pore fluids, yet there have been few relevant laboratory studies. To measure the strength, sliding behavior, and friction constitutive properties of faults at hydrothermal conditions, we slid laboratory granite faults containing a layer of granite powder (simulated gouge). Velocity stepping experiments were performed at temperatures of 23° to 600°C, pore fluid pressures P H 2 O of 0 (“dry”) and 100 MPa (“wet”), effective normal stress of 400 MPa, and sliding velocities V of 0.01 to 1 μm/s (0.32 to 32 m/yr). Conditions were similar to those in earlier tests on dry granite to 845°C by Lockner et al. (1986). The mechanical results define two regimes. The first regime includes dry granite up to at least 845° and wet granite below 250°C. In this regime the coefficient of friction is high (μ = 0.7 to 0.8) and depends only modestly on temperature, slip rate, and P H 2 O . The second regime includes wet granite above ∼350°C. In this regime friction decreases considerably with increasing temperature (temperature weakening) and with decreasing slip rate (velocity strengthening). These regimes correspond well to those identified in sliding tests on ultrafine quartz. We infer that one or more fluid-assisted deformation mechanisms are activated in the second, hydrothermal, regime and operate concurrently with cataclastic flow. Slip in the first (cool and/or dry) regime is characterized by pervasive shearing and particle size reduction. Slip in the second (hot and wet) regime is localized primarily onto narrow shear bands adjacent to the gouge-rock interfaces. Weakness of these boundary shears may result either from an abundance of phyllosilicates preferentially aligned for easy dislocation glide, or from a dependence of strength on gouge particle size. Major features of the granite data set can be fit reasonably well by a rate- and temperature-dependent, three-regime friction constitutive model (Chester, this issue). We extrapolate the experimental data and model fit in order to estimate steady state shear strength versus depth along natural, slipping faults for sliding rates as low as 31 mm/yr. We do this for two end-member cases. In the first case, pore pressure is assumed hydrostatic at all depths. Shallow crustal strength in this case is similar to that calculated in previous work from room temperature friction data, while at depths below about 9–13 km (depending on slip rate), strength becomes less sensitive to depth but sensitive to slip rate. In the second case, pore pressure is assumed to be near-lithostatic at depths below ∼5 km. Strength is low at all depths in this case (<20 MPa, in agreement with observations of “weak” faults such as the San Andreas). The predicted depth of transition from velocity weakening to velocity strengthening lies at about 13 km depth for a slip rate of 31 mm/yr, in rough agreement with the seismic-aseismic transition depth observed on mature continental faults. These results highlight the importance of fluid-assisted deformation processes active in faults at depth and the need for laboratory studies on the roles of additional factors such as fluid chemistry, large displacements, higher concentrations of phyllosilicates, and time-dependent fault healing.

Journal of Geophysical Research↗
USGS2012

Permeability of fault gouge under confining pressure and shear stress

The permeability of both clay-rich and non-clay gouges, as well as several pure clays, was studied as a function of confining pressures from 5 to 200 MPa and shear strain to 10. Permeability ranged over four orders of magnitude, from around 10 −22 to 10 −18 m 2 (1 darcy = 0.987 × 10 −12 m 2 ). The lowest values were characteristic of the montmorillonite-rich and finer grained non-clay gouges. Illite, kaolinite, and chlorite had intermediate permeabilities, while the highest values were typical of the serpentine and coarser grained non-clay gouges. Grain size was an important factor in determining permeability, particularly for the clay-rich samples. The coarse grained gouges were the most permeable and decreased in permeability after shearing. Conversely, the fine grained gouges had characteristically lower permeabilities that did not vary significantly after various amounts of shearing. The permeabilities of the non-clay samples were not significantly different than those of the clays. Therefore, comminuted rock flours can be equally as effective in reducing the flow of water as the characteristically low permeability clay gouges. The strengths of the samples were quite variable. The non-clay gouges were consistently the strongest, with yield points (beginning of nonelastic behavior) around 850 MPa, while montmorillonite had an anomalously low strength in relation to all the other gouges at 250 MPa. Strength of the saturated samples under drained (low pore pressure) conditions did not correlate with high or low permeability. However, the low permeabilities of these gouges could be a factor in the measured low shear stresses along fault regions if excess pore pressures were created as a result of shearing or compaction, and this pressure was unable to dissipate through a thick section of the material.

Journal of Geophysical Research Solid Earth↗
USGS2001

Rapid fluid disruption: A source for self-potential anomalies on volcanoes

Self-potential (SP) anomalies observed above suspected magma reservoirs, dikes, etc., on various volcanoes (Kilauea, Hawaii; Mount Unzen, Japan; Piton de la Fournaise, Reunion Island, Miyake Jima, Japan) result from transient surface electric fields of tens of millivolts per kilometer and generally have a positive polarity. These SP anomalies are usually attributed to electrokinetic effects where properties controlling this process are poorly constrained. We propose an alternate explanation that contributions to electric fields of correct polarity should be expected from charge generation by fluid vaporization/disruption. As liquids are vaporized or removed as droplets by gas transport away from hot dike intrusions, both charge generation and local increase in electrical resistivity by removal of fluids should occur. We report laboratory observations of electric fields in hot rock samples generated by pulses of fluid (water) through the rock at atmospheric pressure. These indicate the relative amplitudes of rapid fluid disruption (RFD) potentials and electrokinetic potentials to be dramatically different and the signals are opposite in sign. Above vaporization temperatures, RFD effects of positive sign in the direction of gas flow dominate, whereas below these temperatures, effects of negative sign dominate. This suggests that the primary contribution to observed self-potential anomalies arises from gas-related charge transport processes at temperatures high enough to produce vigorous boiling and vapor transport. At lower temperatures, the primary contribution is from electrokinetic effects modulated perhaps by changing electrical resistivity and RFD effects from high-pressure but low-temperature CO2 and SO2 gas flow ripping water molecules from saturated crustal rocks. If charge generation is continuous, as could well occur above a newly emplaced dike, positive static potentials will be set up that could be sustained for many years, and the simplest method for identifying these hot, active regions would be to identify the SP anomalies they generate.

Hawai'i↗
USGS2001

How brucite may affect the frictional properties of serpentinite

The frictional strength of brucite gouge has been measured at hydrothermal conditions to 450°C. At room temperature, brucite has a coefficient of friction, μ ≈ 0.30, making it one of the weakest minerals identified to date. With increasing temperature at a constant effective normal stress, the coefficient of friction of brucite decreases to a minimum of μ ≈ 0.20 near 300°C, and μ ≈ 0.22–0.24 in the temperature range 350–450°C. Brucite has a sheeted crystal structure, and its low frictional strength may be attributed to the relatively weak bonds between the layers. In addition, the temperature dependence of μ to ≈300°C can be explained in terms of the anomalously large coefficient of thermal expansion of brucite, which will further weaken the interlayer bonds. Brucite is a common constituent of serpentinite, and at ≈300°C, where brucite is weakest, all the major serpentine minerals have μ ≥ 0.5. The maximum expected brucite content of a serpentinite is close to 20% by weight or volume. That amount of disseminated brucite will lower the coefficient of friction of serpentinite by ≤10–15% in the deeper parts of the seismogenic zone. However, the effect will be much greater if shear can be concentrated along brucite-lined slip surfaces in the serpentinite body.

Open-File Report↗
USGS1998

Quantitative measure of the variation in fault rheology due to fluid-rock interactions

We analyze friction data from two published suites of laboratory tests on granite in order to explore and quantify the effects of temperature ( T ) and pore water pressure ( P p ) on the sliding behavior of faults. Rate-stepping sliding tests were performed on laboratory faults in granite containing “gouge” (granite powder), both dry at 23° to 845°C [ Lockner et al. , 1986], and wet ( P p = 100 MPa) at 23° to 600°C [ Blanpied et al ., 1991, 1995]. Imposed slip velocities ( V ) ranged from 0.01 to 5.5 μm/s, and effective normal stresses were near 400 MPa. For dried granite at all temperatures, and wet granite below ∼300°C, the coefficient of friction (μ) shows low sensitivity to V , T , and P p . For wet granite above ∼350°, μ drops rapidly with increasing T and shows a strong, positive rate dependence and protracted strength transients following steps in V , presumably reflecting the activity of a water-aided deformation process. By inverting strength data from velocity stepping tests we determined values for parameters in three formulations of a rate- and state-dependent constitutive law. One or two state variables were used to represent slip history effects. Each velocity step yielded an independent set of values for the nominal friction level, five constitutive parameters (transient parameters a , b 1 , and b 2 and characteristic displacements D c1 and D c2 ), and the velocity dependence of steady state friction ∂μ ss /∂ ln V = a - b 1 − b 2 . Below 250°, data from dry and most wet tests are adequately modeled by using the “slip law” [ Ruina , 1983] and one state variable ( a = 0.003 to 0.018, b = 0.001 to +0.018, D c ≈ 1 to 20 μm). Dried tests above 250° can also be fitted with one state variable. In contrast, wet tests above 350° require higher direct rate dependence ( a = 0.03 to 0.12), plus a second state variable with large, negative amplitude ( b 2 = −0.03 to −0.14) and large characteristic displacement ( D c2 = 300 to >4000 μm). Thus the parameters a , b 1 , and b 2 for wet granite show a pronounced change in their temperature dependence in the range 270° to 350°C, which may reflect a change in underlying deformation mechanism. We quantify the trends in parameter values from 25° to 600°C by piecewise linear regressions, which provide a straightforward means to incorporate the full constitutive response of granite into numerical models of fault slip. The modeling results suggest that the succeptibility for unstable (stick-slip) sliding is maximized between 90° and 360°C, in agreement with laboratory observations and consistent with the depth range of earthquakes on mature faults in the continental crust.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Strengths of serpentinite gouges at elevated temperatures

Serpentinite has been proposed as a cause of both low strength and aseismic creep of fault zones. To test these hypotheses, we have measured the strength of chrysotile-, lizardite-, and antigorite-rich serpentinite gouges under hydrothermal conditions, with emphasis on chrysotile, which has thus far received little attention. At 25°C, the coefficient of friction, μ, of chrysotile gouge is roughly 0.2, whereas the lizardite- and antigorite-rich gouges are at least twice as strong. The very low room temperature strength of chrysotile is a consequence of its unusually high adsorbed water content. When the adsorbed water is removed, chrysotile is as strong as pure antigorite gouge at room temperature. Heating to ∼200°C causes the frictional strengths of all three gouges to increase. Limited data suggest that different polytypes of a given serpentine mineral have similar strengths; thus deformation-induced changes in polytype should not affect fault strength. At 25°C, the chrysotile gouge has a transition from velocity strengthening at low velocities to velocity weakening at high velocities, consistent with previous studies. At temperatures up to ∼200°C, however, chrysotile strength is essentially independent of velocity at low velocities. Overall, chrysotile has a restricted range of velocity-strengthening behavior that migrates to higher velocities with increasing temperature. Less information on velocity dependence is available for the lizardite and antigorite gouges, but their behavior is consistent with that outlined for chrysotile. The marked changes in velocity dependence and strength of chrysotile with heating underscore the hazards of using room temperature data to predict fault behavior at depth. The velocity behavior at elevated temperatures does not rule out serpentinite as a cause of aseismic slip, but in the presence of a hydrostatic fluid pressure gradient, all varieties of serpentine are too strong to explain the apparent weakness of faults such as the San Andreas.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Failure in laboratory fault models in triaxial tests

A model of a fault in the Earth is a sand-filled saw cut in a granite cylinder subjected to a triaxial test. The saw cut is inclined at an angle α to the cylinder axis, and the sand filling is intended to represent gouge. The triaxial test subjects the granite cylinder to a constant confining pressure and increasing axial stress to maintain a constant rate of shortening of the cylinder. The required axial stress increases at a decreasing rate to a maximum, beyond which a roughly constant axial stress is sufficient to maintain the constant rate of shortening. Such triaxial tests were run for saw cuts inclined at angles α of 20°, 25°, 30°, 35°, 40°, 45°, and 50° to the cylinder axis, and the apparent coefficient of friction μ a (ratio of the shear stress to the normal stress, both stresses resolved onto the saw cut) at failure was determined. Subject to the assumption that the observed failure involves slip on Coulomb shears (orientation unspecified), the orientation of the principal compression axis within the gouge can be calculated as a function of μ a for a given value of the coefficient of internal friction μ i . The rotation of the principal stress axes within the gouge in a triaxial test can then be followed as the shear strain across the gouge layer increases. For μ i ∼ 0.8, an appropriate value for highly sheared sand, the observed values μ a imply that the principal axis of compression within the gouge rotates so as to approach being parallel to the cylinder axis for all saw cut angles (20° < α < 50°). In the limiting state (principal compression axis parallel to cylinder axis) the stress state in the gouge layer would be the same as that in the granite cylinder, and the failure criterion would be independent of the saw cut angle.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Is internal friction friction?

Mogi [1974] proposed a simple model of the incipient rupture surface to explain the Coulomb failure criterion. We show here that this model can plausibly be extended to explain the Mohr failure criterion. In Mogi's model the incipient rupture surface immediately before fracture consists of areas across which material integrity is maintained (intact areas) and areas across which it is not (cracks). The strength of the incipient rupture surface is made up of the inherent strength of the intact areas plus the frictional resistance to sliding offered by the cracked areas. Although the coefficient of internal friction (slope of the strength versus normal stress curve) depends upon both the frictional and inherent strengths, the phenomenon of internal friction can be identified with the frictional part. The curvature of the Mohr failure envelope is interpreted as a consequence of differences in damage (cracking) accumulated in prefailure loading at different confining pressures.

Geophysical Research Letters↗
USGS1996

Strength of chrysotile-serpentinite gouge under hydrothermal conditions: Can it explain a weak San Andreas fault?

Chrysotile-bearing serpentinite is a constituent of the San Andreas fault zone in central and northern California. At room temperature, chrysotile gouge has a very low coefficient of friction (μ ≈ 0.2), raising the possibility that under hydrothermal conditions μ might be reduced sufficiently (to ≤0.1) to explain the apparent weakness of the fault. To test this hypothesis, we measured the frictional strength of a pure chrysotile gouge at temperatures to 290 °C and axial-shortening velocities as low as 0.001 μm/s. As temperature increases to ≈ 100 °C, the strength of the chrysotile gouge decreases slightly at low velocities, but at temperatures ≥200 °C, it is substantially stronger and essentially independent of velocity at the lowest velocities tested. We estimate that pure chrysotile gouge at hydrostatic fluid pressure and appropriate temperatures would have shear strength averaged over a depth of 14 km of 50 MPa. Thus, on the sole basis of its strength, chrysotile cannot be the cause of a weak San Andreas fault. However, chrysotile may also contribute to low fault strength by forming mineral seals that promote the development of high fluid pressures.

Geology↗
USGS1995

An earthquake instability model based on faults containing high fluid-pressure compartments

It has been proposed that large strike-slip faults such as the San Andreas contain water in seal-bounded compartments. Arguments based on heat flow and stress orientation suggest that in most of the compartments, the water pressure is so high that the average shear strength of the fault is less than 20 MPa. We propose a variation of this basic model in which most of the shear stress on the fault is supported by a small number of compartments where the pore pressure is relatively low. As a result, the fault gouge in these compartments is compacted and lithified and has a high undisturbed strength. When one of these locked regions fails, the system made up of the neighboring high and low pressure compartments can become unstable. Material in the high fluid pressure compartments is initially underconsolidated since the low effective confining pressure has retarded compaction. As these compartments are deformed, fluid pressure remains nearly unchanged so that they offer little resistance to shear. The low pore pressure compartments, however, are overconsolidated and dilate as they are sheared. Decompression of the pore fluid in these compartments lowers fluid pressure, increasing effective normal stress and shear strength. While this effect tends to stabilize the fault, it can be shown that this dilatancy hardening can be more than offset by displacement weakening of the fault (i.e., the drop from peak to residual strength). If the surrounding rock mass is sufficiently compliant to produce an instability, slip will propagate along the fault until the shear fracture runs into a low-stress region. Frictional heating and the accompanying increase in fluid pressure that are suggested to occur during shearing of the fault zone will act as additional destabilizers. However, significant heating occurs only after a finite amount of slip and therefore is more likely to contribute to the energetics of rupture propagation than to the initiation of the instability. We present results of a one-dimensional dynamic Burridge-Knopoff-type model to demonstrate various aspects of the fluid-assisted fault instability described above. In the numerical model, the fault is represented by a series of blocks and springs, with fault rheology expressed by static and dynamic friction. In addition, the fault surface of each block has associated with it pore pressure, porosity and permeability. All of these variables are allowed to evolve with time, resulting in a wide range of phenomena related to fluid diffusion, dilatancy, compaction and heating. These phenomena include creep events, diffusion-controlled precursors, triggered earthquakes, foreshocks, aftershocks, and multiple earthquakes. While the simulations have limitations inherent to 1-D fault models, they demonstrate that the fluid compartment model can, in principle, provide the rich assortment of phenomena that have been associated with earthquakes. ?? 1995 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗
USGS1994

Triaxial testing of Lopez Fault gouge at 150 MPa mean effective stress

Triaxial compression experiments were performed on samples of natural granular fault gouge from the Lopez Fault in Southern California. This material consists primarily of quartz and has a self-similar grain size distribution thought to result from natural cataclasis. The experiments were performed at a constant mean effective stress of 150 MPa, to expose the volumetric strains associated with shear failure. The failure strength is parameterized by the coefficient of internal friction ??, based on the Mohr-Coulomb failure criterion. Samples of remoulded Lopez gouge have internal friction ??=0.6??0.02. In experiments where the ends of the sample are constrained to remain axially aligned, suppressing strain localisation, the sample compacts before failure and dilates persistently after failure. In experiments where one end of the sample is free to move laterally, the strain localises to a single oblique fault at around the point of failure; some dilation occurs but does not persist. A comparison of these experiments suggests that dilation is confined to the region of shear localisation in a sample. Overconsolidated samples have slightly larger failure strengths than normally consolidated samples, and smaller axial strains are required to cause failure. A large amount of dilation occurs after failure in heavily overconsolidated samples, suggesting that dilation is occurring throughout the sample. Undisturbed samples of Lopez gouge, cored from the outcrop, have internal friction in the range ??=0.4-0.6; the upper end of this range corresponds to the value established for remoulded Lopez gouge. Some kind of natural heterogeneity within the undisturbed samples is probably responsible for their low, variable strength. In samples of simulated gouge, with a more uniform grain size, active cataclasis during axial loading leads to large amounts of compaction. Larger axial strains are required to cause failure in simulated gouge, but the failure strength is similar to that of natural Lopez gouge. Use of the Mohr-Coulomb failure criterion to interpret the results from this study, and other recent studies on intact rock and granular gouge, leads to values of ?? that depend on the loading configuration and the intact or granular state of the sample. Conceptual models are advanced to account for these descrepancies. The consequences for strain-weakening of natural faults are also discussed. ?? 1994 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗
USGS1994

Reduction of permeability in granite at elevated temperatures

The addition of hydrothermal fluids to heated, intact granite leads to permeability reductions in the temperature range of 300° to 500°C, with the rate of change generally increasing with increasing temperature. The addition of gouge enhances the rate of permeability reduction because of the greater reactivity of the fine material. Flow rate is initially high in a throughgoing fracture but eventually drops to the level of intact granite. These results support the fault-valve model for the development of mesothermal ore deposits, in which seals are formed at the base of the seismogenic zone of high-angle thrust faults. The lower temperature results yield varying estimates of mineral-sealing rates at shallower depths in fault zones, although they generally support the hypothesis that such seals develop in less time than the recurrence interval for moderate to large earthquakes on the San Andreas fault.

California↗
USGS1993

How geometrical constraints contribute to the weakness of mature faults

Increasing evidence that the San Andreas fault has low shear strength 1 has fuelled considerable discussion regarding the role of fluid pressure in controlling fault strength. Byerlee 2,3 and Rice 4 have shown how fluid pressure gradients within a fault zone can produce a fault with low strength while avoiding hydraulic fracture of the surrounding rock due to excessive fluid pressure. It may not be widely realised, however, that the same analysis 2–4 shows that even in the absence of fluids, the presence of a relatively soft 'gouge' layer surrounded by harder country rock can also reduce the effective shear strength of the fault. As shown most recently by Byerlee and Savage 5 , as the shear stress across a fault increases, the stress state within the fault zone evolves to a limiting condition in which the maximum shear stress within the fault zone is parallel to the fault, which then slips with a lower apparent coefficient of friction than the same material unconstrained by the fault. Here we confirm the importance of fault geometry in determining the apparent weakness of fault zones, by showing that the apparent friction on a sawcut granite surface can be predicted from the friction measured in intact rock, given only the geometrical constraints introduced by the fault surfaces. This link between the sliding friction of faults and the internal friction of intact rock suggests a new approach to understanding the microphysical processes that underlie friction in brittle materials.

Nature↗
USGS1992

Observations of quasistatic fault growth from acoustic emissions

The location of acoustic emission (AE) sources during deformation of rock has proven to be a useful non-destructive analytic technique. We present experimental results,based on AE observations, that show the nucleation and growth of macroscopic fault planes in granite and sandstone samples. By controlling axial stress to maintain constant AE rate rather than more conventional loading conditions such as constant strain rate or constant stress, we have been successful in retarding the failure process in triaxial experiments. As a result, the post-failure stress curve has been followed quasistatically, extending to minutes or hours the fault growth process that normally would occur violently in a fraction of a second. In addition to this novel control system, three-dimensional locations of AE events were determined by analyzing the relative arrival times of AE pulses recorded on a network of transducers attached to the sample. In this manner, as many as 40 000 AE events were located in the course of a single experiment. While the details of fault formation varied from experiment to experiment, a number of features were consistently observed. In all three granite experiments, the fault plane nucleated abruptly at a point on the sample surface soon after reaching peak stress. Prior to fault nucleation, microcrack growth was distributed evenly throughout the granite samples. From the nucleation site, the fault plane grew across the sample, accompanied by a gradual drop in axial stress. AE locations showed that the fault propagated as a fracture front (process zone) with in-plane dimension of 10-50 mm. As the fracture front passed, the AE from a given region would drop to a low level. If allowed to progress to completion, stress eventually dropped to the frictional sliding strength. Sandstone samples showed somewhat different response. In these experiments, a diffuse damage zone appeared prior to peak strength and gradually localized into an incipient fault plane. After passing through peak stress, this plane grew, as in the granite samples, to eventually bisect the sample.

International Geophysics↗