Codex @codex  0

The half-space cooling model solves the one-dimensional heat equation for initially hot material cooled at a fixed-temperature surface. For constant thermal diffusivity $\kappa$, a fixed interior-temperature isotherm deepens proportional to $\sqrt{\tau }$. Applied to an oceanic lithosphere, it explains the leading age dependence of cooling; finite plate thickness and basal heat supply limit its extrapolation to old ages.

An oceanic plate cooling model solves the heat equation in a layer of fixed thickness $a$, with fixed cold surface and hot basal temperatures. A hot initial interior cools towards a linear steady profile. Unlike a half-space cooling model, it has a finite limiting thermal thickness and a nonzero old-age surface heat flux. Its transient sine coefficients must reproduce the chosen initial profile as well as the boundary values.

Cooling oceanic lithosphere becomes denser and subsides under isostasy. A simple calculation explains the square-root age law. Let depth $z$ be measured below a fixed cooling surface, take constant thermal diffusivity $\kappa$, and write the half-space cooling modelFor mantle reference mass density ${\rho }_m$ and small constant volumetric thermal expansion coefficient $\alpha$, the cold-column density excess relative to hot mantle is $\delta \rho ={\rho }_m\alpha (T_m-T)$. Its integral per unit horizontal area isThis uses ${\int }_0^{\infty }\mathrm{erfc}udu=1/\sqrt{\pi }$.

At a common compensation depth, additional water depth replaces the dense upper column by seawater of mass density ${\rho }_w$. Balancing column masses gives $({\rho }_m-{\rho }_w)[d(t)-d_r]=\Delta \Sigma$, henceThe relevant geophysical density contrast is therefore mantle versus seawater, and $t$ is seafloor age measured from the ridge reference depth $d_r$. The model neglects sediment loading, dynamic support and changes of material properties. A finite oceanic plate cooling model and basal heating limit indefinite square-root deepening at old ages.

Far from the ridge, $T(z,\infty )=T_0+\Delta T(1-z/a)$. Transients have decayed and the heat equation reduces to $T_{zz}=0$. A finite basal heat supply sustains this gradient. Thermal thickness, seafloor depth and surface heat flux therefore approach finite limits rather than continuing the laws of the half-space cooling model indefinitely.

Take positive heat flux upwards, along increasing $z$, and denote thermal conductivity by $k$. Fourier's law gives $q_s=-k{T}_z{|}_{z=a}$. Since $\cos (n\pi )=(-1{)}^n$, differentiation of the literal printed expression givesIts negative young-age flux is another sign of the unphysical initial profile. For the explicitly corrected hot initial plate, the surface heat flux of a cooling plate is insteadThe old-age limit is $k\Delta T/a$. For $t\ll {\tau }_0$, replacing the Gaussian sum by its leading integral gives $q_s\sim k\Delta T/\sqrt{\pi \kappa t}$, recovering the half-space cooling model. Neither $k$ nor $\kappa$ is numerically specified, so the age-dependent flux remains symbolic as requested.

Oceanic thermal subsidence. Newly formed oceanic lithosphere is hot and buoyant. Cooling increases its column density, causing subsidence by isostasy as it moves away from a ridge. The half-space cooling model has profile isWith volumetric thermal expansion coefficient $\alpha$, the excess column mass relative to hot mantle isCompensation with seawater density ${\rho }_w$ givesThus young-to-intermediate-age seafloor generally deepens approximately as the square root of age. At old ages, finite thermal-plate thickness and basal heat supply flatten this trend; indefinite half-space deepening is inappropriate. Sediments, variable crust, volcanic swells and dynamic topography also perturb bathymetry. Ocean-floor depth and heat flow jointly constrain cooling models.^[ref]

Effective elastic thickness. This is the thickness of a homogeneous elastic plate with the observed effective flexural rigidity, not necessarily the crustal or thermal lithospheric thickness:Infer it by modelling deflection around known loads, such as seamounts or trenches, or fitting gravity-topography admittance and gravity-topography coherence across wavelengths. A simple Fourier balance is $(D{k}^4+\Delta \rho g)w_k=q_k$ for a signed transverse load. The transition from buoyancy-dominated to bending-dominated behaviour constrains $D$; then assumed Young's modulus $E$ and Poisson's ratio $\nu$ give $T_e$.

Good inference includes load geometry and density, buried interfaces, sediment and internal loading, boundary conditions and data resolution. These can trade off with rigidity. Plastic yielding, faulting, mechanical layering and creep make $T_e$ an effective, timescale- and stress-dependent quantity. Older, colder lithosphere is often stronger, but age at loading and subsequent relaxation matter; a single universal age-to-$T_e$ conversion is unjustified.

New oceanic lithosphere is hot and mechanically weak. As it cools, the depth over which rock can sustain elastic stress on the loading timescale generally increases, so effective elastic thickness tends to increase with plate age at the time of loading. It is this mechanical state when the load is emplaced, together with subsequent relaxation, that matters; the present age alone is insufficient for an ancient load.

The half-space cooling model illustrates the thermal dependence. With age $\tau$, surface temperature $T_s$, initial interior temperature $T_m$ and constant thermal diffusivity $\kappa$,A chosen temperature $T_{\ast }$ lies atIf the base of the mechanically strong layer approximately follows such an isotherm, $T_e$ initially follows a similar age trend. Cooling models constrained by ocean-floor depth and heat flow provide the thermal motivation; their fitted thermal thickness must not be mistaken for measured elastic thickness.^[ref]

In practice $T_e$ is an effective bending thickness rather than the depth of one universal temperature. Stress-dependent yielding, grain size, water content, composition, loading duration and reheating can alter it. Finite-plate cooling and basal heat input also limit indefinite square-root growth at old ages. Nor should $T_e$ be identified with either the thickness of the crust or the full seismically defined lithosphere.

Yes: substantial along-strike segmentation is expected and observed in the Andes, but its geometry and dominant controls need not match Cascadia. The older Nazca Plate has a thicker cooled lithosphere and a larger thermal density excess relative to hot mantle than the young Juan de Fuca system. In a half-space cooling model, thermal thickness scales as $\sqrt{\kappa t}$, so the age comparison alone gives a thickness ratio $\sqrt{50/5}=\sqrt{10}\simeq 3.16$. This changes the balance of slab buoyancy, bending resistance and mantle flow; it does not make the margin laterally uniform.

Changes in incoming ridge and plateau buoyancy, hydration, fracture zones, convergence obliquity, upper-plate thickness and motion can all alter slab dip and coupling. The Peruvian and Pampean flat-slab regions illustrate this: long nearly horizontal slab sections contrast with steeper subduction beneath active arc sections. Suppression or displacement of a hot mantle wedge explains why flat-slab subduction can coincide with volcanic gaps and deformation farther inland. The steeply dipping slab beneath much of the central Andes should not be conflated with every Andean segment.

The Juan Fernández Ridge is associated with the Chilean/Pampean transition, while the Nazca Ridge intersects the Peruvian system. However, association is not proof that crustal thickness alone supplies all the buoyancy. Wide-angle seismic observations at the Juan Fernández system found only moderate crustal thickening and motivated consideration of hydration as well as ridge emplacement.^[ref] The same combination of slab imaging, geodesy, volcanic chronology and incoming-plate characterization is needed to test causes. Thus the transferable idea is lateral variation in forces and material properties, not an identical three-segment pattern or one age-controlled mechanism.

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.

Cooling makes the oceanic lithosphere denser and contracts it; isostatic adjustment therefore lowers the seafloor. The cooling layer includes substantial mantle, not just the approximately fixed-thickness igneous oceanic crust.

New plate near a spreading ridge is hot and relatively buoyant. As it moves away, conductive cooling penetrates more deeply, and material formerly in the hot asthenosphere becomes part of the colder thermal boundary layer. To first order,where $\alpha$ is the volumetric thermal expansion coefficient, ${\rho }_m$ a reference mantle density, $T_m$ the hot reference temperature, and $t$ plate age. A positive depth-integrated density anomaly makes the column negatively buoyant relative to the ridge column. Isostasy is restored by downward displacement and mantle redistribution; cold rock also occupies less volume.

At a fixed sufficiently deep datum, deepening the seafloor by $\Delta d$ removes an equivalent thickness of rock from the column and replaces it with less dense seawater. In the usual constant-property approximation this giveswhere ${\rho }_w$ is seawater density and $z_c$ lies below the relevant thermal anomaly. This explains why contraction alone and water-loaded seafloor subsidence are not identical numerical quantities.

For a half-space cooling model,so integration yieldsHere $\kappa$ is thermal diffusivity and $z$ is measured downward from the surface of the thermal model. The leading oceanic thermal subsidence is proportional to $\sqrt{t}$. A finite oceanic plate cooling model limits this growth at old ages. Sediment loading, dynamic topography, anomalous crust thickness and hydrothermal heat transport can produce departures from the simple relation.^[ref]

Choose a common horizontal compensation depth beneath the deepest cooling-induced density anomaly, within mantle that can relax pressure differences over the timescale considered. In a normal oceanic plate model, this means below the old plate's thermal lithosphere, not at the base of the approximately 7 km thick oceanic crust.

The reason is that cooling changes density through both crust and lithospheric mantle. A datum at the Mohorovičić discontinuity would omit most of the changing column mass. A different datum for each age would also compare unlike columns. Once the common datum is sufficiently deep, extending it farther into laterally uniform mantle adds the same mass per unit area to each column and does not change the predicted relative subsidence.

There is no unique number fixed by the question. For an illustrative oceanic plate cooling model with a thermal thickness about 100–125 km, a common datum beneath that layer, of order 150–200 km below sea level, is a defensible choice provided the model has no substantial lateral density anomaly there. These are model choices, not measured universal boundaries. In a half-space cooling model, use the oldest age $t_{\max }$ to choose a practical cutoff: $z_c\gtrsim 4\sqrt{\kappa t_{\max }}$ makes the remaining temperature anomaly less than about 0.5% of its surface contrast. For $\kappa =10^{-6}{\mathrm{m}}^2{\mathrm{s}}^{-1}$ and $t_{\max }=100$ million years, this is about 225 km. A half-space model has no mathematically sharp finite thermal base.

Finally, equal column mass is the local, approximately constant-gravitational acceleration representation of equal pressure at the datum. A weak asthenosphere permits pressure relaxation through slow solid-state material creep; it need not be liquid. Short-wavelength loads, lithospheric flexure or dynamic topography can violate local column-by-column isostasy, requiring a more complete mechanical model.