key: cord-0951238-d6azk9z1
authors: Andreo, Verónica; Porcasi, Ximena; Guzman, Claudio; Lopez, Laura; Scavuzzo, Carlos M.
title: Spatial Distribution of Aedes aegypti Oviposition Temporal Patterns and Their Relationship with Environment and Dengue Incidence
date: 2021-10-09
journal: Insects
DOI: 10.3390/insects12100919
sha: c4f8192dc45e4436cfdae731fb7c5f02f73f3021
doc_id: 951238
cord_uid: d6azk9z1

SIMPLE SUMMARY: Aedes aegypti, the mosquito species that transmits dengue virus among others, is fully adapted to thrive in urban areas. Their activity, however, varies in time and space and this might imply different transmission risk. We hypothesize that the temporal differences in mosquito activity are determined by local environmental conditions. Hence, we explore the existence of groups of temporal patterns in weekly time series of ovitraps records and we associate those patterns to environmental variables derived from remote sensing data and also to dengue incidence. We found three groups of temporal patterns that showed association with land cover diversity, heterogeneity and variability in vegetation and humidity indices estimated over 50-m radius buffer areas surrounding ovitraps. Dengue incidence on a neighborhood basis showed a weak but positive association with the percentage of pixels belonging to one of the patterns detected. The understanding of the spatial distribution of temporal patterns and their environmental determinants might then become relevant to guide the allocation of prevention and monitoring interventions. ABSTRACT: Aedes aegypti, the mosquito species transmitting dengue, zika, chikungunya and yellow fever viruses, is fully adapted to thrive in urban areas. The temporal activity of this mosquito, however, varies within urban areas which might imply different transmission risk. In this work, we hypothesize that temporal differences in mosquito activity patterns are determined by local environmental conditions. Hence, we explore the existence of groups of temporal patterns in weekly time series of Ae. aegypti ovitraps records (2017–2019) by means of time series clustering. Next, with the aim of predicting risk in places with no mosquito field data, we use machine learning classification tools to assess the association of temporal patterns with environmental variables derived from satellite imagery and predict temporal patterns over the city area to finally test the relationship with dengue incidence. We found three groups of temporal patterns that showed association with land cover diversity, variability in vegetation and humidity and, heterogeneity measured by texture indices estimated over buffer areas surrounding ovitraps. Dengue incidence on a neighborhood basis showed a weak but positive association with the percentage of pixels belonging to only one of the temporal patterns detected. The understanding of the spatial distribution of temporal patterns and their environmental determinants might then become highly relevant to guide the allocation of prevention and potential interventions. Further investigation is still needed though to incorporate other determinants not considered here.

Aedes aegypti is the main vector of dengue, zika, chikungunya and yellow fever viruses worldwide. In recent decades, aided by global warming, urbanization, trade and human migration, Ae. aegypti has invaded many temperate areas of the world [1, 2] , reaching latitudes as south as 40°S [3] . This mosquito species is highly anthropophilic and it is fully adapted to thrive in urban areas, where it can fulfill all its ecological needs. Its success is mainly due to its behavior and survival ability, which includes egg laying in natural and man-made water containers associated with domestic and peridomestic settings and eggs' resistance to desiccation [4] . These characteristics, plus a wide distribution and high densities of breeding sites, constitute the key factors that mainly determine the circulation and transmission of dengue and related viruses in urban areas [5] .

Dengue fever causes the greatest human disease burden, with an estimated 10,000 deaths and 100,000 million symptomatic infections per year in over 125 countries (53% of the global population at risk) [6] . The incidence of dengue has grown dramatically in recent decades, with a concomitant increasing frequency of outbreaks, especially in South America during the past 10-12 years. Two of the most important dengue outbreaks in Argentina occurred in 2009 and 2016, affecting more than 25,000 and 40,000 people, respectively, and reaching temperate cities such as Córdoba and Buenos Aires [7] [8] [9] [10] . During 2020, in the midst of COVID-19 pandemic, Argentina experienced the largest dengue outbreak to date, with more than 58,000 confirmed cases, and near 3000 in the city of Córdoba [11] .

Given the lack of a well established and accepted vaccine, dengue control and prevention is traditionally based on vector control and entomological surveillance to estimate the potential risk for virus transmission and disease [12] . Prevention programs are therefore typically focused on removal of Ae. aegypti breeding sites in order to eliminate vector larval stages, treatment of larval habitats and insecticidal spraying to reduce adult density [13] . The latter being usually applied as a reactive measure once cases appear.

In this context, understanding the determinants of mosquito spatial distribution and temporal variations in abundance is key [14, 15] . Indeed, many studies have focused on mapping vector hotspots and disease clusters, as well as understanding the causes of spatial heterogeneity at different scales as key tools for decision making regarding prevention programs and control actions [16] [17] [18] [19] [20] . Others have addressed the prediction of temporal changes in mosquito abundances, either through mechanistic [21] or empirical models [22] [23] [24] with the aim of forecasting risk. In most cases, Earth Observation (EO) data of different spatial and temporal resolution were used as the main source of environmental information to relate to mosquito data and yield predictions [14] . In a different path but with the same aim, various clustering techniques have been used, i.e., spatiotemporal clustering was applied to identify clusters of human cases during outbreaks [8, 9] , coclustering techniques were used to identify favorable space-time conditions triggering outbreaks [25] , spatial clustering was used in the identification of different cover types or groups of neighborhoods [26] .

Time series clustering is a special type of clustering that handles dynamic data and has recently received more attention [27] . Examples in the literature commonly include the use of Dynamic Time Warping (DTW) distance combined with a clustering method to assess patterns in economic time series [28] , crop classification [29] , hydrodynamic behavior of ground water level [30] and phenological regions delineation [31] . Remarkably, only once has clustering been applied to time series of vector data in an attempt to understand if there might be different temporal patterns that could explain the spatial differences usually observed in timing of cases occurrence and transmission risk [32] .

In this contribution, we explore a novel approach that combines time series clustering, EO data, machine learning and dengue incidence, to assess the existence and spatial distribution of Ae. aegypti oviposition temporal patterns and their relationship with environment and dengue incidence in the city of Córdoba. Under the hypothesis that local environmental features might determine differences in female mosquitoes oviposition patterns and that different temporal patterns might be related to different risk levels of dengue transmission, we use time series clustering to group time series of eggs' counts. Furthermore, since our ultimate goal is to enhance operational tools and predict risk in places where we do not have data, we then assess the relationship between temporal patterns and environmental variables extracted from the analysis of high resolution (10 m) EO data in the surroundings of each ovitrap. This allows us to predict types of temporal patterns over space, which in turn facilitates the characterization of neighborhoods according to oviposition activity pattern and evaluate the relationship with dengue occurrence in the last outbreak.

Córdoba is the second largest city in Argentina with a population of 1,330,023 inhabitants in 2010 [33] . It has a surface of 576 km 2 and it is located at 31°24 S, 64°11 W, 450 m.a.s.l. The urban area is surrounded by agricultural fields and small forest patches.

Córdoba city has a temperate climate, with mean annual precipitation of 800 mm. The winter is markedly dry and most precipitation occurs in the summer months. The rainy season spans between October and March, with the highest precipitations from December to February. The mean annual temperature is 21°C (range 12-38°C). Winters are temperate, with several frost days in June and July.

The Suquía River, its tributary La Cañada and numerous additional water channels run through the city. Human activities have resulted in a landscape characterized by a highly developed urban core represented by buildings and green areas in the form of urban parks. Suburban areas are characterized by residential neighborhoods, primarily single-family houses with yards, interspersed with parks and other green spaces.

Entomological data consisted of 300 ovitraps distributed in 150 houses over 5 different areas of the city (Figure 1 ). Houses with ovitraps were at least 150 m apart, with an average distance of 350 m. Two ovitraps were placed in the front yard of each house, usually in shaded places and below or close to bushes or pots with plants. The possibility of trap installation depended on householders' written consent. Each ovitrap consisted of a black 1000 mL plastic container filled with 250 mL of water and a wooden tongue depressor (15 × 2 cm) held vertically to serve as substrate for mosquito oviposition [34] . Eggs laid in each wooden paddle were counted under magnifying glass. All ovitraps were replaced every week, from October 2017 to December 2019. Since the two ovitraps per house are not discriminated in the field, i.e., they are not labeled distinctively, we use the average count of eggs per house per week as the input for further analysis. We performed temporal linear interpolations when data were missing for different reasons. We removed houses (i.e., time series) with more than 10 consecutive missing values to avoid artifacts in temporal interpolations arising from not enough valid data points.

We used Sentinel 2 data for 3 different periods covering the core mosquito seasons in the area during our study period: November 2017 to March 2018, November 2018 to March 2019, and November 2019 to March 2020. A detail of the scenes processed can be found in Table S1 of the Supplementary Materials. Sentinel 2 imagery consists of 12 bands (3 visible and near infrared at 10 m, 3 red-edge at 20 m, 1 near infrared at 20 m, 2 short wave infrared at 60 m, among others). It has a revisit time of around 5 days, given that there are 2 satellites, Sentinel 2A and Sentinel 2B. We selected level 1 scenes with no clouds, one per month, in similar dates, for the 3 periods described above. All bands were atmospherically corrected by means of the Atmospheric and Radiometric Correction of Satellite Imagery (ARCSI) software and imported into GRASS GIS 7.8 [35] where all further remote sensing processing was done.

Common vegetation and water indices such as the Normalized Difference Vegetation Index (NDVI) and the Normalized Difference Water Index (NDWI) were derived for each image. Then, they were temporally averaged to obtain a single summary index image per season. We also estimated a synthetic panchromatic band for each season by temporally aggregating red, green and blue bands using the median (RGB composites for each season are shown in Figure S1 in the Supplementary Materials), and then, averaging them.

Following [20] , a k-means unsupervised classification with 15 different classes was chosen as to identify different spectral covers (see Figure S2 in the Supplementary Materials). Buffers of ≈50and 100-m radius from each ovitrap were overlaid upon bands and different statistics and texture measures were obtained, i.e., number of classes, the most common class, class diversity, mean and standard deviation (sd) NDVI, mean and sd NDWI, contrast, entropy, interspersion, etc. Buffer sizes were established according to the commonly used 100 m flight range of host-seeking female mosquitoes. All variables derived from remote sensing data are defined and explained in Table S2 in Supplementary Materials.

Locations of imported and autochthonous dengue cases for the 2019-2020 season were provided by the Health authorities of Córdoba province. In total, there were 2755 dengue cases in Córdoba city in the 2019-2020 season, 28 were imported and 2727 autochthonous out of which we recovered 2491 coordinates (see Section 3).

We tested different partitional time series clustering algorithms and different distance measures with number of clusters (k) ranging from 3 to 10 (Table 1 ) [27] . We ran 10 repetitions of maximum 100 iterations for each combination of k, distance, and centroid extraction method. For the case of Dynamic Time Warping (DTW) clustering, there were other parameters considered such as vector norm (i.e., Manhattan and Euclidean distances, L1 and L2, respectively) and window size (from 1 to 5, in steps of 1). Since the shape-based distance (SBD) requires data normalization, we applied this pre-processing step in all cases for the sake of consistency, and we tested configurations both with and without centering. Since clustering is an unsupervised method, we used internal cluster validity indices (CVI) to evaluate the results. These indices only consider the partitioned data and try to quantify cluster purity. Each index defines its range of values and whether they are to be minimized or maximized to assign a vote to a certain configuration within a set [36] . Several CVIs were estimated for each combination of settings and majority vote was used to decide on a final result. The CVIs used were: Silhouette (Sil), Dunn (D), COP index (COP), Davies-Bouldin (DB) and Modified Davies-Bouldin (DB*) [36] . We compared all 6 combinations from Table 1 . We then performed majority vote among the 6 configurations per subset to select one best configuration for further comparisons. For time series clustering and evaluation of results, we used the dtwclust package [37] in the R software [38] . The results of the best configurations were then compared using different similarity measures from the clue package [39] .

To understand if local environmental variables derived from remote sensing in 50and 100-m radius areas could explain the clustering of ovitraps time series obtained, we ran random forest (RF) classifications for each period. The dependent variable or outcome was the cluster number. These analyses were run with the package caret [40] in the R software [38] . In all cases, we split the data into training and test sets (70 and 30%, respectively). Given that clusters were unbalanced, we used an up-sample approach before RF training to increase n in smaller groups. Furthermore, we centered and scaled variables in the training sample to avoid that differences in values might affect variable importance. We used repeated cross-validation for hyper-parameter tuning in the training phase, with 5 folds and 10 repetitions. We performed variable selection by means of anova scores filter. This method uses anova p-values as weights to decide which variables or features are to be used in the classification step and we estimated variable importance to infer relationships with clusters. Finally, we evaluated the classification's results with the 30% test dataset that was left aside at the beginning. Overall accuracy was used as performance measure both for training and validation.

Since the ovitrap sampling could not continue given the control activities that had to be performed because of the dengue outbreak starting in January 2020 all over the province, and the COVID-19 pandemic afterwards, we do not have complete ovitrap records for 2019-2020 as to carry out the analysis. Therefore, we used 2018-2019 RF model to predict the time series pattern over centered and scaled remotely sensed data for 2019-2020.

In the attempt of uncovering potential associations among the occurrence of different oviposition temporal patterns and dengue cases, we aggregated data by neighborhood. To estimate dengue incidence over 10,000 people per neighborhood, we first estimated population count by adding up the 100-m resolution data UN adjusted produced by WorldPop (https://www.worldpop.org/, accessed on 30 June 2021) for 2020 in each neighborhood of the city of Córdoba. We then obtained the percentage of pixels belonging to different temporal patterns in each neighborhood. Afterwards, we compared incidences and percentage of each oviposition temporal pattern along all neighborhoods with at least 1 dengue case by means of Spearman rank correlation analysis.

All GRASS and R scripts used to perform data processing, analysis and visualization are available at: https://github.com/veroandreo/mosquito-ts-clust.

Average egg counts varied markedly both in time and space ( From all the clustering configurations tested, we obtained the best distance-centroid combination for each period considered ( Table 2 ). The selection of the best combination of parameters within each distance-centroid combination and period was made according to CVIs and majority vote. We then performed a second round of votes to obtain the best clustering configuration per period. All the best configurations split the 143 ovitraps time series in 3 groups (i.e., the minimum k tested). Two of the selected configurations are based on DTW clustering and the third one corresponds to K-shape clustering (SBD distance + shape extraction centroid). Both cluster sizes and cluster spatial distribution differed according to periods (Figure 3) . Generally, however, the algorithms identified one large group and the other two much less so ( Table 2) . (Table 3) . Median starting date, however, was a bit more variable, determining longer or shorter oviposition seasons in each case. Remarkably, the median start date in 2018-2019 season was around one month earlier than in the previous season (Table 3) . Hence, the median duration of oviposition season was ≈20-30 days longer in 2018-2019. The mean maximum number of eggs is 2 to 3 times higher in 2018-2019 than in the previous season, while the median date of the peak is approximately the same.

The difference among seasons was also evident in the plots showing mean egg counts and 10 and 90 percentiles for each cluster in the three time periods considered. Egg counts were much higher in 2018-2019 than the year before ( Figure 4) . Statistically, however, only in 2017-2018 were there significant differences among clusters in terms of egg counts.

Regarding the shape of the curves, both in 2017-2018 and 2018-2019, there seems to be a temporal pattern with lower abundances and no clear peak (group 2 in 2017-2018 and group 1 in 2018-2019, Figure 4 ). In 2017-2018, among the groups with the highest egg counts, group 1 appears to maintain higher abundances over a longer period than group 3 (Figure 4b) . In 2018-2019, meanwhile, group 2 reaches very high counts earlier than group 3, though the curves are quite similar in their shapes (Figure 4c ). This latter difference is clear when observing the centroids that resulted selected. Centroids of clusters 1, 2 and 3 for each period studied are presented in Figure S3 in Supplementary Materials. Comparisons among clusterings for the different time periods showed that the most similar groups are those from the full series and season 2018-2019 according to several dissimilarity measures (see Table S3 in Supplementary Materials). Indeed, some houses are grouped differently in different years. Figure S4 in Supplementary Materials shows the flow among groups in different periods.

To try to understand if cluster type was related to any particular environmental condition surrounding ovitraps, we ran random forest classifications with the variables derived from Sentinel 2 image analysis and predicted the results over space ( Figure 5 ). In general, results for 2018-2019 were better than those for 2017-2018 both in terms of overall accuracy ( Table 4) Since the season 2019-2020 was rather similar to 2018-2019 in terms of weather and environmental data (see Figures S6 and S7 in Supplementary Materials), we used 2018-2019 RF model to predict oviposition temporal pattern over remotely sensed data for 2019-2020, the dengue outbreak season (Figure 7) . Since the RF model trained with 50-m radius buffer sizes had a higher training and testing OA and seems to relate better to the spatial configuration of the city too, we used its prediction over space to draw further inferences. The spatial predictions for season 2019-2020, similar to that of 2018-2019, showed a matrix mainly composed of the second temporal pattern, with group 3 mostly within the urban fabric and group 1 in the outskirts or border areas of the city as well as green areas within the city like river banks and the park (Figures 4c and 7 for reference regarding the temporal patterns).

The distribution of dengue autochthonous cases and incidence in the city of Córdoba is shown in Figure 8 . Though cases were recorded all over the urban area, there is a clear concentration in the south-west quadrant with some hotspots towards the east, which is also reflected in the raw incidences. In order to relate dengue occurrence to oviposition temporal patterns for the season 2019-2020, we estimated dengue incidence over 10,000 people and we obtained the percentage of pixels belonging to different temporal patterns in each neighborhood. The Spearman correlations among dengue incidences on a neighborhood basis and the proportion of pixels representing the three different temporal patterns found did not show tight associations (Figure 9 ). Indeed, only the temporal pattern 3 (see Figure 4c ) appeared weakly but positively and significantly related to dengue incidence (ρ = 0.15). Higher proportions of this pattern might eventually imply higher dengue incidence. The other temporal patterns, namely 1 and 2 (Figure 4c) , showed negative and barely negative but non-significant correlation with dengue incidences (ρ = −0.13 and −0.03, respectively). In any case, the second temporal pattern is the most commonly found and showed an earlier peak with higher mean egg counts. 

This contribution presents a novel approach that combines time series clustering of vector data, remote sensing data and, machine learning to relate oviposition temporal patterns with occurrence of dengue cases. For the different time periods studied, we found three groups of temporal patterns that showed association with environmental features derived from remote sensing data such as land cover diversity, variability in vegetation and water/humidity indices and, heterogeneity as measured by texture estimated over buffer areas of 50-and 100-m radii surrounding ovitraps. Dengue incidence on a neighborhood basis showed however weak association with the percentage of pixels belonging to different temporal patterns.

Consistently with several other studies, we found that oviposition was highly variable though a seasonal pattern was still evident in all cases even with marked differences among periods regarding egg counts and length of seasons [41] [42] [43] [44] . Within this variability, our time series clustering approach, allowed us to identify 3 different temporal patterns over the city in all 3 periods considered (full period, 2017-2018 and 2018-2019). Most algorithms, however, identified one large group and two smaller ones, similarly to the preliminary results presented in [32] .

Comparing clusters' features among seasons (2017-2018 and 2018-2019), we observed that the main differences laid in length of the mosquito season, given by a 20-30 days earlier onset in 2018-2019, and mean egg counts that was two to three times higher in 2018-2019. In terms of time series curve shapes, both seasons displayed a less abundant temporal pattern of low activity and no marked peak (group 2 and group 1, in 2017-2018 and 2018-2019, respectively), and groups with higher egg counts that mainly differ in peak timing, length and height. These latter patterns deserve some attention since they might indicate higher risk if we consider that more oviposition might imply a higher biting rate [45] . The temporal patterns named group 1 in 2017-2018 and group 2 in 2018-2019 are the most common patterns (i.e., higher frequency) and given their characteristics would represent the riskiest conditions in their respective seasons. Indeed, 2018-2019 group 2 implies 2-3 times more oviposition that group 1 in 2017-2018. Group 1 for 2017-2018 displays the highest peak, with oviposition increasing early and fast and the longest period with high egg counts. Similarly, group 2 for 2018-2019 denotes the highest egg counts and the steepest slope (especially evident from Figure S3 in Supplementary Materials). In both seasons, there is a third group which represented a sort of intermediate state, but still more similar to riskier patterns and with similar egg counts. Importantly, the percentage cover of this pattern was the only one that showed a positive relationship with dengue incidence by neighborhood. The temporal groups might then provide guidance on where and when to prioritize epidemiological surveillance actions and promotional campaigns to eliminate mosquito breeding sites.

The spatial distribution of temporal patterns over the city was also variable, with all 3 patterns appearing in different areas of the city. This suggests that they might indeed be explained by local factors, i.e., the landscape around the ovitrap. Notably, despite the interannual difference in terms of oviposition activity, the predictions denote a spatial dominance of the temporal patterns representing potentially higher or intermediate risk Several studies have attempted to uncover the landscape variables determining breeding hotspots within cities by means of remote sensing data of different spatial resolution. Some have found that heterogeneity and vegetated/urbanized cover are highly relevant to explain Ae. aegypti infestation levels [18] [19] [20] [46] [47] [48] . Importantly, it appears that high environmental variability and heterogeneity in the surroundings of ovitraps might be good predictors of suitability. This can be related to the possibility of finding more potential breeding sites options for females to lay eggs, more places to hide and rest and/or more biting opportunities. Indeed, the best results in terms of train and test overall accuracy, pointed to texture measures and variability in vegetation and water/humidity indices in a radius of 50 m as the most important predictors of oviposition temporal patterns in 2018-2019. In fact, the surroundings of ovitraps belonging to groups 2 and 3 in 2018-2019 were characterized by higher entropy values implying more complex variability, lower IDM and ASM that stand for homogeneity and uniformity, respectively and, higher standard deviation in water and vegetation indices (see Figure S5 ). Furthermore, results for 2017-2018 season pointed to several diversity measures and average conditions in vegetation and water/humidity indices but OA were much lower, especially for the testing set.

Both in terms of weather and environment, seasons 2018-2019 and 2019-2020 were similar; hence, we used 2018-2019 models to predict oviposition temporal patterns for the season of dengue outbreak. Indeed, predictions were very similar. The relationship among dengue incidence and percentage of different temporal patterns per neighborhood resulted inconclusive. There is however a weak and positive association of incidence and proportion of the 3rd oviposition temporal pattern that would require further refinement with other modeling approaches that were out of the scope of this contribution. Indeed, previous studies reviewing the relationship between dengue and vector indices report highly variable results [12, 49] . Still, Cromwell et al. found that entomological information from temporal studies was related to dengue seroconversion while data from cross-sectional studies were not. Hence, sustaining vector surveillance in space and time will likely be beneficial in understanding the vector-disease relationship.

The approach we propose here is aimed at contributing to the understanding of mosquito activity patterns within urban areas in contrast to other studies that aggregate both mosquito and environmental data and consider cities as homogeneous packs and predict mean oviposition or abundance values [22, 23, 43, 50, 51] . In this regard, the only studies assuming differences within cities regarding different temporal patterns are [24, 32] . The former, however, does not use time series clustering but a regular k-means in order to group data before training deep learning LSTM networks and forecasting adult mosquito abundance by group. Meanwhile, in our previous contribution [32] , we raised some potential pitfalls that have been addressed to some extent here. In fact, we used different buffer sizes, we used aggregated remote sensing data instead of just one scene and we performed the analysis with effectively two years of data instead of only one. Furthermore, given the occurrence of a dengue outbreak, we were able to somehow relate the oviposition pattern distribution with dengue cases occurrence, though we were not able to assess any thresholds. This would still require further investigation because given the COVID-19 pandemic and suppression of many activities, we could not have simultaneous oviposition and dengue data. In any case, a similar analysis in places with usual circulation of dengue virus would be relevant as to understand if there are thresholds that might indicate risk and how different groups identified by clustering relate to risk.

It is important to remark that there might be some limitations to our approach that could influence the results. Indeed, since the first step in our workflow can be considered a data driven analysis, we do not have ground truth data to validate the groups obtained. Moreover, because of the limited amount of time series per cluster we used up-sampling before training random forest classification in order to have somehow balanced groups. This might have influenced model's fit as seen in the comparison among training and test accuracy. Similarly, the lower proportion of surface predicted for the lowest activity groups (2 in 2017-2018 and 1 in 2018-2019) might be related to the limited number of time series that were clustered together in these groups. In this sense, the distribution of monitoring ovitraps over the city of Córdoba has changed for season 2020-2021, and those data could be used in further analysis as an independent set to assess predictions. Importantly, because our study focuses on land cover influences, weather variables such as temperature and precipitation were not included. We do not discard, however, that micro-climate differences between ovitraps locations could be relevant to explain the variations in temporal patterns. However, while we might obtain land surface temperature from Landsat at 30 m every 16 days at best, there is not yet a ready-to-use gridded precipitation product with such resolution, and still, better temporal revisit would be needed anyway. Furthermore, our analysis as well as most other studies, cannot account for human behavior in the surrounding of monitoring stations, i.e., plant watering/container re-filling, availability of other containers inside and outside the dwelling, breeding sites removal, pets tipping over containers, etc. However, we pledge that to be able to intervene in a timely manner, we need to study within city variations and micro-habitat mosquito preferences, especially in large cities such as Córdoba and in the face of limited resources. In an attempt to go into this direction, we are currently analyzing the environmental variability at very high spatial resolution (50 cm) in order to suggest a distribution of trapping stations that accounts for such within city variations.

In the present work, with the ultimate goal of predicting risk in places with no mosquito field data, we have explored a novel approach that combines time series clustering of weekly ovitrap records, satellite imagery and machine learning classification tools to assess the association of temporal patterns with environmental variables and finally test the relationship between oviposition and dengue incidence.

The temporal patterns found denoted differences regarding egg counts, peak existence, timing, length and height. The temporal clusters showed association with environmental features, such as land cover diversity, variability in vegetation and water/humidity and texture indices estimated from remote sensing data. Dengue incidence on a neighborhood basis showed a weak but positive association with the percentage of pixels belonging to a temporal pattern occurring only within the urban area and with the highest average of maximum egg counts.

The understanding of the spatial distribution of temporal patterns and their environmental determinants becomes highly relevant to guide the timing and allocation of interventions, especially if we consider that more oviposition might imply higher female biting rates in previous days within a buffer area of a certain ovitrap. Further investigation is still needed though to incorporate other determinants (especially related to human behavior) not considered here. In any case, sustaining long term vector surveillance programs is key for a better understanding of the vector-environment-disease system.

Aedes aegypti, Aedes albopictus, and dengue in Argentina: Current knowledge and future directions. Memórias do Instituto Oswaldo Cruz

Climate change may enable Aedes aegypti infestation in major European cities by 2100

Aedes aegypti spreading in South America: New coldest and southernmost records. Memórias do Instituto Oswaldo Cruz 2020

History of domestication and spread of Aedes aegypti-A Review. Memórias do Instituto Oswaldo Cruz

Urbanization creates diverse aquatic habitats for immature mosquitoes in urban areas

The global burden of dengue: An analysis from the Global Burden of Disease Study

Outbreak of indigenous dengue in the Buenos Aires Metropolitan Area. Experience of the F

Spatio-temporal dynamics of dengue 2009 outbreak in Córdoba City

Analytical report of the 2016 dengue outbreak in Córdoba city

Arbovirus emergence in the temperate city of

Boletín Integrado de Vigilancia

Is Dengue Vector Control Deficient in Effectiveness or Evidence? Systematic Review and Meta-analysis

Characteristics of the spatial pattern of the Dengue vector, Aedes aegypti

Systematic Review: Land Cover, Meteorological, and Socioeconomic Determinants of Aedes Mosquito Habitat for Risk Mapping

A Mapping Review on Urban Landscape Factors of Dengue Retrieved from Earth Observation Data, GIS Techniques, and Survey Questionnaires

An operative dengue risk stratification system in Argentina based on geospatial technology

Modelling the distribution of the vector Aedes aegypti in a central Argentine city: Modelling Aedes aegypti distribution

Could land cover influence Aedes aegypti mosquito populations?

Predicting Aedes aegypti infestation using landscape and thermal features

Towards a workflow for operational mapping of Aedes aegypti at urban scale based on remote sensing

Implementation of a proactive system to monitor Aedes aegypti populations using open access historical and forecasted meteorological data

Environmental Data to Assess Chikungunya, Dengue, and Zika Diseases Through Aedes (Stegomia) aegypti Oviposition Activity Estimation

Modeling Dengue vector population using remotely sensed data and machine learning

Dengue Vector Population Forecasting Using Multisource Earth Observation Products and Recurrent Neural Networks

Identifying Favorable Spatio-Temporal Conditions for West Nile Virus Outbreaks by Co-Clustering of Modis LST Indices Time Series

Urban environmental clustering to assess the spatial dynamics of Aedes aegypti breeding sites

Time-series clustering-A decade review

Intertemporal Similarity of Economic Time Series: An Application of Dynamic Time Warping

Sentinel-2 cropland mapping using pixel-based and object-based time-weighted dynamic time warping analysis

Clustering Groundwater Level Time Series of the Exploited Almonte-Marismas Aquifer in Southwest Spain

The Dynamic-Time-Warping-based k-means++ clustering and its application in phenoregion delineation

Proceedings of the 2019 XVIII Workshop on Information Processing and Control (RPIC)

National Population, Homes and Houses Census; Censo Nacional de Población, Hogares y Viviendas

Operational satellite-based temporal modelling of Aedes population in Argentina

Geographic Resources Analysis Support System (GRASS GIS) Software; Version 7.8; Open Source Geospatial Foundation

An extensive comparative study of cluster validity indices

Time Series Clustering Along with Optimizations for the Dynamic Time Warping Distance

R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing

Cluster Ensembles; R Package Version 0.3-59; 2021

Classification and Regression Training; R Package Version 6.0-88

Oviposition activity and seasonal pattern of a population of Aedes (Stegomyia) aegypti (L.) (Diptera: Culicidae) in subtropical Argentina. Memórias do Instituto Oswaldo Cruz

Prevention of Dengue Outbreaks Through Aedes aegypti Oviposition Activity Forecasting Method. Vector-Borne Zoonotic Dis

Exploring satellite based temporal forecast modelling of Aedes aegypti oviposition from an operational perspective

Understanding the role of temporal variation of environmental variables in predicting Aedes aegypti oviposition activity in a temperate region of Argentina

From eggs to bites: Do ovitrap data provide reliable estimates of Aedes albopictus biting females? PeerJ

Fine scale spatial urban land cover factors associated with adult mosquito abundance and risk in Tucson, Arizona

Spatial pattern evolution of Aedes aegypti breeding sites in an Argentinean city without a dengue vector control programme

An operational machine learning approach to predict mosquito abundance based on socioeconomic and landscape patterns

The relationship between entomological indicators of Aedes aegypti abundance and dengue virus infection

Application of an artificial neural network (ANN) model for predicting mosquito abundances in urban areas

Temporal Dynamics and Spatial Patterns of Aedes aegypti Breeding Sites, in the Context of a Dengue Control Program in Tartagal

The authors would like to thank the European Space Agency (ESA) for the freely available Sentinel 2 imagery and the Health Ministry of Córdoba province for the mosquito data. V.A. wants to especially acknowledge all free and open source software (FOSS) developers, users and enthusiasts whose contributions made this work possible.

The authors declare no conflict of interest.

The following are available online at https://www.mdpi.com/article/ 10.3390/insects12100919/s1: Table S1 : List of Sentinel 2 scenes used. Table S2 : Variables derived from satellite image analysis used for modeling the cluster type in Córdoba (Argentina). Table S3 : Clustering comparisons by means of dissimilarity measures. Figure