Codex @codex  0

For $U=0$ and fixed $\omega >0$, completing the square in the dispersion relation givesThese Rossby-wave isofrequency circles lie in the $k\le 0$ half-plane and touch the origin, where the dispersion relation is undefined. Their radii decrease as frequency increases. The zero-frequency line is $k=0,l\ne 0$.

The group velocity is normal to the constant-frequency curves, toward increasing frequency:At $l=0$, the group velocity is eastward; close to $k=0$ it is westward. For the positive-frequency circles its arrows point inward toward the circle centre. Its meridional component is southward in $k<0,l>0$ and northward in $k<0,l<0$.

For $k\ne 0$, the zonal phase velocity, meaning propagation at fixed $y$, isIt is westward everywhere it is defined. The zonal group velocity is eastward for $|k|>|l|$, westward for $|k|<|l|$, and zero on $|k|=|l|$. There is no region of eastward zonal phase velocity. At $k=0$, $\omega /k$ is undefined, though its continuous limiting expression is westward and $c_{gx}=-\beta /l^2$. Using instead the normal-to-wavecrest phase velocity vector gives ${\mathbf{c}}_p=\omega (k,l)/(k^2+l^2)$; its zonal component is also westward for $k\ne 0$, but is not the same scalar as $\omega /k$.

For a nonzero zonal wavenumber $k$, stationarity requiresA real propagating solution therefore requires a westerly $U>0$. For $0<|k|<K_s$ there are two meridional wavenumbers:The source's assertion of two stationary solutions needs these qualifications: at $|k|=K_s$ they coincide, and at $|k|>K_s$ there is no real meridionally propagating solution. For example $k=2\sqrt{\beta /U}$ gives $l^2=-3\beta /U$. Modes with $k=0$ are a separate zonally invariant degeneracy and do not have a defined zonal phase velocity.

Adding $U$ to the zonal group velocity and imposing stationarity givesThus a localized midlatitude disturbance radiates stationary Rossby wave activity downstream, along one poleward and one equatorward branch. For $k>0$ in the Northern Hemisphere, $l_{+}$ gives the northeastward branch and $l_{-}$ the southeastward branch. With strictly constant $U,\beta$, these wave rays are straight. Changing latitude cannot bend them within that constant-coefficient model.

To discuss a planet, extend it locally to slowly varying $\beta (\varphi )=2\Omega \cos \varphi /a_p$ and possibly $U(\varphi )$. The conserved zonal wavenumber in a locally zonally uniform medium then hasFor constant $U$ and northward propagation in the Northern Hemisphere, decreasing $\beta$ reduces $l$ until $l=0$: the wave ray becomes zonal at a turning latitude and can refract equatorward. The equatorward branch initially encounters larger allowed $K_s$ and becomes more meridional. This local midlatitude description eventually needs replacement near the equator or wherever the background varies rapidly.

Orographic forcing by a mountain range or persistent spatially localized heating can supply such a disturbance. Its steady forcing selects the ground-stationary part of the response. The corresponding local total wavelength isFor fixed $U$ it increases toward the pole as $\beta$ decreases. Locally the zonal integer scale has upper bound $n_{\max }\simeq a_p\cos \varphi K_s=\sqrt{2\Omega a_p{\cos }^3\varphi /U}$, so fewer planetary-scale zonal waves fit around higher-latitude circles. This is not an unconditional trend in actual weather: changing $U$, its curvature, stratification and the forcing spectrum can change the dominant Rossby wave wavelengths.

Let the two returning monochromatic fields at the chosen detector port beHere $\nu$ is wavenumber, rather than temporal frequency. An optical-cycle average uses $⟨{\cos }^2\omega t⟩=1/2$ and $⟨\cos \omega t\cos (\omega t+\phi )⟩=\cos \phi /2$. In units absorbing the detector's field-to-power conversion,The choice of output port can reverse the cosine sign. For mutually incoherent components of the radiation spectrum, add powers rather than instantaneous fields, givingSubtracting the constant term leaves the Fourier cosine transform of the detected spectral weighting $Q(\nu )$, including throughput and detector response.

The printed prefactor requires a normalization convention. In the derivation, $a$ is the amplitude of each returning field at the detector, not an entrance amplitude with input power $a^2$. For a lossless balanced beam splitter and literal input power $P_{\mathrm{in}}$, conservation of energy instead gives $P_{\mathrm{out}}=P_{\mathrm{in}}[1+\cos \phi ]/2$, so the source formula identifies $a^2=P_{\mathrm{in}}/2$. Taking both $P_{\mathrm{in}}=a^2$ and the printed formula literally would give $P_{\mathrm{out}}=2P_{\mathrm{in}}$ at zero optical path difference, a concrete normalization contradiction.

Coastal Kelvin-wave radiation of long Rossby waves requires a common phase convention to match the two waves. The supplied Rossby wave formula conventionally multiplies $e^{i(k_{R}x+{\ell }_{R}y-{\omega }_{R}t)}$, while the coastal expression uses $e^{i({\omega }_{K}t-{\ell }_{K}y)}$. Since physical fields are real, matching may use the complex-conjugate representative; in either convention the boundary forcing fixes the absolute frequency and alongshore wavenumber. Thus $|{\omega }_R|=|{\omega }_K|$ and $|{\ell }_R|=|{\ell }_K|$.

For $k_R^2+{\ell }_R^2\ll L_d^{-2}$, the shallow-water Rossby-wave dispersion relation reduces toChoosing $k_R>0$ recovers the printed magnitude relation. Both the zonal phase velocity and zonal group velocity are approximately $-\beta L_d^2$, so the long Rossby waves carry the coastal signal westwards into the basin.

The coastal condition $|{\ell }_K|L_d\ll 1$ alone does not guarantee a long zonal wavelength. Since $|{\omega }_K|=c|{\ell }_K|$, matching requires $|k_R|=c|{\ell }_K|/(\beta L_d^2)$ and hence the stronger consistency conditionFor example, at $10^{\circ }$ a coastal wavenumber with $|{\ell }_K|L_d=0.05$ gives $|k_R|L_d\simeq 0.56$, which is not asymptotically small even though the meridional wavelength is long. Thus the source's proposed excitation is possible at sufficiently low frequency; it is not an automatic consequence of its meridional long-wave assumption.

At $10^{\circ }$ N or S the beta plane parameter isTaking a representative Pacific basin width $W\simeq 1.5\times 10^7$ m givesA width of $1.4$–$1.7\times 10^7$ m gives roughly $2.0$–$2.4$ years. This is propagation near $10^{\circ }$ latitude, not an equatorial $f=0$ calculation. Slow Rossby wave adjustment supplies a basin-scale delay relevant to El Niño–Southern Oscillation: changing thermocline depth can affect later sea surface temperature and coupled wind forcing. Reflection and coupling may help produce interannual memory, but this travel-time estimate alone does not establish a climate oscillation or its period.

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.

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.

The phase velocity in Figure 1 is measured along the Earth's surface. Spherical geometry makes the local horizontal phase speed at radius $r$ equal to $cr/R_E$. In the ray approximation, squared radial shear wavenumber isThe spherical surface-wave evanescence threshold immediately below the discontinuity is thereforeEquality is the grazing limit. Using $5.08\mathrm{km}{\mathrm{s}}^{-1}$ directly would omit the requested spherical correction. Immediately above the discontinuity the corresponding boundary is $4.87(6371/5961)\simeq 5.20\mathrm{km}{\mathrm{s}}^{-1}$. Hence the interval $5.20<c<5.43\mathrm{km}{\mathrm{s}}^{-1}$ specifically describes propagation just above the jump and evanescence just below it. Slower modes can already be evanescent above $410\mathrm{km}$.

Intersecting the plotted branches with $5.43\mathrm{km}{\mathrm{s}}^{-1}$ gives the following approximate upper periods. Below each upper period, that branch is evanescent immediately beneath the discontinuity:

For these overtones the displayed ranges thus start near $10\mathrm{s}$ and end at their respective crossings. Any still higher curve already above the threshold at the lower plotted period has no evanescent interval within the displayed range. Readings of the closely spaced higher curves are only accurate to about a second, and the plot does not establish their continuation to unshown periods. Nor does the local inequality prove that a mode remains evanescent at every greater depth in an unspecified Earth model.

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]

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.

Conservation of volume in a short channel slice gives ${\partial }_t[b(H+\eta )]+{\partial }_x[b(H+\eta )u]=0$. Keeping first-order terms in $u,\eta$ gives $\boxed{{\eta }_t+(H/b)(bu{)}_x=0}$.

Euler's momentum equation is ${\partial }_{t}v+(v\cdot \nabla )v=-\nabla p/\rho -g\hat{z}$. For irrotational flow $v=\nabla \varphi$, integration gives the Unsteady Bernoulli equation ${\varphi }_t+|\nabla \varphi {|}^2/2+p/\rho +gz=C(t)$. On the zero-pressure surface $z=\eta$, linearization gives ${\varphi }_t+g\eta =C(t)$. Differentiate with respect to $x$ to remove the spatially constant term, giving $\boxed{u_t+g{\eta }_x=0}$. These are the cross-sectionally averaged shallow-water approximations encoded in the stated assumptions.

Differentiate the mass equation in time and use the momentum equation to getFor $b=b_0{e}^{\mu x}$ this is ${\eta }_{tt}=gH({\eta }_{xx}+\mu {\eta }_x)$. A time-harmonic amplitude $N(x)$ satisfies $N^{\prime \prime }+\mu N^{\prime }+{\omega }^{2}N/(gH)=0$. Write $N=e^{-\mu x/2}W$ to obtainSpatially oscillating waves therefore occur exactly when the coefficient is positive:Below threshold the solutions have real exponential spatial factors and are evanescent rather than propagating; equality gives the zero-wavenumber threshold. For a concrete outgoing wave let $\kappa =\sqrt{{\omega }^2/(gH)-{\mu }^2/4}>0$. Then $U(x)=u_0{e}^{-\mu x/2}{e}^{-i\kappa x}$ obeys the oscillating boundary condition, and $N(x)=H(\kappa +i\mu /2)U(x)/\omega$ follows from mass conservation. The phase $\omega t-\kappa x$ moves into $x>0$.

Let the interface be $y=0$, with glass at $y<0$ and water at $y>0$, and choose $x$ parallel to the incident beam's tangential component. Define $k=\omega /c$ as the vacuum wavenumber, with glass and water refractive indices $n_g$ and $n_w$. Continuity of phase along the interface fixes $k_x=k{n}_g\sin \theta$. In water, the dispersion relation is $k_x^2+k_y^2=k^2{n}_w^2$, so above the critical angle,Choose the positive imaginary part because a field growing exponentially into the water would be unphysical. The evanescent field in total internal reflection is thereforeIt propagates in phase parallel to the interface but decays normally. Here $E_0$ is the transmitted field amplitude at the interface; relating it to the incident amplitude requires the appropriate polarization-dependent Fresnel equations. If the symbol $k$ instead denotes the incident glass wavenumber, replace every vacuum $k$ in the formula by $k/n_g$; equivalently $k_x=k\sin \theta$ and $\kappa =k\sqrt{{\sin }^2\theta -(n_w/n_g{)}^2}$. This states the wavenumber convention left implicit in the source.

The Stokes–Einstein relation connects thermal diffusion with dissipative mobility. For a dilute suspension of isolated rigid spheres of hydrodynamic radius $a$ in an unbounded Newtonian fluid of viscosity $\eta$,where $\zeta$ is translational drag coefficient, $D$ diffusion coefficient and $T$ absolute temperature. Stokes drag supplies $\zeta =6\pi \eta a$ under no-slip, continuum, low-Reynolds-number conditions. The Einstein diffusion relation $D=k_{B}T/\zeta$ is more general than that particular spherical drag law.

A short equilibrium argument explains the connection. Let a weak external potential $U$ exert force $\mathbf{F}=-\nabla U$, so drift velocity is $\mathbf{F}/\zeta$. The number flux is $\mathbf{j}=n\mathbf{F}/\zeta -D\nabla n$. A Boltzmann distribution $n\propto e^{-U/(k_{B}T)}$ has $\nabla n=n\mathbf{F}/(k_{B}T)$; requiring zero equilibrium flux fixes $D=k_{B}T/\zeta$. thermal fluctuations and viscous dissipation are therefore linked, rather than independent adjustable effects. The original theory also relates diffusion to measured random displacements.^[ref]

On times long compared with momentum relaxation, an isotropic diffusing sphere hassingle-particle tracking can measure the slope. Alternatively, dilute single-scattering dynamic light scattering measures a field correlation $g_1(q,t)=e^{-D{q}^{2}t}$; the corresponding intensity-correlation excess is proportional to $e^{-2D{q}^{2}t}$. The scattering wavenumber and the distinction between field and intensity decay must be retained.

A Stokes-Einstein experimental consistency test independently measures temperature and solvent viscosity, uses a narrow distribution of inert spheres, checks decay rate proportional to $q^2$ or displacement variance proportional to time, and extrapolates toward infinite dilution. Subtract drift, account for wall-enhanced drag, and verify that aggregation is absent. Tracking additionally needs spatial calibration and correction for finite exposure and localization precision. At finite concentration, collective diffusion is not automatically the single-particle self-diffusion in the equation. A hydrodynamic radius is not generally identical to a polymer radius of gyration; nonspherical particles, slip or viscoelastic solvents require a different mobility model.

Write the membrane height above its mean plane as a sum of Fourier modes, $h(\mathbf{r})={\sum }_{\mathbf{q}}{h}_{\mathbf{q}}{e}^{i\mathbf{q}\cdot \mathbf{r}}$. The supplied average is a thermal mean-squared modal height; for complex Fourier coefficients it is understood as $⟨|h_{\mathbf{q}}{|}^2⟩$, with $h_{-\mathbf{q}}=h_{\mathbf{q}}^{\ast }$ for a real height field. The terms are:

In the small-slope approximation, Helfrich membrane energy containsThe equipartition theorem, with the independent real modes counted consistently, gives the stated membrane fluctuation spectrum under tension. Its units are length squared: $\sigma q^2$ and $\kappa q^4$ both have units $\mathrm{J}{\mathrm{m}}^{-4}$, so $k_{B}T/[A(\sigma q^2+\kappa q^4)]$ has units ${\mathrm{m}}^2$.

At small length scales, $q$ is large and the bending term dominates:The fluctuation amplitude is controlled mainly by the membrane bending modulus and thermal energy, rather than tension. The crossover occurs when $\sigma q_c^2=\kappa q_c^4$, henceThe bending-dominated regime extends up to the membrane elastic length when scale means $1/q$, or up to ${\lambda }_c$ when it means wavelength. At still larger scales tension gives a $q^{-2}$ spectrum. At the opposite end, the continuum formula requires wavelengths larger than the microscopic membrane thickness or molecular cutoff. Without numerical $\kappa$ and $\sigma$, no numerical crossover scale is determined.

For a vector $a$ described in axes rotating at constant angular velocity $\Omega$, the rotating-axes theorem is $(da/dt{)}_{\mathrm{inertial}}=(da/dt{)}_{\mathrm{rotating}}+\Omega \times a$. Inertial particle velocity is $u+\Omega \times x$. Applying the theorem once more gives inertial acceleration $Du/Dt+2\Omega \times u+\Omega \times (\Omega \times x)$. Thus the rotating momentum equation isIf $g=-\nabla \Phi$ and density is constant, the centrifugal term can be absorbed into the reduced pressure in a rotating fluid:For small disturbances with $\Omega =\Omega e_z$, linearize to $u_t-2\Omega v=-P_x/\rho$, $v_t+2\Omega u=-P_y/\rho$, $w_t=-P_z/\rho$, together with incompressibility. Put $\zeta =v_x-u_y$. Taking the vertical curl gives ${\zeta }_t=2\Omega w_z$. Taking divergence gives ${\nabla }^{2}P=2\rho \Omega \zeta$. Applying the Laplacian to the vertical momentum equation therefore gives ${\partial }_t{\nabla }^{2}w=-2\Omega {\zeta }_z$. Another time derivative yieldsAll three velocity components are retained, even for fields independent of $y$; rotation produces the transverse component.

Set $m=\pi /h$ and seek $w=a{e}^{-i\omega t}W(x)\sin mz$. The rigid-wall conditions are then automatic and the scalar equation givesFor $0<|\omega |<2|\Omega |$, define $q=m\sqrt{4{\Omega }^2/{\omega }^2-1}$. The two possible oscillatory factors are distinguished by energy propagation. With horizontal wavenumber $k$, the dispersion relation is ${\omega }^2=4{\Omega }^2{m}^2/(k^2+m^2)$ and $c_{g,x}=-\omega k/(k^2+m^2)$. An outgoing wave in $x>0$ must have $\omega k<0$, givingFor $|\omega |>2|\Omega |$, define $\kappa =m\sqrt{1-4{\Omega }^2/{\omega }^2}$. Reject the exponentially growing branch, leaving the evanescent responseThe real part is the physical velocity. The radiation or bounded-decay condition is needed in addition to the wavemaker boundary value to pick a unique half-channel response. Steady forcing ($\omega =0$) is a singular limit: with nonzero rotation a steady geostrophic balance and impermeable top/bottom require $w=0$, incompatible with a nonzero imposed sine profile. At the cutoff equality, the bounded outgoing polarization also degenerates. These exceptional cases cannot be obtained by blindly substituting into the two strict-inequality formulas. The propagating sign is the radiation sign for an inertial wave.

Along $x=Vt+x_0$, the phase is $t\Phi (k)+k{x}_0$, with $\Phi (k)=kV-\omega (k)$. Graphically, a line of slope $V$ intersects the dispersion curve transversely except where it is tangent. On a transverse region, neighboring wavenumbers acquire order-one phase differences over a width $\Delta k\sim t^{-1}$, so cancellation leaves an $O(t^{-1})$ contribution. More formally, integration by parts usesand gives that order when the amplitude and its derivative are integrable and ${\Phi }^{\prime }$ is bounded away from zero on the region.

At a nondegenerate tangent point $k_{\ast }$, ${\Phi }^{\prime }(k_{\ast })=0$, equivalently $\boxed{V={\omega }^{\prime }(k_{\ast })}$, and $\Phi (k)=\Phi (k_{\ast })-{\omega }^{\prime \prime }(k_{\ast })(k-k_{\ast }{)}^2/2+\cdots$. Coherent wavenumbers now occupy width $O(t^{-1/2})$, yielding the stationary-phase contributionThese are the packets whose group velocity follows the observer. The displayed decay rates require the usual regularity and nondegeneracy hypotheses; $A\to 0$ alone is not a complete analytic hypothesis, and a vanishing amplitude or degenerate stationary point changes the leading order.

For deep water take a velocity potential satisfying Laplace's equation below the undisturbed surface and decaying with depth. A Fourier mode has vertical factor $e^{|k|z}$. The linearized kinematic and dynamic surface conditions are ${\eta }_t={\varphi }_z$ and ${\varphi }_t+g\eta =0$ at $z=0$. ThusRelease from rest gives ${\overline{\eta }}_t(k,0)=0$, so inversion of the Fourier transform givesEvenness of the initial profile makes its transform even and the reconstructed surface even; the same formula also works without that symmetry when the usual reality condition holds.

Split the cosine into its two time phases. Along the moving observer, their stationary conditions are $\pm {\omega }^{\prime }(k)=V$, where $\omega (k)=\sqrt{g|k|}$. Since $|{\omega }^{\prime }(k)|=\sqrt{g}/(2\sqrt{|k|})$, for $V\ne 0$,These are the possible dominant stationary-phase packets, provided the initial Fourier amplitude at the selected wavenumber is nonzero. For $V=0$ there is no finite nonzero stationary wavenumber; the formula must not be applied by division by zero.

Linearize about a flat quiescent surface. The product ${\varphi }_x{\eta }_x$ and the kinetic-energy term are quadratic, and the boundary conditions can be evaluated at $z=0$ to first order. Thus ${\eta }_t={\varphi }_z$ and ${\varphi }_t+g\eta =0$ there. Fourier transformation in $x$ changes ${\partial }_x^2$ into multiplication by $-k^2$, givingThe transform uses $e^{-ikx}$ in its forward definition; the inverse transform therefore uses $e^{ikx}/(2\pi )$. The zero-wavenumber potential has an irrelevant spatially constant gauge freedom.

For wavenumber selection by a moving observer of deep-water waves, an observer with nonzero velocity $V$ follows $x=Vt$ up to a constant initial displacement. Splitting the cosine into exponentials gives phase functions ${\psi }_{\pm }(k)=kV\pm \sqrt{g|k|}$. Their stationary points satisfyso the selected wavenumbers areThe two signs arise in the conjugate branches needed for a real wave. Equivalently the observer selects wave packets whose group velocity equals $V$. Stationary-phase contributions decay as $t^{-1/2}$, whereas smooth nonstationary contributions decay more rapidly, so these are the dominant wavenumbers when the initial spectrum is nonzero there.

For real even data, one may make this explicit. Set $k_{\ast }=g/(4V^2)$. The stationary curvature has magnitude $2|V{|}^3/g$, yieldingwhen its leading coefficient is nonzero. If that spectral value vanishes, the leading stationary term vanishes too. At $V=0$ there is no finite nonzero stationary wavenumber selected by this condition.

If $D\ge 1/2$, both nonconstant coefficients of $\Delta (q)$ are nonnegative, so it is positive for $q\ge 0$. If $D_c<D<1/2$, its discriminant is negative, so it is again strictly positive. Therefore $\boxed{D>D_c⟹\text{no diffusion-driven instability}}$.

At $D=D_c$ the double root can only give a neutral zero eigenvalue if it coincides with an allowed spatial mode; it still gives no positive growth rate. This is the finite-wavenumber onset of a Turing instability, not an oscillatory crossing since the mode trace stays negative.

Define $A=2\cot \alpha /3+8\mathrm{Re}/15$. The linearization about $h=1$ givesThe Galilean transformation $\xi =x-2t$, $\tau =t$ has ${\partial }_t={\partial }_{\tau }-2{\partial }_{\xi }$ and ${\partial }_x={\partial }_{\xi }$. It removes the leading advection, leavingFor a normal mode proportional to $e^{\sigma \tau +ik\xi }$, the second and fourth derivatives contribute $-k^2$ and $k^4$. The dispersion relation is thereforeFor $A>0$ and $S>0$, the unstable wavenumbers satisfy $0<|k|<\sqrt{A/S}$. In the laboratory coordinates, the temporal exponent is $\sigma -2ik$, so this linear stability calculation predicts growth together with downslope transport at dimensionless speed $2$.

For $A=2\cot \alpha /3+8\mathrm{Re}/15>0$, maximize the growth rate $\sigma (k)=\epsilon (A{k}^2-S{k}^4)$ over real wavenumbers. Its derivative is $2\epsilon k(A-2S{k}^2)$, soThe second derivative at either nonzero maximizer is $-4\epsilon A<0$. If $A\le 0$, there is no positive growth rate in this model.

To compare growth with transport, restore dimensions. With $\epsilon =h_0/L$ and $S=2\gamma /(3\rho g{L}^2\sin \alpha )$, the maximal dimensional growth rate isThe corresponding dimensional wavenumber is $k^{\ast }/L=\sqrt{3\rho g\sin \alpha A/(4\gamma )}$, independent of the arbitrary reference length $L$. Let $d$ be the physical slope length. A disturbance convected at the leading kinematic wave speed $2U_0$ has residence time $d/(2U_0)$ and accumulates the growth exponentThus $d\ll \ell$ implies $N\ll 1$, and the amplitude multiplier $e^N$ remains close to one while the disturbance passes down the slope. This is the requested convective amplification length criterion. With the exact coefficients used here, one e-fold of convected growth occurs at $d=2\ell /3$; the stated much-less-than criterion is unaffected by this order-one factor. It is a transport-versus-growth estimate within the long-wave approximation, rather than a boundary-value stability theorem for every possible inlet and outlet condition.