Causal-prior-constrained deep learning for rapid-intensification warning of western North Pacific tropical cyclones
Organizations: Department of Atmospheric Sciences, School of Earth Sciences, Zhejiang University, Hangzhou, 310058, China · School of Mathematical Sciences, Zhejiang University, Hangzhou, 310058, China
Abstract
Rapid intensification (RI) is one of the most dangerous and unpredictable phases of tropical cyclones, and its environmental and structural drivers are difficult to identify and separate in observations. Here we introduce CausalVisionRI, a causal-prior-constrained deep-learning system for 24-h RI warning in the western North Pacific. It integrates infrared remote sensing imagery, atmospheric and oceanic fields, and an observationally discovered causal graph into a unified model. The graph and eight sign-verified effects on 24-h intensity change provide an auditable causal prior for training. Discovered from 522 cyclones in 2000-2017, it identifies causal pathways linking atmospheric and oceanic conditions, the cyclone's dynamical and structural state, and its 24-h intensity change. On a forward holdout covering 2020-2024 and 3,004 records, the five-member ensemble achieves a Brier skill score (BSS) of 0.259 relative to a fixed development climatology, significantly exceeding four reference systems; the three LDA baselines reach at most 0.105. It also attains higher detection rates than a published knowledge-block scheme and a reconstructed contrastive-learning scheme at matched false-detection rates under our USA 30-kt standard and the published RITCF standard. Ablations show that the graph pathway adds 0.032 BSS across the holdout and 0.063 over 2020-2021, with a significant contribution in one ensemble member. A deep-learning model constrained by observationally discovered causal structure thus offers a new route to physically informed RI warning.
Figures & tables
| Model | BSS dev | ROC-AUC | PR-AUC | POD | POFD | FARatio | CSI | BSS [90% CI] |
|---|---|---|---|---|---|---|---|---|
| CVRI-Ens | 0.259 | 0.900 | 0.518 | 0.680 | 0.095 | 0.635 | 0.311 | – |
| SHIPS-12 LDA (COHC) | 0.105 | 0.820 | 0.276 | 0.347 | 0.060 | 0.686 | 0.197 | [ , ] |
| SHIPS-12 LDA (OHCD) | 0.095 | 0.813 | 0.266 | 0.473 | 0.095 | 0.715 | 0.216 | [ , ] |
| LDA (21 panel predictors) | 0.098 | 0.813 | 0.278 | 0.730 | 0.250 | 0.811 | 0.176 | [ , ] |
| LR (uncalibrated) | 0.826 | 0.298 | 0.649 | 0.173 | 0.769 | 0.205 | [ , ] |
| Variable | Meaning | Construction | Source |
| Best-track state | |||
| VMAX | Current intensity | USA 1-min sustained wind (kt) | IBTrACS |
| PER | Recent intensification | VMAX( ) VMAX( h) | IBTrACS |
| SPD | Translation speed | Great-circle displacement rate between best-track positions (kt; nmi h -1 ) | IBTrACS |
| LAT | Latitude | Storm-centre latitude | IBTrACS |
| LAND_DIST | Distance to land | Great-circle distance (km) | GSHHG coastline |
Appendix figures & tables16 assets
Supplementary material from the paper’s appendix.
Appendix
| Treatment | Outcome | Effect | CI low | CI high | E-value | E-value CI limit | Refut. | Adjustment set | |
|---|---|---|---|---|---|---|---|---|---|
| SHRD | DV24 | 4.788 | 5.387 | 4.147 | 1.789 | 1.705 | 0.005 | 3/4 | VORT |
| SHRD | RI | 0.729 | 0.872 | 0.595 | 3.386 | 2.891 | 0.005 | 2/4 | VORT |
| RHLO | DV24 | 2.791 | 2.045 | 3.485 | 1.525 | 1.423 | 0.005 | 2/4 | W500 |
| RHLO | RI | 0.112 | 0.003 | 0.211 | 1.453 | 1.059 | 0.040 | 1/4 | W500 |
| VORT | DV24 | 3.825 | 4.356 | 3.237 | 1.662 | 1.585 | 0.005 | 3/4 | LAT |
| VORT | RI | 0.152 | 0.277 | 0.041 | 1.567 | 1.240 | 0.010 | 2/4 | LAT |
| Rank | Arm | Family | Threshold | BSS 2018 | BSS 2019 | Worst year | Selected |
|---|---|---|---|---|---|---|---|
| 1 | ens_lbB_partial4_sf | lingbot_sf | 0.150 | 0.260 | 0.314 | 0.2604 | yes |
| 2 | ens_lbB_lora16_temp5 | lingbot_mf | 0.250 | 0.383 | 0.247 | 0.2466 | yes |
| 3 | ens_lbB_temp2 | lingbot_mf | 0.225 | 0.322 | 0.231 | 0.2306 | yes |
| 4 | ens_lbB_lora16_sf | lingbot_sf | 0.200 | 0.268 | 0.227 | 0.2275 | yes |
| 5 | ens_lbB_temp5_slim | lingbot_mf | 0.250 | 0.299 | 0.227 | 0.2266 | yes |
| 6 | ens_lbB_lpft4_sf | lingbot_sf | – | 0.256 | 0.219 | 0.2195 | no |
| Year | n | BSS dev | POD | POFD | FARatio | CSI |
|---|---|---|---|---|---|---|
| CVRI-Ens (Top-5) | ||||||
| 2020 | 559 | 0.276 | 0.674 | 0.097 | 0.633 | 0.312 |
| 2021 | 757 | 0.298 | 0.659 | 0.073 | 0.642 | 0.302 |
| 2022 | 571 | 0.159 | 0.639 | 0.099 | 0.697 | 0.258 |
| 2023 | 606 | 0.287 | 0.660 | 0.076 | 0.545 | 0.368 |
| 2024 | 511 | 0.252 | 0.761 | 0.142 | 0.653 | 0.312 |
| Ensemble | K | n | BSS dev | ROC AUC | PR AUC | POD | POFD | FARatio | CSI |
|---|---|---|---|---|---|---|---|---|---|
| K4 (Top-4) | 4 | 3,004 | 0.259 | 0.898 | 0.513 | 0.667 | 0.087 | 0.621 | 0.319 |
| Top-5 (main) | 5 | 3,004 | 0.259 | 0.900 | 0.518 | 0.680 | 0.095 | 0.635 | 0.311 |
| K6 (Top-6) | 6 | 3,004 | 0.279 | 0.902 | 0.527 | 0.662 | 0.076 | 0.589 | 0.339 |
| all-22 | 22 | 3,004 | 0.244 | 0.905 | 0.509 | 0.743 | 0.128 | 0.684 | 0.285 |
| Member | Condition | BSS | 90% CI | Significant | |
|---|---|---|---|---|---|
| Part4 | C1 | 0.093 | 0.160, 0.038 | 0.024 | yes |
| Part4 | C2 | 0.014 | 0.079, 0.053 | 0.681 | no |
| Part4 | C3s1 | 0.076 | 0.148, 0.017 | 0.174 | no |
| Part4 | C3s2 | 0.071 | 0.142, 0.014 | 0.174 | no |
| Part4 | C3s3 | 0.059 | 0.142, 0.007 | 0.292 | no |
| Part4 | C4adj | 0.072 | 0.135, 0.018 | 0.174 | no |
| Member | gate mean | causal share | share 0.5 | gate share | gate RI | gate non-RI | gate |
|---|---|---|---|---|---|---|---|
| Part4 | 0.629 | 0.689 | 0.189 | 0.426 | 0.665 | 0.626 | 0.039 |
| T2 | 0.620 | 0.299 | 0.201 | 0.188 | 0.677 | 0.615 | 0.062 |
| T5 | 0.488 | 0.260 | 0.240 | 0.141 | 0.506 | 0.487 | 0.020 |
| LoRA-T5 | 0.698 | 0.251 | 0.249 | 0.175 | 0.666 | 0.701 | 0.035 |
| LoRA-sf | 0.438 | 0.230 | 0.270 | 0.091 | 0.414 | 0.440 | 0.026 |
| Sample definition | n | Use |
|---|---|---|
| Official test set (main table) | 3,004 | All headline numbers |
| WRPS1 same-sample subset | 2,892 | WRPS1 delta-BSS bootstrap |
| RITCF reconstruction, CMA 30-kt frame | 2,385 | RITCF CMA 30-kt protocol |
| Panel A intersection | 2,178 | Consistent-comparison panel, our label standard |
| 2020–2021 subset | 1,316 | Two-year supplementary contrast |
| RITCF as published | 1,149 | Their reported operating point (external anchor) |
| Panel | n | RI USA | RI CMA | Rate USA | Rate CMA | Jaccard | Agreement |
|---|---|---|---|---|---|---|---|
| A (our standard) | 2,178 | 209 | 120 | 0.096 | 0.055 | 0.530 | 0.954 |
| B (their standard) | 949 | 83 | 66 | 0.087 | 0.070 | – | – |
| # | Predictor | Definition (radius) | Source in this study | Fidelity |
|---|---|---|---|---|
| 1 | VMAX | Maximum sustained wind (kt) | Panel VMAX (IBTrACS USA_WIND, 1-min) | EXACT |
| 2 | INCV | 6-h intensity change, t 6 h to t0 (kt) | Panel PER (6-h IBTrACS track difference) | EXACT |
| 3 | POT | Potential intensity minus VMAX (kt), 0–200 km | tcpyPI 1.4.1 [ 2 ] on ERA5 soundings, 0–200 km; Bister–Emanuel algorithm [ 1 ] | APPROX |
| 4 | OHCD | Ocean heat content (kJ cm -2 ), 0–100 km | GLORYS12V1 daily field, closed window (source_time 24 h t0) | BETTER |
| 5 | SHRD | 850–200 hPa shear magnitude (kt), 200–800 km | ERA5 u/v annular mean difference | EXACT |
| 6 | T150 | 150 hPa temperature ( ∘ C), 200–800 km | Log-p extrapolation from 200/300 hPa (no local 150 hPa level) | SUBSTITUTED |
| Assumption | Handling here | If violated |
|---|---|---|
| Causal sufficiency (no unobserved confounder) | All candidate drivers drawn from reanalysis and best-track fields; E-values quantify the confounding strength needed to explain each effect away (Table S1) | Estimates may be biased in either direction; E-values give the minimum confounding strength required to explain an effect away (DV24: 1.13–2.08; RI: 1.33–3.39) |
| Stationarity of links | Discovery on the 2000–2017 fit window only; four sensitivity configurations re-estimated (Note S2) | Within-window sensitivity does not test stationarity; the configurations probe lag selection, climate indices, latitude and OHC dependence (Note S2), and alternative DAG constructions are audited separately (Note S1) |
| Faithfulness (observed independencies reflect the graph) | Masked tests with mask_type=xyz ; BH control at for lagged links and unadjusted for contemporaneous links | Weak-signal links may be missed (e.g., MFC), limiting detection power rather than demonstrating unfaithfulness |
| Temporal ordering | Six-hourly storm-aligned grid; target never interpolated | Lagged-link orientations rely on the six-hourly ordering; lag-0 links are handled separately (flagged or removed in the strict configuration) |