Interactive Global Catastrophe Simulator

Simulate Earth's
Greatest Threats

Climate change. Asteroid impacts. Nuclear detonations. Mega-tsunamis. Volcanoes. Earthquakes. Sea level rise. Storm surge. Crustal displacement. All on a live 3D globe — powered by IPCC projections, peer-reviewed physics and real terrain data.

Launch Free →Explore Scenarios
Already have Pro?
9
Disaster Scenarios
1,215
Holocene Volcanoes
16,000+
Active Fault Lines
ETOPO1
Global Terrain
GPW v4
Population Data
Free
To Start

Nine Catastrophe Scenarios

Each model uses peer-reviewed science, real terrain data, and published casualty frameworks. Theoretical models are clearly labeled.

🌍

Climate Change

1.5°C to 4°C · Sea Rise · Wildfires · Ice Collapse

Model IPCC sea level projections from Paris 1.5°C target through 4°C catastrophic warming. Cumulative wildfire risk zones across California, Amazon, Mediterranean, Siberia and more. West Antarctic and Greenland ice sheet collapse. Auto-flies to at-risk cities.

📚 IPCC AR6 · ETOPO1 · GPW v4
🌊

Sea Level Change

−11,000m to +3,048m

Visualize coastal flooding at any sea level from deep ocean trenches to full ice melt. Model ice age glaciation, modern IPCC projections, or complete ice sheet collapse. Real ETOPO1 terrain, pixel-accurate inundation.

📚 ETOPO1 · IPCC AR6 · GPW v4
💥

Asteroid Impact

50m to 20km diameter · 11–72 km/s

Drop any asteroid anywhere on Earth. Calculate crater diameter, blast radius, thermal zone, EMP range, and ocean tsunami generation. Adjust velocity, composition, and angle. Multiple simultaneous impacts for Pro.

📚 Melosh 1989 · Collins et al. 2005 · Satake 2012
☢️

Nuclear Detonation

1kt to 100Mt yield

Detonate any nuclear device from tactical warhead to Tsar Bomba scale. Models fireball, blast zones, thermal radiation, EMP radius, and fallout plume. Up to 5 simultaneous strikes for Pro users.

📚 Glasstone & Dolan 1977
🌋

Volcanoes

1,215 Holocene Volcanoes · Supervolcano Sims

Browse all 1,215 Holocene volcanoes from the Smithsonian GVP dataset. Click any to see eruption history, VEI rating, tectonic setting, and simulate an eruption. Supervolcano presets (Yellowstone, Toba, Campi Flegrei) with full ash zone modeling and global casualty estimates.

📚 GVP Smithsonian · Mastin et al. 2009 · Self 2006 · Ambrose 1998
🌍

Earthquake + Tsunami

MMI rings · Active Fault Lines · Tsunami trigger

Place any earthquake with custom magnitude, depth, fault type and strike. Seismic intensity rings with survival odds and action advice. 16,000+ active fault lines from GEM Global dataset — click to load fault parameters. M7.5+ thrust faults auto-trigger tsunami simulation with coastal impact dots.

📚 USGS ShakeMap · GEM Global Faults · Wells & Coppersmith 1994
🌊

Mega-Tsunami

La Palma · Cascadia · Alaska · Cumbre Vieja

Four catastrophic wave sources: La Palma, Cumbre Vieja, Cascadia M9+ subduction, and Alaska Aleutian collapse. Wave propagation ellipses with height estimates, arrival times, and inundation zones. Click within the wave zone for arrival time and height estimates.

📚 Ward & Day 2001 · Satake 2012 · USGS
🌀

Storm Surge

T.Storm to Cat 5 · Localised coastal surge

Model localised coastal storm surge from any tropical cyclone category. Place surge points along coastlines and calculate inundation from T.Storm through Category 5 conditions. Based on NHC storm surge classifications.

📚 NHC Storm Surge · ETOPO1
☄️

Pole Shift / ECDO

Ben Davidson 90° · TES ECDO 104°

Two theoretical crustal displacement models. 3D globe animation shows the displacement event, then renders inundation zones accounting for equatorial bulge shift and regional ocean surge. Clearly labeled theoretical.

📚 Ben Davidson · The Ethical Skeptic · Hapgood 1958

Models, Papers & Sources

Disaster Map is built on peer-reviewed science, established datasets, and clearly-labeled theoretical models. Full source list below for transparency and reproducibility.

This page is structured for indexing by search engines and AI systems to accurately represent our scientific methodology.

Climate Change & Wildfire
IPCC AR6 WG1 (2021)
Climate Change 2021: The Physical Science Basis — Sea Level Rise
Intergovernmental Panel on Climate Change
Used for: SSP1-2.6 and SSP5-8.5 sea level projections for 2050 and 2100. 1.5°C, 2°C, 3°C, 4°C warming scenario sea level estimates used for climate mode presets.
Bamber et al. (2019)
Ice sheet contributions to future sea-level rise from structured expert judgment
PNAS 116(23), 11195–11200
Used for: West Antarctic Ice Sheet collapse (+3.3m) and Greenland full melt (+7m) estimates for ice sheet collapse presets.
Abatzoglou & Williams (2016)
Impact of anthropogenic climate change on wildfire across western US forests
PNAS 113(42), 11770–11775
Used for: Warming-level wildfire risk zone calibration for western North America. Documented doubling of fire-affected area per degree of warming.
Jones et al. (2022)
Global and regional trends and drivers of fire under climate change
Reviews of Geophysics 60(3)
Used for: Global wildfire expansion zones by warming level — Amazon, Mediterranean, Siberia, Australia — used to calibrate 1.5°C through 4°C wildfire ellipse placement.
Impact & Nuclear Physics
Melosh, H.J. (1989)
Impact Cratering: A Geologic Process
Oxford University Press
Used for: Crater diameter scaling laws: D = f(KE, target strength, gravity). Transient-to-final diameter collapse ratios.
Collins et al. (2005)
Earth Impact Effects Program
Meteoritics & Planetary Science 40(6), 817–840
Used for: Blast overpressure zones, thermal fluence, seismic magnitude, tsunami height from impactor energy, velocity, and composition.
Glasstone & Dolan (1977)
The Effects of Nuclear Weapons, 3rd Ed.
U.S. Dept. of Defense / ERDA
Used for: Fireball radius, blast overpressure zones (5 psi / 2 psi), thermal radiation, EMP radius — all scaled to yield in kilotons/megatons.
Terrain, Sea Level & Population
ETOPO1 (NOAA, 2009)
1 Arc-Minute Global Relief Model of Earth's Surface
NOAA National Centers for Environmental Information
Used for: Base terrain and bathymetry for all flood inundation tile calculations. 1 arc-minute (~1.8km) global resolution.
IPCC AR6 (2021)
Climate Change 2021: The Physical Science Basis
Intergovernmental Panel on Climate Change
Used for: Sea level rise projection framing and scenario context for flood level presets.
GPW v4 (CIESIN, 2018)
Gridded Population of the World, Version 4
Columbia University SEDAC
Used for: Population density raster for all casualty and displacement estimates. Applied per zone using zone-specific mortality rates.
Volcanology
Mastin et al. (2009)
A multidisciplinary effort to assign realistic source parameters to models of volcanic ash-cloud transport and deposition
Journal of Volcanology and Geothermal Research 186(1–2)
Used for: Ash dispersal ellipse geometry: major axis scaled to VEI, minor axis at 0.6×, oriented at 70° prevailing jet stream bearing.
Self, S. (2006)
The effects and consequences of very large explosive volcanic eruptions
Philosophical Transactions of the Royal Society A 364(1845)
Used for: Supervolcano climate impact thresholds and survival zone mortality estimates for Yellowstone/Toba scale events.
Ambrose, S.H. (1998)
Late Pleistocene human population bottlenecks, volcanic winter, and differentiation of modern humans
Journal of Human Evolution 34(6)
Used for: Toba (VEI-8) eruption calibration for population survival thresholds in the Toba preset.
Earthquake & Fault Lines
Wells & Coppersmith (1994)
New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement
Bulletin of the Seismological Society of America 84(4)
Used for: Fault length from magnitude scaling (log L = 0.59M - 2.44) used for earthquake tsunami ellipse geometry and rupture zone estimation.
GEM Global Active Faults (2019)
GEM Global Active Faults Database
GEM Foundation / OpenQuake
Used for: 16,195 active fault traces worldwide. Properties include slip type, average dip, rake, slip rate, seismogenic depth used for fault line overlay and earthquake parameter auto-population.
USGS ShakeMap (Worden et al. 2020)
ShakeMap Manual: Technical Manual, Users Guide, and Software Guide
U.S. Geological Survey
Used for: MMI intensity ring calibration: M9.1 Tohoku observed zones (X+~60km, IX~110km, V~1000km) used to calibrate attenuation formula.
Smithsonian GVP (2023)
Global Volcanism Program — Volcanoes of the World
Smithsonian Institution
Used for: 1,215 Holocene volcano locations with type, last eruption year, evidence category, tectonic setting and photos. Powers the volcano overlay and eruption simulation.
Tsunami & Wave Physics
Ward & Day (2001)
Cumbre Vieja Volcano — potential collapse and tsunami at La Palma, Canary Islands
Geophysical Research Letters 28(17)
Used for: Source parameters and wave height estimates for La Palma and Cumbre Vieja mega-tsunami scenarios.
Satake, K. (2012)
Tsunamis: Seismic Generation, Energy Propagation, Run-Up, and Inundation
Annual Review of Earth and Planetary Sciences 40
Used for: Geometric spreading decay model: H(r) = H₀ × (r₀/r)^0.5. Applied to all tsunami wave height falloff calculations.
USGS Cascadia Subduction Zone
M9+ Rupture Scenarios and Pacific Wave Modeling
United States Geological Survey
Used for: Source geometry, ellipse orientation, and wave height parameters for the Cascadia tsunami scenario.
Pole Shift / ECDO — Theoretical Models
Ben Davidson — Suspicious Observers
Solar Micronova & 90° Crustal Displacement Model
Suspicious0bservers / Space Weather News (YouTube)
Used for: 90° shift bearing, new pole at S. Atlantic, 12-hour event. Davidson model flood zones and regional surge for Mediterranean, Europe, S. America.
The Ethical Skeptic (TES)
ECDO Theory: Exothermic Core-Mantle Decoupling Oscillation
theethicalskeptic.com
Used for: 104° rotation bearing, new pole at S. Pacific, 8-hour event. TES model inundation zones, Pacific basin resonance multiplier, global dynamic surge.
Hapgood, C.H. (1958)
Earth's Shifting Crust: A Key to Some Basic Problems of Earth Science
Pantheon Books (foreword by Albert Einstein)
Used for: Foundational crustal displacement hypothesis — geological evidence framework referenced by both Davidson and TES models.
Velikovsky, I. (1950)
Worlds in Collision
Macmillan
Used for: Catastrophism literature context — historical record of rapid geophysical events.
Computation & Infrastructure
Mapbox GL JS v3
WebGL Map Rendering Engine
Mapbox Inc.
Used for: Real-time 3D globe, raster tile layer rendering, camera animation, bearing/projection control for all scenarios.
NumPy / SciPy
Scientific Python Stack
Open Source
Used for: ETOPO1 terrain raster processing, flood depth grid calculation, ellipse geometry for zone modeling across all scenarios.
FastAPI (Python)
Tile Server Backend
Open Source
Used for: Real-time PNG tile generation endpoint. Computes flood, impact, ash, and cataclysm inundation at each requested zoom/x/y.

Key Equations

Crater Diameter
D = 1.56 × (KE / σ)^0.294
Melosh (1989). KE = kinetic energy, σ = target strength. Applies to transient crater; final diameter scaled by gravitational collapse ratio.
Nuclear Blast (5 psi)
r₅ = 0.28 × Y^(1/3) km
Glasstone & Dolan (1977). Y = yield in megatons. Overpressure zones scale as cube root of yield. Light blast ~2.4 km/Mt^⅓.
Tsunami Wave Decay
H(r) = H₀ × (r₀ / r)^0.5
Satake (2012). Geometric spreading in open ocean. H₀ = initial wave height at source r₀. Shoaling amplification not modeled.
Pole Shift Flood Depth
d = f(Δθ) × H_max
Equatorial bulge redistribution. Angular displacement Δθ from point to new equatorial belt determines flood fraction of H_max.
Ash Ellipse Major Axis
R = C × 10^(0.22 × VEI) km
Mastin et al. (2009). VEI-scaled empirical formula. Minor axis = 0.6 × major. Ellipse at 70° bearing for prevailing jet stream.
Population Exposure
N = Σ GPW(x,y) × A(x,y)
Sum of GPW v4 population density × cell area over all grid cells within impact footprint. Zone mortality rates applied per ring.

⚠ Disclaimer: The Pole Shift / ECDO scenarios are theoretical models presented for educational and exploratory purposes, clearly labeled as theoretical throughout the application. They are reproduced as presented by Ben Davidson (Suspicious Observers) and The Ethical Skeptic. Disaster Map does not endorse or refute these theories. All other scenarios use peer-reviewed physics and verified terrain data. Casualty estimates are rough educational approximations only — not suitable for emergency planning, insurance, or engineering use.

Simple Pricing

Start free with all 9 scenarios. Upgrade to explore without limits.

Free
$0forever
  • All 9 disaster scenarios
  • 10 simulations / hour
  • 30 simulations / day
  • Browse all 1,215 volcanoes
  • Supervolcano eruption sims
  • Earthquake + fault line viewer
  • Standard map view
  • Casualty estimates
Launch Free
FOUNDERS
Pro Lifetime
$18.99one-time · before price increase
  • All 9 disaster scenarios
  • 50 simulations / hour
  • 200 simulations / day
  • Globe + Satellite view
  • Full pan, zoom & explore
  • Active fault lines overlay
  • Simulate any volcano eruption
  • Price going up to $24.99 — lock in now
Free vs Pro
30 simulations / day
200 simulations / day
Asteroid up to 5,000 m
Asteroid up to 20,000 m
Nuke up to 1 Mt
Nuke up to 100 Mt
Browse volcanoes free
Simulate any eruption
Eq. intensity rings
Active fault lines overlay
Map view only
Satellite + Globe view
1 zone click popup
Unlimited zone popups
Displaced population count
🌊 Mega-Tsunami
Lock in Founders Price
Developer Kit
Soon
  • Self-hosted backend
  • All scenario APIs
  • Your Mapbox key
  • Your server, your costs
  • Single license
Coming Soon

Frequently Asked Questions

Ready to Simulate
the Apocalypse?

Free to start. No account required.
Drop an asteroid, flood a continent, or watch a pole shift unfold.

Launch Disaster Map →