Codex @codex  0

Jun Korenaga coauthored a comparison of the Earth's thermal evolution with melting conditions inferred from ancient basalts. Reconstructing primary magmas provides estimates of mantle potential temperature that can be compared with heat-budget and mantle convection models.^[ref]

The comparison separates evidence for hotter ambient mantle from exceptionally hot plume-derived melts.

A large, potentially global region of extensive silicate melting in a young or strongly heated rocky body. Its crystallization involves latent heat, crystal settling or flotation, and chemical fractionation. An early magma ocean is distinct from the later, mostly solid mantle convection of a planet. The Moon's proposed example is the lunar magma ocean.

Mantle convection is buoyancy-driven motion in a rocky planetary mantle. Heating and cooling create density differences; chemical and mineral-phase differences can also affect buoyancy. Although the mantle is mostly solid, deformation by material creep permits slow flow over geological time.

A simple model balances viscous resistance with buoyancy and transports heat by both flow and thermal conduction. Its Rayleigh number compares those effects. Temperature-dependent mantle viscosity and the strength of the lithosphere help determine whether the surface moves as plates or remains a stagnant lid. Seismic structure, surface deformation and heat flux constrain different parts of this model.

For typical well-mixed structures, the dominant processes are:

Representative temperatures must specify the level: Earth has about $288\mathrm{K}$ at the surface (about $255\mathrm{K}$ effective emission temperature); Jupiter has about $165\mathrm{K}$ near one bar (about $125\mathrm{K}$ effective temperature); hot Jupiters commonly have photospheric temperatures of order $1000$–$2500\mathrm{K}$; and the Sun's photosphere is about $5800\mathrm{K}$. Upper layers, nightsides, deep interiors and the solar corona have different temperatures. These are characteristic values, not constant temperatures throughout each atmosphere.

Plate tectonics is the movement and recycling of a planet's lithosphere as discrete plates over a deformable mantle. Mantle convection, the negative buoyancy of cool subducting slabs (slab pull), and gravitational sliding from elevated spreading ridges (ridge push) supply driving stresses; deformation and friction resist motion. The mantle mostly deforms by slow solid-state creep, rather than being a global liquid layer.

Its effects include subduction and recycling of crust, creation of new crust, mountain building, earthquakes and volcanism, transport of internal heat, and recycling of water and carbon. The carbonate-silicate cycle can couple volcanic CO2 supply to weathering and long-term climate.

Three major controls on a super-Earth's tectonic mode are:

A larger mass alone does not establish active plate tectonics. The competition between driving stress, yielding and sustained slab buoyancy must be evaluated for the planet's composition and thermal history.

Treat the warm ice shell as a very viscous layer and estimate its Rayleigh numberThe temperature-dependent ice viscosity should be evaluated carefully, because deformation is concentrated near the warm base; a useful first estimate uses a representative basal or depth-averaged viscosity. The onset time is obtained by inserting $h(t)$ and $\Delta T(t)=T_m-T_s(t)$ into this expression and finding when $\mathrm{Ra}$ first exceeds the critical value for the shell's mechanical boundary conditions, usually of order $10^3$.

Above onset, mantle convection within the ice transports heat more efficiently than thermal conduction, increases the basal heat loss for a given thickness, and generally limits further thickening. For otherwise similar bodies, a larger planetary radius normally gives larger gravitational acceleration and hence a larger Rayleigh number. Convection is therefore more likely and begins in a thinner or younger shell, although changes in pressure-dependent melting temperature, shell thickness, and viscosity can alter that comparison.

Homologous temperature and deformation. Homologous temperature is $\Theta =T/T_m$, with both temperatures in kelvin and $T_m$ the appropriate pressure- and composition-dependent melting temperature. It measures proximity to melting more meaningfully than an absolute temperature alone. Pressure, composition and water can alter both $T_m$ and the flow law, so equal $\Theta$ does not imply identical rheology.

A useful thermally activated material creep law iswith stress $\sigma$, grain size $d$, activation energy $E^{\ast }$ and activation volume $V^{\ast }$. Water activity can add another factor. Increasing temperature raises the rate at fixed stress when the activation enthalpy is positive, greatly lowering effective viscosity. Solid material can therefore flow on geological timescales well below melting.

Diffusion creep is approximately linear in stress: bulk-diffusion Nabarro-Herring creep has grain-size scaling $d^{-2}$, while grain-boundary Coble creep has $d^{-3}$ in their ideal regimes. Dislocation creep usually has a stress exponent greater than one and weaker grain-size dependence. Its effective viscosity scales as ${\eta }_{\mathrm{eff}}\propto {\sigma }^{1-n}$ at fixed temperature, so stronger forcing changes the apparent viscosity. At lower homologous temperature, brittle/frictional deformation or high-stress dislocation mechanisms can become important. The dominant process also depends on pressure and confinement; low temperature is not by itself a guarantee of brittle failure deep in a planet.

A deformation mechanism map compares the competing rates or strengths as functions of stress and homologous temperature at specified grain size, pressure and composition. Fine grains favour diffusion creep; larger stress favours nonlinear dislocation creep. There is no universal temperature dividing all these regimes. This sensitivity helps explain a strong outer lithosphere above a slowly flowing solid mantle and why modest thermal or water variations can reorganize deformation.

Determining the actual planform of mantle convection. The task is to infer the geometry of upwellings, downwellings and connecting flow in the present Earth, rather than to draw a preferred laboratory convection cell. Plate velocities and the distribution of subduction zones directly constrain surface motion and likely cold downwellings. Ridges, volcanism and hotspots provide further clues, but a ridge need not mark the centre of a deep plume, nor is every hotspot uniquely explained by one plume geometry.

Seismic tomography supplies three-dimensional velocity structure from body-wave travel times, waveforms and surface-wave dispersion. Fast slab-like structures can often be traced beneath subduction zones into the lower mantle. The 1997 global P-wave study reported long narrow lower-mantle structures connected with convergence and interpreted some slab segments as reaching the deepest mantle; that is evidence for deep circulation, not proof that every slab follows one route.^[ref] Temperature affects seismic speed, but composition, phase changes, water and anisotropy also matter, so velocity is not itself a flow map.

Gravity and the geoid constrain the spatially integrated density field and its induced boundary deformation. Dynamic topography supplies information about vertical stresses, after separating crustal structure, isostatic loads and thermal subsidence. Seismic anisotropy can constrain deformation and fabric orientations, subject to a mineral-physics relation between fabric and flow. Discontinuity topography and the mantle's heat budget provide further tests. Each observation has a different sensitivity to density, temperature and viscosity, making their combination more useful than any one image.

Compare these observations with three-dimensional geodynamic models that solve creeping-flow force balance, mass conservation and thermal transport, including buoyancy from temperature and composition, realistic plate motions, viscosity structure and phase boundaries. Predict seismic observables, gravity, geoid and dynamic topography from the same model, and assess data coverage and resolution before interpreting a feature. Whole-mantle exchange, impeded or deflected slabs near transition zones, and broad deep anomalies can coexist; the evidence does not require identical stacked or single-layer cells everywhere. The actual planform is constrained by joint seismic, gravitational, surface-kinematic and thermal evidence, with substantial nonuniqueness in its detailed flow and history.

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.

Distinguish the depth-dependent mantle temperature from mantle potential temperature $T_p$: the latter is the temperature a parcel would have after reversible decompression in an adiabatic process to a chosen reference pressure. A convecting mantle is approximately adiabatic through much of its interior, with steep thermal boundary layers near the surface and the core-mantle boundary. It is neither isothermal nor wholly molten.

Early accretion, large impacts and gravitational planetary differentiation supplied heat; metal segregation released gravitational energy and extensive melting could produce a magma ocean. Rapid surface heat loss then crystallized much of the silicate mantle. Subsequent evolution reflects competition between radiogenic heating, heat entering from the core and heat carried out through the upper boundary. In a schematic mantle energy balance,Here $C_m$ is mantle heat capacity, $H_m$ mantle radioactive power and $Q_m$ heat leaving the mantle; continental crustal radioactive production is not included in $H_m$. The convective Urey ratio $H_m/Q_m$ compares internal production with mantle heat loss. A value below one does not alone prove cooling, because core heat also enters; cooling requires $H_m+Q_{\mathrm{cmb}}<Q_m$. Radioactive decay reduces the heating supply, while mantle convection, plate tectonics, melt extraction and changing viscosity control the heat loss. A hotter mantle does not automatically imply proportionally faster plate motion: melting, dehydration and boundary strength can alter that relationship.

The robust long-term picture is an energetic early Earth followed by substantial cooling and a presently mostly solid, convecting mantle. The shape of the solid-state history is less certain. Petrological reconstructions commonly place late-Archean ambient $T_p$ above modern values. One reconstruction infers about $1500$–$1600{}^{\circ }\mathrm{C}$ at $2.5$–$3$ Ga versus about $1350{}^{\circ }\mathrm{C}$ today, and permits warming before that late-Archean maximum after the initial magma-ocean cooling.^[ref] This is a model-dependent ambient-mantle history, not proof that the accreting Earth was cooler than the Archean mantle or that every mantle reservoir followed one curve. An early molten transient, subsequent solid-state heating or cooling, and later secular cooling must be distinguished.

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.

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.

A planetary ice shell is a solid water-rich layer overlying an ocean or deeper planetary interior. Its thickness can be controlled by thermal conduction, basal heat flux, surface radiation, and solid-state mantle convection.

The planetary mantle is the predominantly silicate rock region beneath a rocky planet's crust. On geological timescales, solid rock can undergo material creep and support mantle convection without wholesale melting.

A planet whose outer mechanically strong shell does not participate continuously in mantle convection as a connected system of independently moving tectonic plates.

Institutional explanation by the British Geological Survey. Its plate-driving discussion distinguishes mantle convection, ridge push and slab pull and explains the association between subducting edges and plate speed.