Codex @codex  0

The geothermal gradient is the rate at which temperature varies with depth. For steady one-dimensional thermal conduction with constant thermal conductivity $k_T$, the magnitude of upward heat flux is $q=k_{T}dT/dz$.

A heat flow probe combines subsurface temperature measurements with an estimate of thermal conductivity to infer conductive heat flux. Multiple sensor depths help separate a stable geothermal gradient from surface-temperature transients and installation disturbance.

The sign depends on the coordinate convention. Let $z$ increase downward and let $k>0$ be the vertical thermal conductivity. Fourier's law gives downward heat flux $q_z=-kdT/dz$. The upward-positive heat flux is thereforeA subsurface temperature that increases with depth thus corresponds to positive upward conductive loss. The underlying gradient law is the same as in an ordinary conducting slab. ^[ref]

For example, a measured gradient of $25\mathrm{K}\mathrm{k}{\mathrm{m}}^{-1}$ and $k=2\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$ imply $q_{\mathrm{up}}=0.050\mathrm{W}{\mathrm{m}}^{-2}$. This inference requires representative thermal conductivity and an appropriate conduction model. Conductive heat flux can differ from total heat transport where fluid circulation also carries substantial energy.

Use the supplied surface temperature, $420^{\circ }\mathrm{C}\simeq 693\mathrm{K}$, as the premise. With otherwise comparable material and geothermal gradient, this hotter boundary increases thermally activated material creep at shallow depth and places a chosen mechanical cutoff temperature closer to the surface. That effect alone would favour a thinner elastic layer than in a cool terrestrial oceanic plate.

Water weakening of rocks makes dryness act in the opposite direction. Water can enhance material creep; sufficiently dry material can remain strong at temperatures at which its wet equivalent relaxes. Experiments on dry diabase provide a concrete basis for a strong dry crust and possible mechanical coupling between crust and mantle, rather than requiring a weak crustal layer solely because the surface is hot.^[ref]

For an illustrative steady conductive geotherm with thermal conductivity $k_T$ and heat flux $q$, a fixed strength cutoff givesTherefore the comparison would scale aswhere all temperatures are compared consistently and a dry rheology may have a higher effective cutoff $T_{\ast ,V}$. This is a conditional strength-layer estimate, not a universal definition of $T_e$.

Venus can retain appreciable elastic thickness despite its hot surface; temperature favours thinning, while dryness favours strengthening. The supplied facts alone determine neither a numerical $T_e$ nor a definite ordering relative to Earth's oceanic plates. The comparison requires heat flow, subsurface thermal history, material rheology, interior water content and loading timescale. Little atmospheric water is evidence relevant to the hypothesis of dry rocks, not a direct measurement of deep water activity.

I would test whether the ancient farside crust grew locally with an early thermal/compositional gradient or acquired a substantial later accreted layer. The research programme should estimate a globally consistent crustal mass and chronology before tuning either mechanism to a single attractive map.

First, map crustal geometry with independent density constraints. Jointly invert GRAIL gravity and laser-altimeter surface topography, allowing spatially varying density and porosity rather than converting every anomaly to thickness. A new broad-band seismic network should include farside stations, long-lived polar power systems and a communications relay. Natural impacts and moonquakes, supplemented by calibrated active sources near landing sites, would constrain layer velocities, attenuation, crustal interfaces and the crust-mantle boundary. Crossing ray paths and receiver-side conversions are essential to distinguish a thick crust from an incorrect density or velocity assumption. A coherent buried layer or age-correlated discontinuity across the farside would favor addition or stacking; a gradual hemispheric compositional gradient would better support asymmetric local growth. Neither pattern is guaranteed to survive later impacts.

Second, obtain samples that address the initial contrast. Collect geographically distributed ancient highland rocks from both hemispheres, using basin ejecta, central peaks and drilling to sample beneath the impact-gardened regolith. Return samples with mapped geological context. Mineral chemistry and bulk Ca–Al–Mg–Fe abundances would test a primary refractory gradient; crystallization ages would date crustal construction, while shock-reset ages would date overprinting. Isotopic and trace-element measurements could identify an accreted population, but similar isotopes are possible for material derived from the same early debris system. A young deposit draped over older crust would favor accretion; coeval flotation-related crust with a systematic chemical gradient would favor local asymmetric differentiation. Multiple sites are necessary because a single basin's ejecta can mix distant sources.

Third, measure the thermal and volcanic consequences. Place shielded, sufficiently deep heat flow probes in representative highlands and mare terrains. Combine their geothermal gradients and measured thermal conductivity with orbital gamma-ray constraints on heat-producing elements. Map mare volumes and eruption ages, rather than only their exposed areas, to separate melt production from storage and eruption efficiency. A strong correlation between late volcanism, subsurface heat producers and high inferred mantle temperatures would support radiogenic amplification, while old crustal asymmetry without that correlation would show that later heating is not the entire cause.

Finally, compare forward models against withheld observations. Couple irradiation, vapor transport, condensation, lunar magma ocean mixing and crystallization for the thermal scenario; couple orbital evolution, impact deposition, mixing and later viscous relaxation for accretion. Account for later true polar wander before treating present-day hemispheric orientation as the original formation geometry. Require each model to conserve mass and energy and to reproduce the independent thickness, chemistry, age and volcanic data within their uncertainties. Vary unknown initial conditions and compare predicted observables at sites not used for calibration. The discriminating outcome is a conjunction such as a coherent accreted unit with a consistent emplacement age, or an ancient chemical gradient quantitatively consistent with transport under Earth irradiation. A lower gravity misfit alone cannot establish the historical mechanism, because gravity is a nonunique inverse problem.

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.

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.

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.

A depth-dependent estimate of the differential stress a rock column can sustain at a specified strain rate, composition, thermal structure and fluid state. In a simple competing-mechanism model, failure occurs by whichever of frictional sliding and ductile material creep requires the lower stress, giving the displayed lower envelope.

Effective confining pressure generally increases frictional strength with depth. In contrast, a power-law creep model $\dot{ϵ}=A{\sigma }^n\exp [-(Q+P{V}^{\ast })/(RT)]$ gives ${\sigma }_{\mathrm{creep}}=(\dot{ϵ}/A{)}^{1/n}\exp [(Q+P{V}^{\ast })/(nRT)]$. At fixed strain rate, the rise of temperature along a geothermal gradient generally lowers the required creep stress. Their crossing is a brittle-ductile transition, while composition changes at the Moho can produce additional peaks. A rock strength envelope is not uniquely determined by depth or by the labels “continental” and “oceanic”.