Codex @codex  0

At the interface the pressure is continuous. Immediately below it the temperature is $T_s$ and mean molar mass is $M_d(1+ax)$; immediately above it the dry hydrogen has temperature $T_1$. The ideal gas law therefore givesAn infinitesimal interchange is unstable to convection when the lower fluid is less dense than the upper fluid. HenceEquality gives neutral interfacial buoyancy; reversing the inequality gives a stable density jump. This is an interfacial criterion, not a claim that the whole dry upper atmosphere is stable. Its vertical temperature profile has not been specified. Also, $\mathrm{RH}(p_1)=p_1/P<1$, so infinitesimal displacement at the interface initially involves no condensation or release of latent heat.

There are two competing effects. Increasing $T_s$ directly reduces the lower density, favouring convection. Increasing the mole fraction of heavy water vapor increases its mean molar mass, opposing convection. With $x\ll 1$, the criterion is approximately $T_s>T_1$. With $x\to 1$, it tends to $T_s>ϵT_1\simeq 9T_1$: substantial surface heating need not overcome the compositional stabilization in hydrogen.

For a specified $P_d,\Delta p,T_1$, solve $T_s/T_1-1-ax(T_s)=0$ on the admissible branch to locate the transitions. Its derivative is $1/T_1-a{x}^{\prime }(T_s)$, so there is no universal monotone trend: rapid growth of the saturation pressure can suppress an instability that was present at lower $T_s$. Without numerical parameter values, neither the number of crossings nor a unique interval of unstable temperatures can be asserted.

Before condensation, assume a closed parcel with fixed atmospheric mixing ratio, constant specific heat capacities, and negligible irreversible heating. The first law of thermodynamics for a dry adiabatic process isThe parcel remains in approximate pressure equilibrium with a hydrostatic environment. With the standard approximation that parcel and environmental densities are close, $dp/dz\simeq -\rho g$, henceThis is the dry adiabatic lapse rate, about $9.8\mathrm{K}\mathrm{k}{\mathrm{m}}^{-1}$ for air. More precisely the environmental hydrostatic balance gives $dT/dz=-g{\rho }_{\mathrm{env}}/(c_p{\rho }_{\mathrm{parcel}})$; the familiar constant dry adiabatic lapse rate uses their ratio close to one.

For the spin separation at a smooth quantum Hall edge, linearize the confinement energy at a crossing. Compensating a spin energy difference by a potential-energy change givesHere $V^{\prime }$ is interpreted as the specified slope of electron potential energy. If it denotes electrostatic voltage per length instead, the denominator is $e|V^{\prime }|$. The two edges of a neck and the two sides of an antidot have oppositely directed slopes; the formula is a positive local separation on any one edge. For local filling two, only the $n=0$ spin pair is occupied at the neck minimum: $E_{0,+}+U_{\min }<E_F<E_{1,-}+U_{\min }$. Raising both branches with the confinement produces two crossings on each side, as shown in the original neck sketch. This fixed-slope estimate precedes self-consistent electrostatic screening.

An antidot is a depleted potential hill. guiding centres circulate around its closed equipotentials, with chirality opposite to the outer sample edge. A phase-coherent round trip contributes magnetic phase $2\pi BS/{\Phi }_0$, with ${\Phi }_0=h/e$. A localized orbit resonant with the incident energy can transfer an electron between extended channels. If the outgoing channel returns toward the injecting reservoir, the process is resonant reflection and reduces electrical conductance. If coupling instead connects an otherwise blocked incident channel to the far reservoir, it is resonant transmission and adds electrical conductance. Transmission peaks require that forward connectivity; a closed orbit by itself does not guarantee them. Interchannel scattering can change which orbit participates.^[ref]

For the antidot resonance conductance bounds with one active spin-resolved resonance, writeHere ${\Gamma }_L,{\Gamma }_R$ are energy broadenings; the corresponding escape rates are ${\Gamma }_L/\hslash ,{\Gamma }_R/\hslash$. In a reflection-dominated configuration with two otherwise transmitted channels, $G=(e^2/h)(2-\mathcal{L})$: a $2e^2/h$ background can dip to $e^2/h$. In a transmission-dominated configuration with one spectator channel and the second otherwise blocked, $G=(e^2/h)(1+\mathcal{L})$: peaks rise from $e^2/h$ toward $2e^2/h$. These are the two ideal channel-count examples plotted below. More generally, $G=(e^2/h)(N_{\mathrm{open}}-R)$ or $(e^2/h)(N_{\mathrm{dir}}+T)$; a different spectator count changes the integer offset. The constriction filling alone does not specify every tunnel coupling or resonance height. Equal couplings permit unit resonant transfer in the isolated-level model; asymmetry, temperature and dephasing reduce the contrast. Holding the effective orbit area approximately fixed gives $\Delta B\simeq h/(eS)$ for one branch, not a universal charging-corrected period.

In quantum-dot versus antidot resonance topology, a quantum dot is a potential well holding a finite electron island. Its closed guiding-centre circulation is reversed relative to an antidot at the same field. With two barriers in series, current follows source $\to$ dot $\to$ drain; off resonance there need not be a direct conducting path. A gate shifts a level or an addition electrochemical potential into alignment with the reservoirs and produces a electrical conductance peak. For an isolated coherent spin-resolved level, the zero-temperature peak is at most $e^2/h$ and reaches it only for equal couplings; unresolved noninteracting spin degeneracy doubles this bound. In weak sequential tunnelling, peak height depends on rates and temperature rather than being quantized. Charging produces Coulomb blockade peaks with spacing set by addition energy and gate lever arm. An antidot also supports discrete orbit/charging resonances, but can be side-coupled to already open channels and produce dips. Gate voltage changes both orbit area and confinement, whereas a magnetic sweep directly changes enclosed flux.

Coulomb interactions require a self-consistent confinement potential. A gradual decrease in electron density cannot leave arbitrary uncompensated donor charge at no cost. Partly occupied Landau levels have a large charge response: electrostatic screening flattens their effective energy near the Fermi energy over compressible quantum Hall edge strips. Between them, integer-filled incompressible quantum Hall edge strips have a gap, weak density response and most of the potential variation. In an idealized cross-section the effective branch thus alternates between flat, partially occupied portions and rising portions crossing a gap; its carrier density varies oppositely. The bare levels still have Landau level and spin splittings locally, but the common electrostatic potential is no longer an imposed straight slope. exchange interaction can additionally renormalize spin gaps. Widths must follow electrostatics and the density gradient, so the unscreened $E_Z/|V^{\prime }|$ separation is not a quantitative screened-strip width.^[ref]

Hawaii: short-wavelength volcanic loading and long-wavelength mantle support must be distinguished. A dense volcanic edifice replacing seawater adds mass and gives a positive free-air gravity anomaly. Its load bends the lithosphere down, producing a surrounding moat and an outer flexural arch. The associated displaced density interfaces contribute negative and positive gravity signals at neighbouring distances. Observed gravity near the Hawaiian ridge includes strong positive anomalies over the ridge, negative flanks, and a broader positive belt. At much longer wavelengths, the Hawaiian swell also correlates with a positive gravity anomaly, with a smaller gravity-to-topography ratio that requires compensation; the precise source of support is not determined by that correlation alone.^[ref]

The mechanical distinction can be expressed with a simple elastic-plate model. For load height $h$, downward deflection $w$, plate rigidity $D$, loading density contrast $\Delta {\rho }_l$ relative to water, and restoring density contrast $\Delta {\rho }_b$ at an interface,For a Fourier component of horizontal wavenumber $q$, this gives $w_q=\Delta {\rho }_{l}g{h}_q/(D{q}^4+\Delta {\rho }_{b}g)$. The load's direct gravity contribution competes with the negative contribution of a depressed interface at reference depth $t$. In a linear thin-interface approximation,At short wavelengths rigidity suppresses compensation, leaving relatively high positive gravity admittance. At long wavelengths the plate can deflect and compensate more completely, reducing admittance. Fitting the wavelength dependence constrains effective elastic thickness through $D=E_Y{T}_e^3/[12(1-{\nu }^2)]$, subject to density structure, load distribution and subsurface loading assumptions. This explains why local Airy compensation alone misses the spatially extended flexural response.

The broad swell needs a deeper explanation than the narrow volcanic pile: hot, buoyant mantle and mantle flow can support regional topography, while deep density anomalies and surface elevation jointly determine its gravity. A free-air anomaly keeps the attracting mass of the topography; a Bouguer correction removes an assumed direct topographic attraction and can reveal its compensating mass deficit. Gravity plus topography constrain support as a function of wavelength; neither alone uniquely fixes the mantle density distribution.

For the OR alternative, triple junctions are points where three plate boundaries meet. Each boundary can be a spreading ridge, subduction trench, or transform fault. The three plate velocities satisfyThis velocity triangle is necessary but not enough: a single junction velocity ${\mathbf{u}}_J$ must also keep the meeting point on all three moving boundaries. For symmetric spreading between A and B, the ridge constraint is $[{\mathbf{u}}_J-({\mathbf{v}}_A+{\mathbf{v}}_B)/2]\cdot \mathbf{n}=0$. For a transform, ${\mathbf{v}}_A\cdot \mathbf{n}={\mathbf{v}}_B\cdot \mathbf{n}$, and the junction must have that same normal component. For a trench, its normal migration must match the adopted subduction/rollback model; a trench fixed relative to the overriding plate shares that plate's normal velocity.

In velocity space each boundary supplies a line of possible junction velocities. A common intersection means the junction can migrate while preserving its local configuration; incompatible constraints force boundary reorganization. This is kinematic stability, not the claim that the junction is stationary. A symmetric ridge–ridge–ridge example has three arms meeting at $120^{\circ }$ and a common junction velocity; unequal spreading rates do not require those equal angles. Transform and trench junctions depend on the allowed relative motions and boundary geometry. Stable and unstable evolution were formulated through plate-velocity geometry in the original triple-junction analysis.^[ref]

These original diagrams show an illustrative flexural response and a symmetric ridge-junction construction; they are not maps or measured Hawaiian profiles.

Use the intended hot-ridge state and a linear law of thermal expansion about mantle density ${\rho }_m$ at $T_1$. Let $s$ be depth below the old seafloor. Its excess density isand its excess column mass is ${\int }_0^a\delta \rho ds={\rho }_m\alpha \Delta Ta/2$. Isostasy requires equal pressure at a compensation depth below both columns. Lowering the old seafloor by $h$ replaces rock by water, reducing column mass by $({\rho }_m-{\rho }_w)h$. Therefore the isostatic ridge elevation in a plate model isHere $\alpha \Delta T=0.052$; the result is a first-order estimate and neglects sediments, crustal-composition differences and flexural support. If one instead insists on the literal printed zero-age profile, its column is denser than the old plate and the formal relative elevation reverses sign, $-3.9\mathrm{km}$ to this order. That profile cannot describe the elevated hot ridge posited by the question.

Oceanic magnetic anomalies. Newly formed oceanic crust contains magnetic minerals. As it cools through their magnetic blocking temperatures, it acquires thermoremanent magnetization recording the ambient field polarity. Continuing seafloor spreading carries that record away from the ridge. Geomagnetic reversals therefore generate approximately paired marine magnetic stripes parallel to the ridge, with older crust farther from its axis. A magnetometer measures the superposed field above this magnetized material; subtracting a reference field reveals magnetic anomalies. The anomaly sign and shape depend on magnetization direction, latitude, source depth and geometry, so a positive measured anomaly is not universally identical to normal rock polarity.^[ref]

Matching a sequence of stripes to an independently dated reversal chronology gives plate ages and spreading histories. A stripe of width $w$ formed over interval $\Delta t$ implies half-spreading rate $v=w/\Delta t$ under steady symmetric spreading; full separation speed is $2v$. Unequal widths on opposite flanks can reflect asymmetric spreading or later deformation. Alteration, remagnetization and overlapping magnetic sources complicate inversion. The sequence correlation is stronger evidence than an isolated stripe or a field-strength peak.

Rayleigh-Bénard convection. Consider a horizontal fluid layer of thickness $d$, hotter below by $\Delta T>0$. The conductive state is motionless with a linear temperature gradient. A rising warm parcel is buoyant, but viscosity resists its motion and thermal diffusion erases its temperature anomaly. In the Boussinesq approximation, density variations are retained in buoyancy while the flow is otherwise incompressible. The competing effects are measured by the Rayleigh number and Prandtl number,For stress-free impermeable, fixed-temperature boundaries, linear stationary disturbances with horizontal wavenumber $q/d$ giveMinimizing gives $q_c=\pi /\sqrt{2}$ and the critical Rayleigh number for an isothermal layer, ${\mathrm{Ra}}_c=27{\pi }^4/4\simeq 657.5$. Two rigid no-slip, fixed-temperature boundaries instead give approximately $1708$; the critical value is not independent of boundary conditions.^[ref]

Above onset, buoyancy can sustain overturning rolls and heat transfer exceeds pure conduction, with Nusselt number $\mathrm{Nu}>1$. Nonlinear advection eventually limits the disturbance; larger driving can produce time-dependent or turbulent states. The ideal Boussinesq problem does not require hexagonal cells: planform selection depends on boundaries and non-Boussinesq effects. Convection driven by surface tension at a free surface is a distinct mechanism. An adverse temperature gradient supplies driving, but overcoming diffusion and viscosity requires a finite threshold.

Use the flat two-layer model, with thickness $h=30\mathrm{km}$, crustal shear speed ${\beta }_1=3.5\mathrm{km}{\mathrm{s}}^{-1}$ and mantle speed ${\beta }_2=4.5\mathrm{km}{\mathrm{s}}^{-1}$. For SH-wave travel times in a crustal waveguide, unfolding each reflected path gives all reflection orders at once:with $x$ in kilometres and times in seconds. The direct wave begins at the origin; reflected orders have intercepts $60m/3.5=17.14m\mathrm{s}$ and approach the direct-wave slope at large distance.

The critical angle, measured from the vertical, is $i_c=\arcsin (3.5/4.5)\simeq 51.1^{\circ }$. Associated critical refraction and its multiples giveEach head-wave line is tangent to its reflection curve at the critical distance. In particular the primary head wave exists before it becomes the first arrival; it overtakes the direct wave only around $170\mathrm{km}$.

The slope of a reflected travel-time curve is $p=d{t}_m/dx=\sin i/{\beta }_1$. For $x<x_{c,m}$, $p<1/{\beta }_2$, so the mantle supports outgoing shear radiation and successive reverberations leak energy. For $x>x_{c,m}$, $1/{\beta }_2<p<1/{\beta }_1$: the mantle field is evanescent and ideal interface reflection is total. These leaking and trapped SH-wave reverberations explain the annotations in the original travel-time drawing below. A reflected ray is not itself a separate Love mode: constructive interference of the trapped contributions, including reflection phase, selects the normal modes.

For a horizontally propagating mode with ${\beta }_1<c<{\beta }_2$, define vertical wavenumbersThe free upper surface and welded interface impose the Love-wave dispersion conditionwhere ${\mu }_j={\rho }_j{\beta }_j^2$ are shear moduli. Each phase-velocity branch falls towards $3.5\mathrm{km}{\mathrm{s}}^{-1}$ with increasing frequency. Higher branches join $4.5\mathrm{km}{\mathrm{s}}^{-1}$ at their low-frequency cutoffs. At cutoff $s\to 0$, soThus the first three overtones have maximum periods about $10.77$, $5.39$ and $3.59\mathrm{s}$. The fundamental mode has $n=0$ and no positive-frequency cutoff in this flat layer-over-half-space model: it approaches $c={\beta }_2$ as $f\to 0$, with ever-deeper penetration. The phase-velocity sketch uses an explicitly illustrative equal-density choice, because the question gives speeds but no density or rigidity ratio. That missing ratio changes detailed branch curvature, not the cutoff values or the requested branch topology.

For the actual finite spherical Earth, the very-long-period continuation must instead be described by discrete toroidal normal modes, not a flat infinite waveguide at arbitrary wavelength. The flat-model answer therefore does not assert a continuum of arbitrarily low-frequency global nonrigid oscillations. A global lowest oscillation frequency cannot be computed from the supplied crust-and-half-space data alone.^[ref]

Assume the original ice sheet had time to reach local isostatic compensation, melted rapidly relative to subsequent rebound, and left one exponentially relaxing depression mode. Neglect elastic layering and assign the whole present gravity deficit to that unrelaxed depression. Load balance givesThe supplied sheet-gravity approximation gives the residual amplitudeFrom $h(t)=h_0{e}^{-t/\tau }$,For an order-of-magnitude single-mode estimate, choose the ice-sheet diameter as wavelength, $\lambda \simeq 2R=6000\mathrm{km}$. Taking $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$, the more accurate specified relation then givesA circular load has a spectrum, not one uniquely specified wavelength. Choosing $\lambda \sim R$ instead gives about $1.3\times 10^{21}\mathrm{Pa}\mathrm{s}$; another wavelength convention rescales the answer linearly. The robust inference is viscosity of order $10^{21}\mathrm{Pa}\mathrm{s}$ within this simple model. A realistic estimate would also model gradual melting, other gravity anomalies, density interfaces, lithospheric flexure and depth-dependent mantle viscosity.

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.

Assume a stable, lossless elastic body with the usual symmetry of the elasticity tensor and boundary conditions that eliminate the surface terms in integration by parts, such as a traction-free outer surface. Its elastic stiffness operator is then self-adjoint. When density varies spatially, the precise eigenproblem is generalized:Here $\alpha$ abbreviates the pair of mode labels. Equivalently, ${\rho }^{-1}\hat{H}$ is a self-adjoint operator in the mass inner product. Stability makes ${\omega }_{\alpha }^2$ real and nonnegative. Modes at different eigenvalues are orthogonal, and an orthonormal basis can be selected within each degenerate eigenspace. Completeness permits expansion of a suitable displacement field in these mass-weighted elastic normal modes.

Insert the expansion and project onto ${\mathbf{u}}_{\alpha }^{\ast }$. Mass normalization diagonalizes both the acceleration term and the restoring term. For a symmetric stress glut tensor supported away from the boundary,The symmetry of the stress glut tensor removes the antisymmetric part of the displacement gradient, leaving the strain. The result is one independent forced oscillator per mode, with the stated excitation integral. If the modes instead have mass norm $N_{\alpha }$, that integral is divided by $N_{\alpha }$. A source touching the boundary requires the corresponding surface forcing to be retained rather than silently dropping it.

Rigid translations and rotations can be zero-frequency modes of a free body. Their strain vanishes, so an internal symmetric seismic moment tensor does not excite them through this coupling. Attenuation, radiative boundaries or other non-self-adjoint effects require modifications of the simple conservative expansion.

For seismic source inversion, collect calibrated, consistently filtered samples from all available stations and components into a data vector $\mathbf{d}$. For each trial source position and origin time, compute the corresponding elastodynamic Green tensor in an Earth model, including the instrument response used for the data. Write the symmetric seismic moment tensor in six independent symmetric basis tensors:Each column of $\mathbf{G}$ is the synthetic seismogram for one basis tensor. An off-diagonal basis tensor contains both symmetric entries, so the double contraction and its factors of two remain consistent.

A suitable objective iswith noise covariance ${\mathbf{C}}_d$; a justified penalty can be added for regularization. At fixed $\mathbf{q}$, solve for the six entries of $\mathbf{m}$ by weighted least squares, preferably using a QR decomposition or singular value decomposition rather than explicitly inverting the matrix of the normal equations. Then search over the four nonlinear parameters using a coarse location/time grid followed by local refinement. This variable projection reduces a ten-parameter search to a four-parameter nonlinear search with an inner six-parameter linear fit. Merely having more than ten data samples does not ensure identifiability: the design matrix must have adequate rank and the network must resolve the source's radiation patterns.

Potential errors include inaccurate velocity, density, anisotropy and attenuation in the assumed Earth model; neglected finite rupture and source-time effects; timing, instrument calibration or component-orientation errors; correlated noise; uneven station coverage; and tradeoffs between depth, origin time and the moment tensor. Imposing a double-couple seismic source mechanism without evidence can hide genuine non-double-couple components. Regularization also changes the recovered source and should be reported.

Inspect residuals, singular values, recovery of synthetic sources and uncertainty under resampling or noise perturbations. Linearized covariance or a posterior distribution can quantify uncertainty within the assumed forward model, but does not automatically account for errors in that model. A good waveform fit alone is insufficient evidence that all ten parameters are uniquely recovered.

Take upward deflection as positive, use $k>0$, and write $\Delta \rho ={\rho }_m-{\rho }_c>0$. A positive added load depresses the plate, so $w$ is negative beneath it. Because bending moves the original surface while the load adds thickness, the final topography is $e=s+w$. For the cosine mode the horizontal biharmonic operator multiplies $w$ by $k^4$, givingEliminate $s=e-w$ to obtainThe original loading equation contains ${\rho }_m$, whereas the relation after eliminating the added thickness contains $\Delta \rho$. These must not be interchanged prematurely.

At first order in the relief, represent a displaced density interface by a sinusoidal mass sheet. A sheet of surface-density amplitude $\Sigma$ at depth $z$ produces a downward gravity amplitude $2\pi G\Sigma e^{-kz}$ at the datum. This follows by solving the Fourier-transformed gravitational Poisson equation above and below the sheet: the decaying potential has amplitude $-2\pi G\Sigma e^{-k|z-z_0|}/k$, with the normal derivative jump fixed by the sheet mass.

The top interface contributes ${\rho }_{c}e$; the displaced crust–mantle interface contributes $\Delta \rho w$ at depth $t$. ThusThe second term is negative beneath a positive load: depression replaces denser mantle material by lighter crust. Substituting the deflection-to-topography ratio yieldsThis flexural gravity admittance combines mechanical compensation and attenuation of the buried interface's gravitational signal. It assumes small relief, a flat uniform thin plate, uniform density layers and the stated surface loading, without additional internal density anomalies. Here $t$ is the depth of the crust–mantle interface. It is not automatically the effective elastic thickness $T_e$; the latter enters through $D=E{T}_e^3/[12(1-{\nu }^2)]$.

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.

gravity and surface topography. gravity responds to all density structure, while surface topography records both load and support. First specify the reference field and measurement height. latitude and planetary rotation affect normal gravity. At elevation $h$, the reference field decreases approximately by $2gh/R$; a free-air gravity anomaly restores this elevation effect but retains topographic mass. A Bouguer gravity anomaly additionally subtracts an approximate topographic slab contribution $2\pi G{\rho }_{c}h$. terrain correction replaces the unrealistic slab near valleys and ridges, and spherical geometry matters over large areas. A density choice is required for these reductions; a numerical crustal density is not provided here.

In Airy isostasy, surface topography of height $h$ and crustal density ${\rho }_c$ is compensated by a root of thickness $b$ replacing mantle of density ${\rho }_m>{\rho }_c$. Equal column masses give ${\rho }_{c}h=({\rho }_m-{\rho }_c)b$. The positive gravity of the relief is partly cancelled by the negative density contrast of the root. At long wavelengths a compensated mountain can therefore have a modest free-air anomaly and a negative Bouguer anomaly. Pratt isostasy instead varies column density at a roughly common compensation depth. Neither model says that surface topography itself exerts no gravity.

Potential-field depth filtering makes wavelength crucial. In a planar Fourier approximation, a density interface at depth $z$ is attenuated by $e^{-kz}$ at horizontal wavenumber $k$. For a small Airy-compensated height perturbation and a compensation interface at depth $D$, the gravity-topography admittance is approximatelyLong wavelengths allow strong cancellation between load and root. At short wavelengths the root's gravity is filtered out, so the observed near-surface mass can dominate even if the column is compensated. In a flexural isostasy model, an elastic lithosphere supports and spreads the load; for downward deflection $w$ and imposed load $q$, $D_f{\nabla }^{4}w+\Delta \rho gw=q$, with $D_f=E{T}_e^3/[12(1-{\nu }^2)]$. Thus $w_k=q_k/(D_f{k}^4+\Delta \rho g)$, and the effective elastic thickness controls the support scale.

Estimate admittance from cross-spectra of gravity and surface topography, and compare its wavelength dependence and coherence with local-compensation, elastic-flexure and subsurface-load models. A primary Atlantic study sampled ship-track free-air gravity and bathymetry and compared their Fourier-domain relationship with compensation models, illustrating this observational route rather than inferring a root from the sign of one gravity value.^[ref] Check navigation, instrument drift, terrain/density reductions, sediment and water loads, and independent seismic crustal constraints. Density-depth trade-offs and mantle dynamic support make the inverse problem nonunique; surface topography and gravity together constrain support, but do not uniquely map every density anomaly.

mantle viscosity. A second, independent route concerns deformation rates rather than static compensation. glacial isostatic adjustment follows loading and unloading of a deformable Earth: an ice sheet depresses the surface, and subsequent removal leaves time-dependent uplift, horizontal displacement, sea-level and gravity changes. Fit dated relative sea-level indicators and raised shorelines together with geodetic displacement rates and a reconstructed ice/water load history to a viscoelastic Earth model. Relative sea level includes solid-surface motion, changing geoid and redistribution of ocean water; treating it only as local vertical rock motion biases the inference.

For scale, a sinusoidal surface perturbation of wavenumber $k$ on a uniform incompressible viscous half-space relaxes under gravity with ${\tau }_k\sim 2\eta k/(\rho g)$. A Maxwell material separately has a stress-relaxation scale ${\tau }_M=\eta /\mu$. These are different experiments and generally different times: rebound modes also involve gravity, spherical geometry, a layered mantle and elastic lithosphere. Fit a family of relaxation modes or a full load-history convolution, rather than equating every observed rebound time to $\eta /\mu$.

Use independent ice-margin, exposure-age and sea-level constraints to reduce the trade-off between ice history and viscosity; compare horizontal as well as vertical GNSS motions, and separate present ice loss, groundwater changes and tectonic signals from the older loading response. A primary global analysis jointly tested ice-history and mantle-response predictions against space-geodetic observations, rather than treating a local uplift rate as a unique viscosity measurement.^[ref]

postseismic relaxation probes shallower or faster relaxation after a stress perturbation, but must be distinguished from afterslip and poroelastic deformation. Long-wavelength geoid and mantle-flow models provide complementary constraints when density structure is estimated independently from seismic structure and mineral physics; unknown density-to-velocity conversion and boundary tractions remain trade-offs. Laboratory creep experiments constrain constitutive laws as functions of temperature, stress and grain size, but require a large extrapolation to geological conditions. Combine spatially and temporally distinct observations; there is no single directly measured mantle viscosity independent of depth, rheology and loading assumptions.

First-motion source mechanisms. Locate an earthquake from arrival times and use a seismic velocity model to calculate each station's source takeoff direction. Measure the initial direct P wave polarity after checking instrument sign, component orientation, timing and phase identification. An upward trace at an instrument is not automatically a source compression without the propagation and component convention. Plot the inferred source radiation signs on a focal sphere, commonly using a lower-hemisphere equal-area projection, and seek nodal planes separating positive and negative observations.

For a shear dislocation with orthogonal unit slip and fault-normal vectors $\mathbf{s},\mathbf{n}$, the double-couple fault source has $M=M_0(\mathbf{s}{\mathbf{n}}^{\mathsf{T}}+\mathbf{n}{\mathbf{s}}^{\mathsf{T}})$. Its far-field P radiation sign is that ofThe two zeros are orthogonal nodal planes; the four quadrants alternate between outward and inward initial radiation. Their orientations and quadrant signs determine candidate strike, dip and rake. Normal, reverse and strike-slip mechanisms give different projected patterns. The earthquake pressure and tension axes are the negative and positive moment-tensor principal directions, respectively; they do not by themselves measure the complete regional stress tensor.

There is a fundamental fault-plane ambiguity: exchanging $\mathbf{s}$ and $\mathbf{n}$ leaves the radiation tensor unchanged. First motions therefore give a fault plane and an auxiliary plane without saying which actually slipped. Use mapped rupture, aftershock geometry or independent geodetic displacement to select the physical fault where possible. Polarity signs alone do not measure scalar seismic moment, source duration or rupture extent. A poorly sampled focal sphere or many near-nodal observations can admit several mechanisms rather than one defensible beachball.

Propagate uncertainty in event depth, velocity gradients and takeoff angles, allow for mispicked or reversed polarities, and report the range of mechanisms consistent with the observations. A primary method did this by testing alternative locations and velocity models and validated its uncertainty estimates using noisy synthetic data and clusters of events with similar waveforms.^[ref] Broadband waveforms and amplitude ratios add information, especially for non-double-couple sources, but they are additional measurements rather than a way to extract amplitude information from a sign-only data set.

Radial density structure. In an isotropic adiabatic material, the seismic parameter isFor a chemically homogeneous region whose background stratification is adiabatic, $d\rho /dP=\rho /K_S=1/\Phi$. Combine this with spherical hydrostatic balance and mass conservation to obtain the Adams-Williamson equationsThe inward density increase follows compressibility under the increasing pressure. A specified starting density and measured $\Phi (r)$ allow integration within a homogeneous adiabatic layer. Seismic velocities alone, however, do not determine that density normalization, since they contain ratios of moduli to density. Chemical boundaries, phase transitions, nonadiabatic thermal gradients and anisotropy require additional modeling; a composition jump is not obtained by blindly integrating the smooth relation across it.

Use travel times and surface waves to constrain elastic velocities, reflected/converted phases to locate interfaces, and normal-mode frequencies to constrain density as well as elastic structure through their sensitivity kernels. Impose independently known total mass and moment of inertia:Because these weight depth differently, they help distinguish admissible profiles, but two integrals alone do not specify an arbitrary radial function. Laboratory equations of state and mineral physics restrict composition, compressibility and density jumps. pressure and enclosed mass are matched across interfaces, while density may jump. In a liquid region $\beta =0$ and $\Phi ={\alpha }^2$, so the shear-wave absence constrains phase state without itself providing a direct density measurement.

The Preliminary Reference Earth Model is a concrete joint inversion of seismic travel times, normal modes, mass and inertia, including density and attenuation along with elastic structure. The original analysis combined these complementary data rather than deriving density uniquely from P wave speed.^[ref] A defensible radial density model combines hydrostatic/thermodynamic relations with seismic and global integral constraints, and states where homogeneity and adiabaticity fail.

For uniform density $\rho >0$ and stiffness, linear momentum balance without body force is $\rho {\ddot{u}}_i={\partial }_j{\sigma }_{ij}$. The minor elastic symmetries make the two symmetric-strain terms equal, so ${\partial }_j{\sigma }_{ij}=c_{ijkl}{u}_{k,lj}$. Write $s=ct-{\hat{p}}_j{x}_j$ and $u_i=a_{i}f(s)$. Then ${\ddot{u}}_i=c^2{a}_i{f}^{\prime \prime }$ and $u_{k,lj}=a_k{\hat{p}}_l{\hat{p}}_j{f}^{\prime \prime }$. Thus a nontrivial plane-wave profile satisfiesThe density-normalized elastic acoustic tensor is consequentlyThe alternative dummy-index form $B_{ik}={\rho }^{-1}{c}_{ijkl}{\hat{p}}_j{\hat{p}}_l$ is identical. In ordinary conservative elasticity, major symmetry $c_{ikjl}=c_{jlik}$ gives $B_{ji}=B_{ij}$ after interchanging the dummy indices. This symmetry comes from a strain-energy density; it is an assumption of the ordinary elastic tensor, not a consequence of homogeneity alone.

For any such eigenpair the profile $f$ is arbitrary, since both sides carry the same $f^{\prime \prime }$. This nondispersive result uses linear, local, time-independent elasticity without attenuation or spatial dispersion. A profile with $f^{\prime \prime }=0$ is a trivial uniform or affine displacement field and cannot be used by itself to infer the wave-speed eigenvalue condition.

Use geometric rays for body-wave arrivals and a traction-free upper surface. A compact isotropic explosion source has moment tensor proportional to ${\delta }_{ij}$ and produces direct P body radiation in a homogeneous isotropic medium. Thus the basic near-offset arrival is the direct upper-layer P wave, with time $X/{\alpha }_1$. A pure explosion is not a shear fault: it does not automatically emit a direct S wave of comparable strength. Interfaces can nevertheless convert its P energy into SV motion, and the free surface can support Rayleigh waves.

Downward P energy reflects from the interface to give a PP reflection returning to the receiver; for a horizontal interface its ray is symmetric and its time is $\sqrt{X^2+4H^2}/{\alpha }_1$. Other transmitted P energy enters the faster half-space. At $\sin i_c={\alpha }_1/{\alpha }_2$, a critical ray travels along the interface and feeds an upgoing seismic head wave, often labeled Pn in this simplified crustal model. Its ideal travel-time branch isIt requires enough offset for the two inclined upper-layer legs, $X\ge 2H\tan i_c$. Once $X$ exceeds the larger first-arrival crossover distance $X_c=2H\sqrt{({\alpha }_2+{\alpha }_1)/({\alpha }_2-{\alpha }_1)}$, this branch arrives before the direct upper-layer wave because it spends most of the path in faster material. A transmitted ray descending into an otherwise uniform infinite half-space does not return as a separate reflection without an additional reflector.

At oblique incidence, the bonded-interface conditions couple P and SV polarizations. A downgoing P ray can return as a converted SV reflection in the upper layer, with common horizontal slowness $p=\sin i_P/{\alpha }_1=\sin i_S/{\beta }_1$. There are transmitted P/SV branches and, after further interaction with the free surface or interface, mixed-mode reflections and layer reverberations. Conversion coefficients depend on density and impedance as well as wave speeds, and some vanish at normal incidence. The supplied inequalities do not order ${\beta }_2$ against ${\alpha }_1$, so one cannot universally draw a P-incident critical head ray travelling at the lower-layer S speed: that direct conversion construction requires ${\beta }_2>{\alpha }_1$. A shear head-wave phase can be present if the appropriate SV energy and ray conditions are available, but is not an inevitable direct explosion arrival from the stated data alone.

At later times a shallow explosion can excite dispersive Rayleigh-wave trains, with coupled P-SV motion decaying with depth. Short wavelengths mainly sample the slow upper layer, whereas long wavelengths penetrate the lower half-space; group and phase speeds need not be equal. Multiple reflected/converted body waves can overlap this coda. The shear-speed contrast also permits Love-wave modes structurally, but SH motion is decoupled from P-SV in this flat isotropic model. An ideal isotropic explosion does not excite Love waves here; appreciable SH/Love energy would require source asymmetry, scattering or a broken symmetry of the model.

Exact amplitudes, resolved arrival separation and attenuation cannot be calculated without densities, source spectrum, damping and instrument response. The qualitative robust features are a near-offset direct P arrival, reflected and converted phases, a faster-lower-layer P head wave that can become first at long range, reverberations and P-SV surface-wave energy. The ray diagram uses declared illustrative speed ratios solely to show the geometry; they are not additional supplied material constants.

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.

The positive centers in the supplied maps are examples of lunar mascons, especially associated with large impact basins. An empty excavation by itself removes mass and would not explain a positive central anomaly. Basin formation can thin the low-density lunar crust and raise denser lunar mantle beneath the floor; impact melting, collapse and later cooling can also create dense material and lock in a mass excess. Where a basin was later flooded, dense mare basalt adds another contribution. Mare infill is therefore one possible contributor, not a necessary explanation of every mascon.^[ref]

The difference between the maps is the modeled topographic contribution. On a common reference surface, write schematicallyThe free-air gravity anomaly retains relief-related attraction. Rims, excavated floors and central peaks therefore contribute fine structure. The Bouguer gravity anomaly removes the attraction predicted from the measured topography at an assumed crustal density, making internal density and interface structure more apparent. A low basin floor has a negative topographic contribution; subtracting it can enhance an underlying positive central anomaly. Removing a positive rim contribution can weaken a rim high. The supplied lower map is correspondingly smoother at short wavelengths and shows prominent residual basin centers and surrounding rings. This is a change of gravity reduction, not a physical change of the craters.^[ref]

For scale, a Bouguer slab approximation gives $\Delta g\simeq 2\pi G\Delta \rho d$. A $5$ km layer with geophysical density contrast $300\mathrm{kg}{\mathrm{m}}^{-3}$ produces about $63$ mGal above an ideal infinite slab. Actual basin geometry, burial, porosity and mantle uplift require a finite three-dimensional model. Smaller craters can have central lows, damaged porous interiors or more complex rings; the maps do not imply that every crater has a positive anomaly or that positive anomalies measure basalt thickness uniquely.

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.