Asteroid Impact Probability Lab: Orbit Uncertainty, the B-Plane and Why Odds Collapse

Watch an impact probability rise, peak and collapse as observations shrink the uncertainty ellipse. A calm look at how scientists assess risk.

Trending ScienceBlender 5.2 assets
Drag to orbit, wheel / pinch to zoom

The story: how a risk-list entry evolves

Illustrative list: the other entries are made-up background objects with tiny odds.

Scenario

Live readouts

Impact probability (numerical integral)-
Torino scale (simplified)0
Monte-Carlo (20 000 virtual impactors)-
Palermo scale-
Impact energy-
Impact speed-
Effective Earth radius R_eff-
1-sigma semi-axes (major / minor)-
Ellipse area (1-sigma)-
Keyhole probability-
3-D: integrated hits / disc prediction-
Observations (N)-

Probability versus number of observations

Experiments

1. Three ways to get the same number

Circular Gaussian offset from a disc: closed form (Marcum Q) vs the numerical integral vs 1 000 000 Monte-Carlo virtual impactors.

2. Probability evolution curve

Sweeps the observation count for the current scenario and overlays Monte-Carlo points (100 000 samples each) on the graph above.

3. Gravitational focusing check

Integrates real flyby trajectories (RK4) and bisects for the grazing impact parameter, then compares with R_E√(1+v_esc²/v∞²).

The science (equations, assumptions and honest simplifications)

The B-plane

An asteroid's position at a future close approach is not known exactly; the orbit fit leaves a cloud of possible positions. Projected onto the plane through the Earth's centre that is perpendicular to the incoming asymptote (the B-plane, Opik 1976; Valsecchi et al. 2003) the cloud is well described by a 2-D Gaussian whose 1-sigma contour is a long thin ellipse: the long axis follows the "line of variations" (the direction in which the orbit is least constrained, mostly timing along the orbit).

Gravitational focusing

Energy and angular momentum are conserved on the flyby hyperbola: ½v²−μ/r = ½v∞² and L = b v∞. Setting the periapsis equal to the Earth's surface gives the largest impact parameter that still hits:

b_crit = R_eff = R_E * sqrt(1 + v_esc^2 / v_inf^2), v_esc = sqrt(2 mu / R_E) = 11.19 km/s

Slow encounters are focused strongly (the Earth looks much bigger than it is). Experiment 3 checks this formula against numerically integrated trajectories.

Impact probability

P = ∬∬_disc(R_eff) N(x; c, Σ) dA, Σ = R(θ) diag(σa², σb²) R(θ)ᵀ

The page evaluates this semi-analytically: a 1-D Simpson integral along the long axis of the exact chord probability across the short axis (Gaussian error function), and independently by Monte Carlo with sampled virtual impactors, P ≈ hits/n ± √(P(1−P)/n). For a circular Gaussian a closed form exists, P = 1 − Q₁(d/σ, R/σ) (noncentral chi-square, 2 degrees of freedom); Experiment 1 compares all three.

Why the odds rise, then collapse

While the ellipse is much larger than the Earth disc but still covers it, P ≈ A_disc / (2π σa σb) rises as the ellipse area shrinks. Real observation arcs shrink the area roughly like 1/N in this simple model (each semi-axis ∝ q/√N, q = astrometric error). The best-fit centre also moves: it equals the (unknown) truth plus an error that shrinks with the ellipse. If the truth lies off the disc, the ellipse eventually shrinks away from the Earth and P falls to nearly zero; if it lies on the disc, P climbs toward 100 %. The probability is a statement about our knowledge, not about the asteroid, and the fall in P is the good outcome of the same process.

Keyholes

A tiny region of the B-plane (kilometres to hundreds of kilometres across) deflects the asteroid onto an orbit that is in resonance with the Earth and brings it back for a later impact (Chodas 1999; Milani et al. 2005). Its probability is the same Gaussian integral over that small circle; the page uses a single illustrative keyhole and assumes a certain hit if it is threaded.

Energy and the two scales

E = ½ m v_imp², m = ρ π/6 D³, ρ = 2600 kg/m³, v_imp = sqrt(v_inf² + v_esc²) (1 Mt TNT = 4.184e15 J) Palermo: PS = log10( P_i / (f_B T) ), f_B = 0.03 E^−0.8 per year (E in Mt) (Chesley et al. 2002)

The Torino scale (0–10; Binzel 2000) combines probability and energy on a chart. This page uses a simplified classifier: energy bands (<1 Mt none, 1–300 Mt local, 300–1e5 Mt regional, above global) and probability thresholds (1e-6, 1e-4, 1 %, 10 %, 50 %, 99 %). It follows the shape of the official chart but is not an official rating. Real energy also depends on density, strength, angle and atmospheric breakup: the order of magnitude is what matters.

The news hook (reported)

Asteroid 2024 YR4 was reported in early 2025 to have briefly reached Torino level 3, the highest rating assigned to any object, with an Earth-impact probability for 2032 reported at about 3 % at its peak. Further observations then reduced the Earth-impact probability to effectively zero; a small chance of a lunar impact was reported afterwards. The preset "Rise then collapse" is a stylised illustration with invented numbers shaped like that story, not the real orbit solution.

Simplifications

Deflection, in one line

Given years of warning, a tiny velocity change of about 1 cm/s moves the B-plane point by thousands of kilometres; see Asteroid Deflection and Asteroid Trajectory Correction. For what an impact does, see Asteroid Impact Simulation: Planetary Effects.

Sources

Opik, E. J. (1976) Interplanetary Encounters. Binzel, R. P. (2000) The Torino Impact Hazard Scale, Planet. Space Sci. 48. Chesley, Chodas, Milani, Valsecchi & Yeomans (2002) Quantifying the risk posed by potential Earth impacts, Icarus 159. Valsecchi, Milani, Gronchi & Chesley (2003) Resonant returns to close approaches, A&A 408. Milani et al. (2005) Icarus 173. NASA/ESA planetary-defence reporting on 2024 YR4 (2025). NASA DART mission (2022) for deflection.

As of public reporting, no known asteroid poses a significant impact threat to Earth. This is an educational model of how risk assessment works, not a prediction about any real object.

Made with 🟠 Blender · Three.js