Codex @codex  0

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.

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.

Test geological associations and their chronology. The most persuasive hypothetical record would combine rigid-block relative motion, creation of oceanic lithosphere, its destruction at a coherent margin, and repeated assembly and dispersal of continents. I would compile independently dated rock and structural occurrences in age bins, map their associations, and distinguish primary formation from later emplacement. The aim is to evaluate alternative tectonic regimes, not simply to collect the oldest example of a familiar name.

Oceanic lithosphere and sutures. An ophiolite association containing residual peridotite, gabbro, sheeted dikes and submarine basalt could demonstrate construction of an oceanic section. Its case for plate tectonics strengthens if it is associated with an oceanward sedimentary succession, a metamorphic sole, an accretionary prism or mélange, and eventual emplacement along a suture zone between formerly distinct blocks. Date construction, metamorphic heating and emplacement separately. A dismembered collection of ultramafic and mafic rocks need not be a complete ophiolite; an isolated volcanic pile or greenstone belt can also arise in a plateau or vertical-overturn setting. Hotter early mantle may produce oceanic sections unlike a modern template.

Cold burial beside a hot arc. Blueschist and eclogite with well-constrained low temperature/pressure ratios would support unusually deep burial before substantial heating. Coesite or metamorphic diamond in a suitable continental crust protolith could establish ultrahigh-pressure metamorphism, consistent with deep continental burial and return. Reconstruct pressure-temperature-time paths using mineral assemblages, zoning, reaction textures and radiometric dating; a peak pressure without its timing and burial geometry gives a weaker constraint. Eclogite alone only indicates appropriate pressure-temperature-composition conditions and may also form at the base of thick crust.

An approximately coeval, spatially organized pair of low thermal-gradient and high thermal-gradient metamorphic belts would be stronger evidence for juxtaposition of a cold descending plate and a hotter volcanic arc or hinterland. Test the original spatial relation rather than pairing unrelated rocks brought together much later. A statistical study of metamorphic temperature/pressure distributions illustrates how gradual widening and bimodality can be used to evaluate changing tectonic behavior. The ratio is an apparent metamorphic thermal gradient, not necessarily a directly measured steady geothermal gradient.^[ref]

The absence of old blueschist is not a decisive absence of old subduction. A hotter thermal regime may suppress its production, and later heating can transform it. The test should therefore include the distribution of pressure-temperature-time paths predicted by early-planet models, rather than demanding the youngest modern mineral assemblages at all ages.

Kinematic structures and basin evolution. Long, coherent thrust belts, imbricated sediments, folds, regionally consistent shear sense, major strike-slip faults, and suture zones can identify lateral shortening and relative displacement. A reconstructed ocean-facing volcanic arc, forearc basin and accretionary prism with a common polarity is more compelling than isolated shortening. Mountain building alone is not exclusive to plate tectonics: vertical instabilities and impacts can also deform rock. The length scale, sequence, strain partitioning and linkage between structures matter.

For repeated opening and closing, seek a dated progression from continental rifting, through a subsiding passive margin and an intervening oceanic assemblage, to convergence, continental collision and renewed rifting. Matching stratigraphy, sedimentary provenance and structural restoration on opposing blocks can test whether a reconstructed Wilson cycle is geometrically possible. Paleomagnetism from independently dated, demonstrably primary magnetizations could show differential latitudinal motion and rotation of separate blocks. Relative motion between blocks is more useful than motion of one block, which could include true polar wander. Remagnetization, inclination shallowing and longitude ambiguity must be assessed.

Secular changes require an archive model. Compile the incidence, spatial scale and recurrence of these associations through time, normalized where possible by exposed rock area, sampling effort and dating opportunities. Include detrital zircon and sedimentary provenance where they preserve information from eroded source rocks, while recognizing zircon fertility and sedimentary sorting biases. Increased abundance of plate-like associations in several independent regions could support their transition from local episodes to a widespread regime. A sharp break is only one possible outcome; a progressive increase in persistence or geographical coverage would support gradual establishment.

A useful conceptual accounting iswith additional redistribution between age bins if ages are reset or misassigned. This is not a calibrated dating law: it shows why an observed first appearance need not equal a true first formation. Subduction consumes oceanic lithosphere; erosion, metamorphism, tectonic reworking and burial selectively remove or obscure evidence. The oldest archive is especially incomplete, and some apparently diagnostic assemblages only become likely as the planet cools. An age peak may measure enhanced preservation during continental collision rather than an increased production rate.

State the chronological inference precisely. If an unambiguous assemblage forms at time $t_e$, then the required process was operating by $t_e$; an onset time measured forward from planetary formation obeys $t_{\mathrm{onset}}\le t_e$. If both are expressed as ages before present, the inequality reverses: the onset age is at least as old as the assemblage age. Constraining how much earlier onset could have occurred requires well-sampled older records and a credible prediction that the sought evidence would have formed and survived. Several coeval, independent records consistent with sustained horizontal recycling give a stronger inference of widespread plate tectonics than the oldest ambiguous occurrence of any single rock type. Combined with geodynamic modelling and geochemistry, the defensible answer would report a range for first local activity and a potentially different range for establishment of a persistent global system.

Three-dimensional geodynamic modelling tests mantle plume-driven subduction initiation and sensitivity to plate structure and weakening. Nature 527, 221–225; DOI 10.1038/nature15752.

Research author on numerical geodynamic modelling.