Abstract
Malaria risk is highly heterogeneous. Understanding village and household-level spatial heterogeneity of malaria risk can support a transition to spatially targeted interventions for malaria elimination. This analysis uses data from cross-sectional prevalence surveys conducted in 2014 and 2016 in two villages (Megiar and Mirap) in Papua New Guinea. Generalised additive modelling was used to characterise spatial heterogeneity of malaria risk and investigate the contribution of individual, household and environmental-level risk factors. Following a period of declining malaria prevalence, the prevalence of P. falciparum increased from 11.4 to 19.1% in Megiar and 12.3 to 28.3% in Mirap between 2014 and 2016, with focal hotspots observed in these villages in 2014 and expanding in 2016. Prevalence of P. vivax was similar in both years (20.6% and 18.3% in Megiar, 22.1% and 23.4% in Mirap) and spatial risk heterogeneity was less apparent compared to P. falciparum. Within-village hotspots varied by Plasmodium species across time and between villages. In Megiar, the adjusted odds ratio (AOR) of infection could be partially explained by household factors that increase risk of vector exposure, such as collecting outdoor surface water as a main source of water. In Mirap, increased AOR overlapped with proximity to densely vegetated areas of the village. The identification of household and environmental factors associated with increased spatial risk may serve as useful indicators of transmission hotspots and inform the development of tailored approaches for malaria control.