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 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.
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.
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.
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 | Source | Resolution | Feeds |
|---|---|---|---|
| Rainfall & PET | TerraClimate | monthly climatology, ~4 km → admin1 | legacy · A · B |
| Vector species mix | Vector Atlas (arabiensis / funestus / gambiae) | admin1 fractions | B — species weighting |
| NDVI (greenness) | MODIS MOD13Q1 | 16-day, 250 m | B — τ calibration target |
| Population & prevalence | WorldPop 2020 · MAP PfPR2–10 | admin1 | decision-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.
τ 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.
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.
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.