This clock has not been scientifically validated. The latest application and JSON dataset underwent an AI-assisted technical and methodological audit using Astra ULTRA. This audit is not independent peer review or predictive validation.
Mission statement

Apocalypse Clock is a transparent, scenario-based systemic-risk model. Its purpose is to make interacting global risks comparable within a single analytical framework: it examines 23 global threats across civilizational, biospheric, and technological domains, not as isolated dangers, but as coupled components of an interconnected global system.

Its parameter data are drawn primarily from official governmental reports, intergovernmental and international-organization assessments, scientific literature, and evidence-based institutional datasets. Each parameter is documented through source references, uncertainty ranges, evidence strength, growth calibration, and explanatory notes.

The model estimates a probabilistic critical horizon through a multi-stage analytical pipeline: evidence-graded parameter scoring, Weighted MCDA, dependency amplification, domain-weight normalization, process-specific horizon functions, Monte Carlo uncertainty sampling, bootstrap uncertainty estimation, structural ensemble aggregation, and Dynamic cascade propagation.

The ensemble compares compensatory aggregation, non-compensatory max-rule aggregation, graph-weighted aggregation, and Dynamic cascade logic. The paired highlighted years represent the P50 median and P90 later edge of the model’s Dynamic cascade horizon. Each displayed year is a model-derived threshold horizon under explicit assumptions: a stress-oriented critical threshold, not a deterministic forecast, a prophecy, or an empirically measured probability of collapse.

The optional Model Diagnostics adds deeper diagnostic tools, including Weibull survival analysis, network eigenvector centrality, Poisson-binomial convergence-tail analysis, Shannon entropy, Fast OAT sensitivity analysis, Sobol/Jansen sensitivity indices, SMAA weight robustness, non-compensatory veto diagnostics, tail-dependence stress testing, and a scientific audit summary. These diagnostics explain and stress-test the baseline result; they do not overwrite the paired headline horizons.

The expanded view below breaks that horizon into the model steps behind it: threat scoring, dependency links, process-specific horizons, uncertainty sampling, aggregation rules, and cascade diagnostics.

1. Threat scoring
Each of the 23 threats is encoded through six MCDA dimensions, an effective growth rate, and a destabilization threshold, producing a comparable systemic-risk profile.
2. Interaction & horizons
The model applies dependency amplification and maps threats into process-specific horizons, distinguishing continuous degradation, event-like shocks, and regime-shift dynamics.
3. Uncertainty engine
Monte Carlo simulation, Beta/log-normal sampling, bootstrap intervals, and P10/P50/P90 percentiles describe uncertainty instead of relying on a single deterministic date.
4. Cascade & diagnostics
Dynamic cascade testing and structural ensemble aggregation determine the paired headline horizons. The Model Diagnostics then adds sensitivity, robustness, stress-test, and distribution diagnostics as an interpretive audit layer, without changing those horizons.

Results remain sensitive to thresholds, dependency strength, domain weights, growth assumptions, and aggregation rules. Optional diagnostics make that sensitivity visible; they are interpretive tests, not a second headline forecast.

Live calculation console
Idle
Waiting for a full model pass. Each listed stage receives a check mark and a calculation-type badge when it finishes.
Projected horizon
Dynamic cascade · paired P50/P90 horizon clocks · projected pathway
Model horizon markers
Compensatory P50 threshold horizon:
Compensatory P10 horizon
Compensatory P50 horizon
Compensatory P90 horizon
Max-rule P50
Graph heuristic P50
Dynamic cascade P50 horizon
Dynamic cascade P90 horizon
Top threat P50
Scenario-conditioned probability of threshold crossing under the compensatory aggregation rule
P ≤ 2035
P ≤ 2050
Horizon-coded outcome SD
P50 rule range
Structured heuristic under deep uncertainty
Not a validated forecast of collapse, extinction or catastrophe
▼
The clock now runs against the embedded data_v1_9_0.json data map as its default and primary parameter source. AI presets and All-AI Average are secondary comparison/audit inputs unless deliberately selected. You can upload either a full source map with {mu, lo, hi} ranges or an evidence overlay with {obs, weight} observations to update the current parameter ranges client-side.
Bundled source map active
AI Preset
Selected AI
JSON Open JSON dataset
Active file
data_v1_9_0.json
Parameter entries
Threat coverage
Mode
Bundled
Bundled data_v1_9_0.json parameter map is embedded directly in this widget and remains the default primary source unless a custom source map, evidence overlay, or AI preset is deliberately selected.
Inspect active parameter map
Loading source map…
On-demand diagnostics
Diagnostic calculations
Optional checks of the finished main run: sensitivity, weight robustness, stress tests and an audit. They run only on request and do not change the clocks.
Not run yet.
Additional Sensitivity Graphs OAT + Sobol ▼ Hide
Sensitivity Analysis OAT + Sobol/Jansen
2 Independent Sensitivity Screens
Compare one-at-a-time sensitivity with interaction-aware sensitivity. These views show how individual threats and connected threat combinations move the current Apocalypse Clock result. Read them as diagnostic tests of model behavior, not direct rankings of real-world danger.
FAST SCREEN
One threat 20% weaker or stronger
INTERACTION SCREEN
Direct bars plus interaction bars
Fast OAT Diagnostic
One-at-a-Time Sensitivity
This chart shows how much the Apocalypse Clock result changes when one threat is made 20% weaker or 20% stronger, while all other threats stay the same. Longer bars mean the model is more sensitive to that threat. Blue bars show a decrease, red bars show an increase. This is a diagnostic test of model behavior, not a direct ranking of real-world danger.
◀ Left bars reduce the response
Right bars ▶ increase it
Sobol Jansen Indices Interaction-Aware
Interaction-Aware Sensitivity
This chart shows which threats contribute most to uncertainty in the Apocalypse Clock result when many possible combinations are tested. The first-order bar shows a threat's direct effect. The total-order bar includes both its direct effect and its interaction with other threats. This helps reveal which risks matter not only alone, but also as part of a connected systemic-risk network.
Exploratory only: Negative indices appeared due to estimator noise.
Experimental Diagnostic Graphs SMAA Veto Tail Shock Distribution Audit ▼ Hide
On-demand experimental layer robustness stress tests distributions audit
Experimental scientific diagnostics
These diagnostics are intentionally excluded from the main clock run. They test weight robustness, non-compensatory domain failure, correlated polycrisis tail shocks, and the shape of the already-computed Monte Carlo distribution. Read them as diagnostic and stress-test views, not as the headline forecast.
RUN MODE
Idle use the Model Diagnostics button
Run the main simulation first, then use the Run diagnostics button above to populate this experimental layer.
Veto stress test
Veto Stress Test
This chart checks whether one major risk domain can cross a high internal stress boundary even when the overall average remains below the main threshold. It prevents severe pressure in one area from being hidden by lower stress elsewhere. Blue markers show the median crossing year. Gold markers show a later high-end estimate. The test is diagnostic, not predictive.
Tail-dependence shock
Tail-Dependence Shock
This chart tests what happens when several connected threats worsen together. It compares a baseline cascade with scenarios where linked risks receive a shared stress shock. If the markers move left, the modeled crossing happens earlier. This helps show whether the Clock is sensitive to clustered systemic pressure rather than only isolated threat changes.
SMAA weight robustness
Weight Robustness
This chart tests whether the main driver remains stable when the model weights are slightly changed. A high percentage means that the same threat remains the top-ranked driver across many alternative weighting runs. It does not mean probability of collapse. It means the ranking is robust under weight uncertainty.
Scientific audit
Run identity and module status
Compact audit of which optional scientific modules were executed and with what sample sizes.
No experimental audit yet.
Probabilistic diagnostics histogram boxplot heatmap CDF
Distribution layer for the current Monte Carlo result
These graphs visualize the already-computed baseline Monte Carlo result: crossing-year density, structural-rule boxplots, domain percentile heatmap, and ensemble CDF curves.
SOURCE
Completed baseline simulation
Distribution Histogram
Crossing-Year Histogram
This chart counts crossings observed by 2100. Runs with no crossing are reported separately as censored, not assigned an event year. Reference lines show full-run P10, P50, and P90 only when identifiable within the horizon.
Structural Boxplot
Aggregator Uncertainty
This chart compares aggregation rules using crossing times conditional on crossing by 2100. Each rule's censored share is reported separately; boxes do not describe unidentified post-horizon times or unconditional quantiles when censoring is present.
Domain Heatmap
Domain Percentile Map
This chart shows P10, P50, and P90 first crossings of the within-domain functional-loss trigger after propagation through the full threat network. Each domain uses the headline model’s criticality, overlap, essential-service and threshold rules, including each run’s sampled structure, with a domain-specific denominator. The domains are not additive, and the years are model percentiles rather than completed-collapse dates or deterministic forecasts.
CDF
Cumulative Crossing Curves
This chart shows how the probability of model crossing accumulates over time. Curves that rise earlier indicate earlier crossing under that aggregation rule. Steeper curves mean outcomes are more concentrated. Flatter curves mean greater timing uncertainty.
Model Diagnostics How the model works

How the Apocalypse Clock model works

This panel separates the model into two computational layers. The core run produces the displayed horizon from evidence-informed scoring, dependency amplification, process-specific horizon functions, Monte Carlo uncertainty sampling, structural ensemble aggregation, and dynamic cascade testing. The Model Diagnostics adds an optional diagnostic layer after the baseline result exists. It visualizes sensitivity, weight robustness, non-compensatory stress rules, correlated tail shocks, and the distributional structure of the current Monte Carlo result. These diagnostics are for auditability and interpretation under deep uncertainty; they do not overwrite the headline clock date.

Interpretation rule All years are model-derived critical horizons under explicit assumptions. The Model Diagnostics audits and stress-tests the baseline result; it does not change the paired headline horizons.
01
Evidence inputs
Evidence-graded estimates define central values, uncertainty ranges, source strength and notes for each threat dimension.
02
Weighted MCDA
Six dimensions are weighted into a baseline systemic score: scale, urgency, acceleration, interdependence, irreversibility and governance failure.
03
Dependency amplification
Threats are treated as coupled system components. Network links amplify priorities when threats reinforce one another.
04
Horizon functions
Different threat types use different time-to-threshold functions for continuous, event-like and regime-shift risks.
05
Monte Carlo
Input uncertainty is sampled repeatedly to produce distributions of possible crossing years rather than one fixed date.
06
Structural ensemble
Multiple aggregation rules are compared: compensatory, max-rule, graph-weighted and dynamic cascade.
07
Dynamic cascade
Directed functional losses propagate when pressure plus loss of supporting capacity reaches a target threshold. A single seed may spread; the aggregate uses global or essential-service criticality shares.
08
Diagnostic trigger
After the core result is complete, the Model Diagnostics can run sensitivity, robustness, veto, tail-shock and distribution diagnostics on demand.
09
Visual audit layer
Histogram, boxplot, heatmap and CDF views expose the distributional shape of the already-computed Monte Carlo and structural ensemble results.
10
Headline horizon
The paired headline clocks remain the P50 median and P90 later edge of the core Dynamic cascade distribution. Diagnostics may show alternative stress views, but they do not replace this rule.
Technical interpretation of each layer
Model core

Weighted MCDA

Each threat receives six dimension scores. The model combines them into the same weighted baseline score used by the ranking code, then applies scenario, domain weighting and dependency effects.

baseᵢ = wSc·Scᵢ + wUr·Urᵢ + wAcc·Accᵢ + wInt·Intᵢ + wIrr·Irrᵢ + wGov·Govᵢ
Network layer

Dependency amplification

Threats are not treated as independent. A threat can amplify another through declared dependency links, producing the same adjusted priority value used for ranking and horizon calculation.

Pᵢ = baseᵢ × depFactorᵢ × domainMultiplierᵢ
Time layer

Process-specific horizons

The model does not convert all risks into time the same way. Continuous, event-like and regime-shift threats use the process-specific horizon functions implemented in the engine.

continuous: Tᵢ = NOW + ln(θᵢ/Pᵢ)/ln(1+gᵢ)
event median: λᵢ(y)=clamp(2gᵢPᵢ/θᵢ·(1+gᵢ)ʸ⁻ᴺᴼᵂ,1e-5,5), Σλᵢ ≥ ln2
regime median: pᵢ=clamp(logistic(2(Pᵢ-θᵢ)),1e-4,.999), Tᵢ=NOW+ln(.5)/ln(1-pᵢ)
Uncertainty layer

Monte Carlo simulation

Uncertain inputs are sampled repeatedly. The output is summarized through P10, P50 and P90 rather than a single deterministic year.

P10 P50 P90 = percentiles of simulated crossings
Aggregation layer

Structural ensemble

Different aggregation rules answer different questions: weighted mass crossing, earliest single threat, graph-linked joint stress and dynamic cascade propagation.

ensemble = {compensatory, max-rule, graph-weighted, dynamic cascade}
Cascade layer

Dynamic cascade test

A single initial functional failure can spread through directed dependencies. The trigger uses criticality weights globally or within an essential-service basket, and each run samples the threshold and the other structural choices; neither extinction of all organisms nor initial failure in every domain is required.

target: qᵢ(t) + κᵢΣⱼaᵢⱼAⱼ(t−lagᵢⱼ) ≥ 1; trigger: max(global share, essential-service shares) ≥ cascadeThreshold
On-demand diagnostics

Sensitivity + weight robustness

OAT shows local one-at-a-time leverage. Sobol/Jansen estimates direct and total variance contribution. SMAA perturbs weights to test whether top-ranked drivers remain stable.

S1 = direct effect ST = total effect SMAA = rank robustness
Stress-test layer

Stress tests + distribution views

Veto and tail-dependence diagnostics test non-compensatory and correlated-stress alternatives. Histogram, boxplot, heatmap and CDF views make the baseline uncertainty structure visible.

diagnostic result ≠ headline result
Top threat cards will update after the model run.
Global Stress Index
0–100 composite
P(threshold ≤ 2035)
Near-term scenario summary
P(threshold ≤ 2050)
Mid-century scenario summary
Central 80% interval width
Model P90 − P10 width
Dominant driver
Current priority summary
Risk horizons
Domain indices & top 5 leading threats
Horizon source
Distant / watch / critical horizon band
Monte Carlo probability band
Model-generated interval
■ Civilizational
■ Biospheric
■ Technological
Causal dependency network
How the 23 threats worsen each other, from the research register of 25 September 2026 (101 links, plus one added link). Click a threat to see what it worsens directly (step 1) and in turn (step 2). Node size is the adjusted priority of the selected scenario. The clocks are still calculated with the dataset's own declared dependencies.
Threat contribution ranking ▼
Full threat register (23 threats)
Displayed P10, P50 and P90 years are model-derived Monte Carlo horizons computed from evidence-informed scenario assumptions. They are not direct empirical forecasts, observed dates, or literature-reported collapse years. Scientific evidence is read from the explicit threat-level assessment in the active JSON. Input calibration summarizes the eight parameter-level grades; U means not supplied.
# Threat Sc Ur Acc Int Irr Gov Base score Adj. score Model Ti Model P10–P90 Effective risk-growth proxy Model threshold Scientific evidence Input calibration Mechanism