The so-called “brittle-ductile transition” is thought to be the strongest part of the lithosphere, and defines the lower limit of the seismogenic zone. It is characterized not only by a transition from localized to distributed (ductile) deformation, but also by a gradual change in microscale deformation mechanism, from microcracking to crystal plasticity. These two transitions can occur separately under different conditions. The threshold conditions bounding the transitions are expected to control how deformation is partitioned between localized fault slip and bulk ductile deformation. Here, we report results from triaxial deformation experiments on pre-faulted cores of Carrara marble over a range of confining pressures, and determine the relative partitioning of the total deformation between bulk strain and on-fault slip. We find that the transition initiates when fault strength (σf) exceeds the yield stress (σy) of the bulk rock, and terminates when it exceeds its ductile flow stress (σflow). In this domain, yield in the bulk rock occurs first, and fault slip is reactivated as a result of bulk strain hardening. The contribution of fault slip to the total deformation is proportional to the ratio (σf − σy)/(σflow − σy). We propose an updated crustal strength profile extending the localized-ductile transition toward shallower regions where the strength of the crust would be limited by fault friction, but significant proportions of tectonic deformation could be accommodated simultaneously by distributed ductile flow.

Under the low pressure and temperature conditions of the upper crust, rocks generally deform by grain-scale microcracking, and crustal-scale deformation is accommodated by slip on discrete fault planes. In this regime, the overall strength of the crust is limited by fault friction (Scholz, 2002; Paterson and Wong, 2005). Deeper in the crust, at higher pressure and temperature, rock deformation becomes more diffuse, and may be driven by crystal plastic phenomena such as dislocation creep. Here, the overall strength of rocks can generally be described by a steady-state flow law sensitive primarily to temperature and strain rate (e.g., Goetze and Brace, 1972; Evans and Kohlstedt, 1995). The transition between these two rheological domains, the so-called “brittle-ductile transition”, occurs over a pressure and temperature range where rocks deform by an interplay of cracking and crystal plasticity. The brittle-ductile transition commonly loosely refers to the progressive change in crustal rheology with increasing depth; here we will use the term “ductile” in the sense described by Rutter (1986), whereby it refers to macroscale distributed flow, regardless of the nature of the deformation mechanism, and will use “brittle” to describe fracturing processes at all scales.

In nature, the brittle-ductile transition zone has been identified in exhumed shear zones showing markers of crystal plasticity (e.g., mylonites) overprinted by slip planes and pseudotachylytes that are inherent to the brittle regime (e.g., Sibson, 1980; Passchier, 1982; Hobbs et al., 1986). Such field evidence suggests that the transition in deformation mechanism is associated with a change in the degree of strain localization, from narrow frictional slip zones to wider plastic shear zones.

Laboratory experiments have shown that the transition from localized fracture to ductile flow generally occurs when the frictional strength of the fault, σf, equals the bulk flow stress of the rock, σflow (Byerlee, 1968; Kohlstedt et al., 1995). However, distributed deformation at the macroscopic scale may still be dominated by brittle microscale processes, and only further increases in pressure and temperature lead to fully crystal-plastic flow. This shows that the macroscale transition in strain localization (localized-ductile transition) is not necessarily the same as the microscale transition in deformation mechanism (brittle-plastic transition) and that the two transitions can occur under different pressure and temperature conditions. The resulting complex interplay between brittle and plastic mechanisms makes the flow stress σflow sensitive to a large number of parameters in the ductile regime (see Evans and Kohlstedt, 1995, and references therein), notably the imposed strain rate and the accumulated strain.

Furthermore, the criterion σflow > σf for the onset of ductile deformation was originally established from studies on initially intact materials undergoing a simple monotonic loading history, and describes deformation regimes in a binary manner (localized or distributed) without emphasizing the potential for coexistence of both fault slip and bulk ductile flow. The applicability of this criterion to the crust might therefore be limited, because crustal-scale deformation is controlled by preexisting structures (faults and shear zones; see, e.g., Goetze and Evans, 1979; Brace and Kohlstedt, 1980). Thus, it remains unclear if and how faults are reactivated across the brittle-ductile transition. Previous experimental studies have commonly used sample geometries that enforce sliding on narrow shear zones between essentially rigid blocks under increasing pressure and temperature conditions (e.g., Shimamoto, 1986; Pec et al., 2016), which do not allow for quantification of partitioning between fault slip and bulk strain.

Here, we conducted rock deformation experiments on pre-faulted samples of Carrara marble and monitored strain partitioning and fault reactivation across the localized-ductile transition. Our experiments were performed at room temperature and confining pressures (Pc) from 5 to 80 MPa. We determined partitioning of the total shortening between fault slip and off-fault matrix strain by subtracting the matrix strain (measured with strain gauges) from the total shortening (measured with external displacement transducers).

Experiments were conducted in two stages. During the first stage, samples were pre-faulted by loading at Pc = 5 MPa until localized brittle failure occurred. Following failure, an additional increment of shortening graphic = ΔL/L (L—length) of either 0.1% or 1% was allowed to accumulate before proceeding to the second stage, in order to test any effect of accumulated fault slip on the transition. In the second stage, Pc was increased stepwise from 5 to 80 MPa in 5 or 10 MPa increments. At each pressure step, the samples were reloaded at an axial shortening rate of graphic = 10−5 s−1 until 0.1% of irrecoverable axial shortening was accumulated, and then unloaded before proceeding to the next pressure step (see Section DR1 and Fig. DR2 in the GSA Data Repository1 for an extended methodology, and Table DR3 for a summary of experimental conditions).

During the first stage (Fig. 1), the sample behaves in a manner typical of the brittle regime, and the stress drop (accompanied by partial relaxation of the off-fault elastic strain) marks the formation of the macroscopic shear fault. During the second stage, at each confining pressure step, the stress-shortening relationship is initially linear, but deviates from linearity at some threshold stress σy, and then tends to plateau (Fig. 1A). This “plateau” stress increases significantly with increasing Pc. At low Pc (10 and 20 MPa), the matrix strain (εmatrix) initially increases at the same rate as the total shortening (ΔL/L), then deviates toward a constant value. The deviation point occurs at a stress denoted σf, and marks the onset of fault slip (triangles in Fig. 1B). At intermediate Pc (30–60 MPa), the same deviation is observed to occur, but εmatrix continues to increase beyond this point, albeit at a lower rate, indicating contributions from both matrix strain and fault slip to the total shortening. This observation appears to be independent of shortening, as demonstrated in an additional experiment where a single, second-stage deformation cycle was performed at Pc = 35 MPa, which shows no further deviation in matrix strain for a total shortening of up to a further 2% (Fig. DR4). Finally, at the highest Pc (70 MPa and above), εmatrix remains equal to ΔL/L throughout the deformation cycle, which implies that the fault is fully locked.

To assess the extent of microcracking in the matrix, we measured the horizontal P-wave speed across the fault during each deformation cycle (Fig. DR5). The wave speed at the start of each cycle increased with confining pressure. During deformation, the wave speed changed very little for cycles at Pc <30 MPa, but decreased progressively for all cycles at higher pressures. The magnitude of the decrease in P-wave speed increased with increasing Pc from 30 to 60 MPa but then decreased at higher confinement.

At Pc = 10 and 20 MPa, the yield stress and the fault strength are equal, and the calculated slip contributes close to 100% of the total shortening (Figs. 2A and 2D). Between Pc = 30 MPa and Pc = 60 MPa, σf increases linearly with Pc, whereas σy remains approximately constant at ∼115 MPa. Over this pressure range, the slip contribution progressively decreases from ∼80% at Pc = 30 MPa down to ∼15% at Pc = 60 MPa. At Pc = 70 MPa and above, the fault is fully locked, σf becomes inaccessible, and the slip contribution drops to zero. During the experiment where more slip is accumulated on the fault (1% rather than 0.1%) prior to stage 2 (Figs. 2B and 2E), σf and σy behave in a comparable manner to that described above, but σf increases with increasing Pc at a slightly higher rate. As a result, the deviation between the two initiates at Pc = 20 MPa and the fault becomes fully locked around Pc = 55 MPa. Similarly, the slip contribution decreases from >60% at Pc = 20 MPa to 20% at Pc = 45 MPa. During the experiment at the higher shortening rate of 10−4 s−1 (Figs. 2C and 2F), the trend remains the same, but the Pc domain over which σf = σy extends up to 35 MPa. From Pc = 40 MPa to Pc = 60 MPa, σf continues to increase linearly with increasing Pc, and σy remains approximately constant at 135 MPa. The slip contribution decreases from ∼80% at Pc = 40 MPa to 0% at Pc = 60 MPa. At the lower shortening rate of 10−6 s−1, the stress at the onset of fault slip σf does not differ significantly from that at higher shortening rates. By contrast with the test performed at 10−4 s−1, where the decrease in slip contribution initiates at Pc ≈ 40 MPa, at the lower rate of 10−6 s−1, that decrease initiates at Pc ≈ 15 MPa.

Our results show that with increasing confining pressure, faulted Carrara marble samples gradually shift from purely localized behavior where most of the deformation is accommodated by slip on the fault, to ductile behavior where strain is homogeneously distributed throughout the sample and the fault is locked. The transition commences at the confining pressure where fault strength becomes larger than matrix yield stress (σf > σy), and terminates when fault strength becomes equal to matrix flow stress (σf = σflow) (Figs. 1A, 1B, and 2). Thus, a transitional behavior where both fault slip and matrix deformation coexist occurs over a range of conditions delimited by σy < σf < σflow.

When σf = σy, no matrix strain is recorded (confirmed by the absence of significant variations in P-wave speed), and the yield stress of the rock is controlled by fault friction alone. This can be explained by the fact that, at low Pc, the fault frictional strength is likely lower than the yield stress of the off-fault matrix material (Fredrich et al., 1989). However, when σf > σy, the rock initially yields in the matrix and deformation is entirely ductile. The associated decrease in P-wave speed indicates that this ductility is driven mostly by diffuse microcracking. However, upon further loading, strain hardening eventually leads to reactivation of the fault when the applied stress reaches σf (confirmed by the existence of a single fault plane in post-mortem samples; Figs. DR6 and DR7). After reactivation, both ductile matrix strain and fault slip operate simultaneously, and partitioning of the total shortening between them is proportional to the ratio (σf − σy)/(σflow − σy), regardless of shortening rate and initial fault slip (Fig. 3). When σf ≥ σflow, the fault is locked and the deformation is fully ductile. The decrease in magnitude of the drop in P-wave speed under these conditions suggests that the contribution of microcracking to the overall deformation decreases with respect to that of crystal plasticity (Fredrich et al., 1989).

Our observations highlight the key role of the yield stress in the partitioning between localized fault slip and bulk deformation of the matrix. In Carrara marble, the control on yield stress switches from microcracking to crystal plasticity at low Pc (∼50 MPa; Fredrich et al., 1989; Fig. 4). This is corroborated by the pressure-insensitive behavior exhibited by our yield stress data at Pc >40 MPa (Figs. 2A–2C). Remarkably, the impact of strain rate on the partitioning of deformation is well captured, to first order, by the rate dependency of yield stress only (Fig. 3).

Our results are compatible with those of previous studies on silicate rocks using initially intact samples, where a similar progression from initial ductile yielding to strain localization and faulting with increasing deformation has been reported for conditions approaching the brittle-ductile transition (Hirth and Tullis, 1994). Additionally, the coexistence of ductile flow and localized shear zones has been observed in granite and feldspar aggregates (Tullis and Yund, 1977, 1992).

The existence of a zone of transitional behavior delimited by the yield stress can be integrated into a crustal-strength profile model (e.g., Kirby, 1980; Brace and Kohlstedt, 1980; Sibson, 1983; Fig. 4). Because yield stress is systematically lower than the flow stress, it appears that the transitional regime where ductile and localized strain coexist extends toward shallower depths compared to previous models of the brittle-ductile transition, into a depth range usually considered to be fully localized. In this zone, crustal strength is still controlled by fault friction, but with increasing depth, a growing proportion of the strain can be accommodated off-fault as the yield stress diverges from the frictional strength. This would suggest an overall widening of the shear zone, which is consistent with geological (e.g., Sibson, 1977; Scholz, 1988; Shimamoto, 1989; Cooper et al., 2010, 2017) and geophysical (e.g., Cowie et al., 2013) observations. Furthermore, high strain rates during seismic and post-seismic slip would increase both yield and flow stresses, therefore shifting the transition zone to greater depth. This is consistent with the existence of a zone of alternating behavior as discussed by Scholz (1988) and the formation of complex overprinted brittle and ductile structures observed in nature (e.g., Sibson, 1980; Melosh et al., 2014). Conversely, lower strain rates during the interseismic period (10−12 s−1 to 10−15 s−1) would reduce yield and flow stresses, which would in turn promote ductile deformation by shifting the transition zone to shallower depths. In this region of the crust, fault reactivation is dependent on the ability of the crust to harden with increasing strain. If recovery mechanisms are active, it is possible that large amounts of tectonic strain can be accommodated off-fault during transient deformation episodes, and if recovery is predominant, fault reactivation never occurs. Therefore, the gray area in Figure 4 represents all possible stress states in the crust. This rheology could explain the abnormally low stresses recorded around major faults (e.g., Behr and Platt, 2014), but the mechanisms responsible for low-temperature strain hardening and recovery are, to date, mostly unknown.

Unfortunately, there is a paucity of systematic data on low-temperature yield stress in crustal rocks. However, laboratory studies on wet quartz single crystals (e.g., Balderman, 1974; Doukhan and Trépied, 1985) suggest low-temperature yield stresses on the order of 50–100 MPa, which would imply a transition zone depth of only a few kilometers in continental crust.

We thank Whitney Behr, Matej Pec, and Nicola De Paola for their constructive reviews. This project received funding from the European Union’s Horizon 2020 research and innovation program under the Marie Sklodowska-Curie grant agreement 642029–ITN CREEP, and under the European Research Council grant agreement 804685–RockDEaF.

1GSA Data Repository item 2019391, extended methodology, supplemental figures, and data, is available online at http://www.geosociety.org/datarepository/2019/, or on request from editing@geosociety.org.
Gold Open Access: This paper is published under the terms of the CC-BY license.