Phase 2 · Probability calibration · February initialization · FMAM and its months

Phase 2 · February initialization · FMAM and its months

Binary calibrators (one-vs-rest, normalized)

Platt, beta and isotonic calibration of each category against the rest, divided by their sum inside the pipeline, against the incumbent (their own family: Holm across these methods and the targets). Confirmed better: none.

Tercile RPSS, outer leave-one-year-out 1993–2016
MethodFebMarAprMayFMAM
Climatology blend (incumbent)+0.126+0.011+0.009+0.008+0.021
Platt (sigmoid) calibration, one-vs-rest + normalization+0.135+0.008+0.011+0.007+0.020
Beta calibration, one-vs-rest + normalization+0.143+0.015+0.008+0.007+0.017
Isotonic regression, one-vs-rest + normalization+0.061+0.017+0.002+0.010+0.018
SEL (chosen each year)+0.137+0.014+0.009+0.008+0.021
Tercile RPSS, fit 1993–2016, scored 2017–2025 (exploratory)
MethodFebMarAprMayFMAM
Climatology blend (incumbent)+0.056+0.029+0.012+0.026+0.059
Platt (sigmoid) calibration, one-vs-rest + normalization+0.037+0.021+0.011+0.014+0.054
Beta calibration, one-vs-rest + normalization+0.060+0.033+0.015+0.020+0.049
Isotonic regression, one-vs-rest + normalization+0.016+0.045+0.014+0.029+0.067
SEL (chosen each year)+0.060+0.029+0.012+0.026+0.059

Rainfall-distribution models (Pathway B)

Each model predicts a full rainfall distribution from the raw ensemble (pooled over cells in units relative to each cell's climatology); terciles at the training CHIRPS boundaries, quantiles and intervals come from it. References, in the same folds: the published probabilities (tercile RPS) and the MVS-corrected ensemble (CRPS, quantiles). ✓ = better than the reference with Holm-adjusted p < 0.05 across the models and targets, and not worse over 2017–2025. Confirmed better — tercile probabilities: none; rainfall amounts (CRPS): none.

Tercile probabilities

Tercile RPSS, outer leave-one-year-out 1993–2016
MethodFebMarAprMayFMAM
Published probabilities (incumbent)+0.137+0.011+0.009+0.008+0.021
EMOS, censored shifted gamma+0.103-0.043-0.003-0.008+0.028
Heteroscedastic extended logistic regression+0.168-0.038+0.004-0.007+0.021
Regularized quantile regression+0.099-0.021-0.006-0.008+0.023

Rainfall amounts

Amount CRPSS against the CHIRPS climatology, 1993–2016
MethodFebMarAprMayFMAM
MVS-corrected ensemble+0.250-0.034-0.001-0.005-0.007
EMOS, censored shifted gamma+0.248-0.013+0.003-0.003+0.026
Heteroscedastic extended logistic regression+0.234-0.067-0.008-0.023+0.017
Regularized quantile regression+0.259-0.010+0.000-0.004+0.020

Intervals and quantiles

Share of observations inside the 10–90% interval, 1993–2016 (80% is calibrated)
MethodFebMarAprMayFMAM
MVS-corrected ensemble55%59%68%67%70%
EMOS, censored shifted gamma83%81%76%78%79%
Heteroscedastic extended logistic regression90%91%84%84%84%
Regularized quantile regression78%77%77%77%77%
Pinball loss of the median (mm), 1993–2016 (lower is better)
MethodFebMarAprMayFMAM
MVS-corrected ensemble3.9413.8520.0120.4234.86
EMOS, censored shifted gamma4.0414.1120.1720.6234.08
Heteroscedastic extended logistic regression3.9814.2820.1520.6634.29
Regularized quantile regression4.0514.0220.1820.5234.14

Test years

Fit 1993–2016, scored 2017–2025 (exploratory)
MethodFebMarAprMayFMAM
Published probabilities / MVS ensemble+0.059 / +0.000+0.029 / +0.065+0.012 / +0.004+0.026 / +0.025+0.059 / +0.071
EMOS, censored shifted gamma (RPSS / CRPSS)+0.028 / +0.021+0.038 / +0.053+0.004 / +0.007+0.036 / +0.022+0.070 / +0.059
Heteroscedastic extended logistic regression (RPSS / CRPSS)+0.026 / +0.033+0.038 / +0.033+0.012 / -0.004+0.044 / +0.007+0.068 / +0.061
Regularized quantile regression (RPSS / CRPSS)+0.000 / +0.023+0.038 / +0.037+0.014 / +0.008+0.046 / +0.025+0.068 / +0.057

Fits and model definitions

  • EMOS, censored shifted gamma: converged in every fold.
  • Heteroscedastic extended logistic regression: converged in every fold.
  • Regularized quantile regression: converged in every fold; quantile rearrangement needed for at most 2.2% of cells.
Models
ModelEstimatesPredictorsFittingOutputs
EMOS, censored shifted gammarainfall distribution with a probability of zeroensemble mean m' and mean absolute difference md' (relative)maximum likelihood, pooled over cellsCDF, quantiles, tercile probabilities
Heteroscedastic extended logistic regressioncoherent CDF across all thresholdssqrt of the ensemble mean m', log spread s' (relative)log loss over each cell’s CHIRPS deciles, pooledCDF, quantiles, tercile probabilities
Regularized quantile regressionconditional quantiles at 19 levels (5–95%), CDF by interpolationensemble mean m' and spread s' (relative)pinball loss per level (Huber-smoothed), L2 1e-4, pooledCDF, quantiles, tercile probabilities
Tercile RPSS and amount CRPSS of the Phase 2 distribution models against the published pipeline.
Tercile RPSS and amount CRPSS of the Phase 2 distribution models against the published pipeline.

← Phase 1Phase 3 →