Codex @codex  0

An empirical relationship summarizing the high friction measured for many rock surfaces as a function of effective normal load. Weak fault gouge or high pore pressure can produce different effective fault strengths. Named for James Byerlee.

Subduction zone segmentation is the observation that slab structure, fault slip and upper-plate behavior change along a margin; their boundaries need not coincide exactly. In Cascadia, it is documented by several independent measurements.

Onshore GNSS velocities, coastal leveling and tide gauges constrain along-strike changes in megathrust locking and interseismic uplift. Episodic tremor and slip has distinct recurrence and migration patterns, and its occurrence is lower in parts of central Oregon than farther north and south. Coastal geological records constrain differences in coseismic subsidence and the extents of past ruptures. Those observations sample different depths and time windows, so a present-day locking boundary should not automatically be declared a permanent earthquake rupture barrier. A 2014 analysis found persistent central-margin partial creep compatible with reduced uplift, reduced paleoseismic subsidence and lower cumulative slow slip.^[ref]

A physically specific hypothesis links these variations to subslab mantle buoyancy. An onshore–offshore seismic tomography study imaged low-velocity asthenosphere beneath northern and southern Cascadia, spatially associated with stronger plate locking and more tremor, with different central-margin structure. Warmer or partly molten asthenosphere can be less dense. Its buoyant traction can push upward against the slab and change contact stresses on the megathrust. Hotspot-influenced material is one proposed northern source; three-dimensional deformation near the Gorda region offers a southern upwelling mechanism.^[ref]

There are additional controls. The Siletzia basaltic basement beneath the overriding plate varies in extent and thickness. Its permeability and stiffness can affect fluid escape and fault stresses: poorly drained material can preserve high pore pressure, lower effective normal stress and promote creep. Along-strike changes in incoming plate age, slab dip, fracture-zone hydration, sediment supply and plate-edge stresses can also modify the thermal and frictional environment. Central-margin creep has specifically been associated with the distribution of Siletzia and inferred fluid trapping.^[ref]

How do we distinguish these causes? Offshore and onshore seismic arrays test the geometry of the anomalies; resolution and synthetic-recovery tests determine whether apparent segmentation could result from uneven station coverage. Seismic attenuation, anisotropy, magnetotelluric sounding, rock physics and heat-flow data help separate temperature, melt, water and composition. Gravity anomalies and seismic interfaces constrain basement and slab geometry; geodetic inversions test whether the inferred locking pattern persists across plausible fault geometries. Forward models then ask whether the measured density or pore-pressure contrasts produce stresses and slip patterns of the required size and sign.

Segmentation is observed robustly; a unique cause is not established by spatial correlation alone. Seismic velocity is not a direct density measurement, central creep and deep tremor are not the same phenomenon, and multiple coupled structures can reinforce one another. The buoyant-mantle and upper-plate fluid hypotheses are testable explanations rather than mutually exclusive certainties.

Use a rock strength envelope: at each depth plot the smaller differential stress needed for frictional failure and for sustained ductile material creep at a specified strain rate. Shallow rock generally strengthens with depth because effective confining pressure increases; sufficiently deep hot rock weakens because thermally activated material creep becomes efficient.

At shallow depth, Byerlee's law or a Mohr-Coulomb failure criterion gives a relation of the formwhere $C$ is cohesion, $\mu$ a coefficient of friction, ${\sigma }_n$ the compressive normal stress and $P_f$ the pore pressure. Increasing overburden generally increases effective stress, hence the frictional strength. Conversion from fault shear stress to differential stress depends on the stress regime and fault orientation; a sketch should not label a universal numerical slope. High pore pressure or weak fault gouge can lower the frictional strength substantially.

A simplified dislocation creep relation isHere $Q$ is activation energy, $V^{\ast }$ activation volume, $R$ the molar gas constant and $T$ absolute temperature; $A$ and $n$ depend on composition and deformation mechanism. At fixed strain rate, increasing temperature usually decreases the required stress strongly. The competing effects of pressure, water weakening of rocks, grain size and composition also matter. The intersection with frictional failure gives a brittle-ductile transition, not a change from solid to liquid.

For the Indian continental interior, a useful first sketch has an upper-crustal strength peak, a weaker lower continental crust, and renewed strength in lithospheric mantle beneath the Moho. Quartz-bearing rocks can undergo ductile material creep at conditions where colder, drier olivine-rich mantle retains appreciable strength. This creates a possible two-peak, or “jelly sandwich”, rock strength envelope. With further depth the mantle also weakens. A cold craton can retain strong lithospheric mantle to much greater depths than a tectonically heated continental region.

For the central Pacific, show thin mafic oceanic crust above a relatively thick layer of strong lithospheric mantle. There need not be a conspicuous weak lower-crustal layer. Much of the integrated strength can reside below the Moho, with a broad strength peak followed by ductile weakening towards the asthenosphere. Cooling and age control the thickness: young hot oceanic lithosphere would give a much shallower and weaker profile than an old central-Pacific plate.

The plotted Moho depths of about 38 km and 7 km, and all plotted stresses, are illustrative assumptions. “India” alone does not specify a unique geothermal gradient, water content or strain rate, and the Pacific includes plates of different ages. The robust comparison is the possible crust–mantle layering in a continent versus the thin oceanic crust and mantle-dominated strength of an old oceanic plate. The Moho is a compositional boundary; neither it nor the deepest small plotted strength is automatically the seismic lithosphere-asthenosphere boundary.

Left: ordinary frictional strengthening with depth followed by thermally activated material creep. Right: a very weak shallow faulted layer above a narrow strong region near 10 km, followed by essentially the same weakening material creep branch. These are the shapes actually plotted in Figure 2. The right curve is not a strong upper crust over a uniformly weak lower crust; neither diagram displays a second mantle-strength peak.

For the left curve, increasing overburden increases the effective normal load on faults. A schematic Coulomb friction law iswhere $C$ is cohesion, $\mu$ the coefficient of friction, ${\sigma }_n$ the compressive normal stress and $P_f$ the pore pressure. With ordinary rock friction and approximately hydrostatic pressure, ${\sigma }_n-P_f$ increases roughly linearly with depth. The stress needed for frictional sliding therefore increases, as in the straight shallow branch. Byerlee's law summarizes the relatively high friction found in many laboratory rocks.^[ref]

At higher temperature, dislocation creep or another ductile mechanism accommodates the imposed deformation at lower stress. For example,Here $\dot{ϵ}$ is strain rate, $Q$ an activation energy, $R$ the molar gas constant and $T$ absolute temperature. At fixed strain rate, increasing $T$ lowers the required stress. The operative envelope is approximately $\min ({\sigma }_f,{\sigma }_d)$, expressed in a consistent stress convention; the peak marks the brittle-ductile transition. This accounts for the left curve's peak near 100 MPa and its decline toward small strength by about 20 km.

For the right curve, normal high-friction sliding with depth-independent fluid-pressure ratio cannot explain the nearly vertical weak shallow branch. Plausible weakening includes high pore pressure tracking the overburden, low-friction fault gouge, or a connected network of pre-existing weak faults. Roughly constant small ${\sigma }_n-P_f$ yields roughly constant small frictional strength. The abrupt increase near the base of this layer requires a change of conditions: for example, the weak fluid-rich fault network terminates or becomes discontinuous, leaving stronger intact or less altered rock below. A change in mineralogy or deformation localization can also give a contrast. Simply saying “there are fluids” does not explain the return to high strength. Below the strong region, thermal material creep again lowers the strength.

The figure alone does not identify the weakening mechanism or a particular lithological boundary. In particular, a 10 km strength maximum is not evidence for a 10 km deep Mohorovičić discontinuity. The strong-fault/weak-fault evidence and the distinction between fault strength and surrounding crustal strength are discussed in Stress orientations at intermediate angles to the San Andreas Fault, California by Jeanne L. Hardebeck and Andrew J. Michael.

Support for the left curve comes from high laboratory friction and substantial in situ differential stress; support for the right curve comes from low inferred fault traction, weak fault materials or high pore pressure. No single observation directly measures the whole plotted envelope.

For the left, stronger shallow-crust model:

For the right, weak shallow-layer model:

To discriminate, measure the strength-controlling quantities through depth and across both faults and their surroundings. Obtain borehole stress, pore pressure and temperature profiles where feasible; use borehole breakouts and hydraulic tests to constrain effective stress. Recover fault and host rocks and test their friction and material creep at relevant fluid conditions. The right curve specifically predicts an end to the shallow weakening near its sharp rise: investigate changes in permeability of a porous medium, alteration and fault gouge with depth, rather than testing only a surface sample.

Combine those measurements with heat flux transects, magnetotelluric sounding for conductive fluids, accurate earthquake depths and focal mechanisms, and GNSS or InSAR observations of deformation. The latter can help constrain the depth distribution of locking and postseismic relaxation, but neither locking nor seismicity alone gives absolute strength. Earthquake stress drop is a change in stress, not the pre-event stress. Low heat flux can also reflect heat transport or transient weakening; high electrical conductivity can have causes other than free water. A joint explanation of all these data is needed to choose between the models locally.

Probably yes in connected fractures or transiently opened pore space, especially where a transform fault admits circulating seawater; a uniform water-filled porous layer several kilometres thick is not guaranteed. Porosity and permeability of a porous medium generally decrease as cracks close or seal, while ongoing faulting can reopen pathways. Fluid can occur at high pressure and temperature, including a supercritical state; it need not resemble cold surface water.

Several complementary tests address presence and origin:

No one anomaly uniquely identifies water. Seismic velocity can also reflect cracks, composition and temperature; electrical conductivity can reflect melt or conductive minerals. Serpentinization demonstrates incorporation of water into minerals, which must be distinguished from fluid currently occupying pore space. At an oceanic transform fault, an external seawater supply is the natural first hypothesis, but isotope ratios and chemical budgets can test contributions from deep fluids rather than presuming the source. A research example of coupled deformation and hydrothermal alteration is Mass transfer and fluid flow during detachment faulting and development of an oceanic core complex, Atlantis Massif by Chiara Boschi and others; a detachment system illustrates the methods without proving identical flow at every transform fault.

At a convergent margin, the incoming oceanic lithosphere is a chemically and mechanically layered package: sediments, oceanic crust and lithospheric mantle. Subduction does not deliver this whole package unchanged to one destination. Solids can be accreted, underplated, exhumed or recycled deeply, while released fluids and melts carry selected components to the overriding plate. The outcome depends on incoming composition and hydration, slab age, convergence rate, thermal structure, interface coupling and subsequent tectonic history.

The incoming material already contains a history. Sediments include terrigenous debris, biogenic carbonates, siliceous deposits, clay minerals, pore water and organic matter. Oceanic crust consists mainly of basalt, sheeted dikes and gabbro; hydrothermal alteration can add water, carbonates and mobile elements before subduction. Beneath it, residual peridotite forms the lithospheric mantle. Faulting at a spreading environment, a transform boundary or an outer-rise bending region can let seawater enter and cause serpentinization. The lithospheric mantle is not uniformly hydrated to its full thickness. Consequently, sediment thickness alone does not measure the incoming volatile inventory.

Sediments: separation starts near the trench. Compaction and increasing effective pressure expel pore fluids, which can escape through the forearc, transport dissolved constituents and influence pore pressure on the plate interface. Some sediments are scraped off into an accretionary prism, shortened by thrusts and retained in the upper plate. Others pass beneath the leading edge, then detach and undergo tectonic underplating beneath the forearc. Deeper underplated material may return to the surface during tectonic exhumation, preserving a record of burial and deformation. At a sediment-starved erosive margin, little prism may form, and tectonic erosion can entrain material from the overriding plate as well. Thus net continental gain and net continental loss are both possible; they can even alternate at one margin.

Sediments that remain attached to the slab undergo metamorphism, dehydration, and potentially decarbonation or carbonate dissolution. The released fluids can transport boron, barium, strontium, lead and other soluble constituents. Carbon can leave as dissolved species, remain in carbonate, or occur as reduced carbon depending on redox conditions and reaction history. Along sufficiently hot paths, sediment-derived melts or mixtures of sediment and altered rock can transfer material into the mantle wedge; buoyant mixed material may rise as a mélange diapir. Cold paths can instead retain more solid material and some volatile-bearing phases to greater depths. These are alternative reaction and transport pathways, not a single compulsory sequence. A carbon-budget study explicitly includes fluid dissolution, sediment diapirism and storage in the upper plate, showing why volcanic carbon dioxide alone cannot close the whole subduction carbon budget.^[ref]

Oceanic crust: metamorphic transformation and chemical filtration. Altered basalt and gabbro are progressively buried. Depending on their pressure-temperature paths and reaction kinetics, they can pass through blueschist facies or warmer assemblages and eventually form eclogite, dominated by garnet and omphacite. Eclogitization increases density and can strengthen negative slab buoyancy, but a kinetic delay can leave metastable rock. It is incorrect to assign every slab the same sequence or the same fixed depths for these changes.

During burial, breakdown of hydrous phases releases water-rich fluids; different mineral assemblages retain different quantities of water. The fluids migrate through or along the interface and into the mantle wedge, where they cause metasomatism and lower the solidus of peridotite. Flux melting of the mantle wedge generates much of the primary magma feeding a volcanic arc. This is different from saying that the entire descending oceanic crust melts beneath the arc. Most of its refractory mass can remain solid. Partial melting of oceanic crust is more plausible for unusually hot or young slabs, or other favorable thermal conditions, and produces a distinct crustal contribution.

The depth of dehydration is controlled by thermal structure and composition. Global thermal-petrological calculations, for example, predict more efficient deep water transport in old, fast, cold slabs than in hot ones, with substantial differences between volcanic crust, gabbro and hydrated peridotite. Such a calculation supports conditional deep retention; its modeled proportions should not be imposed on every margin.^[ref]

Mechanical detachment adds more branches. Upper-crustal slices, seamounts or oceanic plateaus may be accreted or underplated instead of continuing down the slab. A complete or partial ophiolite can be emplaced onto another plate by obduction in suitable convergence histories. Most intact oceanic crust is not saved this way. High-pressure fragments may be entrained in an interface channel or mélange and exhumed; their small preserved volume should not be mistaken for the dominant global mass flux. Magma generated after slab-derived transfer may also pond and crystallize beneath the upper crust by magmatic underplating, which is a different process from mechanical tectonic underplating of solid slices.

Lithospheric mantle: residual rock and a variable water carrier. Most incoming lithospheric mantle is depleted peridotite left after extraction of ridge melts. Its cold, dense, mechanically coherent portion contributes to slab pull and generally sinks with the plate. Hydrated portions contain serpentine and possibly other hydrous minerals. Breakdown of these phases releases water at conditions that depend on the thermal path, feeding further metasomatism of the mantle wedge and sometimes dehydration embrittlement. The latter is a proposed contributor to intermediate-depth earthquakes, not an exclusive explanation for all such events. Very cold paths may transfer some water onward in high-pressure hydrous phases or as hydrogen incorporated in nominally anhydrous minerals.

Do not confuse the incoming hydrated lithospheric mantle with newly serpentinized forearc mantle in the overriding plate. The latter receives slab-derived fluids and can temporarily store their water. Hydrated, relatively buoyant interface material can facilitate return flow and tectonic exhumation, while the colder deeper bulk of the lithospheric mantle continues down. Some peridotite may also be accreted or obducted with an ophiolite. These exchanges mean that even an exposed serpentinite requires structural and chemical context before it can be assigned to its original plate.

Deep recycling and eventual return. The residual eclogite, remaining sediments and lithospheric mantle can enter the deeper mantle. A slab can temporarily stagnate around the mantle transition zone or penetrate into the lower mantle, depending on buoyancy, phase transformations, geometry and mantle viscosity. Crustal material can remain as compositionally distinct domains, mix, or react with surrounding peridotite; it need not retain the near-surface eclogite assemblage indefinitely as pressure rises. Subsequent mantle convection can return recycled components toward spreading regions or mantle plumes, where partial melting incorporates them into new basalt. A primary study of Hawaiian basalt illustrates a specific model in which melts derived from recycled crust react with peridotite to form a pyroxenite source. It is an interpretation of combined chemical evidence, not a unique label attached to every enriched basalt.^[ref]

How the pathways can be observed. Seismic reflection profiles, drilling, wedge geometry and balanced structural sections constrain accretion, fluid escape and tectonic erosion. High-pressure metamorphism, mineral zoning, radiometric dating and pressure-temperature-time paths in exhumed rocks reveal burial, reaction and return. Their reconstructable histories are especially useful when deformation and fluid inclusions are linked to the dated mineral growth. Earthquake distributions, seismic tomography, seismic attenuation and electrical conductivity constrain the geometry and physical state of the slab and mantle wedge, although fluids, melts and temperature can produce nonunique signals.

Arc magmas carry complementary chemical evidence: relatively mobile barium, boron, lead or strontium can trace fluid-mediated transfer, whereas sediment-sensitive thorium and isotope ratios can help identify other contributions. Element mobility changes with fluid composition and pressure-temperature conditions, so no one element is an infallible tracer. Boron isotope ratios, radiogenic strontium isotope ratios and lead isotope ratios should be interpreted with measured incoming sediments and altered oceanic crust. Arc gases and forearc fluid discharge constrain volatile return; ocean-island basalts, mantle xenoliths and diamond inclusions provide possible evidence of deeper recycling. Each observation samples only part of the system.

For any element $X$, a useful control-volume statement isHere permanent accretion is counted outside the chosen reactive control volume, return includes fluid, melt and exhumed-solid exports, and transient underplating or reservoir growth contributes to storage until exported. For a different boundary choice, the terms must be regrouped consistently. The three incoming layers can have very different transfer efficiencies for bulk rock, water, carbon and a particular trace element. A small mass of escaping fluid can dominate an element's return while most of the solid slab continues into the deep mantle.

Specify what “onset” means. I would distinguish the first local, perhaps short-lived subduction event from persistent, self-sustaining recycling in a connected global network of tectonic plates. Modern plate tectonics combines relatively rigid moving lithosphere, localized boundary deformation, seafloor spreading and sustained convergence. A sinking drip or a single imposed descending sheet is not sufficient evidence of that global regime. Different processes may become established at different times, so a transition interval may be a better hypothesis than an instantaneous worldwide switch.

Construct competing physical models. Solve conservation of mass, momentum and energy for evolving mantle convection and lithosphere, including temperature- and stress-dependent rheology, brittle yielding, plastic deformation, melting, compositional buoyancy and relevant phase changes. Use geodynamic modelling to compare a stagnant lid, episodic overturn, drip-dominated recycling and a sustained mobile lid. In a simplified incompressible creeping-flow framework,Here $\mathbf{u}$ is velocity, $p$ pressure, $\eta$ effective viscosity, $\dot{\mathbf{ϵ}}$ the strain-rate tensor, $\rho$ density, $\mathbf{g}$ gravitational acceleration, $c_p$ specific heat capacity, $k$ thermal conductivity, $H$ radiogenic power per unit mass, and $\Phi$ dissipative heating per unit volume. These equations are a framework: compressibility, latent heat and chemical transport may require additional terms. A simple yield stress law ${\tau }_y=C+\mu (\rho gz-P_f)$ represents cohesion $C$, friction coefficient $\mu$ and pore pressure $P_f$, illustrating why hydration and interface weakening can change the outcome.

Vary the early-planet controls rather than assuming the answer. Explore plausible initial mantle potential temperature, radiogenic heating, crust thickness and composition, water inventory, grain-size evolution, surface temperature, and boundary strength. Hotter mantle lowers viscosity but can produce thicker buoyant basaltic crust; it does not automatically make sustained subduction easier. A cold strong plate can transmit slab pull yet resist bending and initiation. Water weakening of rocks, melt-induced weakening, inherited faults and mechanically weak sediments can change this balance. Test proposed triggers such as mantle plume impingement or gravitational instability without presupposing an established global plate-force system.

Useful diagnostics include the Rayleigh number,with thermal expansion coefficient $\alpha$, layer thickness $D$, thermal diffusivity $\kappa$ and temperature contrast $\Delta T$, together with surface mobility, the localization of deformation, plate longevity, melt production and material recycling. A large Rayleigh number establishes vigorous convection under suitable conditions, not necessarily a mobile lid. A convincing model should sustain one-sided subduction, connected boundary evolution and seafloor spreading after any artificial trigger is removed. Compare initiation thresholds with continuation thresholds: hysteresis can allow an established regime to persist under conditions in which it would not spontaneously start.

As an example of the type of experiment, three-dimensional models have tested mantle plume-induced subduction initiation, finding sensitivity to plate buoyancy, magmatic weakening and hydration of the interface. The scientific use of such results here is to identify a feasible mechanism and its predicted signatures, not to adopt their success as a direct measurement of an ancient onset date.^[ref]

Connect a model regime to a date through independent evidence. Couple the experiments to candidate planetary cooling histories, then predict observable pressure-temperature-time paths, magmatic compositions, deformation patterns and recycling rates for each regime. If independently dated rocks constrain a thermal or compositional state, ask which regimes can reproduce that state and the observed geological associations. Repeat over uncertain initial conditions and constitutive parameters, check numerical resolution, and compare two- and three-dimensional results where geometry matters. Models that prescribe a descending slab, fix plate velocity, or insert a weak boundary may be excellent forward models of a working subduction zone but do not demonstrate spontaneous initiation.

The result is a conditional statement: persistent plate tectonics is physically feasible under a range of conditions, and dated geological observations may favor that range at a particular time. Uncertain cooling histories, several feasible regimes and uncertain rheology prevent numerical feasibility alone from dating onset. The strongest constraint comes when the same parameter range explains multiple independently dated observations and alternative regimes fail those tests.

Effective stress and frictional strength generally increase with depth through the shallow crust undergoing brittle fracture. At $1\mathrm{km}$, small overburden and often weak, fractured near-surface materials limit the shear stress and elastic energy that can accumulate. Some shallow materials also undergo stable sliding rather than unstable earthquake nucleation. Near $15\mathrm{km}$, a locked fault prone to brittle fracture can sustain a much larger differential stress before frictional instability releases stored elastic energy.

This depth is plausibly near the lower, strong part of the seismogenic layer: below it, increasing temperature eventually promotes material creep across the brittle-ductile transition. Once nucleated at depth, a dynamic finite-source rupture can propagate into shallow regions that could not independently nucleate the same event. Depth alone does not determine the nucleation point: local pore pressure, frictional stability and stress heterogeneity also matter.

Interpret the vertical, apparently unoffset fractures as opening geological joints. Their normals are east–west, so the least principal stress ${\sigma }_3$ must be east–west and sufficiently tensile in effective stress. The other two principal stresses lie within the vertical north–south geological joint plane; their ordering cannot be recovered from this fracture orientation alone.

One example, with compression positive and coordinates east, north, up, isHere ${\sigma }_1$ is vertical, ${\sigma }_2$ north–south, and ${\sigma }_3=-3\mathrm{MPa}$ east–west. The geological joint opens east–west when its tensile strength is $3\mathrm{MPa}$. An isotropic pore pressure can make the east–west effective stress tensile even while the total normal stress remains compressive. Draw opening arrows pointing away from the vertical fracture, rather than sliding arrows along it.

The reverse fault requires strong east–west compression; the vertical opening geological joints require the least principal stress in the effective stress tensor east–west. They need not have formed at the same time or under the same local stress field. Compression can form the reverse fault during burial and crustal shortening; subsequent tectonic exhumation, unloading, relaxation, cooling or another tectonic episode can rotate the principal stress ordering and permit opening geological joints.

Furthermore, local stress near faults, layer interfaces and topography can differ from the regional stress. A high pore pressure shifts every Mohr circle left by $p$, helping tensile failure, but isotropic pore pressure alone cannot exchange the east–west and vertical principal stress ordering. Determine relative timing using cross-cutting relations, offsets of geological joints, and any generations of fracture-filling minerals. Spatial coexistence is therefore evidence of a deformation history, rather than a requirement that the two failure diagrams describe one simultaneous homogeneous stress state.

Potential leakage paths are permeable or reactivated faults and connected fractures cutting the sealing sequence; lateral migration through connected reservoir layers into an unsealed outcrop or spill point; and natural or artificial penetrations such as poorly sealed abandoned wells or a compromised injection well. Injection raises pore pressure, potentially overcoming the capillary entry pressure of a seal, hydraulically opening fractures, or causing fault reactivation through reduced effective stress.

Seal continuity and porous permeability matter more than the rock name alone: a fractured or porous limestone need not be an effective caprock. Dissolution, chemical alteration or well-cement degradation can change leakage pathways over time. The underlying salt may be an effective barrier beneath the reservoir, but it does not seal the roof. Conversely, dissolved or residually trapped carbon dioxide is less mobile than a connected buoyant free-fluid phase, so connected porosity and fluid state also affect leakage.