This shows exactly what the proxy returned for the selected CME before any enrichment.
Use it to check which DONKI fields are available: time21_5,
enlilList, isMostAccurate,
source signatures, multiple analyses, etc.
Select a CME first, then open this inspector.
This tool is an advanced educational and situational-awareness simulator combining NASA DONKI observations, NOAA real-time solar-wind data, WSA–ENLIL model comparison, and a locally computed numerical propagation model with six physics improvements added in v3.0. It is not an official warning product — defer to NOAA SWPC for operational alerts.
The simulator separates measured inputs, model-derived propagation, and scenario estimates. v3.4 corrects geometry and shock-physics errors found in an external mathematical audit of the v3.3 debug report; every quantity in the report is now reproducible from the inputs printed beside it.
The propagation backbone is a Drag-Based Model integrated with 4th-order Runge-Kutta (RK4). CMEs faster than the ambient solar wind decelerate; slower disturbances accelerate toward the background flow. The equation of motion is:
a = −γ(r) · (v − w(r)) · |v − w(r)|
v is CME speed, w(r) is the distance-dependent ambient solar-wind speed, and γ(r) is a drag parameter that varies with heliocentric distance. Three specific improvements over standard constant-γ DBM:
Dynamic drag γ(r,t) — v3.3 computes drag from first principles at each RK4 step using CME cross-sectional area, live-calibrated solar-wind density, estimated CME mass, virtual mass loading, and a speed-dependent drag coefficient. Wider low-mass CMEs couple harder to the solar wind; massive/narrow CMEs coast longer.
Density-calibrated reference (v3.0) — the 1 AU density reference in the drag calculation is now calibrated against the live NOAA SOLAR-1 L1 measurement rather than a fixed 5 cm⁻³ assumption. During CME sheath passages when density reaches 14–20 p/cc, drag increases by up to 3× relative to quiet-sun conditions. Clamped to [0.3×, 3×] to avoid pathological transient spikes.
Slow-CME injection delay (v3.0) — gradual eruptions (filament-associated, v < 600 km/s) continue accelerating beyond 21.5 R☉ (Zhang & Dere 2006). Using 21.5 R☉ as the universal injection point with the coronagraph-measured speed overestimates transit for slow CMEs. Correction: add 0–3 hours to the injection start time, linearly scaled as (600−v)/600 × 3h. Zero for v ≥ 600 km/s.
Standard DBM predicts when the CME nose reaches 1 AU. When Earth sits at angular offset θ from the CME axis it intercepts the side of the expanding front, which arrives later and moves more slowly.
What changed in v3.4. Versions up to v3.3 stacked a Vršnak cosine projection t_nose(1/cosθ − 1) underneath an ellipsoid shape term and added the two together. That sum had no closed-form justification, could not be reproduced from the printed report, and systematically under-predicted the delay. It has been replaced by the standard self-similar circular ("ice-cream cone") front, solved exactly.
R_apex / d = (1 + sin λ) / (cos Δ + √(sin²λ − sin²Δ))
A sphere of radius r = h·sinλ centred at distance h on the cone axis, apex at R = h(1+sinλ). Solving |d·û − h·n̂| = h sinλ gives the relation above: the radial distance the apex has reached at the moment the flank first touches a target at distance d and off-axis angle Δ.
Checks: Δ=0 gives R_apex = d (direct hit, zero delay). Δ=λ gives R_apex = d(1+sinλ)/cosλ — finite, no divergence at the cone edge. Δ>λ has no solution: Earth is outside the cone.
The delay is then the DBM travel time for the apex to cover R_apex − 1 AU, using the same drag parameter and wind as the main integration — not a geometric multiplier applied to the nose transit time.
v_front(Δ) = v_apex · d / R_apex
Because the front is self-similar in h, the local radial speed where it touches the target follows directly. A flank encounter is always slower than the apex. The simulator now uses this local front speed for shock formation, impact speed and geomagnetic coupling; earlier versions reused the nose speed, which was not geometrically self-consistent.
Not modelled: front flattening / pancaking. An axis ratio is computed and printed in the debug report as a diagnostic, but it is not applied to the arrival time. Treat wide, fast halos as having additional unmodelled delay uncertainty.
When a faster trailing CME launches within 18 hours and is within 35° of a leading CME, it will overtake and cannibalize the slower leader. The merged structure arrives earlier than either CME would independently. Based on Lugaz et al. (2012), the merged arrival time is approximated as a kinetic-energy-weighted mean:
t_merged = (v₁² · t₁ + v₂² · t₂) / (v₁² + v₂²)
where v₁, t₁ are the leading CME's speed and arrival time, and v₂, t₂ are the trailing CME's. The v² weighting reflects the kinetic energy of each ejecta. The leading CME's arrival is shifted earlier by the difference; the correction is only applied when this produces an earlier arrival (physically: the faster CME must actually be catching up).
The original interaction risk warning remains displayed in the event card regardless of whether the arrival shift is applied. This improvement complements the existing flag by making a quantitative timing correction rather than just displaying a warning.
Kp and Dst estimates are scenario outputs, not deterministic forecasts. The simulator uses CME speed, density, shock/sheath structure, impact geometry, and live solar-wind context to estimate plausible storm ranges via a 200-trial Monte Carlo ensemble.
Rotation-sense Bz bias (v3.0) — flux rope magnetic rotation sense depends on which flank Earth intercepts. From Bothmer & Schwenn (1998) and Mulligan (1998):
Confidence in the rotation-sense prediction degrades linearly from the disk center (θ=0°, full weighting) to the cone edge (θ=45°, 50/50 — no geometric information remaining). This shifts the Bz probability distribution asymmetrically rather than using the symmetric bimodal sampling of prior versions, improving Kp/Dst scenario accuracy for off-axis encounters.
The proxy backend continuously scans live NOAA SOLAR-1 solar-wind data for interplanetary shock signatures and matches them against model predictions, building a persistent forecast-accuracy log.
Shock detection criteria (all three within 45 minutes):
Detected shocks are matched to the nearest ENLIL-predicted CME arrival (within 36 hours). The predicted-vs-actual delta in hours is logged persistently. Peak Kp and Dst observed in the 72 hours following detection are tracked and updated continuously, enabling retrospective forecast accuracy assessment. A background cron job runs detection every 15 minutes independently of site visitors.
When a confirmed arrival matches the currently selected CME, a green badge appears in the event card and a toast notification fires. The full log is visible in the Arrival Verification section below the event list.
Different model components have very different confidence levels. The table below reflects v3.0 improvements — arrival time confidence is improved by the flank correction and density calibration, geomagnetic confidence is modestly improved by the rotation-sense bias.
| Component | Confidence | Reason |
|---|---|---|
| CME speed | ★★★★★ | Directly supplied by DONKI analysis. Projection corrected by the v3.4 cone deprojection; a GCS or stereoscopic reconstruction would be better where available. |
| CME direction | ★★★★☆ | Good when latitude/longitude and source region are supplied. PFSS deflection applied. |
| Arrival time — direct hit | ★★★☆☆ | Density-calibrated RK4 + injection delay. Improves with ENLIL comparison. |
| Arrival time — flank | ★★★☆☆ | Flank correction (v3.0) reduces error significantly. Validated: 3.6h → 1.3h for AR 4479 CME. |
| Earth-hit geometry | ★★★☆☆ | Cone geometry is well-constrained; flank boundary uncertainty remains. |
| CME–CME merged arrival | ★★☆☆☆ | Interaction shift (v3.0) improves cases where cannibalization is clearly detected. |
| Kp scenario — direct hit | ★★☆☆☆ | Depends on Bz, which remains unknown before L1 arrival. |
| Kp scenario — flank | ★★☆☆☆ | Reported as four explicit Bz scenarios with priors. Rotation sense informs storm timing only. |
| Dst scenario | ★★☆☆☆ | Ring-current response sensitive to sustained southward Bz duration. |
| Bz orientation | ★☆☆☆☆ | Cannot be predicted from coronagraph geometry. This is a fundamental physics limit, not a model deficiency. |
CME catalog, source location, active region, all analyses, 21.5 R☉ timing, linked events (flares, SEPs, IPS, GST), and ENLIL impact records.
Real-time solar wind speed, density, IMF Bz/Bt from the new SWFO-L1 spacecraft that replaced DSCOVR as primary L1 monitor in 2026. Data delivered within 5 minutes.
Live geomagnetic index comparison for forecast-vs-observed validation as the event arrives at Earth.
Official 3D MHD heliospheric model comparison. Predicted shock arrival, Rmin, Kp at three IMF orientations, and per-spacecraft impact list when available.
Solar flare classification and energetic proton flux for contextual event linking and sheath hazard assessment.
Precise planetary ephemeris for accurate real-time positions of Earth, Mercury, Venus, and deep-space spacecraft in the heliocentric views.
The simulator is strongest at showing geometry, timing envelopes, and uncertainty ranges. It is weakest at predicting internal magnetic orientation (Bz) before L1 arrival — this is a fundamental physical limit shared by all coronagraph-based models including WSA–ENLIL.
The density calibration uses instantaneous L1 density, not the upstream density profile the CME actually encountered during transit. Ensemble spread is within-model only: it varies speed, wind and half-angle but does not represent structural model error, so it should not be read as a calibrated probability. Full operational forecasts require expert interpretation and models that include the structured background heliosphere.
v_r = v_POS(1+sin ω)/(sin ρ + sin ω).