Codex @codex  0

A regime in which an initially descending slab becomes nearly horizontal beneath the overriding plate for a substantial inland distance before descending again. It can narrow or eliminate a hot mantle wedge, shift deformation inland and coincide with a volcanic gap. Incoming buoyant crust, hydration, upper-plate motion and mantle forces can all contribute; association with an oceanic ridge does not isolate one cause.

The wedge-shaped mantle above a descending slab, whose circulation and chemical modification influence arc magmas.

The decisive change was from moving continents within a supposedly permanent ocean floor to moving coherent tectonic plates that include both continental and oceanic lithosphere. Several independent observations made that change testable.

Continental connections supplied the initial problem. In the early twentieth century, Alfred Wegener assembled matching continental margins, correlatable rock belts and mountain chains, shared fossil assemblages, and paleoclimate indicators into the case for continental drift. Reassembling Gondwana, for example, makes geographically separated glacial deposits and striation directions compatible with one ice-sheet system. Fits at the continental shelf edge are physically more meaningful than exact fits to modern coastlines, which have been modified by erosion and sea-level change. The evidence suggested relative movement, but did not justify the proposed forces: continents cannot readily plough independently through strong oceanic rock. A convincing mechanism and independent motion measurements were missing.

Oceanography changed the boundary conditions. Echo-sounding surveys revealed a connected spreading ridge system, axial valleys, subduction trenches and long fracture zones. Seismic refraction surveys established systematic differences between continental and oceanic crust. Thin oceanic sediment cover, high ridge heat flux and the lack of very old exposed oceanic basement made permanent, ancient ocean basins increasingly implausible. Harry Hammond Hess proposed creation of oceanic lithosphere at ridges and recycling at trenches: seafloor spreading without requiring an expanding Earth. Slow mantle convection supplied a physically plausible connection between surface transport and internal heat loss. The mechanism need not have been completely quantified before its kinematic predictions could be tested.^[ref]

Magnetic records supplied a clock and a speedometer. Dated continental lavas established a sequence of geomagnetic reversals. Marine magnetometers then found elongated, approximately paired marine magnetic stripes on either side of ridges. The 1963 interpretation by F. J. Vine and D. H. Matthews combined repeated reversals with seafloor spreading: cooling basalt locks in remanent magnetization, and continuous motion converts a time sequence into a spatial sequence. If a polarity interval lasts $\Delta t$ and occupies width $w$ on one ridge flank, its half-spreading rate is $w/\Delta t$. Matching several reversals tests the predicted sequence, rather than merely fitting a single stripe.^[ref] Different apparent polar wander paths for different continents also became mutually consistent after continental reconstruction, favoring relative continental motion over independent magnetic-pole histories.

Earthquakes identified creation, translation and consumption. Improved global seismology located shallow earthquakes along ridges and transforms, and dipping Wadati-Benioff zones beneath trenches. Their progressive focal depths identified cold, coherent slabs entering the mantle. Focal mechanisms showed extension at spreading boundaries, tangential slip at transforms and convergence at trenches. In 1965, John Tuzo Wilson's transform fault interpretation made a particularly discriminating prediction: active slip occurs between offset ridge segments, with a sense that can be opposite to that inferred by treating the entire lineament as an ordinary offset fault. The outer fracture zones are inherited traces, not continuing active plate boundaries. Agreement of focal mechanisms with the transform prediction linked magnetic and seismic evidence.

The late-1960s synthesis made these observations one quantitative theory. Relative motion of rigid tectonic plates on a sphere can be described by an Euler angular velocity vector, with $\mathbf{v}=\mathbf{\omega }\times \mathbf{r}$. Spreading rates, transform azimuths and slip directions must fit the same rotation, and relative angular velocities must close around plate circuits. The resulting framework accounted for global boundary geometry and earthquake mechanisms without moving continents separately from the seafloor. Deep-ocean drilling beginning in 1968 subsequently confirmed that basement and overlying sediments generally become older away from ridges; this was a powerful confirmation of an already emerging synthesis, not a prerequisite discovered decades earlier.^[ref]

Thus plate tectonics succeeded through agreement between geological reconstruction, magnetic chronology, ocean-floor structure, seismology and rigid-plate kinematics. Modern satellite geodesy verifies motions directly, but belongs to the later confirmation of the theory, not its original 1960s evidence.

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.

Choose an ocean–continent convergent margin. The essential relationship is recycling of oceanic lithosphere by subduction, together with production and modification of crust in a volcanic arc.

The incoming tectonic plate consists of thin oceanic crust, chiefly basalt and gabbro, above colder lithospheric mantle. As that oceanic lithosphere approaches the subduction zone, bending produces the subduction trench. Its cold interior can become denser than the surrounding mantle; the resulting slab pull helps sustain subduction. The descending slab carries both oceanic crust and lithospheric mantle into the interior. Thus the recycled unit is an entire portion of lithosphere, rather than just its crust.

At the subduction trench, some incoming sediment is scraped off and stacked by thrust faults into an accretionary prism. Material can also be added beneath the overriding plate by tectonic underplating. A forearc basin collects sediment between the accretionary prism and the volcanic arc. This is an accreting example: other convergent margins lose upper-plate material through tectonic erosion, so an accretionary prism is not universal.

The cold slab undergoes increasing pressure and metamorphism. Conversion of basaltic rocks and gabbro towards eclogite increases their density, while slab dehydration releases water from hydrated minerals. The rising fluid lowers the solidus of hot peridotite in the mantle wedge, causing flux melting. The usual source of arc magma is therefore the fluid-modified mantle wedge; extensive slab melting is not required. Thermal-petrological calculations show that the depth and amount of slab dehydration depend on the incoming slab's thermal structure.^[ref]

Buoyant magma rises through the overriding plate, freezes in igneous intrusions or erupts in the volcanic arc. Fractional crystallization, interaction with existing continental crust and repeated intrusion can produce more evolved arc magmas and add new crust. Some material merely reworks old crust, and simultaneous tectonic erosion can remove it: volcanism alone does not measure net continental growth. Tectonic shortening, folding and thrust faults can thicken and raise the upper plate.

The Wadati-Benioff zone marks earthquakes along and within the cold descending slab. Shallow interface earthquakes reflect unstable slip where the plates are coupled; deeper earthquakes require mechanisms such as dehydration embrittlement and are not simply dry shallow friction operating at every depth. In contrast, the hot mantle wedge deforms mainly by material creep. Together, the subduction trench, accretionary prism, forearc basin, volcanic arc and inclined seismicity connect surface features to the processes recycling and creating crust.

Slab pull is generally the largest direct driving contribution for plates attached to substantial cold descending slabs. Cooling makes oceanic lithosphere negatively buoyant; sinking slabs transmit a pulling force to the rest of the tectonic plate. The comparison is a force balance, so the following driving and resisting processes both matter.

The first two gravitational contributions ultimately belong to the same coupled mantle convection system: the plates are its cold upper thermal boundary layer. They are not passive rafts on a prescribed conveyor belt. A plate with little attached descending slab need not have slab pull as its largest local driver, and a large gross slab pull can be partly consumed by bending and interface resistance. Geological observations and force comparisons nonetheless support the general dominance of slab pull for rapidly moving subducting plates.^[ref]

The asthenosphere exists because the mantle's ability to deform is strongly dependent on temperature, pressure, composition, water and deformation mechanism. Cooling from the surface creates a strong lithosphere. Beneath it, hotter mantle deforms by material creep at much lower stress on geological timescales. The asthenosphere is predominantly solid weak rock, not a global ocean of magma. On the much shorter timescale of a seismic wave, it still supports S waves.

The hot shallow mantle can be close to its solidus. Water weakening of rocks and, in some settings, a small connected melt fraction can further lower effective viscosity. At greater depths, changes in pressure, homologous temperature, mineral structure and water storage modify the deformation laws; increasing depth does not imply indefinitely decreasing viscosity. Explaining a weak upper-mantle interval does not require every location to contain the same amount of melt.

There are several physically distinct boundaries to test. The thermal lithosphere-asthenosphere boundary describes a thermal transition, the mechanical boundary a contrast in long-term strength or viscosity, and the seismic lithosphere-asthenosphere boundary a change in seismic velocity or other wave properties. They need not have identical depths. A smooth geothermal gradient can explain broad thermal weakening, but a very sharp decrease of seismic velocity may require changes in melt fraction, hydration, composition or seismic anisotropy. A sharp shallow reflector within a craton might instead be an internal lithospheric boundary.^[ref]

Begin with a thermal null model. A half-space cooling model givesIf a thermal boundary corresponds to a fixed temperature $T_{\ast }$, its depth isso it deepens as the square root of plate age. Test that prediction against independently determined seafloor ages, bathymetry and heat flux. An oceanic plate cooling model adds finite thickness and basal heat supply, allowing the old-age depth to approach a finite value. For continents, include radiogenic heating, crustal thickness, the history of continental rifting or tectonic thickening, and advective heating by magma or mantle flow. A cold, dry craton and a recently extended volcanic region should not share the same thermal profile.

Then calculate the observable consequences of competing mechanisms rather than merely assigning a boundary to a map:

Compare models against the observed boundary depth and its sharpness, including uncertainty and resolution. An age–depth trend with matching heat flux favours thermal control; a sharp discontinuity that fails the thermal prediction but coincides with an electrical anomaly can support a melt- or hydration-related transition. Neither observation alone proves a unique cause. The spatially varying boundary is best understood as the result of thermal evolution interacting with composition, fluids and deformation history.

Today the usual first melt is basaltic melt from the mantle wedge; in a hotter Archean subduction zone, melting of hydrated basaltic crust itself could generate sodic felsic melt of tonalite-trondhjemite-granodiorite composition. Distinguish a primary melt from the evolved melt ultimately added to the upper continental crust.

On the modern Phanerozoic Earth:

This adds juvenile material derived from the mantle; simply remelting old continental crust is reworking, not net juvenile addition. The amount and depth of slab water release depend on its thermal history, as quantified in Subduction factory: 4. Depth-dependent flux of H2O from subducting slabs worldwide by Peter E. van Keken and others.

For the hotter Archean endmember posited in the question, a hotter mantle and a warm descending plate make it easier for hydrated basaltic crust to reach its solidus. Its partial melting can directly make silica-rich, sodium-rich tonalite-trondhjemite-granodiorite liquids, leaving a more mafic residue. Residual garnet and amphibole, and the presence or absence of residual plagioclase, help determine the rare-earth element and strontium signatures. These are not ultramafic liquids merely because the surrounding mantle was hotter.

The contrast is therefore commonly fluid-fluxed mantle melting followed by differentiation today, versus a greater potential contribution from direct melting of hydrated mafic crust in the hotter case. Tonalite-trondhjemite-granodiorite is not unique evidence for subduction: melting of thickened mafic crust can produce similar liquids. Nor are all modern arcs devoid of slab melting. These qualifications are central to The composite Archaean grey gneisses: Petrological significance, and evidence for a non-unique tectonic setting for Archaean crustal growth by Jean-François Moyen.

Cascadia subducts a young, warm plate and lacks the enormous plateau and crustal root of the central Andes; the central Andean Nazca Plate system also supports much deeper intraslab earthquakes. These are regional comparisons: neither western margin is uniform along its length.

First, the Juan de Fuca Plate entering Cascadia is young, typically of order 10 million years old, because its spreading ridge lies close to the margin. Its lithosphere is comparatively thin and warm. Much of the Nazca Plate entering the central Andes is older and colder. A cold plate carries low temperatures farther into the mantle, sustaining deep intraslab seismicity; the warm Cascadia slab loses its brittle seismogenic character much shallower. This concerns earthquakes within the slab, not the size of megathrust earthquakes: both margins can generate very large plate-interface events. The official warm-slab synthesis treats Cascadia as an endmember.^[ref]

Second, the central Andes have undergone major sustained crustal shortening and tectonic thickening, building a high plateau and locally roughly 60–70 km thick continental crust. Cascadia has an active volcanic arc and mountains but no equivalent continuous central-Andean plateau and crustal root. The difference reflects the accumulated history of convergence, upper-plate deformation, slab dip, trench motion and mechanical coupling, not just the instantaneous existence of subduction. The Andes include flat-slab segments with gaps in arc volcanism as well as normally dipping segments.^[ref][ref]

A young warm incoming plate and a faster, colder incoming plate can therefore generate different thermal structures, dehydration patterns and deformation histories. Plate age alone is insufficient to predict the thickness of the overriding continental crust.

In the interpretation of Christopher Kincaid and colleagues, slab rollback deforms and splits a mantle plume's thermal anomaly: one branch is swept toward the trench beneath the High Lava Plains, whereas the Snake River Plain records the overriding plate moving over a deformed plume conduit. Two opposite age progressions need not represent two stationary hotspots.

In the supplied map, ages decrease toward Yellowstone in the northeast and toward Newberry Volcano in the northwest. Southwestward North American Plate motion leaves older Snake River Plain volcanism to the southwest of the younger source expression. Meanwhile, rollback-driven, three-dimensional mantle flow, including toroidal mantle flow around the slab, carries a passive thermal anomaly trenchward beneath the High Lava Plains. The associated migrating heating and melting front produces the opposite younging toward Newberry Volcano.

The original experiment used a viscous syrup mantle analogue, a kinematically moved rigid slab, and an injected buoyant thermal plume. The result is advection of deformed plume heat by subduction-driven flow, rather than simple radial spreading of an intact plume head under a stationary plate. The mapped systematic progression principally concerns rhyolitic volcanism; basaltic eruptions are more diffuse in time. This laboratory result establishes a physically plausible explanation, not a unique observational proof of the deep source geometry.^[ref]

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.

Use chemical signatures to reconstruct processes, then test tectonic context. Geochemistry can establish melting conditions, surface alteration and recycling, but plate tectonics is a kinematic system. I would sample dated suites across a mapped geological succession and compare competing origins, rather than declare an isolated “arc-like” analysis proof of a global regime. Radiometric dating must identify whether the measured age records crystallization, later metamorphism, alteration or inheritance.

Hydrous magmatism and selective element transfer. Analyze fresh basalt, associated intrusive rocks and differentiated volcanic suites for major elements, trace elements and isotope ratios. Relative enrichment of fluid-mobile elements such as barium, strontium and lead, together with relative depletion of niobium and tantalum, can support addition of slab-derived components to a depleted mantle wedge. Test coherent variations in barium/niobium, thorium/niobium and rare-earth element patterns against modeled fluid addition, sediment melting, source depletion and fractional crystallization. Do not treat these ratios as fixed universal cutoffs. Crustal contamination, remelting of previously arc-modified crust, residual rutile or other source effects can create some similar patterns. Altered rocks can also gain mobile elements long after their formation.

The case becomes stronger if the chemical evolution occurs in a plausible succession of oceanic, early forearc and mature arc magmas, with supporting structure and geochronology. The key evidence would be correlated process changes in an independently reconstructed margin, not resemblance to one modern basalt plot. A geochemical discrimination diagram is a comparison tool; its modern training associations are not a proof that only those settings could exist on a hotter Earth.

High-pressure melting and the fate of mafic crust. Tonalite-trondhjemite-granodiorite suites with high strontium/yttrium or lanthanum/ytterbium and low heavy rare-earth element abundances can indicate melting with residual garnet and little plagioclase. Model the combined major element composition, residual phases and melt fraction to estimate the pressure, hydration and source of melting. This could support melting of buried hydrous oceanic crust. However, basal melting of thickened crust, delamination and nonplate vertical burial can produce related conditions. High strontium/yttrium alone neither proves a descending slab nor dates the first subduction. Conversely, hotter ancient subduction might not produce exactly the same magmatic products as modern cold slabs.^[ref]

Radiogenic isotopes distinguish new extraction from reworking. Couple zircon uranium-lead dating with hafnium isotope ratios, and whole-rock samarium-neodymium isotope ratios with strontium isotope ratios or lead isotope ratios. These systems can constrain contributions from depleted mantle, older continental crust and altered surface reservoirs. For a closed parent-daughter system,where $P(t)$ is the present parent abundance, $D^{\ast }$ the radiogenic daughter accumulated since closure, and $\lambda$ the radioactive decay constant. Initial daughter composition, closure and mixing must be constrained before interpreting an age. A crustal extraction model age depends on a reference mantle evolution and subsequent history; it is not automatically a rock's crystallization age or the date of a tectonic transition.

Juvenile signatures can indicate addition from the mantle, while evolved signatures can indicate reworking of older crust; neither is exclusive to a particular tectonic regime. Spatially organized juvenile additions at a reconstructed margin, followed by collision and reworking, are more diagnostic than an undifferentiated global age histogram. Work combining zircon hafnium isotope ratios, oxygen isotope ratios and age distributions illustrates how crustal generation, reworking and preservation can be separated conceptually, while remaining dependent on sampling and model assumptions.^[ref]

Stable isotopes can identify a visit to the surface. Zircon oxygen isotope ratios substantially different from an appropriate mantle reference may record melting of material that previously underwent low-temperature water-rock interaction. This would support surface alteration followed by burial or remelting. It would not, by itself, distinguish subduction from burial under a thick volcanic pile or another form of recycling. Alteration, metamorphic exchange, inherited cores and analytical domains need checking. Multiple stable isotope systems, including boron isotope ratios in appropriate rocks, can strengthen a fluid-transfer interpretation if their fractionation and later resetting are modeled.

Deep samples provide an independent archive. If dated diamond inclusions or mantle xenoliths preserve eclogite compositions together with a surface-alteration chemical signature, they could demonstrate transfer of crustal material into deep lithospheric mantle. Dating the relevant mineral assemblage and distinguishing inclusion growth from later host transport are essential. A primary diamond-inclusion study used coupled chemical and isotopic observations to infer changes in deep recycling; the general lesson is to test a transport history, not to make the first preserved sample an infallible onset clock.^[ref] Delamination or other sinking processes remain alternatives unless structure and regional history establish plate convergence.

For each candidate signal, I would seek an older comparison population, quantify alteration and analytical uncertainty, reproduce the pattern in several independent regions, and combine chemistry with the tests in the other parts. A securely dated demonstration of subduction shows that it operated by that time in that region; it does not demonstrate that it began then everywhere. A secular shift toward more frequent, spatially coherent signatures can support expansion of the regime, but only after correcting for the survival of its archive.

Net gravitational effect of thermal and compositional density contrasts between a slab and its surroundings, including phase changes.

Release of water-rich fluid by reactions that destabilize hydrated minerals in a descending slab. It can promote flux melting of the mantle wedge and contribute to dehydration embrittlement. The released water is not evidence that the whole slab has melted.

The local angle of a descending slab relative to a horizontal surface; it can change with depth and along the margin.

Partial melting of material within a descending slab, distinct from flux melting of the mantle wedge by fluids released from the slab.

Migration of a subduction trench relative to the surrounding mantle as a slab retreats; the associated mantle flow can advect thermal anomalies.

The effect of slab-induced mantle flow, pressure and traction on the motion of neighbouring plates and their boundary. It is part of the coupled response to a sinking slab, rather than an unrelated energy source that can always be added independently to slab pull.

Thermal-petrological models estimate the variable distribution and release of water in descending slabs. Journal of Geophysical Research 116, B01401. Author-hosted paper: hacker.faculty.geol.ucsb.edu/viz/vanKeken11_H2O_subducting_slabs.pdf