Modeling the association between Aedes aegypti ovitrap egg counts, multi-scale remotely sensed environmental data and arboviral cases at Puntarenas, Costa Rica (2017–2018)

Abstract

Problems with vector surveillance are a major barrier for the effective control of vector-borne disease transmission through Latin America. Here, we present results from a 80-week longitudinal study where Aedes aegypti (L.) (Diptera: Culicidae) ovitraps were monitored weekly at 92 locations in Puntarenas, a coastal city in Costa Rica with syndemic Zika, chikungunya and dengue transmission. We used separate models to investigate the association of either Ae. aegypti-borne arboviral cases or Ae. aegypti egg counts with remotely sensed environmental variables. We also evaluated whether Ae. aegypti-borne arboviral cases were associated with Ae. aegypti egg counts. Using cross-correlation and time series modeling, we found that arboviral cases were not significantly associated with Ae. aegypti egg counts. Through model selection we found that cases had a non-linear response to multi-scale (1-km and 30-m resolution) measurements of temperature standard deviation (SD) with a lag of up to 4 weeks, while simultaneously increasing with finely-grained NDVI (30-m resolution). Meanwhile, median ovitrap Ae. aegypti egg counts increased, and respectively decreased, with temperature SD (1-km resolution) and EVI (30-m resolution) with a lag of 6 weeks. A synchrony analysis showed that egg counts had a travelling wave pattern, with synchrony showing cyclic changes with distance, a pattern not observed in remotely sensed data with 30-m and 10-m resolution. Spatially, using generalized additive models, we found that eggs were more abundant at locations with higher temperatures and where EVI was leptokurtic during the study period. Our results suggest that, in Puntarenas, remotely sensed environmental variables are associated with both Ae. aegypti-borne arbovirus transmission and Ae. aegypti egg counts from ovitraps.

Keywords

• Oviposition
• MODIS
• Landsat 8
• Sentinel 2
• Schmalhausen’s law
• Model selection