Yarkovsky Effect

Yarkovsky effect is a non-gravitational thermal recoil perturbation acting on rotating small Solar System bodies—primarily asteroids and meteoroids ranging from ~10 centimeters to ~40 kilometers in diameter—caused by the anisotropic emission of thermal infrared photons from a surface with finite thermal inertia. Because a rotating asteroid's surface requires a finite time to heat up under solar illumination and cool down after sunset, its warmest surface point is displaced away from the subsolar noon meridian toward the afternoon hemisphere. Consequently, the net momentum carried away by outbound thermal infrared photons (λ ≈ 5–30 μm) is not directed radially back toward the Sun, but includes a small transverse (along-track) force component that continuously increases or decreases the asteroid's orbital mechanical energy E = −G M_☉ m / (2a), driving a cumulative secular drift in its orbital semi-major axis (da / dt).

Originally proposed in a Russian engineering pamphlet around 1900 by the Polish-Russian civil engineer Ivan Osipovich Yarkovsky (1844–1902), the effect would have been completely forgotten had the Estonian astronomer Ernst J. Öpik not read Yarkovsky's pamphlet around 1909 and recalled it from memory forty years later in 1951. Modern radar and spacecraft tracking (Steven R. Chesley et al., 2003 on asteroid 6489 Golevka; OSIRIS-REx on 101955 Bennu) has established the Yarkovsky effect as the primary physical mechanism that feeds main-belt asteroids and meteorites into chaotic orbital resonances with Jupiter and Saturn, delivering both Earth-crossing asteroids and prebiotic carbonaceous chondrite volatiles to the terrestrial planets, while simultaneously governing long-term planetary defense impact hazard predictions for near-Earth asteroids such as 99942 Apophis and 101955 Bennu.

Every second, a spherical asteroid of radius R and Bond albedo A orbiting at heliocentric distance r intercepts solar radiation power P_abs = (1 − A) π R² Φ_☉(r) and re-radiates that exact same thermal power in the mid-infrared at equilibrium temperature T_eq = [ (1 − A) Φ_☉ / (4 ε σ_SB) ]^1/4. Since each emitted infrared photon of energy E_γ carries linear momentum p_γ = E_γ / c, an element of surface area dS at temperature T with outward unit normal n̂ experiences a Lambertian thermal recoil force:

• dF_th = −(2 / 3c) ε σ_SB T⁴ n̂ dS (where the factor of 2/3 arises from integrating cos² θ over the outward hemisphere).

If the asteroid had zero thermal inertia (Γ = 0), its surface temperature would track incident sunlight instantaneously, peaking exactly at local solar noon (12:00) and creating a thermal recoil force directed 100% radially outward along the Sun–asteroid line (which merely weakens solar gravity slightly without changing the semi-major axis a). Conversely, if the asteroid had infinite thermal inertia (Γ → ∞) or rotated infinitely fast, its temperature would be uniform across all longitudes, yielding zero net recoil force.

In a real regolith-covered asteroid, finite thermal inertia (Γ = √(κ ρ C_p), in units of J m⁻² K⁻¹ s^−1/2, where κ is thermal conductivity, ρ is bulk density, and C_p is specific heat capacity) delays peak surface temperature by a thermal phase lag angle φ past local solar noon into the afternoon (~14:00 to 15:00 local solar time). This produces two distinct components of the Yarkovsky effect:

1. The Diurnal Yarkovsky Effect (Spin-Sense Dependent):
   • Prograde Rotators (Obliquity γ ≈ 0°, spinning in the same counter-clockwise sense as their orbit around the Sun): As the surface rotates from noon (facing the Sun) toward dusk, the hot afternoon hemisphere faces backward against the direction of orbital motion. Therefore, the excess thermal infrared photons are beamed backward like a miniature rocket exhaust, imparting a positive along-track thrust (F_trans > 0) in the direction of orbital motion. By Gauss's planetary perturbation equation (da/dt = [2 a^3/2 / √(G M_☉)] · (F_trans / m)), positive along-track thrust increases orbital mechanical energy and causes the asteroid to spiral outward away from the Sun (da / dt > 0, proportional to cos γ).
   • Retrograde Rotators (Obliquity γ ≈ 180°, spinning opposite to their orbital motion): Because the surface rotates from noon toward the leading orbital hemisphere, the hot afternoon side faces forward along the orbit. Outbound infrared photons are beamed forward, creating a braking force (F_trans < 0) that drains orbital energy and causes the asteroid to spiral inward toward the Sun (da / dt < 0).
2. The Seasonal Yarkovsky Effect (Always Inward Drag, da / dt < 0):
   • Discovered by David P. Rubincam (1995, 1998), the seasonal component arises when an asteroid's spin axis lies near its orbital plane (obliquity γ ≈ 90°). During one half of the asteroid's heliocentric year, the northern pole faces the Sun, reaches maximum temperature after summer solstice due to deep seasonal thermal conduction (penetration depth l_s = √(κ / [ρ C_p n]) ~ 1–10 meters), and radiates excess infrared heat opposite to the velocity vector as it moves toward the autumn equinox. During the other half of the year, the southern pole does the exact same thing. Consequently, the seasonal Yarkovsky effect always acts as an orbital drag (da / dt ∝ −sin² γ ≤ 0), causing both prograde and retrograde bodies with tilted spin axes to drift inward.

In linear heat-diffusion theory (David Vokrouhlický, 1998–1999), the strength of the diurnal Yarkovsky force is governed by the dimensionless thermal parameter Θ, which compares the thermal relaxation rate of the surface to the asteroid's rotation angular velocity ω = 2π / P_rot:

• Θ = (Γ √ω) / (ε σ_SB T_sub³).

The transverse Yarkovsky acceleration reaches its maximum when Θ ≈ 1 (yielding a thermal phase lag angle φ ≈ 20°–30°):

• a_Yark,trans ≈ −(4 / 9) [ (1 − A) Φ_☉ / (c ρ R) ] · [ Θ / (1 + 2Θ + 2Θ²) ] · cos γ.

Notice the 1 / (ρ R) scaling in the denominator: because thermal thrust scales with cross-sectional area (∝ R²) while mass scales with volume (∝ ρ R³), smaller, lower-density rubble-pile asteroids experience much stronger Yarkovsky drift than large monolithic bodies—until diameter drops below ~10 cm, where the thermal penetration wave l_d = √(κ / [ρ C_p ω]) heats the entire pebble isothermally and switches off the temperature gradient.

1. V-Shaped Asteroid Families and Age Dating: When a catastrophic collision shatters a parent asteroid in the Main Belt (producing families such as the Karin, Koronis, Eos, or Massalia families), smaller fragments (D = 1–5 km) experience faster Yarkovsky semi-major axis drift (da/dt ∝ 1/D) than large fragments (D = 20–50 km), spreading outward (prograde) and inward (retrograde) from the parent semi-major axis a₀. Plotting inverse diameter 1 / D versus semi-major axis a reveals a sharp "Yarkovsky V-shape" whose boundary slope directly measures the collisional age of the asteroid family (ranging from 5.8 Myr for the young Karin cluster to > 2 Gyr for ancient families).
2. Delivery of Meteorites and Near-Earth Asteroids via Kirkwood Gaps: For decades, dynamicists could not explain why freshly chipped meteorites (whose cosmic-ray exposure ages are typically 5 to 50 million years for stony chondrites) continuously arrive on Earth long after the Kirkwood mean-motion resonances (3:1 at 2.50 AU, 5:2 at 2.82 AU) and the ν₆ secular resonance should have been emptied of bodies. William F. Bottke, David Vokrouhlický, David Rubincam, and David Nesvorný (2000–2006) proved that the Yarkovsky effect is the conveyor belt: collisions in the Main Belt produce meter-to-kilometer fragments that slowly drift in semi-major axis via Yarkovsky forces over 10–100 Myr until they walk across the edge of the 3:1 or ν₆ resonance, where eccentricities are pumped up to Earth-crossing values (q < 1.0 AU) in < 2 Myr.
3. High-Precision Radar & Spacecraft Measurements (Golevka, Bennu, and Apophis):
   • 6489 Golevka (2003): Using Arecibo and Goldstone delay-Doppler radar ranging between 1991 and 2003, Steven Chesley and collaborators achieved the first direct detection of the Yarkovsky effect on an asteroid, measuring a 15-kilometer orbital displacement over 12 years (da/dt = −6.39 × 10⁻⁴ AU/Myr) and deriving Golevka's bulk density (ρ ≈ 2.7 g/cm³).
   • 101955 Bennu (OSIRIS-REx, 2019–2021): NASA's OSIRIS-REx spacecraft measured the retrograde carbonaceous asteroid Bennu (D ≈ 492 m, γ = 177.6°) drifting inward at da / dt = −284 ± 0.2 meters per year (−1.899 × 10⁻³ AU/Myr), allowing mission navigators (Davide Farnocchia et al., 2021) to determine Bennu's mass (M = 7.329 × 10¹⁰ kg, bulk density ρ = 1.190 g/cm³, ~50% macroporosity) and constrain its cumulative Earth impact probability between 2178 and 2290 (1 in 1,750, peaking on September 24, 2182) to meter-level accuracy through the 2135 gravitational keyhole.
   • 99942 Apophis (2020–2021): Optical and Goldstone/Green Bank radar observations proved that Apophis (D ≈ 340 m) experiences an inward Yarkovsky drift of da / dt ≈ −170 meters per year, which shifted its trajectory relative to the April 13, 2029 Earth close-approach keyhole and definitively ruled out the previously feared 2068 Earth impact scenario.

1. Öpik, E. J. (1951). Collision probabilities with the planets and the distribution of interplanetary matter. Proceedings of the Royal Irish Academy. Section A: Mathematical and Physical Sciences, 54, 165–199. https://www.jstor.org/stable/20488532
2. Chesley, S. R., Ostro, S. J., Vokrouhlický, D., Čapek, D., Giorgini, J. D., Nolan, M. C., ... & Benner, L. A. M. (2003). Direct detection of the Yarkovsky effect by radar ranging to asteroid 6489 Golevka. Science, 302(5651), 1739–1742. https://doi.org/10.1126/science.1091452
3. Bottke, W. F., Jr., Vokrouhlický, D., Rubincam, D. P., & Nesvorný, D. (2006). The Yarkovsky and YORP effects: Implications for asteroid dynamics. Annual Review of Earth and Planetary Sciences, 34(1), 157–191. https://doi.org/10.1146/annurev.earth.34.031405.125154
4. Farnocchia, D., Chesley, S. R., Takahashi, Y., Rozitis, B., Vokrouhlický, D., Rush, B. P., ... & Lauretta, D. S. (2021). Ephemeris and hazard assessment for near-Earth asteroid (101955) Bennu based on OSIRIS-REx data. Icarus, 369, 114594. https://doi.org/10.1016/j.icarus.2021.114594