The emergence of crack-like behavior of frictional rupture: The origin of stress drops

Rapid slip along interfaces separating bodies in frictional contact is mediated by the spatiotemporal dynamics of frictional rupture [1, 2]. Frictional rupture is a fundamental process of prime importance for a broad range of physical systems, e.g. it is responsible for squealing in car brake pads [3], for bowing on a violin string [4], and for earthquakes along geological faults [5, 6, 7], to name just a few well-known examples. The essence of frictional rupture propagation is that a state of relatively high slip rate (the rate of interfacial shear displacement discontinuity) behind the rupture edge propagates into a low/vanishing slip rate state ahead of it, cf. Fig. 1. As such, frictional rupture appears to be essentially similar to ordinary tensile (opening) cracks, where a finite tensile displacement discontinuity (broken material) state behind the crack edge propagates into a zero tensile displacement discontinuity (intact material) state ahead of it [8].

There is, however, an important and fundamental difference between frictional rupture and ordinary tensile cracks that manifests itself in the stress states associated with these two processes. A tensile crack, i.e. a crack subjected to opening forces, is composed of surfaces that cannot support stress, so the stress behind its edge vanishes. Consequently, tensile (opening) crack propagation is a process in which the applied stress ahead of the crack edge drops to zero behind it. This stress drop, which accompanies tensile crack propagation, has dramatic implications. Most notably, the loss of stress-bearing capacity along the crack surfaces is compensated by large concentration of deformation and stress near the crack edge, oftentimes in a way that mimics a mathematical singularity, whose intensity increases with increasing stress drops [8]. Frictional rupture is different from tensile cracks because the finite frictional interaction between the two bodies in contact behind the rupture edge generically implies that the stress there cannot drop to zero, but rather remains finite.

The close relations between frictional rupture and tensile cracks can be maintained if, as is widely assumed, the stress behind the frictional rupture edge — the residual stress $\tau_{\rm res}$ — is well-defined and is generically smaller than the far-field stress $\tau_{\rm d}$ that is required to drive rupture. Moreover, the residual stress $\tau_{\rm res}$ is generally assumed to be an intrinsic interfacial property of the slipping contact interface, typically related to the kinetic friction coefficient. Under these assumptions, a finite stress drop $\Delta\tau\!\equiv\!\tau_{\rm d}-\tau_{\rm res}\!>\!0$ exists and effective crack-like properties of frictional rupture, e.g. edge singularity, are expected to emerge. These assumptions have been adopted in an extremely broad range of theoretical and numerical studies [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27], and their implications have been consistent with geophysical observations [28, 26] and have been confirmed in some recent laboratory experiments [29, 30, 31, 32, 33, 34, 35, 36]. In fact, stress drops are among the few remotely observable parameters in earthquake science, providing a key link to the frictional properties of faults, which can be only indirectly inferred.

Yet, to the best of our knowledge, currently there is no basic understanding of how and under what conditions the effective crack-like behavior of frictional rupture emerges from fundamental physics. More specifically, there is a need to understand what the physical origin of a finite stress drop $\Delta\tau\!>\!0$ is and under what conditions it emerges. Here we address these basic questions; first, we show that for rapid rupture a finite and well-defined stress drop is not an interfacial property, as is widely assumed, but rather it is directly related to wave radiation from the frictional interface to the bodies surrounding it (the so-called radiation damping effect [37, 38, 39, 40]) and to long-range bulk elastodynamic interaction effects. Second, we show that the emergence of a stress drop is a finite time effect, limited by the wave travel time in finite systems. Third, we show that for slow rupture, i.e. rupture that is significantly slower than the elastic wave-speeds [41, 42, 43, 44], stress drops are transiently controlled by the long-range quasi-static elasticity of the bodies surrounding the frictional interface. Reanalysis of very recent experimental results, reported by two different experimental groups, provides strong support to our theoretical predictions, for both rapid and slow rupture. All in all, our findings elucidate the origin of stress drops in frictional rupture, i.e. they offer a comprehensive and fundamental understanding of why, how and to what extent frictional rupture might be viewed as an ordinary fracture process.

The starting point for our discussion is a physically-motivated interfacial constitutive law, i.e. a relation between the dynamical and structural variables that characterize a frictional interface and the frictional resistance stress $\tau$ [45]. A frictional interface is formed when two bodies come into contact. Each of them satisfies its own continuum momentum balance equation $\rho\ddot{\bm{u}}(\bm{r},t)\!=\!\nabla\!\cdot{\bm{\sigma}}(\bm{r},t)$, where $\rho$ is the mass density, $\bm{u}$ and $\bm{r}\!=\!(x,y)$ (in two-dimensions) are the displacement and position vector fields respectively, and $\bm{\sigma}$ is the stress tensor field (a superposed dot represents a time derivative). $\bm{\sigma}$ in each body is related to $\bm{u}$ through a bulk constitutive law (i.e. a constitutive law that characterizes the bodies forming the interface), oftentimes Hooke’s law of linear elasticity, to be adopted below as well. Note that body forces are neglected in the momentum balance equation.

The interfacial constitutive law involves three bulk quantities evaluated at the interface located at $y\!=\!0$: (i) the slip rate/velocity $v(x,t)\!\equiv\!\dot{u}_{x}(x,y\!=\!0^{+},t)\!-\!\dot{u}_{x}(x,y\!=\!0^{-},t)$, where $+/-$ correspond to the upper/lower bodies, respectively (ii) the shear stress $\sigma_{xy}(x,y\!=\!0,t)$ that is balanced by the frictional stress, $\tau(x,t)\!=\!\sigma_{xy}(x,y\!=\!0,t)$ and (iii) the normal stress $\sigma(x,t)\!\equiv\!-\sigma_{yy}(x,y\!=\!0,t)$. A large body of evidence accumulated in the last few decades indicates that the interfacial constitutive law must also involve a set of non-equilibrium order parameters $\{\phi_{i}\}$, sometimes termed internal-state fields, that represent the structural state of the interface and encode its history [46, 47, 48, 49, 45, 50, 51, 52]. In a minimal formulation, adopted in numerous studies [53, 38, 54, 55, 39, 56, 57, 58], a single internal-state field $\phi(x,t)$ is used. This assumption is adopted here, without loss of generality.

The interfacial constitutive law, at any position $x$ along the interface and at any time $t$, is described by the following local relation

which must be supplemented with a dynamical equation for the evolution of $\phi$. Extensive evidence indicates that $\phi$ physically represents the age/maturity of the contact [47, 48, 49, 45, 50, 51, 52] and that its evolution takes the form

with $g(1)\!=\!0$ and where $\phi$ is of time dimension. The characteristic slip displacement $D$ controls the transition from a stick state $v\!\approx\!0$, with a characteristic structural state $\phi\!=\!\phi_{0}$, to a steadily slipping/sliding state $v\!>\!0$, with $\phi_{\rm ss}\!=\!D/v$. The precise functional form of $g(\cdot)$ (with $g(1)\!=\!0$) plays no role in what follows.

The function $f(|v|,\phi_{\rm ss}\!=\!D/v)\!=\!\tau_{\rm ss}(v)/\sigma$, under steady-state sliding conditions and a controlled normal stress $\sigma$, has been measured over a broad range of slip rates $v$ for many materials [45]. Together with general theoretical considerations [59], it is now established that the steady-state frictional stress $\tau_{\rm ss}(v)$ is generically $N$-shaped, as shown in Fig. 2a. Consider then a frictional system driven by a shear stress $\tau_{\rm d}$, which is larger than the minimum of the $\tau_{\rm ss}(v)$ curve, cf. Fig. 2a. What are the generic properties of frictional rupture that might emerge under these conditions?

As explained in Sect. I, frictional rupture is a propagating spatiotemporal object that features a relatively high slip state $v\!>\!0$ behind its edge and a stick (no/slow slip) $v\!\approx\!0$ state ahead of it, as shown in Fig. 1. This spatiotemporal dynamical process can be directly related to Fig. 2a, where the driving stress $\tau_{\rm d}$ intersects the $\tau_{\rm ss}(v)$ curve at three points. The leftmost intersection point features an extremely small slip velocity, $v\!\approx\!0$, which corresponds to the state ahead of the rupture edge. We stress that in general the state of the interface ahead of the rupture edge may be far from steady-state and the results to follow are largely independent of its detailed properties. The rightmost intersection point features a relatively large slip velocity, $v\!>\!0$, which corresponds to the state behind the rupture edge. The transition between these two states takes place in the edge region, and is controlled by $\dot{\phi}$ in Eq. (2) and by spatiotemporal bulk dynamics. In this transition region, the slip velocity $v$ also goes through the intermediate intersection point, which is not a stable fixed point as the other two. The crucial observation is that the stress behind steadily propagating frictional rupture is $\tau_{\rm d}$, i.e. the residual stress equals to the driving stress, $\tau_{\rm res}\!=\!\tau_{\rm d}$. This implies that we expect no stress drop to emerge at all, $\Delta\tau\!=\!\tau_{\rm d}-\tau_{\rm res}\!=\!0$, and consequently no crack-like behavior.

In many studies available in the literature, a steady-state friction curve $\tau_{\rm ss}(v)$ that does not feature a minimum is adopted [38, 39, 60, 61, 62]. We consequently discuss here such a no-minimum friction curve and plot an example of it in Fig. 2a. In this case, the driving stress $\tau_{\rm d}$ behind the rupture edge cannot be balanced by the friction stress and the slip velocity in this region is expected to continuously accelerate. As such, we cannot expect a well-defined steady-state stress drop to emerge, though the stress will definitely drop below $\tau_{\rm d}$ behind the rupture edge.

The discussion above, both for the $N$-shaped and the no-minimum steady-state friction curves, seems to lead to the quite remarkable conclusion that based on basic physics considerations we expect no finite and well-defined steady-state stress drops to emerge at all in the context of friction rupture, and hence no crack-like behavior as well. This appears to be in sharp contrast to ample evidence indicating the existence of finite and well-defined steady-state stress drops in various frictional systems [29, 30, 32, 33, 63, 64, 35, 1]. How can one reconcile the two apparently conflicting conclusions?

To address this question, let us write down in more detail the general expression for interfacial shear stress $\sigma_{xy}(x,y\!=\!0,t)$, valid also out of steady-state, when a constant driving stress $\tau_{\rm d}$ is applied at the far boundaries of the systems, say at $y\!=\!\pm H$ ($H$ is the height of each of the two bodies in contact). For bulk linear elastodynamics, we have $\sigma_{xy}(x,y\!=\!0,t)\!=\!\tau_{\rm d}+\tilde{s}(x,t)$, where $\tilde{s}(x,t)$ is a spatiotemporal integral that quantifies the long-range (in both space and time) elastodynamic interaction between different parts of the interface [65, 66, 67]. Under strict homogeneous steady-state conditions, we have $\tilde{s}(x,t)\!\to\!0$ and consequently $\sigma_{xy}(x,y\!=\!0,t)\to\tau_{\rm d}$, which corresponds to the rightmost intersection point Fig. 2a attained far behind the rupture edge (see discussion above). At finite times, before strict steady-state conditions are attained, the spatiotemporal integral term $\tilde{s}(x,t)$ makes a finite contribution to $\sigma_{xy}(x,y\!=\!0,t)$, which quantifies the deviation from steady-state.

Under these conditions, and in particular for times in which information regarding the evolution of the slip velocity $v(x,t)$ relative to some initial/reference slip velocity $v_{0}$ does not have enough time to propagate to the boundaries at $y\!=\!\pm H$ and back to the interface, the spatiotemporal integral term $\tilde{s}(x,t)$ can be decomposed into two contributions, one is a local contribution of the form $\frac{\mu}{2c_{s}}\!\left(v(x,t)-v_{0}\right)$, where $\mu$ is the linear elastic shear modulus and $c_{s}$ is the shear wave-speed, and the other is a non-local (in space and time) contribution $s(x,t)$ [65, 66, 67]. This decomposition is valid for times shorter than ${\cal O}(H/c_{s})$, before wave interaction with the boundaries is possible, and for these times the interfacial shear stress takes the form [65, 68, 55, 66, 67, 57]

In many studies available in the literature the idealized infinite system limit $H\!\to\!\infty$ is considered, for which Eq. (3) is valid at all times. The term $\frac{\mu}{2c_{s}}\!\left(v(x,t)-v_{0}\right)$ physically represents wave radiation from the interface to the bodies that form it and is therefore known as the radiation damping term [37, 38, 39, 40]. It is associated with “damping” because from the perspective of the interface it acts as a viscous stress with $\mu/2c_{s}$ being the effective viscosity. This term makes an important contribution to stress drops in frictional rupture, as is shown next.

Consider a point along a frictional interface that is initially located ahead of a propagating frictional rupture and whose slip velocity is $v\!\approx\!0$, which represents $v_{0}$ in Eq. (3). When the frictional rupture goes through this point, the stress and slip velocities vary significantly. Suppose then that the system height $H$ is sufficiently large such that the spatiotemporal integral decays once the rupture went sufficiently far ahead, $s(x,t)\!\to\!0$, but the radiation damping contribution is still valid (i.e. shear waves did not have enough time to propagate to the far boundaries and back). Under these conditions, the slip velocity $v_{\rm res}$ at the spatial point under consideration, now far behind the frictional rupture, is determined by the shear stress balance $\sigma_{xy}\!=\!\tau$ and satisfies

where $v^{0}_{\rm res}\!\gg\!v_{0}$ and $s(x,t)\!\to\!0$ have been used. Note that the superscript ’$0$’ in $v^{0}_{\rm res}$ represents the fact that this is the theoretically predicted residual slip velocity under the assumption that $s(x,t)\!=\!0$ far behind the propagating rupture front. The residual stress at this point takes the form $\tau_{\rm res}\!=\!\tau_{\rm d}-\frac{\mu}{2c_{s}}v^{0}_{\rm res}$, and consequently a finite stress drop of magnitude

is expected to emerge on times shorter than ${\cal O}(H/c_{s})$.

A geometric representation of Eqs. (4)-(5) is shown in Fig. 2b for the $N$-shaped $\tau_{\rm ss}(v)$ and in Fig. 2c for the no-minimum $\tau_{\rm ss}(v)$, both shown in Fig. 2a. In Figs. 2b-c, the left-hand-side of Eq. (4) $\tau_{\rm ss}(v)+\frac{\mu}{2c_{s}}v$ is regarded as an effective steady-state curve and is plotted by a dashed line. The radiation damping contribution $\frac{\mu}{2c_{s}}v$ in Fig. 2b shifts the location of the effective (finite time) steady-state slip rate to lower rates (compared to the strict steady-state represented by the rightmost intersection point in Fig. 2a). In Fig. 2c it gives rise to an effective (finite time) steady-state slip rate, which simply does not exist for the no-minimum curve in Fig. 2a. This shows that Eq. (5) is valid independently of the properties of $\tau_{\rm ss}(v)$. Note that relevant solutions $v^{0}_{\rm res}$ to Eq. (4) exist only if $\tau_{\rm d}$ is larger than the effective (finite time) minimum of the steady-state friction curve, as is highlighted in Figs. 2b-c. The very same condition played a central role in the analysis of [39], where the conditions for rupture mode selection (self-healing pulses vs. crack-like) have been extensively discussed. Somewhat related issues have also been discussed in [71].

We would like to note that once $\Delta\tau$ is assumed to exist, the relation in Eq. (5) between $\Delta\tau$ and $v_{\rm res}$ as given quantities, is a known elastodynamic relation [72], which has previously received some experimental support (e.g. see Fig. 8 in [73]) and is strongly supported by recent experimental data we have extracted from recent experimental work [35, 36], see Sect. IV. To the best of our knowledge, however, none of these works (except for [39], as mentioned above) predicted the residual slip velocity $v_{\rm res}$ to emerge from the combined effect of the steady-state friction curve $\tau_{\rm ss}(v)$ and of the radiation damping term $\frac{\mu}{2c_{s}}v$ as predicted in Eq. (4), and none of them interpreted Eq. (5) as a major contribution to stress drops $\Delta\tau$ in frictional rupture. In the context of a quasi-static fault model of [74], the radiation damping term $\frac{\mu}{2c_{s}}v$ has been added in an ad hoc manner in order to avoid unbounded slip velocities to emerge during frictional instabilities in the quasi-static formulation. It has been noted in this context [74] that the effective viscosity $\mu/2c_{s}$ affects the magnitude of the stress drop $\Delta\tau$ when $\tau_{\rm ss}(v)$ has no-minimum (cf. Fig. 2c); nevertheless, the general physical picture in which the radiation damping term significantly contributes to stress drops in frictional rupture, independently of the properties of $\tau_{\rm ss}(v)$, has not been discussed.

When the radiation damping term $\tfrac{\mu}{2c_{s}}v$ in Eq. (4) does not faithfully represent the physics of a given system, no well-defined finite stress drop is expected to accompany rapid frictional rupture. This can happen in two generic cases; first, in the limit of thin bodies $H\!\to\!0$, where essentially there is no bulk to radiate energy into, the radiation damping term $\tfrac{\mu}{2c_{s}}v$ simply does not exist to begin with. In this case, frictional rupture exists, but it is not accompanied by any stress drop, as shown previously in [75] and here in Fig. A1. Second, the radiation damping term $\tfrac{\mu}{2c_{s}}v$, which exists at relatively short times (shorter than ${\cal O}(H/c_{s})$), is expected to vanish in the long time limit $t\!\gg\!H/c_{s}$. This limit can be probed by performing experiments or simulations with long enough systems for long enough times. Indeed, simulations of effectively long systems yield rupture fronts with no stress drop, see Fig. 6a in [76] and details therein. Finally, note that the radiation damping term $\tfrac{\mu}{2c_{s}}v$ is expected to decrease to zero in discrete steps corresponding to each wave reflection from the system’s boundaries. Experimental evidence for the stepwise nature of the decrease in the radiation damping term (associated with discrete wave reflections) will be discussed below.

At this point, we first set out to test the predictions in Eqs. (4)-(8) against extensive numerical simulations. To that aim, we consider two semi-infinite bodies in frictional contact. The advantage of considering infinite-height bodies, i.e. the $H\!\to\!\infty$ limit, is that the interfacial relation in Eq. (3) becomes exact at all times, unlike for finite bodies. We also employ periodic boundary conditions, with periodicity $W$, in the sliding direction. We performed spectral boundary integral method [65, 66, 67] calculations under anti-plane shear (mode-III symmetry) deformation conditions, which are similar to — yet somewhat simpler than — the in-plane shear (mode-II symmetry) deformation conditions considered up to now [8]. The main simplification is that the displacement field in the mode-III problem ${\bm{u}}(x,y,t)\!=\!u_{z}(x,y,t)\hat{\bm{z}}$ (the unit vectors satisfy $\hat{\bm{z}}\,\bot\,\hat{\bm{x}},\hat{\bm{y}}$) is essentially scalar. The basic field $u_{z}(x,y,t)$ satisfies the bulk elastodynamic equation $\mu\nabla^{2}u_{z}\!=\!\rho\,\ddot{u}_{z}$, together with $v(x,t)\!\equiv\!\dot{u}_{z}(x,y\!=\!0^{+},t)\!-\!\dot{u}_{z}(x,y\!=\!0^{-},t)$ and $\tau(x,t)\!\equiv\!\sigma_{yz}(x,y\!=\!0,t)\!=\!\mu\,\partial_{y}u_{z}(x,y\!=\!0,t)$. Equation (3) remains valid, where $\sigma_{xy}(x,y\!=\!0,t)$ is replaced by $\sigma_{yz}(x,y\!=\!0,t)$ and the integral term $s(x,t)$ corresponds to mode-III, see Appendices for more details.

The employed interfacial constitutive law features the generic properties discussed above, with $f(|v|,\phi)$ of Eq. (1) that reduces under steady-state conditions to either the $N$-shaped or the no-minimum curves of Fig. 2a (the exact expressions for $f(|v|,\phi)$ can be found in Appendix A2), and with $g(\cdot)$ of Eq. (2) that is given by $g\!=\!1-|v|\phi/D$ [46, 45, 52, 5, 49]. The bodies are loaded by a constant driving stress $\tau_{\rm d}$, as depicted schematically in Fig. 2a, and frictional rupture is nucleated by introducing Gaussian perturbations of proper amplitude into a homogeneous state of very low slip velocity $v_{0}$ (that corresponds to the leftmost intersection point in Fig. 2a), following the theoretical framework of [80] (the details of the nucleation procedure are described in the Appendices).

An example of the stress distribution of an emerging frictional rupture is shown in Fig. 3a. The figure reveals two rapid rupture fronts propagating in opposite directions (at $84\%$ of the shear wave-speed $c_{s}$), where the stress ahead of the two fronts is the applied stress $\tau_{\rm d}$. As the emerging rupture is rapid, i.e. propagating at a speed comparable to the elastic wave-speed, the relevant prediction is given by Eq. (5). As predicted, the observed stress left behind the two rapid fronts, $\tau_{\rm res}$, is constant and smaller than the driving stress $\tau_{\rm d}$, giving rise to a finite stress drop $\Delta\tau$. Following the discussion above, since in these calculations $H\!\to\!\infty$, the finite stress drop persists indefinitely (while in finite size systems it persists for times $\sim{\cal O}(H/c_{s})$, cf. the experiments of [35], to be discussed later). The stress drop $\Delta\tau$ observed in Fig. 3a quantitatively agrees with the prediction in Eq. (5), as stated in the figure legend.

In Fig. 3b we present the slip velocity distribution that corresponds to the snapshot shown in Fig. 3a. As predicted, the slip velocity attains a plateau level $v_{\rm res}$ behind the propagating rupture fronts, and the residual velocity $v_{\rm res}$ is used to normalize the slip velocity distribution. The simulationally-measured residual velocity $v_{\rm res}$ deviates by $\sim\!35\%$ from the theoretical prediction in Eq. (5), in quantitative agreement with the estimate provided in Sect. II.2 for the spatiotemporal integral contribution to the stress drop in this case (see also Appendix A7). Consequently, while the upper equation in (6) implies that $s$ is in fact not very small, the prediction for $\Delta\tau$ works perfectly fine due to the $\epsilon\!\ll\!1$ factor in the lower equation in (6).

In order to quantitatively test the latter prediction for rapid rupture over a range of physical conditions, we first solved Eq. (4) to obtain $v^{0}_{\rm res}(\tau_{\rm d})$ and then plugged it in Eq. (5) to obtain $\Delta\tau(\tau_{\rm d})$, where the latter is plotted in solid lines in Fig. 3c (the two solid lines correspond to both the $N$-shaped or the no-minimum steady-state friction laws shown in Fig. 2a). We then numerically calculated $\Delta\tau$ (as demonstrated in Fig. 3a) for various driving stresses $\tau_{\rm d}$, for both the $N$-shaped or the no-minimum steady-state friction curves shown in Fig. 2a. The numerical results (discrete symbols) are superimposed on the theoretical prediction in Fig. 3c (the lowest numerical data point on the lower curve corresponds to panels (a)-(b)). The agreement between the theoretical prediction and the numerical results for $\Delta\tau$ is very good, where it is better for the $N$-shaped steady-state friction law (lower curve) than for the no-minimum steady-state friction law (upper curve). This difference is fully accounted for by the magnitude and sign of $\epsilon$, directly affected by the $\frac{\partial{f}_{\rm ss}}{\partial\log(v^{0}_{\rm res})}$ term (cf. Eq. (6)).

The contribution of the spatiotemporal integral $s(x,t)$, which is smaller than the radiation damping contribution for rapid rupture, may be a dominant effect for slow rupture, as discussed around Eqs. (7)-(8). To test this possibility, we generated slow rupture by changing the frictional parameters, the loading level and the nucleation procedure, as explained in Appendix A5. An example of such slow rupture is shown in Fig. 4 (solid line), exhibiting a rupture front propagating at a velocity $2$ orders of magnitude smaller than $c_{s}$ and leaving behind it a stress drop $\Delta\tau$. We first verified that $\frac{2c_{s}\Delta\tau}{\mu\,v_{\rm res}}\!\gg\!1$, as predicted by Eq. (7) (and as stated in the figure legend). To test whether indeed $\Delta\tau\!\simeq\!-s(t)$, cf. Eq. (7), we used the slip displacement gradient $\partial_{x}\delta(x,t)$ obtained from the dynamical simulation and used it to calculate the quasi-static integral in Eq. (8) for any point $x$ along the interface at $t$ corresponding to the snapshot in Fig. 4. The result for $\tau_{\rm d}+s(x,t)$ is superimposed on $\tau(x,t)$ of Fig. 4 (dashed line), demonstrating excellent agreement with the fully dynamic result, and in particular validating $\Delta\tau\!\simeq\!-s(t)$ for this $t$. Consequently, our simulations strongly support the theoretical predictions in Eqs. (4)-(8), for both rapid and slow frictional rupture. Next, we discuss direct experimental support for these predictions.

In the previous section, we provided strong simulational support to the theoretical predictions in Eqs. (4)-(8). Our goal here is to test these predictions against experimental data. Equation (4) predicts the slip velocity $v_{\rm res}$ behind the frictional rupture once the steady-state friction curve $\tau_{\rm ss}(v)$ is known. The latter is not always known a priori over a sufficiently wide range of steady-state slip velocities. In fact, measuring both the frictional stress $\tau_{\rm res}$ and the slip velocity $v_{\rm res}$ behind rupture fronts allows to extract $\tau_{\rm ss}(v)$. In this case, any triplet $(\tau_{\rm d},\tau_{\rm res},v_{\rm res})$ is predicted to follow either Eq. (5) for rapid rupture or Eq. (7) for slow rupture, independently of the steady-state friction law $\tau_{\rm ss}(v)$.

Measurements of both $\tau_{\rm res}$ and $v_{\rm res}$ behind rupture fronts for various $\tau_{\rm d}$ have been recently performed by two independent experimental groups using two different experimental systems and techniques [35, 36]. The first focused on the frictional dynamics along the interface between two blocks of Homalite probed through a novel ultrahigh full-field imaging technique [35]. The second focused on the frictional dynamics along the interface between two blocks of poly(methylmethacrylate) (PMMA) probed through a combination of high speed interfacial imaging (via a method of total internal reflection) and simultaneous measurements of the deformation fields slightly above the interface [36]. Here we use the data reported in these works to quantitatively test our predictions.

We start with [35], where the observed rupture fronts were all in the rapid (fully inertial) regime, and hence the relevant prediction is given in Eq. (5). To test this prediction, we extracted the relevant experimental data from [35] and plotted in Fig. 5 the resulting $\Delta\tau$ vs. $v_{\rm res}$ against the linear $\Delta\tau(v_{\rm res})$ relation of Eq. (5), with a slope that corresponds to $\mu/2c_{s}$ of Homalite (see figure caption for additional details). The results reveal excellent agreement between the experimental data and the theoretical prediction, without any free parameter. The theoretical relation in Eq. (5) predicts $v^{0}_{\rm res}$ by assuming that the long-range spatiotemporal contribution $s(x,t)$ is negligibly small, while in the experimental data we just use the measured $v_{\rm res}$. The fact that there are no significant horizontal deviations of the experimental data from the theoretical line of Eq. (5), i.e. that $v_{\rm res}$ is close to $v^{0}_{\rm res}$, indicates that $s(x,t)$ is smaller in these experiments compared to the estimate provided in Sect. II.2. The origin of the latter discrepancy is not yet clear.

The experimental data included in Fig. 5 have been obtained in the short time regime, before any wave reflection from the system’s boundaries, for which Eqs. (4)-(5) are valid. Yet, for a single case, measurements are reported after the first wave reflection (but before the second one), cf. Fig. 4c in [35]. Under these conditions, we expect the ordinary radiation damping term $\tfrac{\mu}{2c_{s}}v$ to be replaced by a smaller term. Consequently, Eq. (4) predicts that $v_{\rm res}$ should increase, exactly as is experimentally observed in Fig. 4c in [35] (denoted as “Reflected rupture” in the figure), thus providing direct evidence for the stepwise reduction of the radiation damping term with discrete wave reflections.

We now turn to the experiments of [36], which report on extensive measurements of $\tau_{\rm res}$ and $v_{\rm res}$, cf. Fig. 3a there. In order to test our predictions we need also measurements of the driving stress $\tau_{\rm d}$, which will allow us to extract the stress drop $\Delta\tau\!=\!\tau_{\rm d}-\tau_{\rm res}$. The corresponding extensive measurements of $\tau_{\rm d}$ are not presented in [36], though $\tau_{\rm d}$ is presented for two highly relevant examples in Fig. 2b there (where $\tau_{\rm d}$ is denoted by $\tau_{0}$). One example corresponds to very rapid rupture (in fact, it is supershear rupture, propagating at nearly the dilatational wave-speed) and the other to slow rupture (propagating at about $10\%$ of the Rayleigh wave-speed). These data are exactly what is needed in order to test our predictions, in particular regarding the change in behavior for slow rupture. For rapid rupture, we extracted from Fig. 2b of [36] (blue data) $\Delta\tau\!=\!1.17$ MPa and from Fig. 1d (bottom) $v_{\rm res}\!\simeq\!1.5$ m/s (see additional details in Appendix A6). Using the reported values $\rho\!=\!1170$ kg/m³ and $c_{s}\!=\!1345$ m/s [36], together with $\mu\!=\!\rho c_{s}^{2}$, we obtain $\frac{2c_{s}\Delta\tau}{\mu\,v_{\rm res}}\!=\!0.99$. This is in great quantitative agreement with the prediction for rapid rupture in Eq. (5), and is fully consistent with the independent experimental data presented in Fig. 5. This result indicates that the long-range spatiotemporal contribution $s(x,t)$ is indeed small also in the experiments of [36].

We then extracted the corresponding data for the slow rupture, obtaining from Fig. 2b of [36] (red data) $\Delta\tau\!=\!0.242$ MPa and from Fig. 1d (top) $v_{\rm res}\!\simeq\!0.02$ m/s. Using these data we obtain $\frac{2c_{s}\Delta\tau}{\mu\,v_{\rm res}}\!\simeq\!15\!\gg\!1$, in agreement with the theoretical prediction for slow rupture in Eq. (7). The results presented in the last two sections provide strong simulational and experimental support for our theoretical predictions and hence for the proposed picture of the physical origin of stress drops in frictional rupture. As stress drops are important for the possible emergence of crack-like behaviour in frictional rupture, future studies of the latter should be based on the results presented in this paper.

The possible deep relations between frictional rupture and ordinary fracture provide a powerful conceptual and quantitative framework to understand frictional dynamics in a wide variety of physical contexts. This framework is extensively used to interpret and quantify geophysical observations [28, 26], as well as a broad spectrum of laboratory phenomena [29, 30, 31, 32, 33, 34, 35]. For example, a recent series of careful laboratory experiments [32, 33, 34] demonstrated that when the analogy between frictional rupture and ordinary fracture holds, the dynamic propagation of laboratory earthquakes and their arrest can be quantitatively understood to an unprecedented degree [25]. Yet, the fundamental physical origin and range of validity of the analogy between frictional rupture and ordinary fracture are not yet fully understood. In this paper, we developed a comprehensive and fundamental understanding of why, how and to what extent frictional rupture might be viewed as an ordinary fracture process.

A important ingredient in the analogy is the emergence of a finite and well-defined stress drop $\Delta\tau\!=\!\tau_{\rm d}-\tau_{\rm res}$, the difference between the applied driving stress $\tau_{\rm d}$ and the residual stress $\tau_{\rm res}$, in frictional rupture. In the first part of the paper we showed that, contrary to widely adopted assumptions, the residual stress $\tau_{\rm res}$ is not a characteristic property of frictional interfaces. Rather, for rapid rupture $\tau_{\rm res}$ is shown to crucially depend on elastodynamic bulk effects — mainly wave radiation from the frictional interface to the bodies surrounding it, but also long-range elastodynamic bulk interactions (encapsulated in the integral term $s(x,t)$ in Eq. (3)) — and the applied driving stress $\tau_{\rm d}$ itself, in addition to the contribution of the slip rate dependence of the constitutive friction law. Notably, we showed that for rapid rupture the deviation of $\tau_{\rm res}$ from $\tau_{\rm d}$, i.e. the existence of a finite stress drop $\Delta\tau$, is a finite time effect, mainly limited by the wave travel time in finite systems. For slow rupture, it is shown that if a stress drop exists, it is intimately related to the long-range quasi-static elasticity of the bodies surrounding the interface, again not exclusively to interfacial physics. Our theoretical predictions are supported by extensive computations and existing experimental data from two independent laboratory experiments.

Our findings have important implications that go beyond their basic nature; first, the results show that the widely used slip-weakening models [9, 10], in which the existence of a residual stress $\tau_{\rm res}\!<\!\tau_{\rm d}$ is a priori assumed (as a fixed interfacial property), should be employed with care. In particular, as $\tau_{\rm res}$ has been shown to depend on the externally applied stress $\tau_{\rm d}$, on the properties of the bodies surrounding the interface and on the rate-dependence of the frictional constitutive behavior, $\tau_{\rm res}$ cannot be assumed to be fixed. Rather, it should be self-consistently calculated from the coupled interface-bulk problem. Somewhat related conclusions in relation to slip-weakening models, based on measurements of evolving local friction during spontaneously developing laboratory earthquakes, have been drawn in [35]. We also note that the dependence of $\tau_{\rm res}$ on the bulk constitutive relation and on long-range bulk-mediated interactions may give rise to interesting effects for more complicated bulk constitutive relations (e.g. viscoelasticity) and in the presence of repeated rupture events [82]. Such effects should be explored in future work.

The existence of a finite stress drop $\Delta\tau$ generically leads to accelerating frictional rupture under stress-controlled far-field loading conditions $\tau_{\rm d}$ if $\Delta\tau$ is independent of the rupture size $L$. In these situations, inertia-limited rapid rupture is expected to emerge on time scales for which wave interaction with the outer boundaries does not exist, and $\Delta\tau$ is controlled by elastodynamic bulk effects. These conditions are typically realized in many geological and laboratory earthquakes [28, 29, 30, 31, 32, 33, 26, 34, 35]. On the other hand, slow rupture propagation — a widely observed, yet highly debated and elusive phenomenon [41, 42, 43, 44] — is expected to feature smaller stress drops $\Delta\tau$ that may decrease with increasing rupture size $L$, such that rupture acceleration is limited.

The possible $L$ dependence of $\Delta\tau$ and its possible relations to the emergence of slow rupture should be further explored in the future. In addition, as the whole discussion in this paper is valid for stress-controlled far-field loading conditions (characterized by $\tau_{\rm d}$), future work should also consider velocity-controlled far-field loading conditions, where finite stress drops might emerge from different physical considerations [83]. Finally, as the existence of a finite stress drop $\Delta\tau$ does not in itself guarantee the crack-like behavior of frictional rupture, future work should clarify to what extent the analogy to ordinary cracks can in fact quantitatively account for the dynamics of frictional rupture. All in all, as stress drops are key quantities in frictional failure dynamics, we expect our results to provide a conceptual and quantitative framework to address various fundamental and applied problems in relation to the rupture dynamics of frictional interfaces, with implications for both laboratory and geophysical-scale phenomena.

Acknowledgements E. B. and J.-F.M. acknowledge support from the Rothschild Caesarea Foundation. E. B. acknowledges support from the Israel Science Foundation (Grant No. 295/16). J.-F.M., F. B. and T. R. acknowledge support from the Swiss National Science Foundation (Grant No. 162569). This research is made possible in part by the historic generosity of the Harold Perlman Family. We are grateful to Ilya Svetlizky for pointing out the differences between rapid and slow rupture in the context of our work, and the relevance of the experimental data of [36] to it. His comment has stimulated the discussion of slow rupture in this paper. In addition, we are grateful to him for pushing us to look at Broberg’s self-similar solutions, eventually leading to the extraction of the estimates presented in Sect. II.2.

The simulations discussed above relied on a spectral boundary integral formulation of the elastodynamic equations [65, 66, 67]. The latter relates the traction stresses acting along the interface between two linearly elastic half-spaces and the resulting displacements. For the mode-III (anti-plane shear) elastodynamic problem studied in the manuscript, the interface is initially uniformly pre-stressed by $\tau_{\rm d}$ and is set to slide at an extremely small steady velocity $v_{0}$, such that the shear tractions at the interface take the form

The second right-hand-side term represents the instantaneous response to changes in the sliding velocity, the so-called radiation damping term. As discussed in the manuscript, this term can be understood as the damping of interfacial energy due to elastic waves radiated into the infinite domain. The third term $s(x,t)$ accounts for the history and spatial distribution of interfacial displacements $u(x,t)$. Both $s(x,t)$ and $u(x,t)$ are related in the spectral domain via a convolution integral, whose expression can be found in [67]. Due to the spectral nature of the formulation, the simulated domain is taken to be periodic in the lateral direction, with periodicity $W$. The latter is chosen to be large enough to prevent any effect of the periodicity on the results reported in the manuscript.

Rupture is nucleated at the center of the domain by introducing a Gaussian perturbation of the slip velocity into an initial steady sliding state at $v_{0}$. The sliding velocity is then computed by combining Eq. (A1) and the rate and state-dependent friction law $\tau=\sigma\mathrm{sgn}(v)f(|v|,\phi)$ (see Appendix A2 for more details). $u(x,t)$ is then integrated in time using an explicit time-stepping scheme

Note that the factor $0.5$ on the right-hand-side of Eq. (A2) ensures that $v(x,t)$ is indeed the slip velocity. In order to guarantee the stability and the convergence of the numerical scheme, $\Delta t$ is defined as the time needed for a shear wave to travel a fraction $0.2$ of one grid spacing, i.e. $\Delta t\!=\!0.2\,\Delta x/c_{s}$. Additional information about the numerical scheme and the nucleation procedure can be found in [80], together with videos of similar rupture events.

The friction laws used in this work, and whose steady-state behaviors are plotted in Fig. 2, are related to the one used previously in [58, 80]. The friction law is defined by the relation between the shear stress $\tau\!\equiv\!\sigma_{xy}$ and the compressive normal stress $\sigma\!\equiv\!-\sigma_{yy}$ at the interface, $\tau\!=\!\sigma\sgn(v)f\left(\left|v\right|,\phi\right)$, and by the evolution equation for state variable $\phi$, $\dot{\phi}\!=\!g\left(\left|v\right|,\phi\right)$. The constitutive functions $f\left(\left|v\right|,\phi\right)$ and $g\left(\left|v\right|,\phi\right)$ used in this work take the form

where $\phi$ represents the typical age/maturity of contact asperities that compose the interface at a microscopic scale [45]. In Eq. (A3), $f_{0}$ sets the scale of the dimensionless frictional resistance (friction coefficient), $b$ is the aging coefficient and $a$ is related to the thermally-activated rheology of contact asperities [45]. The function $\sqrt{1+\left(v_{*}/v\right)^{2}}$ that appears also in $g\left(\left|v\right|,\phi\right)$ ensures that for vanishingly small steady-state velocities, $\phi$ saturates after extremely long times to a finite value of $D/v_{*}$, rather than diverges. As discussed in [80], this regularization makes no significant difference in the results discussed above. For the sake of notational simplicity, the regularization is hence omitted in the main text (though it is included in the calculations). Eqs. (A3) and  (A4) lead to the steady-state friction curve of Fig. 2 with a minimum at an intermediate $v$ (brown solid line), while the no-minimum steady-state friction curve (dash-dotted orange line) is obtained after neglecting the “+1” in the $b$ term.

The reader is referred to [58, 80] for additional discussions about the formulations of Eqs. (A3)-(A4), which go beyond the conventional rate-and-state friction laws. Nevertheless, the results and conclusions discussed above are independent of the choice of the rate-and-state formulation.

In this Appendix we describe propagating steady-state rupture fronts in thin (quasi-1D) systems, where no stress drops emerge. We consider two long and thin linear elastic bodies of height $H$ in frictional contact, such that the momentum balance equation $\rho\ddot{\bm{u}}\!=\!\nabla\!\cdot{\bm{\sigma}}$ reduces to [75, 84, 80]

where $u\!\equiv\!u_{x}$, $\bar{\mu}$ and $c_{\rm 1D}$ are the effective shear modulus and wave-speed [75, 84], respectively, and $\tau_{\rm d}$ is a constant driving stress (see Fig. 2).

Propagating 1D steady-state solutions then satisfy [80]

where we defined a co-moving coordinate $\xi\!\equiv\!x\!-\!\beta c_{\rm 1D}\,t$, integrated out $u$ and eliminated partial time-derivatives.

Steady-state rupture propagation is a dynamical process in which a homogeneous $V$ state invades a homogeneous $v_{0}\!\ll\!V$ state [75, 32, 35, 80], both shown in Fig. 2 as the intersections of the velocity strengthening branches of the friction law with the driving stress $\tau_{\rm d}$. We found these solutions for the friction law described in Sect. Appendix A2: The friction laws, using a shooting method [85] (similar to that used in [75, 84]). The solution is shown in Fig. A1. To normalize the stress fields we used the definition $\tau_{\rm m}$ is the maximal value $\tau$ attains in the profile, and $\ell\!\equiv\!\frac{1-\beta^{2}}{\beta}\frac{HV\bar{\mu}}{c_{\rm 1D}\left(\tau_{\rm m}-\tau_{\rm d}\right)}$ is the lengthscale over which the fields change, which can be calculated by a scaling analysis of Eq. (A6).

As seen in Fig. A1, the stress both ahead and behind the rupture front equal $\tau_{\rm d}$, i.e. there exists no stress drop. A corollary is that no singularity is observed in Fig. A1 (compare to Figs. 3a and 4).

The parameters used for all the calculations described in this work, except for the slow rupture to be discussed in Appendix A5, are given in Table 1.

Note that the values of the listed parameters are characteristic of some laboratory experiments (see [59] for details). However, the generic properties of the derived results are independent of the exact numbers, and are relevant to a broad range of materials and physical situations. For example, $v_{*}$ that controls the velocity scale below which the system is in the stick phase, can be taken to be significantly smaller.

The very same constitutive framework can give rise to slow rupture fronts ($c_{\rm r}\!\ll\!c_{s}$), as demonstrated in Fig. 4. While the emergence of slow rupture is of great interest in general, in the present context we are just interested in generating slow rupture and studying its properties in relation to the theoretical prediction in Eqs. (7)-(8). One way to generate slow rupture within our constitutive framework is to use friction parameters that shift the steady-state friction curve to smaller slip velocities and to employ a different nucleation procedure.

In particular, using the different set of parameters listed in Table 2, we obtain the dash-dotted orange steady-state curve in Fig. A2a (the solid brown line is identical to the one shown in Fig. 2a). For these parameters, and for the same value of the normalized driving stress $\tau_{\rm d}/\sigma$, the effective steady-state friction curve shown in Fig. A2b (dashed orange line, obtained by adding the radiation damping term $\frac{\mu}{2c_{s}}v$), is practically indistinguishable from the steady-state friction curve in the relevant slip velocities range. In addition, rupture is nucleated by introducing a perturbation to the internal state field $\phi$ of the form

with $k\!=\!2\pi/W$ and $\varepsilon\!=\!10^{-4}$, into an interface that slides homogeneously at a velocity $v_{\rm vw}$ that corresponds to a fixed-point on the velocity-weakening branch (it is marked by the black diamond in Fig. A2b). This nucleation procedure is different from the one used elsewhere in the paper, where Gaussian perturbations are introduced into an essentially locked-in interface, as described in detail in [80].

These modifications are sufficient to generate the slow rupture shown in Fig. 4. The physics behind the emergence of slow rupture, which is very interesting in itself, is not thoroughly discussed here. It deserves an investigation of its own, which we hope to pursue in the future.

The data reported in Fig. 5 are obtained from the experimental measurements [35] of ruptures propagating along a frictional interface formed by two plates of Homalite. Figure 8 of [35] reports the steady-state friction coefficient versus slip rate measured at the interface in the wake of the propagating rupture front. In terms of the notation used in this manuscript, the former is the ratio $\tau_{\rm res}/\sigma$, while the latter corresponds to $v_{\rm res}$. In the experiments of [35], the frictional interface is pre-cut at an angle $\alpha$ from the principal direction of the imposed compressive stress $P$, such that

Rubino et al. [35] distinguished data measured with respectively high and low levels of accuracy. From the data sets reported in Fig. 8 of [35] and their associated boundary conditions listed in Table 3, we compute the triplets $(\tau_{\rm d},\tau_{\rm res},v_{\rm res})$, which are then used to construct Fig. 5.

In Sect. IV, we also analyzed experimental data extracted from [36] in order to test the theoretical predictions in Eqs. (6)-(7). The bulk parameters $\rho\!=\!1170$ kg/m³ and $c_{s}\!=\!1345$ m/s, reported on in [36], have been used together with $\mu\!=\!\rho c_{s}^{2}$. The stress drops $\Delta\tau$ can be read off Fig. 2b of [36] (all figure indices in this and the next paragraphs refer to [36]), where the shear stress distribution near the edge of both slow (red) and rapid (here supershear, blue) rupture is presented. For rapid rupture, we extracted from Fig. 2b (blue data) $\Delta\tau\!=\!1.17$ MPa. The corresponding particle velocity distribution $\dot{u}_{x}(x,y\!=\!3.5\,\hbox{mm})$, measured $3.5$ mm above the interface, is presented in Fig. 1d (bottom). We used the leftmost value behind the edge, $\dot{u}_{x}(x,y\!=\!3.5\,\hbox{mm})\!\simeq\!0.75$ m/s as an estimate for the tail particle velocity at the interface $\dot{u}_{x}(y\!=\!0^{+})$, from which we estimate the residual slip velocity to be $v_{\rm res}\!=\!2\dot{u}_{x}(y\!=\!0^{+})\!\simeq\!1.5$ m/s. Note that the latter estimate is in very good agreement with the slip velocity reported for rapid (supershear) rupture in Fig. 3 (blue circles). Moreover, it is also in very good agreement with the normalized real area of contact $A_{r}/A_{0}\!\simeq\!0.6$, reported on in Fig. 1c, which according to Fig. 3b indeed corresponds to $v_{\rm res}$ slightly larger than $1$ m/s.

Using these estimates, we obtain

in great agreement with the theoretical prediction Eq. (5) in the manuscript. Repeating this procedure for slow rupture, we extracted from Fig. 2b (red data) $\Delta\tau\!=\!0.242$ MPa and from Fig. 1d (top) $v_{\rm res}\!=\!2\dot{u}_{x}(y\!=\!0^{+})\!\simeq\!2\times 0.01=0.02$ m/s. Using these estimates, we obtain

in great agreement with the theoretical prediction Eq. (7).

We provide here details about the estimates of the spatiotemporal integral contribution to the stress drop presented in Sect. IIB. As explained in Sect. IIA, the relative magnitude of the latter, denoted by $s$, is given by $\delta v_{\rm res}/v^{0}_{\rm res}\!=\!(v_{\rm res}-v^{0}_{\rm res})/v^{0}_{\rm res}$. Once $\Delta\tau$ is known, $v^{0}_{\rm res}$ is obtained through Eq. (5). Our goal is to show how $v_{\rm res}$ can be estimated using benchmark crack solutions [78]. Consider first anti-plane shear (mode-III symmetry) self-similar crack solutions, where the crack face displacement is given in Eq. (6.9.145) in [78] and takes the form

where the notation has been adapted to be consistent with the present manuscript and ${\bm{E}}(\cdot)$ is the complete elliptic integral of the second kind [86]. Note that the center of the crack, whose length at time $t$ is $2c_{\rm r}t$, is located at $x\!=\!0$.

In order to estimate $v_{\rm res}$ based on Eq. (A13), follow these steps: (i) replace the far-field applied shear stress $\tau_{\rm d}$ by the stress drop $\Delta\tau$ to account for a finite frictional resistance, not included in the crack solution (ii) take the limit $c_{\rm r}\!\to\!c_{s}$ ($c_{s}$ is the limiting propagation speed in mode-III) because only in this limit the frictional rupture we consider mimics self-similar propagation (iii) focus on a region far behind the propagating tip, i.e. set $x\!\to\!0$, which is relevant for the definition of a stress drop (iv) take the time derivative to obtain $\dot{u}_{z}(x\!=\!0,y\!=\!0^{+},t)$ and finally (v) set $v_{\rm res}\!=\!\dot{u}_{z}(x\!=\!0,y\!=\!0^{+},t)-\dot{u}_{z}(x\!=\!0,y\!=\!0^{-},t)\!=\!2\dot{u}_{z}(x\!=\!0,y\!=\!0^{+},t)$. The result reads

Using Eq. (5), which implies $v^{0}_{\rm res}\!=\!2c_{s}\Delta\tau/\mu$, together with Eq. (A14), we obtain (as reported on in Sect. IIB)

where ${\bm{E}}(0)\!=\!\pi/2$ has been used.

In order to repeat this procedure for in-plane shear (mode-II symmetry) self-similar crack solutions, we start with the crack face displacement given in Eq. (6.9.85) in [78]

for $|x|\!<\!c_{\rm r}t$, where the function $g_{2}(\cdot)$ is given in Eq. (6.9.87) in [78] and $c_{d}$ is the dilatational wave speed. Considering the limit $c_{\rm r}\!\to\!c_{R}$ ($c_{R}$ is the Rayleigh wave-speed, the limiting speed for mode-II) and following the procedure described above, we obtain

which leads to

The latter is a function of $c_{s}/c_{d}$ (note that $c_{R}/c_{s}$ is also determined by $c_{s}/c_{d}$). Evaluating $g_{2}(c_{R}/c_{s},c_{s}/c_{d})$ [88] for the whole range of admissible $c_{s}/c_{d}$ values (which in turn depends on Poisson’s ratio) and plugging in Eq. (A18) shows that $\delta{v}_{\rm res}/v_{\rm res}^{0}$ varies between $-0.613$ and $-0.417$, as reported on in Sect. IIB.