---
title: Off-Fault Damage Zone in Rupture Mechanics
url: https://www.emergentmind.com/topics/off-fault-damage-zone
type: topic
---

# Off-Fault Damage Zone in Rupture Mechanics

Searching arXiv for recent and foundational papers on off-fault damage zones in dynamic rupture, laboratory characterization, and related modeling.
An off-fault damage zone is the region of rock surrounding a principal fault trace or slipping fault core that undergoes distributed brittle deformation or inelastic deformation during rupture, rather than slip being confined to the principal surface itself. Across laboratory, field, seismological, and numerical studies, it includes microfractures and secondary cracks in wall rock, distributed strain around complex fault junctions, and, in some fault zones, multiscale fracture networks that can themselves host dynamic rupture. Its significance follows from energy partitioning: part of the elastic strain energy released by shear failure is dissipated as off-fault damage, altering rupture dynamics, slip-weakening behavior, seismic radiation, and near-fault elastic and hydraulic properties [2007.09687] [1804.03705] [2307.14229].

## 1. Definition, geometry, and structural expression

The off-fault damage zone is distinct from the principal slip surface and from the much narrower on-fault process zone. In 2D plane-strain rupture mechanics, on-fault slip is accommodated within a process zone of size  
$$
R_0=\frac{9\pi}{32(1-\nu)}\frac{D_c\mu}{(\mu_s-\mu_d)(-\sigma_{yy})},
$$
whereas the damage zone width $W_d$ is defined as the distance from the fault trace to the locus where a chosen damage threshold is first exceeded [2212.12709]. Outside $W_d$ the rock remains essentially elastic; inside $W_d$ irreversible deformation has consumed part of the elastic strain energy radiated by the rupture front [2212.12709].

Its geometric expression depends on scale and setting. In continental surface rupture, off-fault damage zones can span hundreds to thousands of meters, as in the Kaikoura triple-junction area, where the Papatea fault shows a distributed damage zone width of $\sim2$ km, in contrast to tens of meters on the Jordan and Kekerengu faults [1804.03705]. In laboratory granite experiments, the same phenomenon appears as a millimetre-scale zone of microfracturing: after quasi-static failure the damage zone width is around $10$ mm, and after dynamic failure around $20$ mm [2007.09687]. At seismogenic depth, aftershock clustering reveals finite-thickness planar patches with mean thickness $T=4\sigma_3\approx290$ m, suggesting that deformation is volumetric rather than confined to an infinitely thin core [1006.0885].

In explicit multiscale fault-zone models, the damage zone is volumetric and networked rather than a uniform halo. One 3D representation uses a $10$ km (strike) $\times1$ km (normal) $\times6$ km (depth) region around an $8$ km $\times4$ km listric main fault, containing $854$ explicit fractures with radii $R$ between $100$ m and $500$ m, following a power-law size distribution $\propto R^{-2}$ [2412.15416]. In the related multiscale study, two conjugate fracture families strike on average N120°E and N20°E, each with a $\pm10^\circ$ Fisher-distributed scatter in strike and dip; fracture lengths follow $N(R)\propto R^{-2}$ between $R_{\min}\approx100$ m and $R_{\max}\approx500$ m; and fracture density decays away from the fault core as
$$
P_{10}(d)=C\,d^{-m},
$$
with $C\approx2.5$ and $m\approx0.8$ [2307.14229]. About $78\%$ of fractures connect to more than one neighbor, a structural property that directly conditions rupture transfer [2307.14229].

Near-surface field characterization also shows that “damage zone” is not internally homogeneous. In the Vado di Corno Fault Zone, a high-strain damage zone forms a very slow layer with $V_p\le1.0$–$1.5$ km/s and thickness $\sim10$–$15$ m in shallow seismic tomography, while field mapping places high-strain damage-zone lenses $\sim50$–$160$ m outward from the principal fault surface [2601.04717]. This suggests that the off-fault damage zone is better regarded as a mechanically heterogeneous envelope composed of units with different fracture intensity, sealing state, and compliance, rather than as a single material class.

## 2. Energetics and constitutive mechanics

In rupture mechanics, the strain energy released by shear failure is partitioned into radiated seismic waves, frictional work on the slip surface, and breakdown work. Breakdown work contains the fracture energy $\Gamma$ dissipated in the tip process zone, of which a fraction $\Gamma_{\mathrm{off}}$ is spent in creating off-fault damage [2007.09687]. A corresponding dynamic balance used in multiscale rupture models is
$$
G_{\mathrm{release}} \ge G_c(L)+G_{\mathrm{radiated}}+G_{\mathrm{off\text{-}fault}},
$$
where $G_{\mathrm{release}}$ is the elastic energy release rate [2307.14229].

A recurrent result is that off-fault damage is an active energy sink rather than a passive elastic buffer. Laboratory measurements in Lanhélan granite give $\Gamma_{\mathrm{off}}\simeq3$ kJ/m$^2$ for quasi-static rupture and $\simeq5.5$ kJ/m$^2$ for dynamic rupture [2007.09687]. In the quasi-static run, total fracture energy on the principal slip zone was $\Gamma\approx27$ kJ/m$^2$, so that $\Gamma_{\mathrm{off}}/\Gamma\approx10\%$; in the dynamic runs, the additional $\sim2.5$ kJ/m$^2$ of off-fault energy implies a $\sim10\%$ increase in total $\Gamma$ with rupture speed [2007.09687]. In 2D FDEM simulations, the energy fraction associated with off-fault fracture energy rises with rupture length and depth: at $20R_0$, $E_G^{\rm off}\approx10\%$ of total released energy at $2$ km, $\approx15\%$ at $6$ km, and $\approx20\%$ at $10$ km [1901.01771].

Several constitutive descriptions are used to represent this energy sink. In multiscale fracture-network rupture, both main fault and fractures obey a laboratory-based rate-and-state friction law with rapid velocity weakening, with effective normal stress
$$
\sigma_n'=\sigma_n-P_f(z),
$$
and relative pre-stress ratio
$$
\mathcal{R}=\frac{\tau_0-f_w\sigma_n'}{(f_p-f_w)\sigma_n'}.
$$
An optimally oriented fracture has $\mathcal{R}\to1$; sustained slip generally requires $\mathcal{R}\gtrsim0.3$, while truly favorable fractures have $\mathcal{R}\gtrsim0.6$ [2307.14229]. In those simulations, the state-evolution distance is tied to fracture radius by $L\approx2\times10^{-5}R$, corresponding to $G_c(R)\approx400\,\mathrm{Pa}\times R$ [2307.14229].

Other frameworks distribute damage in the bulk. Micromechanical models define damage by growth of tensile wing cracks from pre-existing penny-shaped cracks, with
$$
D=\frac{4\pi}{3}N_v\,(a\cos\Phi+l)^3,
$$
and evolving elastic moduli $\mu(D)$ and $\lambda(D)$ [2406.18408]. Continuum damage–breakage models use a scalar damage $\alpha$ and breakage variable $B$, causing reductions in shear modulus $\mu$ and bulk modulus $\lambda$ and producing either smooth distributed damage or localized shear bands once brittle failure is reached [2503.21260]. Phase-field models instead use a continuous damage variable $d(\mathbf{x})\in[0,1]$ and a regularization length $L$, so that fault propagation and off-fault damage follow from a single energy functional [2206.07933] [2503.06576].

## 3. Observation and quantification

Off-fault damage zones are quantified by different observables in different regimes, but several studies converge on measurable width, density, and property-reduction metrics.

In laboratory triaxial experiments on intact Lanhélan granite at $100$ MPa confining pressure, time-resolved 3D $P$-wave velocity tomography and post-mortem microstructures quantify off-fault damage during and after failure [2007.09687]. Tomography reveals a narrow low-velocity band around the slip surface with up to $25\%$ drop in $V_P$, and SEM fracture maps show fracture densities decaying to background over $\sim10$ mm after quasi-static rupture and $\sim20$ mm after dynamic rupture [2007.09687]. The conversion from $V_P(\theta)$ to crack density uses an effective-medium model for dry rock with penny-shaped cracks, fitting $V_{PH}$ and $V_{PV}$ in each voxel to recover the crack-density tensor components $\alpha_{11}$ and $\alpha_{33}$ [2007.09687].

Field-scale geodetic mapping resolves much broader off-fault deformation. In the Kaikoura earthquake, optical image correlation combines SPOT6 and Pleiades imagery, with a high-resolution correlation at $1.8$ m pixel and detection threshold $\sim20$ cm [1804.03705]. Swath profiles show maximum right-lateral slip of $\sim11$ m on the Kekerengu fault and total left-lateral slip of $\sim6$ m distributed over $\sim2$ km width on the Papatea fault [1804.03705]. There, damage zone width is defined where the displacement gradient is diffuse rather than localized [1804.03705].

The supershear-transition literature uses the off-fault damage zone itself as a diagnostic. For the 2001 Kunlun earthquake, optical image correlation yields a mean fault-zone width of $238\pm80$ m for most of the rupture, but over an $\sim11$ km stretch the width collapses to $127\pm39$ m, coincident with the inferred sub-to-supershear transition [1909.12672]. Early aftershock catalogs for the Izmit, Denali, and Craig earthquakes show corresponding along-strike gaps in cumulative seismic-moment density over a few-kilometer segment at the transition location [1909.12672].

Seismological and near-surface structural methods add further constraints. A Gaussian-mixture clustering analysis of the Mount Lewis aftershocks reconstructs anisotropic planar patches with mean length $L\approx2.3$ km, mean width $W\approx0.9$ km, and mean thickness $T\approx290$ m, a thickness an order of magnitude larger than relative location errors [1006.0885]. In the Vado di Corno Fault Zone, ultrasonic measurements at $\sim1$ MHz and field tomography at $4.5$ Hz together show that the fault core is slower than some damage-zone units at laboratory scale, yet tomography resolves a very slow intensely fractured high-strain damage zone at shallow depth [2601.04717].

| Context | Width or thickness | Diagnostic |
|---|---:|---|
| Laboratory granite rupture | $\sim10$ mm quasi-static; $\sim20$ mm dynamic | Microfracture density and $V_P$ drop [2007.09687] |
| Kaikoura Papatea fault | $\sim2$ km | Diffuse displacement gradient [1804.03705] |
| Kunlun supershear study | $238\pm80$ m to $127\pm39$ m | Optical image correlation [1909.12672] |
| Mount Lewis aftershocks | $\approx290$ m thickness | Gaussian-kernel seismic clustering [1006.0885] |
| Vado di Corno HSDZ | $\sim10$–$15$ m tomography; $\sim50$–$160$ m field lenses | Seismic tomography and field mapping [2601.04717] |

These measurements are not interchangeable. They sample different materials, frequencies, depths, and definitions of damage. A plausible implication is that “width” is operational rather than unique: it depends on whether the threshold is based on crack density, velocity reduction, inelastic strain, scalar damage, or diffuse displacement.

## 4. Rupture interaction, scaling, and dynamic regimes

Off-fault damage zones interact dynamically with rupture through branching, jumping, stress transfer, and energy competition. In multiscale fracture networks, rupture on one fracture imposes dynamic shear stress changes $\Delta\tau$ and normal stress changes $\Delta\sigma_n$ on neighbors; if $\Delta\tau$ raises the local shear above the peak strength while $\Delta\sigma_n$ unclamps neighboring fractures, neighboring fractures may nucleate [2307.14229]. Two transfer modes are distinguished: branching at direct intersections, and jumping by dynamic overstress at distance $\le80$ m [2307.14229].

Successful cascading rupture within a fracture-network damage zone requires a specific combination of structure, prestress, and frictional scaling. In the 3D listric-fault models, cascades arise if at least one fracture family is favorably oriented with $\mathcal{R}\gtrsim0.6$ under ambient $SH_{\max}$ azimuth $\Psi\approx50^\circ$–$80^\circ$, while the conjugate set has $0.3\lesssim\mathcal{R}\lesssim0.6$; fracture connectivity and spacing $<80$ m ensure rupture transfer; and $L$ scaling with fracture size provides sufficient energy for self-sustained rupture [2307.14229]. Under these conditions, cascades involving hundreds of fractures can release moment magnitudes up to $M_w\approx5.6$ without engaging the main fault, and cascading rupture within the fracture network discourages rupture on the main fault [2307.14229]. In “pure cascade” cases near $\Psi\approx65^\circ$, the main fault remains dynamically quiescent except for small self-arresting patches, whereas at $\Psi\approx100^\circ$–$120^\circ$ runaway rupture occurs on the main fault with limited off-fault fracture slip [2307.14229].

Rupture speed leaves a direct signature in damage-zone width. Laboratory experiments show that for slips up to a few mm, the spatial extent of off-fault peak stresses, and hence damage width, is controlled by rupture speed $v$, not by cumulative slip; the observed doubling from $\sim10$ to $\sim20$ mm implies $v\simeq0.8\,c_S$ in the dynamic runs [2007.09687]. Fracture-mechanics analysis of supershear transition predicts that damage-zone half-width scales approximately as
$$
w(v)\sim\left[\frac{K_{II}^{\rm dyn}(v)}{\sigma_y}\right]^2,
$$
so that width first broadens and then collapses as $v\to c_R$ because the dynamic stress intensity factor contracts [1909.12672]. FDEM simulations of planar supershear rupture show a corresponding “damage gap” around the transition point, with minimal off-fault failure during the jump [1911.03468].

Depth dependence is robust, but model-specific. In one micromechanical formulation,
$$
w(z)\approx c\,R_0(z),
\qquad
c\approx
\begin{cases}
3.4 & D_{\rm thres}=0.2,\\
2.3 & D_{\rm thres}=0.9,
\end{cases}
$$
with $R_0(z)\propto[-\sigma_{yy}^0(z)]^{-2}$, implying narrower yet denser damage at depth [2406.18408]. In 3D continuum damage–breakage simulations, the $1\%$ $\Delta\mu$ contour extends $\sim2.5$ km from the fault at $2.5$ km depth but only $\sim1.3$ km at $7.5$ km depth, summarized by a simple scaling argument $W\propto d^{-1}$ [2503.21260]. This suggests that narrowing with depth is consistent across models, while the precise exponent is formulation-dependent.

## 5. Numerical representations and methodological issues

A large part of the literature concerns how to represent off-fault damage numerically without collapsing it into an imposed plastic halo or ignoring it entirely.

Combined finite–discrete element methods represent off-fault fractures explicitly on element interfaces. In HOSSedu, each potential interface can undergo tensile or shear cohesion softening and slip-weakening friction, allowing spontaneous off-fault crack networks to emerge from local stress concentrations [1911.03468]. This framework has been used for planar faults, kinks, step-overs, and rough faults; with off-fault damage included, a rupture can arrest at a compressional kink, jump a dilational step-over through relay-zone fracturing, or arrest early on a rough fault because cracking amplifies stress heterogeneity [1911.03468]. The Kaikoura rupture reconstruction uses the same 2D FDEM approach to compare alternative rupture paths through the Jordan–Kekerengu–Papatea triple junction, with the preferred scenario reproducing both narrow localization on Jordan/Kekerengu and broad distributed deformation on Papatea [1804.03705].

Micromechanical continuum models resolve damage through internal variables rather than discrete cracks. In one formulation, wing-crack growth updates a scalar damage field $D$ and therefore the homogenized elastic stiffness tensor at each time step, permitting three local regimes: purely elastic, sliding along microcracks, and tensile opening of penny and wing cracks [2406.18408]. A distinct 3D continuum damage–breakage model uses damage $\alpha$ and breakage $B$ to distinguish smooth distributed damage from localized, mesh-independent shear bands, with observed band orientations in agreement with a linearized bifurcation analysis [2503.21260].

Phase-field models unify on-fault rupture and off-fault damage with a single diffuse variable. In the quasi-dynamic rate-and-state formulation, off-fault damage is the region where $0<d(\mathbf{x})<1$, and a practical thickness measure is obtained by choosing a threshold $d_{\rm th}$ and defining the extent normal to the main rupture where $d\ge d_{\rm th}$ [2206.07933]. In the in-plane extension, damage zone thickness $w_d$ is naturally of order the phase-field width $O(L)$, while the model also incorporates pressure-dependent shear strength and a revised propagation criterion that accounts explicitly for the coupling between shear strength and normal stress [2503.06576].

Diffuse-interface hyperbolic formulations take a broader continuum view. The GPR model introduces a damage concentration $\xi\in[0,1]$ and a solid volume fraction $\alpha\in[0,1]$, allowing secondary cracks and complex fault geometries to be represented without mesh alignment or remeshing [2007.01026]. In the dynamic rupture examples, the diffuse fault core is imposed with width $W_{fc}=100$ m and spontaneous off-fault shear cracks develop in the extensional quadrants, with a typical off-fault damage zone width $\sim500$–$1000$ m for $\xi>0.1$ [2007.01026].

A central methodological issue is stress initialization. For 2D strike-slip rupture under plane strain, neglecting the out-of-plane stress or choosing $\sigma_{zz}$ arbitrarily can move the model into a forbidden reverse-fault regime and lead to a six-fold underestimation of damage-zone width: in the studied cases, correct 3D stress ordering gives $W_d\approx1200$ m, whereas the incorrect setup gives $W_d\approx200$ m [2212.12709]. This is not a minor technicality, because the same factor-of-six appears whether damage is diagnosed by $\Delta CFF>0.3\Delta\tau$, plastic strain $\gamma_t>10^{-4}$, or scalar damage $D>0.3$ [2212.12709].

## 6. Seismic radiation, hydrology, and hazard implications

Off-fault damage zones affect not only rupture mechanics but also source spectra, near-field shaking, groundwater response, and rupture reconstruction.

Explicit 3D fracture-network simulations show that cascading ruptures in the damage zone generate distinct ground-motion signatures with high-frequency content, resembling source coda-wave-like signatures [2412.15416]. In the modeled cascading case, near-field corner frequency is $f_c\simeq1.1$ Hz with standard deviation $\simeq0.5$ Hz, versus $\lesssim0.5$ Hz with standard deviation $\lesssim0.1$ Hz in non-cascading cases; the seismologically inferred stress drop is $\sim75$ MPa, an order of magnitude larger than the directly computed average off-fault $\Delta\sigma\simeq9.6$ MPa [2412.15416]. Peak ground acceleration and peak ground velocity are roughly twice those of the non-cascading cases despite the smaller moment magnitude, and pseudo-spectral accelerations at short periods $T=0.1$–$0.2$ s are elevated [2412.15416].

Other dynamic-rupture frameworks yield the same qualitative conclusion by different mechanisms. FDEM simulations show enhanced high-frequency radiation in the near field due to coseismic off-fault damage, with $>20$ Hz signals arriving after the main rupture and a reduction of seismic efficiency by up to $\sim20\%$ at $10$ km depth when damage is included [1901.01771]. Continuum damage–breakage simulations show that rapidly evolving bulk moduli generate broadband energy above $1$ Hz and nearly isotropic fault-normal to fault-parallel amplitude ratios in the $1$–$4$ Hz band, while also reducing rupture speed and increasing supershear transition distance [2503.21260].

The hydraulic consequences are likewise non-negligible. A 3D groundwater model coupled to dynamic rupture with inelastic off-fault deformation shows that inelastic dilation causes notable depressurization within $1$ to $2$ km off the fault at shallow depths $<3$ km, with modeled coseismic pore-pressure drops of approximately $-0.1$ to $-0.2$ MPa within $\pm1$ km of the fault [2512.23912]. Near-fault water-level drawdowns reach $20$–$90$ m at $x=\pm30$ m, $30$–$40$ m at $1$ km, and $5$–$10$ m at $2$ km, whereas elastic-only poroelasticity gives only a few meters and opposite-sign changes across the fault [2512.23912]. This indicates that off-fault damage can act as a transient pressure sink during the coseismic phase and accelerate postseismic recharge [2512.23912].

In porous-media flow modeling, the damage zone can also be represented as a distinct hydrogeologic layer rather than collapsed into the core. A mixed-dimensional Darcy model treats the damage zone as a lower-dimensional manifold $\mu$ of thickness $\epsilon_\mu$, with its own tangential permeability and Robin-type coupling to both rock matrix and fault core [1903.01117]. This captures pressure-drop localization at damage-zone interfaces, as well as flow focusing or diversion around low-permeability patches [1903.01117].

For seismic hazard and rupture interpretation, the principal implication is that off-fault damage is diagnostically informative and dynamically consequential. Off-fault deformation “fingerprints” can distinguish competing rupture paths in complex systems such as Kaikoura [1804.03705]. Highly connected, favorably pre-stressed fracture networks can shadow or inhibit main-fault rupture while releasing off-fault seismic moment up to $M_w\approx5.6$ [2307.14229]. Standard single-plane source models and conventional ground-motion prediction equations may therefore underpredict near-fault PGA and high-frequency content when multiscale fault-zone complexity is ignored [2412.15416].

Source: https://www.emergentmind.com/topics/off-fault-damage-zone