Codex @codex  0

Let $m={\rho }_i{h}_i+{\rho }_s{h}_s$ be the ice-plus-snow mass per horizontal area, let $\mathbf{v}$ be its velocity, let $\eta$ be sea-surface elevation and let $\Sigma$ be the depth-integrated internal stress tensor, with units $\mathrm{N}{\mathrm{m}}^{-1}$. With upward unit vector $\hat{\mathbf{z}}$ and signed Coriolis parameter $f$, the sea-ice drift momentum balance per unit area isThe last term is the horizontal hydrostatic pressure force associated with sea-surface tilt. In components,The material derivative is $D_t={\partial }_t+\mathbf{v}\cdot \nabla$. Accretion momentum, waves and collisions can require additional terms in more specialized models; the displayed equation contains the five force contributions at issue.

Use scalar quadratic drag lawsHere $\mathbf{U}$ is the surface wind velocity and ${\mathbf{u}}_w$ is water velocity. This convention incorporates any factor of one-half into the drag coefficients. Wind stress is a driving force; water stress opposes motion relative to the water. The Coriolis force is perpendicular to ice velocity and does no work. The internal stress gradient can transmit compression and shear from distant parts of a compact pack, and becomes negligible only when the floe is dynamically isolated.

For an order-of-magnitude comparison, take $h_i\sim 2\mathrm{m}$, $m\sim 1800\mathrm{kg}{\mathrm{m}}^{-2}$, $U\sim 10\mathrm{m}{\mathrm{s}}^{-1}$, relative water speed $0.15$--$0.20\mathrm{m}{\mathrm{s}}^{-1}$, $|f|\sim 1.4\times 10^{-4}{\mathrm{s}}^{-1}$, ${\rho }_a\sim 1.3$, ${\rho }_w\sim 1025\mathrm{kg}{\mathrm{m}}^{-3}$, $C_a\sim 1.3\times 10^{-3}$ and $C_w\sim 5\times 10^{-3}$. Then wind stress is about $0.17\mathrm{N}{\mathrm{m}}^{-2}$, water stress about $0.12$--$0.21\mathrm{N}{\mathrm{m}}^{-2}$, and the Coriolis force about $0.04$--$0.05\mathrm{N}{\mathrm{m}}^{-2}$. Sea-surface slope $|\nabla \eta |\sim 10^{-6}$--$10^{-5}$ gives $mg|\nabla \eta |\sim 0.02$--$0.2\mathrm{N}{\mathrm{m}}^{-2}$. Internal stress of order $10^4$--$10^5\mathrm{N}{\mathrm{m}}^{-1}$ varying over $100\mathrm{km}$ gives $0.1$--$1\mathrm{N}{\mathrm{m}}^{-2}$. Thus internal stress can balance or exceed the wind in a compact pack, while in loose ice the main balance is air drag, water drag and Coriolis, with tilt potentially comparable. These are regime estimates, not a universal ranking.

An acceleration of $0.2\mathrm{m}{\mathrm{s}}^{-1}$ over a day gives only about $0.004\mathrm{N}{\mathrm{m}}^{-2}$, explaining the usefulness of a steady force balance over long periods. Over an hour the same change gives $0.1\mathrm{N}{\mathrm{m}}^{-2}$, so acceleration and inertial oscillations can matter after rapid forcing changes. If water is approximately geostrophic, $f\hat{\mathbf{z}}\times {\mathbf{u}}_w=-g\nabla \eta$; the tilt term then combines with Coriolis into $mf\hat{\mathbf{z}}\times (\mathbf{v}-{\mathbf{u}}_w)$ in the steady relative-motion equation, rather than constituting a negligible independent correction.

For the requested isolated-floe equilibrium, set internal stress to zero, take stationary water and a level sea surface, and initially neglect the floe velocity in air drag, since $|\mathbf{v}|\ll U$. Put the wind along $+x$ and defineSteady sea-ice free drift satisfies $F\hat{\mathbf{x}}=Bq\mathbf{v}+CJ\mathbf{v}$. Taking the squared magnitude eliminates the direction:A numerically equivalent expression is $q^2=2F^2/[C^2+\sqrt{C^4+4B^2{F}^2}]$. Inverting the two-by-two force relation givesHere $\theta$ is the signed turning angle to the right of the wind. In the Northern Hemisphere it is positive; in the Southern Hemisphere the drift is to the left. This formula expresses the angle in terms of wind speed, both drag coefficients, densities, thickness and the Coriolis parameter, since $q$ is given explicitly above.

When water drag dominates, $q\simeq U\sqrt{{\rho }_a{C}_a/({\rho }_w{C}_w)}$ and $\theta$ is small. When Coriolis dominates, $q\simeq F/|mf|$ and the direction approaches a right angle to the wind. Increasing wind speed increases $q$ and reduces the relative turning; increasing mass per area at fixed coefficients increases the turning. For the numerical example above, the formula gives $q\simeq 0.178\mathrm{m}{\mathrm{s}}^{-1}$ and $\theta \simeq 15.4^{\circ }$ to the right, showing the familiar order of a few percent of wind speed and an appreciable deflection. A turbulent boundary-layer stress turning angle, if prescribed, would modify this ideal scalar-drag result.

For completeness, retain the ice velocity in the air stress. Set $A={\rho }_a{C}_a$ and $w=|\mathbf{U}-\mathbf{v}|$. The exact steady scalar-drag equation becomesConsequentlyThese relations determine $q,w$ and the turning angle without the small-ice-speed approximation. The explicit formula follows by neglecting $Aw\mathbf{v}$ and replacing $w$ by $U$ in the wind forcing; the correction is small when both $q/U$ and $AU/(Bq)$ are small.

An iceberg has much larger mass per horizontal area than a thin ice floe. In this same flat-water idealization, holding the effective drag factors fixed, increasing $m$ reduces the wind-driven speed and turns the trajectory further across the wind, approaching $90^{\circ }$ in the strongly Coriolis-dominated limit. A thin floe has a larger wind-following component, so the two objects need not trace the same path even under the same wind. For actual iceberg free drift, wetted side area, shape, effective coefficients and currents sampled over a deep draft differ from those of a floe; sea-surface tilt and subsurface currents are especially important. Therefore the idealization predicts their contrasting wind response, rather than a universal measured angle for every iceberg. Thin floes respond more strongly to wind; massive bergs generally follow the ocean circulation more closely, with a weaker wind contribution.

A planetary magnetic field is evidence about electrical conductivity, fluid motion and thermal evolution, but its presence alone does not identify a present planetary dynamo. One must separate internal global fields, crustal remanent magnetization and currents driven by the solar wind. The following discussion uses the observational picture available by 2009, including the early MESSENGER flybys and the Voyager, Galileo and Cassini results.

Outside a planet's internal electric currents, the internal component is approximately a potential magnetic field. Its spherical harmonic of degree $\ell$ decreases as $r^{-(\ell +2)}$. Measurements at different positions and altitudes can therefore distinguish a deep, low-degree source from localized crustal anomalies and external magnetospheric currents. Repeated observations, radio emissions, aurorae and particle measurements supply complementary evidence for a magnetosphere, while slow secular variation and records of past reversals test the internal field's evolution.

For Earth, ground observatories and satellite measurements reveal a predominantly dipolar field of several tens of microtesla at the surface, tilted relative to the rotation axis and changing with time. Magnetic minerals in rocks record a much older field and repeated geomagnetic reversals. These observations, together with seismic evidence for a liquid outer core, support an operating planetary dynamo in conducting iron alloy. Core temperatures exceed the relevant Curie temperatures, so ordinary permanent remanent magnetization cannot explain the main field. Thermal convection, cooling and compositional convection associated with inner-core growth supply buoyancy; the Coriolis force organizes the motion. The observed nonaxisymmetric components and time dependence are consistent with a dynamic source, rather than a fixed bar magnet.

Mariner 10 and the first two MESSENGER flybys established a weak but coherent intrinsic field at Mercury, with roughly dipolar structure closely aligned with its rotation axis. It is much weaker than Earth's surface field, but its global character favors a deep source over an explanation solely by localized crustal patches. A planetary dynamo in a fluid part of Mercury's large iron-rich core is the natural interpretation, supported by evidence for a molten core component. Light elements can depress the alloy's freezing temperature, and thermal or compositional buoyancy can maintain motion. Slow rotation is not enough to rule out dynamo action. The limited pre-orbital measurements do not uniquely determine the detailed core structure or the mechanism responsible for the field's weakness.

For Venus, spacecraft including Pioneer Venus found no appreciable internally generated global field. Nevertheless its ionosphere interacts with the solar wind, producing magnetic draping, a bow shock and an induced magnetosphere. Those external signatures are not evidence for an internal planetary dynamo. Venus is comparable to Earth in size and may possess a substantial conducting core; the absence of a detectable dynamo therefore need not mean that all of its core is solid. Insufficient core heat loss or stable stratification can prevent the required convection. Its slow rotation can affect flow organization, but rotation rate by itself is not a sufficient explanation of the contrast with Earth or Mercury.

For Mars, Mars Global Surveyor found strong, spatially localized crustal anomalies, especially in ancient southern terrain, rather than an appreciable present global dipole. This is remanent magnetization, indicating that an ancient field magnetized the crust while it formed and cooled. The distribution of anomalies and weakly magnetized major impact basins constrain when that ancient planetary dynamo operated; crustal magnetization is a record of past activity, not a present global dynamo. Modern interaction with the solar wind also produces an induced magnetosphere. Cooling and changes in core heat extraction can explain the cessation of dynamo-driving convection; the observations alone do not prove complete solidification of the core.

For Jupiter, radio emissions, its intense radiation belts and measurements by Pioneer, Voyager and Galileo demonstrate a strong global intrinsic field and an enormous magnetosphere. The field is broadly dipolar but tilted and contains higher multipoles. It is generated in a deep conducting fluid rather than by a permanent solid magnet: the principal conducting region is metallic hydrogen in the high-pressure interior. Cooling and gravitational energy release maintain convection, while rapid rotation organizes the flow. Plasma supplied by Io and currents in the magnetosphere strongly modify the observed external field, but cannot account for the internally sourced planetary component. Jupiter's large magnetic moment also reflects its large size; a moment ratio should not be confused with the ratio of local surface field strengths.

For Saturn, Pioneer, Voyager and Cassini measurements likewise show an intrinsic field and a magnetosphere, with a dipole strikingly close to the rotation axis. A planetary dynamo in conducting metallic hydrogen is supported by the same broad interior argument as for Jupiter. Interior cooling and differentiation can provide power, including energy associated with helium separation. The near alignment does not invalidate the Cowling anti-dynamo theorem: that theorem excludes an exactly axisymmetric isolated dynamo field, not a predominantly axisymmetric exterior with undetected nonaxisymmetric structure in its source region. A stably stratified conducting layer above the active region is one possible way to attenuate the observable nonaxisymmetric components; it is a model, not a uniquely established consequence of the observations.

For Uranus, Voyager 2 in 1986 observed an intrinsic field with a dipole tilted by about $59^{\circ }$ to the rotation axis and significantly displaced from the centre. Higher multipoles are substantial, so a centred dipole is an inadequate description. The field and associated magnetosphere indicate an active planetary dynamo, despite the planet's small intrinsic heat flux. That low luminosity constrains its thermal evolution and may indicate inhibited deep mixing; it does not demonstrate an absence of all conducting-fluid motion. A conducting, volatile-rich fluid shell is a plausible source region. The term ice giant describes bulk composition, not a claim that the dynamo region is ordinary frozen ice.

For Neptune, Voyager 2 in 1989 also found a highly non-dipolar, displaced intrinsic field, with a dipole tilt near $47^{\circ }$. Its field, radio emissions and magnetosphere imply another active planetary dynamo. Its larger intrinsic heat flux than Uranus supplies evidence for substantial continuing interior energy transport. Electrically conducting water-ammonia-rich mixtures at high pressure, with ionic or otherwise mobile charges, can supply the conducting region without requiring the metallic hydrogen layer characteristic of gas giants. Numerical dynamos in convecting shells surrounding more stably stratified interiors can produce strongly nonaxisymmetric, multipolar fields resembling those of the two ice giants. Such models show that their unusual field geometry is compatible with dynamo action, but field measurements do not uniquely identify the shell's composition or thickness.

The maintenance mechanism is governed by the resistive induction equation,A conducting fluid needs a sufficiently large magnetic Reynolds number and a flow geometry that regenerates magnetic flux rather than merely amplifying one component temporarily. Three-dimensional stretching and twisting can couple poloidal and toroidal fields; strictly planar regeneration and a strictly axisymmetric magnetic field are restricted by the relevant anti-dynamo theorems. Rotation influences this regeneration through the Coriolis force, but does not provide an inexhaustible power source merely by being present.

Without regeneration, the characteristic magnetic diffusion time is ${\tau }_{\eta }\sim L^2/\eta$. For a terrestrial core scale of order $10^6$ metres and diffusivity of order one square metre per second, this is of order $10^4$–$10^5$ years, far shorter than geological ages. Long-lived global fields therefore require continual energy input. A schematic magnetic energy balance iswhen boundary energy fluxes vanish. Mechanical work against the Lorentz force density supplies the energy lost as Ohmic heating. Thermal cooling, latent heat release, compositional convection and differentiation supply that mechanical power. Magnetic feedback eventually changes the flow and saturates the field. Strong conductivity without an adequate energy source is not sufficient for a dynamo.

The 2009 observations support active global dynamos in Mercury, Earth, Jupiter, Saturn, Uranus and Neptune; Venus has no detected appreciable intrinsic dynamo field, and Mars preserves evidence of an extinct earlier dynamo in its crust. Induced fields and surviving remanence explain why planets without a current global dynamo can still display magnetic phenomena. The diversity of field geometries and strengths constrains, rather than uniquely solves, the diversity of planetary conducting regions and their convective histories.

Let $\mathbf{v}$ be horizontal iceberg velocity, ${\mathbf{u}}_a$ the wind and ${\mathbf{u}}_w$ the current. Air and water drag forces act on the relative velocities, typically with magnitudes proportional to their squares. Their exposed and immersed areas and drag coefficients differ. For iceberg mass $m$, the horizontal balance can be written schematically aswhere $f$ is the Coriolis parameter and $\zeta$ the sea-surface elevation. The last term represents the large-scale pressure gradient; contact with pack ice and grounding require additional forces. Iceberg free drift neglects such contact constraints, not the water drag or Coriolis force.

In still water with a fixed wind force $\mathbf{F}$, a transparent transient model linearizes water drag as $-m\gamma \mathbf{v}$. Put $J\mathbf{v}=\hat{\mathbf{z}}\times \mathbf{v}$. ThenThe damped free drift of an iceberg isThus an initially resting berg accelerates under the wind, turns under the Coriolis force, and approaches a balance of wind, water drag and rotation through damped inertial oscillations. Its equilibrium drift is to the right of the wind in the Northern Hemisphere and to the left in the Southern Hemisphere. With quadratic drag, the drag coefficient in the equilibrium balance depends on speed; the linear formula illustrates the turning and relaxation rather than prescribing an exact percentage of wind speed.

For a moving ocean, the background pressure gradient associated with geostrophic balance must also be included. Large tabular icebergs have substantial submerged area and respond strongly to currents; smaller bergs can have a larger wind-driven velocity relative to the water. Near Antarctica, the mainly westward Antarctic Coastal Current can carry bergs along the continental margin. Export into offshore gyres and then the eastward Antarctic Circumpolar Current changes their trajectories; bathymetry, winds, fronts and evolving size all matter. Gladstone and Bigg's Weddell Sea tracking illustrates the contrast between coastal and offshore drift. The Antarctic Convergence separates colder Antarctic surface water from warmer water farther north, where an exported iceberg can encounter enhanced melting. No single wind direction determines the complete route.

Let $L$ be a core-scale length, $B_{\ast }$ a representative strength of the prescribed poloidal magnetic field, and $v_A=B_{\ast }/\sqrt{{\mu }_0\rho }$ the Alfvén speed. A consistent weakly damped, small-amplitude ordering isThe first ratio is the Lehnert number. It makes the wave slow relative to rotation and allows almost geostrophic flow, with the azimuthal velocity approximately constant along a rotation-axis cylinder. The small wave amplitude gives $|B_{\varphi }|\sim \sqrt{{\mu }_0\rho }|U|\ll B_{\ast }$. Poloidal advection then changes neither the background field nor the wave substantially over one wave period. If there are no-slip boundary layers, also require their damping rate $\sqrt{\nu \Omega }/L\ll \sigma$, or otherwise retain their torque; a small bulk viscosity alone does not remove that boundary torque.

In the rotating magnetohydrodynamic momentum equation, axisymmetry removes the azimuthal pressure derivative. The neglected nonlinear azimuthal advection has size $|{\mathbf{u}}_p|U/L$, and the neglected viscous term has size $\nu U/L^2$. The retained azimuthal Coriolis force is $2\Omega u_r$. A typical compatible local balance has $|u_r|\sim \sigma U/\Omega$, small compared with $U$ when the Lehnert number is small. The azimuthal Lorentz force density can be written using $\nabla \cdot {\mathbf{B}}_p=0$:The physical prefactor is $1/({\mu }_0\rho )$: the extra $i$ printed inside the first reduced momentum equation is a typographical error, also contradicted by the subsequent real-valued system.

In the resistive induction equation, neglecting poloidal advection of the perturbation and magnetic diffusion leaves the winding of the prescribed poloidal magnetic field by differential rotation. ThusThis is the small-amplitude torsional Alfvén wave approximation, with a time-independent background and no external azimuthal forcing.

We now average over the full height of a rotation-axis cylinder. The sphere is impenetrable. Integrating incompressible flow over the volume inside any such cylinder gives zero flux through its cylindrical side, since there is no flux through its spherical caps. ThereforeThis is why the Coriolis force disappears from the cylinder average; it need not vanish locally. An axisymmetric field in the insulating exterior is a potential field with $B_{\varphi }=0$, so continuity gives $B_{\varphi }=0$ on the sphere. Consequently both the end flux terms and the moving-endpoint terms in the magnetic-stress integral vanish. Expanding the divergence givesSince $U$ is independent of $z$, its integrated time derivative is $2z_b{U}_t$. Defining $T={\int }_{-z_b}^{z_b}{B}_r{B}_{\varphi }dz$ givesThe stationary background and the fact that $U/r$ depends only on $r,t$ also givesoTogether these form the cylinder-averaged torsional Alfvén wave system.

Assume $H>0$ on the cylinders carrying the wave, with regular endpoint limits. If $H$ vanishes on an interior cylinder, then $B_r$ and $T$ vanish there and that cylinder has no radial magnetic coupling; the displayed $H^{-1}$ energy term must be interpreted on the coupled region rather than literally dividing by zero. Since $H$ is time independent, differentiation of the proposed conserved energy of a cylinder-averaged torsional wave givesRegularity at the rotation axis makes the lower endpoint vanish, while vanishing magnetic torque at the collapsed outer cylinder makes the upper endpoint vanish. This provesThe magnetic term has a direct interpretation. Decompose $B_{\varphi }=(T/H)B_r+B_{\varphi }^{\perp }$, with $\int B_r{B}_{\varphi }^{\perp }dz=0$. Then $\int B_{\varphi }^{2}dz=T^2/H+\int (B_{\varphi }^{\perp }{)}^{2}dz$, and the induction equation gives $(B_{\varphi }^{\perp }{)}_t=0$. Thus, apart from a common dimensional factor and a stationary orthogonal magnetic component, the conserved quantity is the physical kinetic plus magnetic energy.

For the frequency estimate, put $V=U/r$ and eliminate $T$:Locally this is a wave equation with speedIt is the cylinder-averaged radial field component that matters, not simply the dipole measured at the Earth's surface. Taking a cylinder height of order $6\times 10^6\mathrm{m}$, $\rho \sim 10^4\mathrm{kg}{\mathrm{m}}^{-3}$, and an assumed internal cylindrical radial field $B_{r,\mathrm{rms}}\sim 0.2$--$2\mathrm{mT}$ givesFor a standing wave of radial length $L_c\sim 2$--$3.5\times 10^6\mathrm{m}$, the cyclic frequency is of order $c_T/(2L_c)$ and the angular frequency is of order $\pi c_T/L_c$. HenceGeometry and the assumed internal field can shift the estimate substantially. As a concrete calibration, the 2010 torsional-wave study inferred a cylindrical radial rms field near $2\mathrm{mT}$ and a roughly four-year crossing time, with an observed six-year recurrence. A much weaker assumed interior field instead gives the older multi-decadal estimate. The paper supplies no numerical interior field strength, so a unique numerical frequency cannot be deduced without such an assumption.

Apply the rotating-frame derivative formula first to $r$, then to $v^{\prime }=[dr/dt{]}_{S^{\prime }}$. Since the angular velocity is constant in either frame,The two apparent forces in the rotating equation are the Coriolis force and the centrifugal acceleration multiplied by mass. ThusSet ${\varphi }_{\mathrm{eff}}=\varphi -\frac{m}{2}|\omega \times r{|}^2$. The supplied gradient identity gives $\nabla {\varphi }_{\mathrm{eff}}=\nabla \varphi +m\omega \times (\omega \times r)$. Taking the scalar product with $v^{\prime }$ removes $G$, since it is orthogonal to $v^{\prime }$, and removes the Coriolis force. Therefore the energy balance in a uniformly rotating frame isThe displayed quantity is conserved when the potential is stationary in the rotating coordinates. This includes a static inertial potential invariant under rotations about the chosen axis, such as vertical gravity in part (b). If the word conservative only means a time-independent inertial potential, that extra rotational invariance is needed; a conservative force alone does not imply the claimed conservation.

For an explicit counterexample to that unrestricted interpretation, take unit mass, an inertial force $F=e_x$, $\varphi =-x$, no $G$, rotation $\omega =\omega e_z$, and inertial trajectory $r(t)=(t^2/2,1,0)$ with initial velocity zero. It satisfies Newton's second law. Substituting $v^{\prime }=\dot{r}-\omega \times r$ into the proposed energy gives $E=\omega t$, which is not constant for nonzero rotation. The qualified conservation law above is the one used for the hoop.

Let $M=M_1+M_2$ and $\alpha =M_2/M$, with $0<\alpha <1$. In the centre of mass frame the two primary distances are $r_1=\alpha R$ and $r_2=(1-\alpha )R$. Applying Newton's law of universal gravitation and circular centripetal acceleration to the first primary gives $M_1{\omega }^2\alpha R=G{M}_1{M}_2/R^2$, and henceThe same result follows from the second primary or from the relative two-body equation of motion.

A rotating-frame equilibrium has zero rotating-frame velocity, so its Coriolis force vanishes. Above the orbital plane both primaries have downward gravitational components; below it both have upward components. The centrifugal force has no $z$ component. More explicitly,which vanishes only at $z=0$ for finite nonsingular points. Thus all equilibria lie in the orbital plane.

The equation of motion in a rotating frame for the dwarf isOn the orbital plane the centrifugal contribution is $m{\omega }^2\mathbf{r}$. Thus the rotating-frame effective potential is $E=U-m{\omega }^2(x^2+y^2)/2$; equilibria are exactly its stationary points. The Coriolis force affects motion about those equilibria but not their positions. This is the circular restricted three-body problem, with the dwarf treated as a test mass.

The printed energy unit has the wrong dimensions. The correct total energy unit is $mGM/R$, or $GM/R$ for specific energy; $GM/R^2$ is an acceleration. Dividing $E$ by $mGM/R$ and measuring lengths in units of $R$ gives the dimensionless effective potentialThe primary coordinates are $(-\alpha ,0)$ and $(1-\alpha ,0)$. Restricting to $y=0$ gives $F(x)=\mathcal{E}(x,0)$ andThe absolute values are essential; dropping them would reverse the attraction on the left of a primary. On each of the intervals separated by the two primaries, $F^{\prime }$ decreases strictly from $+\infty$ to $-\infty$. There is exactly one Collinear Lagrange point on each interval, a maximum along the $x$ direction.

For the inner Lagrange point, write $x=1-\alpha -\delta$ with $\delta >0$. For the L2 Lagrange point, write $x=1-\alpha +\delta$. The respective equations areBalancing the leading terms gives $\delta =(\alpha /3{)}^{1/3}+O({\alpha }^{2/3})$. Therefore the small-mass-ratio collinear Lagrange points obeyThe omitted $-\alpha$ displacement of the secondary is smaller than the stated error. The leading distance from the secondary is the dimensionless Hill radius.

For the L3 Lagrange point put $x=-1+c\alpha +O({\alpha }^2)$. It lies to the left of both primaries, soExpanding the two positive inverse-square factors gives $(1-\alpha )/(x+\alpha {)}^2=1+(2c+1)\alpha +O({\alpha }^2)$ and $\alpha /(x-1+\alpha {)}^2=\alpha /4+O({\alpha }^2)$. Thus $F^{\prime }=-(3c+5/4)\alpha +O({\alpha }^2)$, soAll three have $y=z=0$. The small-$\alpha$ equalities in the paper are asymptotic expressions, not exact coordinates for finite mass ratio.

For the remaining Lagrange points, introduce plane polar coordinates about the first primary: $x=-\alpha +s\cos \theta$, $y=s\sin \theta$. Its distance to the second primary is $d=(s^2+1-2s\cos \theta {)}^{1/2}$, andThe angular stationary condition is ${\partial }_{\theta }\mathcal{E}=\alpha s\sin \theta (d^{-3}-1)=0$. For an off-axis point $s>0$ and $\sin \theta \ne 0$, so $d=1$. The radial stationary condition becomesHence $s=1$, and $d=1$ implies $\cos \theta =1/2$. The Triangular Lagrange points are exactlyBoth triangles are equilateral. The collinear and off-axis calculations together exhaust the stationary points.

For the requested sketch, $\mathcal{E}$ tends to $-\infty$ at either primary and at large in-plane radius. Each Collinear Lagrange point is a saddle of the planar effective potential: its $x$ curvature is negative and its $y$ curvature is positive. To see the latter, let $d_1=|x+\alpha |$, $d_2=|x-1+\alpha |$ and $D=(1-\alpha )/d_1^3+\alpha /d_2^3$. At an inner point both distances are below one, so $D>1$. At an exterior stationary point, use barycentric positions $a_1=-\alpha$, $a_2=1-\alpha$ and weights $1-\alpha ,\alpha$: the equation implies $x(D-1)=\sum w_i{a}_i/|x-a_i{|}^3$. To the right, the weighted sum is positive because the positive primary is closer; to the left it is negative because the negative primary is closer. Thus $D>1$ there as well, and ${\mathcal{E}}_{yy}=D-1>0$.

At either Triangular Lagrange point, the planar Hessian isIts trace is $-3$ and determinant $27\alpha (1-\alpha )/4>0$, so these are local maxima. Their height is $\mathcal{E}=-3/2+\alpha (1-\alpha )/2$. The sketch below uses a finite mass ratio to separate the five points clearly; its singular wells are clipped only for display. These curvature labels concern the effective potential, not a stability test that ignores the Coriolis force.

Use rotating barycentric coordinates $\mathbf{R}=(x,y,z)$, rotating velocity $\mathbf{w}$, and $\mathbf{\Omega }=\hat{\mathbf{z}}$. Multiplying the three equations by the corresponding velocity components and adding cancels the Coriolis force terms:Thus the Jacobi constant is $C=2U-|\mathbf{w}{|}^2$. The inertial barycentric velocity is $\mathbf{V}=\mathbf{w}+\mathbf{\Omega }\times \mathbf{R}$. If $H_z=(\mathbf{R}\times \mathbf{V}{)}_z$ and $E_b=|\mathbf{V}{|}^2/2-{\mu }_1/r_1-{\mu }_2/r_2$, thenThe inertial energy $E_b$ and inertial angular momentum $H_z$ need not separately be constant. Their combination is constant because the two gravitational sources rotate steadily at unit angular speed.

For the Tisserand parameter take ${\mu }_2\ll {\mu }_1\simeq 1$, and compare the osculating orbital elements well outside a close encounter, where ${\mu }_2/r_2$ and barycentre-to-primary corrections are negligible. ThenThis requires a circular secondary orbit, the restricted three-body problem approximation, and no dissipative force or other perturber. $T$ is conserved to leading order between encounter episodes; the osculating $T$ during a close encounter need not be constant. The exact invariant is $C$.

For clarity the exact two-body energy relation is $a={\mu }_1/(2{\mu }_1/r_1-v_1^2)$, rather than the printed formula without the numerator ${\mu }_1$. At the leading order ${\mu }_1=1$ used for $T$, the two agree; retaining finite ${\mu }_1$ in just part of that formula is inconsistent.

At a close encounter put the secondary radius and speed equal to one at leading order, and denote the particle's incoming relative speed outside the secondary's strong-deflection region by $u$. Conservation of relative speed in the short gravitational encounter givesAn incoming orbit with $a=1$ has $|{\mathbf{v}}_1|=1$ at $r_1=1$ by the vis-viva equation; escape there requires an outgoing speed exceeding $\sqrt{2}$. The triangle inequality therefore gives the necessary single-encounter escape velocity boundEquality only gives a marginal parabolic limit. This argument deliberately gives a weak necessary bound, not a sufficient scattering criterion.

The planet-encounter relative velocity is related to the Tisserand parameter byHere $v_{1,\mathrm{tangential}}=h_z$ at the encounter radius, and the formula also applies to inclined passages at the planet's orbital plane. For the intended prograde coplanar incoming orbit, $I=0$ and $a=1$, so $u^2=2-2\sqrt{1-e^2}$. Applying the preceding necessary bound givesCoplanar alone also allows $I=\pi$. For that retrograde case $u^2=2+2\sqrt{1-e^2}$, and even a circular incoming orbit can be scattered onto an escaping orbit; the orbital eccentricity restriction therefore presumes prograde motion. Also, retaining the outgoing relative-velocity geometry gives the sharper necessary condition $|{\mathbf{v}}_{1,\mathrm{out}}|\le 1+u$, hence $u>\sqrt{2}-1$ and $e>\sqrt{4\sqrt{2}-5}/2$ for the same prograde incoming orbit. This stronger bound is consistent with, and implies, the weaker bound above.

For the Tisserand periapsis bound, first writeA bound orbit capable of another close encounter must intersect the secondary's circular radius, so $q\le 1\le Q$. For fixed $q<1$, the largest possible $T$ occurs at $I=0$. The remaining expression decreases with $Q\ge 1$, sinceThus every such orbit obeysOn $0<q<1$, $g(q)$ increases monotonically from $2$ to $3$. For $2<T<3$ the least accessible pericentre is the unique root of $T=g(q)$: the limiting orbit is prograde and coplanar, with apocentre just at the planet. Introduce $s=\sqrt{3-T}$. Its tangential speed at this apocentre is $1-s$, and its vis-viva equation gives the particularly well-behaved answerEquivalently, setting $w=\sqrt{2q/(1+q)}$ gives $T=2-w^2+2w$ and $w=1-s$, with $q=w^2/(2-w^2)$. Rationalizing the previous result yieldsThe apparent singularity at $T=\sqrt{8}$ is removable; the first form gives $q_{\min }=(\sqrt{2}-1)/2$ there. The limits are $q_{\min }\to 0$ as $T\downarrow 2$ and $q_{\min }\to 1$ as $T\uparrow 3$. The same lower envelope also excludes a smaller pericentre on an escaping trajectory when $2<T<3$. To see this without introducing an apocentre, use the planet-frame velocity sphere: $u=\sqrt{3-T}<1$ forces the tangential component at an encounter to be at least $w=1-u>0$, so total specific angular momentum obeys $h\ge w$. At pericentre the stellar energy obeys $h^2=2q+2E{q}^2$, while $E=h_z-T/2\le h-T/2$. ThusIf $q<q_{\min }=w^2/(2-w^2)$, then $q^2<w\le h$, and the left side is increasing with $h$. Hence $w^2\le 2q+(2w-T)q^2$. But $2w-T=w^2-2$, and the right side is smaller than $w^2$ for $0\le q<q_{\min }$, a contradiction. This establishes the barrier for all Kepler trajectories that reach the encounter radius, including unbound ones.

For $0<T\le 2$, bound encounter-crossing orbits can approach $q=0$ (take $Q\to 2/T$), so the positive root extrapolated from the displayed expression is not a universal barrier. For $T>3$ there is no trajectory reaching a close encounter in this approximation, since $u^2=3-T$ would be negative. Multiple encounters can approach the derived lower envelope by changing energy, angular momentum and inclination while preserving the Tisserand parameter; the invariant alone does not guarantee a particular encounter history reaches it.

Let $\mathbf{u}$ be ice velocity, $\mathbf{W}$ air velocity and ${\mathbf{U}}_o$ water velocity in an Earth-fixed frame. Put $J\mathbf{v}=\hat{\mathbf{z}}\times \mathbf{v}$ and $f=2\Omega \sin \varphi$. For effective air and water drag areas $A_a,A_w$, a quadratic drag convention commonly used for geophysical stresses givesThe Coriolis force magnitude is $m|f||\mathbf{u}|$, while its sign is retained in the vector equation. The coefficients here incorporate any conventional factor $1/2$; if one uses $\frac{1}{2}\rho C_{D}A|\mathbf{v}|\mathbf{v}$, then $C_D=2C$ describes the same resistance. Windage and water drag areas must be declared: an iceberg can have pressure drag on its exposed and submerged faces as well as basal skin drag. They cannot generally be inferred from ice mass alone.

To obtain an explicit quadratic-drag turning angle, first take still water and no sea-surface slope, and suppose $|\mathbf{u}|\ll W$. Then the wind force is approximately $D{\mathbf{e}}_x$ with $D=K_a{W}^2$, while water resistance remains quadratic. Steady free drift satisfiesIn the northern hemisphere write $\mathbf{u}=U(\cos \theta ,-\sin \theta )$ with $0\le \theta <\pi /2$. Resolving parallel to the drift and toward its left givesSquaring and adding eliminates $\theta$, rather than replacing water resistance by a fixed linear coefficient:The positive root and angle are thereforeThe stable root displayed avoids subtracting nearly equal large quantities. The drift is to the right of the wind for $f>0$ and to the left for $f<0$. Define the mass per drag area $\mu =m/A_w$, the area ratio $r=A_a/A_w$, the wind stress scale $a={\rho }_a{C}_{a}r{W}^2$ and the water-drag coefficient $b={\rho }_w{C}_w$. Dividing the force balance by $A_w$ givesIn particular, $\mu$ has units $\mathrm{kg}{\mathrm{m}}^{-2}$, not metres. A convenient dimensionless form isThe cosine follows by substituting $U=D\sin \theta /B$ into the along-drift equation, giving $1-{\cos }^2\theta =\chi \cos \theta$. Thus wind speed, drag coefficients and latitude are insufficient without mass and area normalization.

For $\chi \ll 1$, water drag dominates the turning correction:For $\chi \gg 1$, Coriolis force dominates the cross-wind balance:The cross-wind force does no work. Multiplying the original equilibrium equation by $\mathbf{u}$ gives $DU\cos \theta =K_w{U}^3$, which independently verifies the along-drift relation.