Research
Which methods are better
29 methods have been compared with the approach the platform publishes, for 4 seasons and their months (20 targets), over All Ethiopia: 12 bias corrections, 8 probability calibrations, 6 rainfall-distribution models, 1 two-stage calibration and 2 statistical downscaling methods. Each was tested on its own against the published approach. So far no method is recommended in place of the published approach; the published forecasts keep it.
How a method is judged
Tercile probabilities
Scored with the RPS against the published MVS + climatology blend.
Rainfall amounts
Scored with the CRPS against the published MVS-corrected ensemble.
- Confirmed better (historical evidence): better over 1993–2016 when each year is forecast from the other years, significant after correcting for the number of methods and targets tested (Benjamini–Hochberg-adjusted p < 0.05, which controls the share of false findings; Holm until 17 October 2026), and not worse over the test years from 2017 on. A confirmed method becomes a candidate.
- Validation: a candidate is then scored on independent seasons as they are observed. It is recommended after 3 seasons that are not worse than the published approach, and not supported if worse.
- One change at a time: each method was tested against the published approach, not combined with another; the published forecasts change only after a recommendation and a decision.
Which method to use in each season
| Season | Tercile probabilities | Rainfall amounts |
|---|---|---|
| FMAM (January start) | Published (MVS + climatology blend) | Published (MVS-corrected ensemble)
Candidate: EMOS, censored shifted gamma for FMAM · validation 1 of 3 seasons |
| FMAM (February start) | Published (MVS + climatology blend)
Candidate: Square-root-space Mean-Variance Scaling for Feb · validation 1 of 3 seasons | Published (MVS-corrected ensemble) |
| JJAS (May start) | Published (MVS + climatology blend) | Published (MVS-corrected ensemble)
Candidate: Gamma Quantile Mapping for Jun · validation 1 of 3 seasons Candidate: Variance Scaling (unbounded) for Jul · validation 1 of 3 seasons Candidate: Principal component regression (PCR) for Aug, JJAS · validation 1 of 3 seasons Candidate: Regularized quantile regression for Sep · validation 1 of 3 seasons |
| ONDJ (September start) | Published (MVS + climatology blend) | Published (MVS-corrected ensemble)
Candidate: Square-root-space Mean-Variance Scaling for Jan · validation 0 of 3 seasons |
A candidate is confirmed on the historical evidence and waiting for validation seasons. The evidence for each target is on Method selection.
Findings
- Tercile probabilities. Confirmed better for at least one target: Square-root-space Mean-Variance Scaling (FMAM (Feb) Feb). No probability calibration, rainfall-distribution model, two-stage calibration or statistical downscaling method is confirmed better for any target. Most consistent: Square-root-space Mean-Variance Scaling, better on average for 14 of 20 targets over 1993–2016 and 16 of 20 over 2017+.
- Rainfall amounts. Confirmed better for at least one target: Square-root-space Mean-Variance Scaling (5), Gamma Quantile Mapping (4), Variance Scaling (unbounded) (4), Principal component regression (PCR) (3), EMOS, censored shifted gamma (2), Power Transformation (2), Regularized quantile regression (1), Leung–Dirichlet calibration (1), Canonical correlation analysis (CCA) (1). Most consistent: Square-root-space Mean-Variance Scaling, better on average for 20 of 20 targets over 1993–2016 and 18 of 20 over 2017+.
Comparison: tercile probabilities
Against the published MVS + climatology blend (RPS). Counts are targets; “Better in 2026” counts the targets of the independent 2026 seasons observed so far.
| Method | Category | Confirmed better | Better 1993–2016 | Better 2017+ | Better in 2026 |
|---|---|---|---|---|---|
| Square-root-space Mean-Variance Scaling | Bias correction | 1: FMAM (Feb) Feb | 14 of 20 | 16 of 20 | 11 of 15 |
| Variance Scaling (unbounded) | Bias correction | 0 | 16 of 20 | 12 of 20 | 11 of 15 |
| Gamma Quantile Mapping | Bias correction | 0 | 14 of 20 | 16 of 20 | 10 of 15 |
| Local Intensity Scaling | Bias correction | 0 | 14 of 20 | 6 of 20 | 2 of 15 |
| Linear Scaling | Bias correction | 0 | 13 of 20 | 5 of 20 | 3 of 15 |
| Power Transformation | Bias correction | 0 | 9 of 20 | 17 of 20 | 10 of 15 |
The other 23 methods (29 in all)
| Method | Category | Confirmed better | Better 1993–2016 | Better 2017+ | Better in 2026 |
|---|---|---|---|---|---|
| Temperature scaling | Probability calibration | 0 | 9 of 20 | 8 of 20 | — |
| Platt (sigmoid) calibration, one-vs-rest + normalization | Probability calibration | 0 | 9 of 20 | 6 of 20 | — |
| Beta calibration, one-vs-rest + normalization | Probability calibration | 0 | 8 of 20 | 9 of 20 | — |
| Regularized Dirichlet calibration | Probability calibration | 0 | 7 of 20 | 7 of 20 | — |
| Isotonic regression, one-vs-rest + normalization | Probability calibration | 0 | 5 of 20 | 8 of 20 | — |
| Bayesian model averaging for precipitation | Rainfall-distribution model | 0 | 5 of 20 | 3 of 20 | — |
| Quantile Delta Mapping | Bias correction | 0 | 4 of 20 | 7 of 20 | 5 of 15 |
| Regularized multinomial logistic regression | Probability calibration | 0 | 4 of 20 | 6 of 20 | — |
| Delta Change | Bias correction | 0 | 4 of 20 | 5 of 20 | 7 of 15 |
| CDF-transform | Bias correction | 0 | 4 of 20 | 5 of 20 | 5 of 15 |
| Ordinal (cumulative-logit) logistic regression | Probability calibration | 0 | 4 of 20 | 5 of 20 | — |
| Principal component regression (PCR) | Statistical downscaling | 0 | 3 of 20 | 8 of 20 | 11 of 15 |
| Detrended Quantile Mapping | Bias correction | 0 | 3 of 20 | 7 of 20 | 6 of 15 |
| EMOS, censored shifted gamma | Rainfall-distribution model | 0 | 3 of 20 | 6 of 20 | 9 of 15 |
| Empirical Quantile Mapping | Bias correction | 0 | 3 of 20 | 5 of 20 | 7 of 15 |
| Regularized quantile regression | Rainfall-distribution model | 0 | 3 of 20 | 5 of 20 | 3 of 5 |
| Leung–Dirichlet calibration | Two-stage calibration | 0 | 2 of 15 | 3 of 15 | 1 of 10 |
| Heteroscedastic extended logistic regression | Rainfall-distribution model | 0 | 2 of 20 | 8 of 20 | — |
| Scaled Distribution Mapping | Bias correction | 0 | 2 of 20 | 3 of 20 | 5 of 15 |
| Pseudocount smoothing (no blend) | Probability calibration | 0 | 2 of 20 | 2 of 20 | — |
| Canonical correlation analysis (CCA) | Statistical downscaling | 0 | 1 of 20 | 13 of 20 | 9 of 15 |
| Quantile regression forest | Rainfall-distribution model | 0 | 1 of 20 | 2 of 20 | — |
| Bayesian joint probability | Rainfall-distribution model | 0 | 0 of 20 | 1 of 20 | — |
Comparison: rainfall amounts
Against the published MVS-corrected ensemble (CRPS). Counts are targets; “Better in 2026” counts the targets of the independent 2026 seasons observed so far.
| Method | Category | Confirmed better | Better 1993–2016 | Better 2017+ | Better in 2026 |
|---|---|---|---|---|---|
| Square-root-space Mean-Variance Scaling | Bias correction | 5: JJAS (May) Jun, JJAS (May) Jul, JJAS (May) Aug, JJAS (May) Sep, ONDJ (Sep) Jan | 20 of 20 | 18 of 20 | 10 of 15 |
| Gamma Quantile Mapping | Bias correction | 4: JJAS (May) Jun, JJAS (May) Aug, JJAS (May) Sep, ONDJ (Sep) Jan | 20 of 20 | 16 of 20 | 9 of 15 |
| Variance Scaling (unbounded) | Bias correction | 4: JJAS (May) Jul, JJAS (May) Aug, JJAS (May) Sep, ONDJ (Sep) Jan | 14 of 20 | 15 of 20 | 6 of 15 |
| Principal component regression (PCR) | Statistical downscaling | 3: JJAS (May) Aug, JJAS (May) Sep, JJAS (May) JJAS | 7 of 20 | 12 of 20 | 6 of 15 |
| EMOS, censored shifted gamma | Rainfall-distribution model | 2: FMAM (Jan) FMAM, JJAS (May) Sep | 15 of 20 | 12 of 20 | 5 of 15 |
| Power Transformation | Bias correction | 2: JJAS (May) Jul, JJAS (May) Aug | 12 of 20 | 16 of 20 | 9 of 15 |
| Regularized quantile regression | Rainfall-distribution model | 1: JJAS (May) Sep | 17 of 20 | 14 of 20 | 0 of 5 |
| Leung–Dirichlet calibration | Two-stage calibration | 1: ONDJ (Sep) Jan | 11 of 15 | 9 of 15 | 5 of 10 |
| Canonical correlation analysis (CCA) | Statistical downscaling | 1: JJAS (May) Sep | 7 of 20 | 12 of 20 | 6 of 15 |
The other 12 methods (21 in all)
| Method | Category | Confirmed better | Better 1993–2016 | Better 2017+ | Better in 2026 |
|---|---|---|---|---|---|
| Bayesian model averaging for precipitation | Rainfall-distribution model | 0 | 13 of 20 | 7 of 20 | — |
| Heteroscedastic extended logistic regression | Rainfall-distribution model | 0 | 9 of 20 | 9 of 20 | — |
| Delta Change | Bias correction | 0 | 9 of 20 | 8 of 20 | 6 of 15 |
| Empirical Quantile Mapping | Bias correction | 0 | 6 of 20 | 7 of 20 | 7 of 15 |
| Quantile regression forest | Rainfall-distribution model | 0 | 5 of 20 | 3 of 20 | — |
| Quantile Delta Mapping | Bias correction | 0 | 4 of 20 | 5 of 20 | 3 of 15 |
| Bayesian joint probability | Rainfall-distribution model | 0 | 3 of 20 | 4 of 20 | — |
| Detrended Quantile Mapping | Bias correction | 0 | 2 of 20 | 4 of 20 | 5 of 15 |
| CDF-transform | Bias correction | 0 | 2 of 20 | 3 of 20 | 4 of 15 |
| Linear Scaling | Bias correction | 0 | 1 of 20 | 2 of 20 | 2 of 15 |
| Local Intensity Scaling | Bias correction | 0 | 0 of 20 | 2 of 20 | 2 of 15 |
| Scaled Distribution Mapping | Bias correction | 0 | 0 of 20 | 0 of 20 | 1 of 15 |
More statistical power
The same out-of-sample scores, no refitting, read three ways that give the tests more power:
- 33 years: each target's 1993–2016 leave-one-year-out years joined with its test years (2017 onward, fitted on 1993–2016), so every year 1993–2025 counts, scored by a model that did not see it.
- Pooled: for each method and year, its skill difference from the published approach averaged over all its targets of all seasons (as a share of the climatology score, so targets of different rainfall add up); one test over the years.
- Two corrections for testing many methods: Holm (controls the chance of any false finding) and Benjamini–Hochberg (controls the share of false findings; less conservative).
The selection rule above applies Benjamini–Hochberg within each benchmark family (it used Holm until 17 October 2026; changed by decision after these analyses were seen). The pooled and 33-year tests below are sensitivity analyses, not part of the rule.
Tercile probabilities
Significantly better than the published approach when pooled: Holm none; Benjamini–Hochberg none.
| Method | Category | Pooled skill gain | Years better | p | Holm p | BH p | Targets confirmed over 33 years (Holm / BH) |
|---|---|---|---|---|---|---|---|
| Square-root-space Mean-Variance Scaling | Bias correction | +0.63% | 25 of 33 | 0.004 | 0.124 | 0.082 | 0 / 0 of 20 |
| Power Transformation | Bias correction | +0.38% | 19 of 33 | 0.047 | 1.000 | 0.339 | 0 / 0 of 20 |
| Gamma Quantile Mapping | Bias correction | +0.33% | 21 of 33 | 0.006 | 0.159 | 0.082 | 0 / 0 of 20 |
| Regularized Dirichlet calibration | Probability calibration | +0.18% | 19 of 33 | 0.199 | 1.000 | 1.000 | 0 / 0 of 20 |
| Beta calibration, one-vs-rest + normalization | Probability calibration | +0.08% | 16 of 33 | 0.265 | 1.000 | 1.000 | 0 / 0 of 20 |
| Variance Scaling (unbounded) | Bias correction | +0.07% | 19 of 33 | 0.033 | 0.897 | 0.321 | 0 / 0 of 20 |
| Platt (sigmoid) calibration, one-vs-rest + normalization | Probability calibration | −0.06% | 13 of 33 | 0.679 | 1.000 | 1.000 | 0 / 0 of 20 |
| CDF-transform | Bias correction | −0.20% | 17 of 33 | 0.744 | 1.000 | 1.000 | 0 / 0 of 20 |
The other 21 methods
| Method | Category | Pooled skill gain | Years better | p | Holm p | BH p | Targets confirmed over 33 years (Holm / BH) |
|---|---|---|---|---|---|---|---|
| Linear Scaling | Bias correction | −0.25% | 14 of 33 | 0.723 | 1.000 | 1.000 | 0 / 0 of 20 |
| Local Intensity Scaling | Bias correction | −0.26% | 14 of 33 | 0.739 | 1.000 | 1.000 | 0 / 0 of 20 |
| Empirical Quantile Mapping | Bias correction | −0.34% | 15 of 33 | 0.938 | 1.000 | 1.000 | 0 / 0 of 20 |
| Detrended Quantile Mapping | Bias correction | −0.35% | 16 of 33 | 0.885 | 1.000 | 1.000 | 0 / 0 of 20 |
| Ordinal (cumulative-logit) logistic regression | Probability calibration | −0.37% | 9 of 33 | 0.949 | 1.000 | 1.000 | 0 / 0 of 20 |
| Regularized multinomial logistic regression | Probability calibration | −0.38% | 9 of 33 | 0.949 | 1.000 | 1.000 | 0 / 0 of 20 |
| Quantile Delta Mapping | Bias correction | −0.39% | 15 of 33 | 0.918 | 1.000 | 1.000 | 0 / 0 of 20 |
| Delta Change | Bias correction | −0.40% | 14 of 33 | 0.904 | 1.000 | 1.000 | 0 / 0 of 20 |
| Temperature scaling | Probability calibration | −0.46% | 11 of 33 | 0.990 | 1.000 | 1.000 | 0 / 0 of 20 |
| Scaled Distribution Mapping | Bias correction | −0.72% | 13 of 33 | 0.979 | 1.000 | 1.000 | 0 / 0 of 20 |
| Leung–Dirichlet calibration | Two-stage calibration | −0.78% | 12 of 33 | 0.992 | 1.000 | 1.000 | 0 / 0 of 15 |
| Isotonic regression, one-vs-rest + normalization | Probability calibration | −1.65% | 4 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| Heteroscedastic extended logistic regression | Rainfall-distribution model | −2.03% | 8 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| EMOS, censored shifted gamma | Rainfall-distribution model | −2.51% | 8 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| Principal component regression (PCR) | Statistical downscaling | −2.91% | 8 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| Canonical correlation analysis (CCA) | Statistical downscaling | −3.17% | 12 of 33 | 0.998 | 1.000 | 1.000 | 0 / 0 of 20 |
| Regularized quantile regression | Rainfall-distribution model | −3.19% | 2 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| Bayesian model averaging for precipitation | Rainfall-distribution model | −3.19% | 5 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| Bayesian joint probability | Rainfall-distribution model | −3.98% | 2 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| Pseudocount smoothing (no blend) | Probability calibration | −5.31% | 3 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| Quantile regression forest | Rainfall-distribution model | −6.49% | 1 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
Rainfall amounts
Significantly better than the published approach when pooled: Holm Regularized quantile regression, Square-root-space Mean-Variance Scaling, Gamma Quantile Mapping, Variance Scaling (unbounded); Benjamini–Hochberg Regularized quantile regression, Square-root-space Mean-Variance Scaling, EMOS, censored shifted gamma, Gamma Quantile Mapping, Variance Scaling (unbounded).
| Method | Category | Pooled skill gain | Years better | p | Holm p | BH p | Targets confirmed over 33 years (Holm / BH) |
|---|---|---|---|---|---|---|---|
| Regularized quantile regression | Rainfall-distribution model | +2.39% | 28 of 33 | < 0.001 | < 0.001 | < 0.001 | 0 / 1 of 20 |
| Square-root-space Mean-Variance Scaling | Bias correction | +1.58% | 29 of 33 | < 0.001 | < 0.001 | < 0.001 | 3 / 12 of 20 |
| EMOS, censored shifted gamma | Rainfall-distribution model | +1.47% | 25 of 33 | 0.003 | 0.053 | 0.013 | 1 / 2 of 20 |
| Gamma Quantile Mapping | Bias correction | +0.97% | 25 of 33 | < 0.001 | < 0.001 | < 0.001 | 2 / 8 of 20 |
| Bayesian model averaging for precipitation | Rainfall-distribution model | +0.68% | 17 of 33 | 0.105 | 1.000 | 0.316 | 0 / 0 of 20 |
| Variance Scaling (unbounded) | Bias correction | +0.59% | 28 of 33 | < 0.001 | < 0.001 | < 0.001 | 3 / 7 of 20 |
| Empirical Quantile Mapping | Bias correction | +0.54% | 20 of 33 | 0.040 | 0.635 | 0.139 | 0 / 1 of 20 |
| Leung–Dirichlet calibration | Two-stage calibration | −0.09% | 18 of 33 | 0.830 | 1.000 | 1.000 | 1 / 2 of 15 |
The other 13 methods
| Method | Category | Pooled skill gain | Years better | p | Holm p | BH p | Targets confirmed over 33 years (Holm / BH) |
|---|---|---|---|---|---|---|---|
| Power Transformation | Bias correction | −0.36% | 18 of 33 | 0.850 | 1.000 | 1.000 | 1 / 3 of 20 |
| Heteroscedastic extended logistic regression | Rainfall-distribution model | −1.11% | 10 of 33 | 0.978 | 1.000 | 1.000 | 0 / 0 of 20 |
| Delta Change | Bias correction | −1.71% | 10 of 33 | 0.997 | 1.000 | 1.000 | 0 / 0 of 20 |
| Bayesian joint probability | Rainfall-distribution model | −1.90% | 5 of 33 | 0.996 | 1.000 | 1.000 | 0 / 0 of 20 |
| Principal component regression (PCR) | Statistical downscaling | −1.93% | 16 of 33 | 0.626 | 1.000 | 1.000 | 0 / 2 of 20 |
| Canonical correlation analysis (CCA) | Statistical downscaling | −2.37% | 15 of 33 | 0.783 | 1.000 | 1.000 | 0 / 1 of 20 |
| Quantile regression forest | Rainfall-distribution model | −3.29% | 6 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| Quantile Delta Mapping | Bias correction | −3.40% | 9 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| Detrended Quantile Mapping | Bias correction | −4.47% | 5 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| Local Intensity Scaling | Bias correction | −4.85% | 5 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| Linear Scaling | Bias correction | −6.30% | 2 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| Scaled Distribution Mapping | Bias correction | −6.40% | 3 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
| CDF-transform | Bias correction | −7.39% | 4 of 33 | 1.000 | 1.000 | 1.000 | 0 / 0 of 20 |
Pooled skill gain: mean over years of the method's improvement over the published approach, as a share of the climatology score (+1% = 1% of the climatology score better). Years better: years where the pooled difference favours the method.
Research pages
- Method selection
The approach the evidence supports for each season, target and product, with the evidence for each change. - Bias correction
Thirteen rainfall-amount corrections against the operational Mean-Variance Scaling. - Probability calibration
Probability-stage methods and rainfall-distribution models (Phases 1–3) against the published probabilities. - Amount and range shadow
EMOS medians and 10–90% ranges next to the MVS ensemble, scored season by season. - Statistical downscaling
Model output statistics (CCA, PCR), the CPT approach of ICPAC and IRI, against the published pipeline. - Leung–Dirichlet calibration
Member-wise Leung correction and identity-anchored Dirichlet, each stage against the published pipeline.