IPM → OpenMalaria handoff

Malaria seasonality — from a stale rainfall proxy to climate-driven emergence

Replacing the 2016-vintage country-level rainfall proxy with an admin1-resolved, climate-driven mosquito-emergence seasonality — built from TerraClimate + Vector Atlas data already on hand. 309 admin1 units across 18 countries. Three methods compared: legacy (rainfall+1mo), A water-balance, B species kernel. Kernel persistence τ (days): arabiensis=80, funestus=100, gambiae=95 · deployment lead: 1 month.

The idea: two decoupled objects

Emergence lags rainfall — deployment anticipates it

The emergence / transmission shape feeds OpenMalaria’s EIR and should lag rainfall (aquatic development + habitat filling). The intervention onset drives product timing and sits at or before rainfall onset — chemoprevention needs drug onboard before biting ramps, and roads wash out once rains arrive.

They share rainfall but diverge on purpose, so the emergence-model lag can never corrupt deployment timing.

Why it matters for the sim. We hand OpenMalaria the lagged emergence shape (per admin1); we keep deployment timing on raw-rainfall onset, nudged earlier by a tunable lead. Both come out of the same build, consistently, for every admin.

How each method works

What data goes in, and how it is used

Legacy — rainfall proxy

status quo · country-level

In: one monthly rainfall series per country (TerraClimate, the 2016-vintage table).

How: normalise to the annual max and shift +1 month as a fixed stand-in for the mosquito lag. No evaporation, no habitat persistence, no species, no sub-national variation — every admin in a country gets the same curve.

A — water-balance

zero-tuning physical baseline · admin1

In: admin1 monthly rainfall + PET (potential evapotranspiration).

How: a Thornthwaite–Mather soil-water bucket — rain recharges, PET drains. The lag and the dry-season tail fall out of the physics with no free parameters. Honest but coarse: soil moisture over-smooths and cannot tell species apart.

B — species kernel

mechanistic · admin1 · NDVI-calibrated

In: admin1 monthly rainfall + PET and Vector Atlas species composition (arabiensis / gambiae / funestus).

How: standing water accumulates with rain and decays at a rate set by PET and by each vector’s habitat persistence τ. Per-species habitat is blended by the local species mix. τ is calibrated to MODIS NDVI phase (~90 days). Captures the lag, the dry-season tail, and species structure.

Data sources — which dataset feeds which part

DataSourceResolutionFeeds
Rainfall & PETTerraClimatemonthly climatology, ~4 km → admin1legacy · A · B
Vector species mixVector Atlas (arabiensis / funestus / gambiae)admin1 fractionsB — species weighting
NDVI (greenness)MODIS MOD13Q116-day, 250 mB — τ calibration target
Population & prevalenceWorldPop 2020 · MAP PfPR2–10admin1decision-impact section

Rainfall and PET (potential evapotranspiration) are the only climate drivers — the same two fields feed all three methods; the models differ only in how they turn them into a mosquito curve. Species composition enters B alone; NDVI is used only to tune B, never as a direct input; population and prevalence are used only to weight the impact analysis.

How the habitat-persistence τ is set

τ is a per-species persistence time-constant (in days) — roughly how long a larval habitat lingers after rain before it dries out. It is a model parameter, not a directly measured quantity, and it enters the kernel through a monthly water-balance recursion:

habitatm = habitatm−1 × decaym + rainm
decaym = exp( −(30 / τ) × PETm / PETref )

Each month, standing water is topped up by rain and drained by evaporative demand. A hot, high-PET month drains it faster; a larger τ makes it persist — carrying transmission into the dry season. So τ is exactly what separates a flashy, rain-tracking vector from one that breeds in water that lasts.

Set in two steps — ordering from biology, magnitude from data.

1 · prior (ordering): arabiensis breeds in transient rain pools (short τ), funestus in permanent / semi-permanent water (long τ), gambiae in between.

2 · calibration (magnitude): fit τ so the kernel’s monthly shape best matches MODIS NDVI phase across 54 admins. Best-fit effective persistence is ~90 days; NDVI confirms the ordering but cannot separate the magnitudes, so we set arabiensis 80 / gambiae 95 / funestus 100 days.

Caveat: NDVI greenness persists longer than larval habitat, so these τ are an upper bound — entomological EIR would refine them downward.

Geospatial: where and when the season concentrates

Peak transmission month — country status quo → admin1 resolution → emergence model

The first panel is the actual status quo: one peak month per country. Panel 2 keeps the same rainfall proxy but at admin1 resolution — sub-national structure appears for free. Panels 3–4 add the water-balance and species-kernel models. Two distinct gains: resolution, then model.

Seasonality amplitude — perennial (dark) to sharply seasonal (bright)

Monthly emergence snapshots (kernel)

Why the decision matters: people & burden

Who is mistimed by the country-level status quo

Across 838M people in these admins, 60% live where the legacy country method mistimes the peak by ≥1 month versus the admin1 kernel (59% of estimated infections). For ≥2-month errors it is 15% of people but 20% of cases — prevalence concentrates exactly where the country average is most wrong.

Peak-transmission month under the legacy country method vs the admin1 kernel, per admin, weighted two ways: by population and by estimated infections (population × MAP PfPR2–10). The decision is not academic — it moves the recommended timing for a majority of people, and the burden-weighted share is higher still where disagreement is largest.

Admin1: profiles across all countries

Country-mean seasonality shape

Admin1: deployment timing

Rainfall onset month (drives deployment = onset − lead)

Calibration against MODIS NDVI

Independent validation — NDVI greenness as the phase target

We fit the kernel persistence τ to MODIS NDVI phase (min–max normalised + correlation, since NDVI’s timing is trustworthy but its amplitude is not). Result: effective persistence ~90 days; the calibrated kernel beats legacy and water-balance on phase correlation, shape error, and peak-month error (~2.0 → ~1.6 mo). NDVI confirms the species ordering but not the magnitudes. Caveat: NDVI over-persists vs larval habitat, so τ is an upper bound — entomological EIR would refine it.

Country explorer

Per-admin1 profiles — solid = rainfall onset · dashed = deployment · dotted = kernel peak