Every wildfire prediction model deployed today shares a dangerous property: none of these methods provides formal guarantees on how much fire spread is missed. Despite extensive work on wildfire spread prediction using deep learning, no prior study has applied distribution-free safety guarantees to this domain, leaving evacuation planners reliant on probability thresholds with no formal assurance. We address this gap by presenting, to our knowledge, the first application of conformal risk control (CRC) to wildfire spread prediction, providing finite-sample guarantees on false negative rate (FNR <= 0.05). We expose a stark failure: across three model families of increasing complexity (tabular: LightGBM, AUROC 0.854; convolutional: Tiny U-Net, AUROC 0.969; and graph-based: Hybrid ResGNN-UNet, AUROC 0.964), standard thresholds capture only 7-72% of true fire spread. CRC eliminates this failure uniformly. Our central finding is that model architecture determines evacuation efficiency, while CRC determines safety: both spatial models with CRC achieve approximately 95% fire coverage while flagging only approximately 15% of total pixels, making them 4.2x more efficient than LightGBM, while the graph model's additional complexity over a simple U-Net yields no meaningful efficiency gain. We propose a shift-aware three-way CRC framework that assigns SAFE/MONITOR/EVACUATE zones for operational triage, and characterize a fundamental limitation of prevalence-weighted bounds under extreme class imbalance (approximately 5% fire prevalence). All models, calibration code, and evaluation pipelines are released for reproducibility.
Figures & tables
Figure 1: Why standard thresholds fail. A cross-section through a U-Net fire prediction shows probability peaking at ∼0.8 and decaying smoothly. The standard threshold ( p^≥0.5 , gray) misses the majority of the fire region; the CRC threshold ( λ^=0.002 , green) captures it entirely.
Figure 2: FNR as a function of decision threshold for all three models. The red dashed line marks α=0.05 . Both spatial models maintain safe FNR over a much wider threshold range than LightGBM, explaining their 4.2 × smaller evacuation zones under CRC. The near-overlapping U-Net and ResGNN-UNet curves illustrate diminishing returns from architectural complexity. Markers show each model’s CRC λ^ .
Category
Feature
Variable
Unit
Topography
elevation
Terrain height
m
Weather
th
Wind direction
∘
vs
Wind speed
m/s
tmmn
Min temperature
K
tmmx
Max temperature
K
sph
Specific humidity
kg/kg
Table 1: NDWS input features (12 channels per 64×64 pixel patch). All features are from remote-sensing and reanalysis sources ( Huot et al., 2022 ) .
Method
Params
λ^
Cov. ↑
FNR ↓
Set ↓
AUC ↑
LGBM ( p≥.5 )
—
.500
.072 ± .02
.928 ± .02
.001
.854
LGBM+CRC
—
.003
.941 ± .01
.059 ± .01
.626 ± .01
.854
LGBM+3way
—
[0,.03]
1.00
.000
.040
.854
UNet ( p≥.5 )
470K
.500
.390 ± .02
.610 ± .02
.007
.969
UNet+CRC
470K
.002
.947 ± .01
.053 ± .01
.149 ± .01
.969
UNet+3way
470K
[0,.02]
1.00
.000
.051
.969
Table 2: Test-set results across three models and three decision frameworks. CRC targets FNR≤0.05 . Three-way CRC uses cost-weighted αcw=0.50 , cfn=5 , cfp=1 , ρ∈[0.9,1.1] . All models trained from scratch on identical 70/15/15 splits. Values show ± 95% bootstrap CI (10K image resamples).
Figure 3: Same U-Net, different thresholds. Standard thresholding ( p^≥0.5 ) catches 56% of fires; CRC ( λ^=0.002 ) catches 100%. Bottom: probability cross-section showing how CRC captures low-probability fire pixels. Green shading marks pixels saved by CRC.
Figure 4: Three-model comparison. Left: AUROC. Center: without CRC, no model reaches 95% coverage. Right: with CRC, all three meet the safety target; both spatial models achieve ∼4× smaller evacuation zones than LightGBM.
Figure 5: The safety gap across three test samples (0.5%, 2.7%, 7.8% fire). Col. 1: ground truth. Col. 2: standard threshold, with missed fire in red (100%, 73%, 34% missed). Col. 3: CRC threshold, with improved detection (green). Col. 4: three-way zones (SAFE/MONITOR/EVACUATE).
Figure 6: Per-pixel fire risk gradient. Col. 1: ground truth. Col. 2: continuous probability heatmap. Col. 3: probability bands with CRC contour (green dashed) and EVACUATE boundary (red). Col. 4: calibrated risk gradient for operational mapping.
Model
SAFE
MONITOR
EVACUATE
Bpw
LightGBM
0.0%
95.96%
4.04%
0.988
U-Net
0.0%
94.90%
5.10%
0.987
ResGNN-UNet
0.0%
85.15%
14.85%
0.987
Table 3: Three-way zone fractions on the test set. The SAFE zone collapses because Bpw=max(cfnπ1,cfpπ0)=0.95 is dominated by cfpπ0 under ∼5% fire prevalence.
Model
Params
Epochs
Optim
Val Loss
AUC
LightGBM
—
100 iter
GBM
—
.854
Tiny U-Net
470K
50
AdamW
0.048
.969
ResGNN-UNet
229K
23 †
AdamW
0.131
.964
† Early-stopped from 30 (patience=7, best epoch 16).
Table 4: Model training summary. All models trained from scratch on identical NDWS splits.
Figure 7: Tiny U-Net training curves: BCE loss and cosine annealing LR over 50 epochs. Smooth convergence with minimal train–val gap.
Appendix figures & tables3 assets
Supplementary material from the paper’s appendix.
Appendix
Figure 8: Safety vs. efficiency for all nine configurations. Top: coverage (safety), where only CRC and three-way methods achieve ≥95% . Bottom: set size (efficiency), where both spatial models with CRC achieve the best trade-off.
Figure 9: Bare model vs. CRC across four test samples. Fire catch rate improves from 8% → 98%, 0% → 83%, 56% → 100%, 48% → 100%.
Figure 10: Qualitative prediction maps: ground truth, probabilities, standard prediction, CRC prediction, three-way zones, and error analysis.
Urban Resilience.AI Lab, Zachry Department of Civil and Environmental Engineering, Texas A&M University, College Station, TX 77840, United States of America