Cometary Non-Gravitational Acceleration: The Outgassing Jets That Bend a Comet's Orbit

What It Is: A Comet as a Leaky Rocket

A comet nucleus is a few-kilometre lump of ice, dust, and frozen gases — a dirty snowball. When it approaches the Sun, solar heating drives sublimation: water ice (and more volatile species like CO and CO₂) turns directly to gas and streams off the surface, dragging dust with it to form the coma and tails.

Crucially, this outflow is not symmetric. Sublimation is strongest on the sunward, daytime hemisphere, and rotation plus a thermal lag mean the peak jet activity is offset from the sub-solar point. By Newton's third law, gas leaving in one direction pushes the nucleus in the opposite direction. The net momentum flux is a tiny thrust — a non-gravitational force — superimposed on solar gravity.

• The force scales with the mass-loss rate, which rises steeply as the comet nears the Sun.
• Its direction depends on the spin axis, active-region geography, and thermal lag — so it has radial, along-track, and out-of-plane components.
• The effect is small per second but cumulative: integrated over an orbit it shifts perihelion time by hours to days.

The Mechanism and the Marsden Model

The thrust equals the momentum carried away by escaping gas: F = Q · m · v_th · f, where Q is the molecular production rate (molecules/s), m the molecular mass, v_th the outflow speed (~0.5–1 km/s for water near 1 AU), and f a dimensionless collimation efficiency (0 for isotropic loss, up to ~1 for a perfect jet). Real comets have f of order 0.01–0.5.

Because tracking Q directly is hard, Brian Marsden, Zdenek Sekanina, and Donald Yeomans (1973) introduced an empirical model still used by JPL and the Minor Planet Center. The acceleration is written in the comet's orbital frame:

• a_ng = [A1·R̂ + A2·T̂ + A3·N̂] · g(r), with R̂ radial (Sun→comet), T̂ transverse (along-track), N̂ normal to the orbit plane.
• A1, A2, A3 are constants fit to astrometry alongside the six orbital elements.
• The dimensionless g(r) encodes how water sublimation ramps up sunward: g(r) = α (r/r₀)⁻²·¹⁵ [1 + (r/r₀)⁵·⁰⁹]⁻⁴·⁶¹, normalized so g = 1 at r₀ = 2.808 AU (α ≈ 0.1113). Inside ~2.8 AU it climbs steeply; far out it falls off roughly as r⁻² like insolation.

The transverse term A2 is the one that changes orbital energy and period — it is what makes Encke arrive early.

Characteristic Numbers and a Worked Example

Take a nucleus of radius R ≈ 1 km, density ρ ≈ 500 kg/m³, so mass M ≈ (4/3)π R³ ρ ≈ 2 × 10¹² kg. Near 1 AU an active comet might lose water at Q ≈ 10²⁸ molecules/s (m ≈ 3 × 10⁻²⁶ kg), i.e. a mass-loss rate of ~300 kg/s, ejected at v_th ≈ 0.5 km/s.

• Thrust: F ≈ (300 kg/s)(500 m/s) · f ≈ 1.5 × 10⁵ · f newtons.
• Acceleration: a = F/M ≈ (1.5 × 10⁵ · f) / (2 × 10¹²) ≈ 7.5 × 10⁻⁸ · f m/s². With f ≈ 0.05, a ≈ 4 × 10⁻⁹ m/s² — comparable to the peak radial value Rosetta measured at 67P (1.28 × 10⁻⁸ m/s²).

Compare to solar gravity at 1 AU (5.9 × 10⁻³ m/s²): the non-gravitational term is ~10⁻⁶ of gravity here, larger nearer perihelion. Over one 3.3-year orbit of Encke, the transverse push advances perihelion by ~2–3 hours; two centuries ago the effect was stronger and it drifted faster. Because a = F/M and F is measurable astrometrically, fitting A2 gives a lower limit on nucleus mass — a rare remote way to weigh a comet.

How It Is Observed and Detected

Non-gravitational acceleration reveals itself as a residual: after fitting a gravity-only orbit (Sun plus planetary perturbations) to precise astrometry, the observed positions still drift systematically. The signature is a smooth departure that peaks around perihelion and correlates with the comet's activity.

• Timing: historically, a comet returning consistently early (Encke) or late signals a persistent A2 term.
• Astrometric fitting: modern orbit solutions (JPL Horizons, MPC) solve for A1, A2, A3 whenever the data span and precision allow — thousands of comets carry these parameters.
• Spacecraft radio tracking: Rosetta measured 67P's acceleration directly via Doppler, resolving radial and normal components peaking (1.28 ± 0.17) × 10⁻⁸ and (0.52 ± 0.20) × 10⁻⁸ m/s² about 15–24 days after perihelion — a clear thermal-lag asymmetry.
• Extreme cases: 'Oumuamua's (4.92 ± 0.16) × 10⁻⁶ m/s² excess was detected at 30σ despite no visible coma.

The steep heliocentric dependence — for 67P the force scaled roughly as r⁻⁶ near perihelion — is itself diagnostic of sublimation-driven thrust rather than, say, radiation pressure (which follows a clean r⁻²).

How It Differs From Related Effects

Several small forces perturb small bodies; non-gravitational cometary acceleration is specifically the sublimation rocket and should not be conflated with its cousins:

• Radiation pressure: sunlight's photon momentum pushes purely radially, ∝ r⁻², and dominates for dust grains and any putative solar sail. It cannot change orbital energy the way a transverse outgassing term can.
• Yarkovsky effect: for asteroids, anisotropic thermal re-emission of absorbed sunlight yields a tiny along-track force (~10⁻¹⁵ m/s²-class semimajor-axis drift). It needs no volatiles and is far weaker; it dominates asteroid orbital evolution over Myr, not per-orbit.
• YORP: the rotational analogue of Yarkovsky, torquing spin rather than translating the orbit.
• Poynting–Robertson drag: makes small dust spiral inward; irrelevant to a km-scale nucleus.

The key discriminants are the steep, non-r⁻² heliocentric scaling, the correlation with observed coma/gas production, and a transverse (energy-changing) component — all hallmarks of outgassing rather than photon or thermal-lag forces on inert rock.

Significance, Famous Cases, and Open Questions

Non-gravitational forces are not a nuisance to be subtracted away — they are a probe of cometary physics and a hazard-assessment necessity. Accurate impact predictions for Earth-crossing comets require them; ignoring A2 can move a predicted perihelion by days.

• Comet Encke: the archetype. Its secular acceleration peaked around 1820 and has smoothly declined for two centuries as volatiles deplete; it has recently even flipped to a slight deceleration — a comet aging toward dormancy in real time.
• 67P/Churyumov–Gerasimenko: Rosetta tied the orbit-scale force directly to mapped active regions and measured ~10.5 × 10⁹ kg (≈0.1%) of mass lost across the 2015 perihelion.
• 1I/'Oumuamua: its strong radial acceleration with no detectable coma ignited debate — hyper-active tiny outgassing (H₂ or N₂ ice?), an ultra-low-density fractal body, or (more speculatively) a solar sail. Its true nature remains unresolved.

Open questions include what sets the collimation factor f, how to model asymmetric pre/post-perihelion activity from first principles, and whether non-gravitational fits can reliably yield nucleus masses for comets we never visit.

Representative non-gravitational acceleration magnitudes (radial component, near perihelion) across bodies

Body	NG acceleration (m/s²)	As fraction of solar gravity	How measured
Comet Encke (short-period)	~10⁻⁷ near 0.34 AU	~10⁻⁴	Perihelion timing drift over centuries
67P/Churyumov–Gerasimenko	1.28 × 10⁻⁸ (radial peak)	~10⁻⁵	Rosetta radio tracking, 2014–2016
1I/'Oumuamua (interstellar)	4.92 × 10⁻⁶ at ~1.4 AU	~2 × 10⁻³	30σ astrometric fit (Micheli 2018)
Typical long-period comet	10⁻⁸ to 10⁻⁶	10⁻⁶ to 10⁻⁴	Astrometry + Marsden A1/A2/A3 fit
Weak/dormant comet	< 10⁻⁹	< 10⁻⁶	Often undetectable in the orbit