# Spatiotemporal modeling of Aedes aegypti risk: enhancing dengue virus control through meteorological and remote sensing data in French Guiana

Sarah Bailly, Vanessa Machault, Samuel Beneteau, Philippe Palany, Camille Fritzell, Romain Girod, Jean-Pierre Lacaux, Philippe Quénel, Claude Flamand

► To cite this version:

Sarah Bailly, Vanessa Machault, Samuel Beneteau, Philippe Palany, Camille Fritzell, et al. Spatiotemporal modeling of Aedes aegypti risk: enhancing dengue virus control through meteorological and remote sensing data in French Guiana. Pathogens, 2024, 13 (9), pp.738. 10.3390/pathogens13090738. hal-04686317

HAL Id: hal-04686317

https://hal.science/hal-04686317v1

Submitted on 4 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Spatiotemporal Modeling of *Aedes aegypti* Risk: Enhancing Dengue Virus Control through Meteorological and Remote Sensing Data in French Guiana

Sarah Bailly <fnref n="1" />,† ID, Vanessa Machault <fnref n="2" />,‡, Samuel Beneteau <fnref n="1" /> ID, Philippe Palany <fnref n="3" />, Camille Fritzell <fnref n="1" />,§, Romain Girod <fnref n="4" />,||, Jean-Pierre Lacaux <fnref n="2" />, Philippe Quénel <fnref n="1" />,<fnref n="5" /> ID and Claude Flamand <fnref n="1" />,<fnref n="6" />,<fnref n="7" />,* ID

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Bailly, S.; Machault, V.; Beneteau, S.; Palany, P.; Fritzell, C.; Girod, R.; Lacaux, J.-P.; Quénel, P.; Flamand, C. Spatiotemporal Modeling of *Aedes aegypti* Risk: Enhancing Dengue Virus Control through Meteorological and Remote Sensing Data in French Guiana. *Pathogens* **2024**, *13*, 738. https://doi.org/10.3390/pathogens13090738

Academic Editor: Yuying Liang

Received: 15 July 2024

Revised: 16 August 2024

Accepted: 26 August 2024

Published: 29 August 2024

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2024 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** French Guiana lacks a dedicated model for developing an early warning system tailored to its entomological contexts. We employed a spatiotemporal modeling approach to predict the risk of *Aedes aegypti* larvae presence in local households in French Guiana. The model integrated field data on larvae, environmental data obtained from very high-spatial-resolution Pleiades imagery, and meteorological data collected from September 2011 to February 2013 in an urban area of French Guiana. The identified environmental and meteorological factors were used to generate dynamic maps with high spatial and temporal resolution. The study collected larval data from 261 different surveyed houses, with each house being surveyed between one and three times. Of the observations, 41% were positive for the presence of *Aedes aegypti* larvae. We modeled the *Aedes* larvae risk within a radius of 50 to 200 m around houses using six explanatory variables and extrapolated the findings to other urban municipalities during the 2020 dengue epidemic in French Guiana. This study highlights the potential of spatiotemporal modeling approaches to predict and monitor the evolution of vector-borne disease transmission risk, representing a major opportunity to monitor the evolution of vector risk and provide valuable information for public health authorities.

**Keywords:** dengue virus; *Aedes aegypti*; spatiotemporal modeling; remote sensing; vector control; French Guiana

1. Introduction

Arboviral diseases, transmitted by infected arthropods, have emerged as a significant public health concern [1–4]. French Guiana, an overseas French department with a population of 300,000 in Northeast South America, has faced multiple outbreaks of dengue, chikungunya, and Zika viruses in recent years, with Ae. aegypti being the sole known

<footnote n="1">¹ Epidemiology Unit, Institut Pasteur in French Guiana, Cayenne 97306, French Guiana; bailly.sarah03@gmail.com (S.B.); samuel.beneteau@ntymail.com (S.B.); camillefritzell@gmail.com (C.F.); philippe.quenel@ehesp.fr (P.Q.)</footnote>

<footnote n="2">² Aerology Laboratory, Observatoire Midi-Pyrénées (OMP), Université Paul Sabatier, 31062 Toulouse, France; vanessamachault@yahoo.com.br (V.M.)</footnote>

<footnote n="3">³ Météo-France, Direction Antilles-Guyane, Fort-de-France 97262, Martinique; philippe.palany@gmail.com</footnote>

<footnote n="4">⁴ Medical Entomology Unit, Institut Pasteur in French Guiana, Cayenne 97306, French Guiana; romain.girod@pasteur.fr</footnote>

<footnote n="5">5 University Rennes, Inserm, EHESP, Irset (Institut de Recherche En Santé, Environnement et Travail)—UMR-S 1085, 35000 Rennes, France</footnote>

<footnote n="6">⁶ Epidemiology and Public Health Unit, Institut Pasteur in Cambodia, Phnom Penh 12201, Cambodia</footnote>

<footnote n="7">⁷ Mathematical Modelling of Infectious Diseases Unit, Institut Pasteur, UMR2000, CNRS, 75015 Paris, France</footnote>

<footnote>* Correspondence: cflamand@pasteur-kh.org; Tel.: +855-12333650</footnote>

<footnote>† Current address: Regional Health Observatory of French Guiana, Cayenne 97300, French Guiana.</footnote>

<footnote>‡ Current address: National Institute for Health and Medical Research (INSERM), UMR 219, Research Institute for Sustainable Development (IRD) EMR 271, Bordeaux Population Health Centre, Bordeaux University, 33600 Bordeaux, France.</footnote>

<footnote>§ Current address: ALIMA—The Alliance for International Medical Action, 75000 Paris, France.</footnote>

<footnote>|| Current address: Scientific Direction, Institut Pasteur, 75015 Paris, France.</footnote>
vector responsible for transmission in the country [5–11]. After a period of no dengue virus circulation between 2015 and 2018, dengue virus activity gradually intensified, culminating in a new epidemic in 2020 involving two serotypes [12,13], resulting in 12,300 clinical cases and four deaths [14].

While vaccines have been developed for some arboviral disease prevention, vector control remains the cornerstone of prevention efforts. In French Guiana, control efforts involve indoor and outdoor deltamethrin spraying to target adult mosquitoes in homes with human cases and within a 100 m radius. Additionally, mechanical and chemical methods are used to eliminate larval breeding sites [15,16]. However, these interventions can be limited by delayed responses and the challenge of accurately identifying contamination sites. Effective mitigation of high transmission risks necessitates the continuous collection of entomological and epidemiological data at the local level. Unfortunately, the intermittent collection of entomological data, due to the substantial resources required for continuous monitoring, often results in imprecise information, hindering the accurate identification of high-risk areas. Consequently, detailed risk maps are essential for pinpointing areas with a high probability of transmission from larval breeding sites, thereby enabling precision public health interventions to optimize control efforts. The continuous collection of fine-scale data, combined with epidemiological information, is crucial for enhancing the effectiveness of control measures and improving our understanding of vector-borne disease transmission.

The correlation between spatial and temporal patterns of *Ae. aegypti* larval distribution and dengue risk has been extensively studied and debated over the decades [17,18]. While some studies have identified significant links [19,20], traditional larval indices, such as the Breteau Index and House Index, have shown inconsistent correlations with dengue transmission [18,21]. This inconsistency raises questions about their reliability as universal predictors of outbreaks, as they frequently fail to establish clear, quantifiable associations for predicting dengue epidemics. The variability in these correlations across different geographic regions and temporal contexts underscores the need for more precise and context-specific models. Such models are crucial for accurately forecasting dengue risk and enabling more effective and targeted public health interventions.

Weather and climate conditions such as rainfall, relative humidity, and temperature [5,6,22], as well as environmental conditions including vegetation indices, land use, and urbanization rate [23,24] play pivotal roles in the evolution of dengue epidemics. These conditions exhibit spatial and temporal heterogeneity. Remote sensing tools, as a valuable resource for collecting environmental and meteorological data with various spatial, spectral, and temporal resolutions, are employed to monitor these conditions at global, regional, and local scales. These tools, combined with different machine learning algorithms, are used to detect associations with entomological indices. For instance, in a study conducted in 2014 in Martinique, Machault et al. used logistic regression to model the presence of houses potentially hosting *Ae. aegypti* larvae [25]. The presence of larvae was significantly associated with maximum humidity and asphalt surface predictors. In Thailand, Sarfraz et al. used a decision tree to predict the probability of vector reproductive habitats using temperature, rainfall, human density, land cover, and elevation factors [26]. Others researchers have utilized time-series [27] or spatial analysis [28] to map abundances and explore correlations with reported cases. Previous studies have demonstrated that the effects of climatic parameters on dengue epidemic incidence can vary significantly from one site to another, depending on the local epidemiological context. Despite the importance of dengue [5,7,9,12,29,30] and the past epidemics transmitted by *Ae. aegypti* in French Guiana [8,31,32], there is currently no model designed to develop an early warning system tailored to the entomological context. This gap is particularly critical given the need for precise and timely interventions to control the spread of dengue. By utilizing historical data, we aim to demonstrate the potential of a straightforward predictive model to enhance surveillance, target high-risk zones, and improve intervention strategies for arbovirus transmission in French Guiana. To address this need, our study focuses on
understanding and modeling the fine spatial and temporal dynamics of *Ae. aegypti* larvae presence based on meteorological and environmental factors. Our research is centered in Matoury, a municipality on the French Guiana coast near Cayenne, to capture the intricate interactions between climatic conditions and mosquito breeding patterns. Additionally, we applied this larval model to all houses in the municipalities of Cayenne, Rémire-Montjoly, and Matoury during the 2020-2021 dengue epidemic [33]. This extrapolation aimed to evaluate the model's ability to identify areas at high risk of dengue transmission across different municipalities.

**2. Materials and Methods**

**2.1. Study Area**

French Guiana, located in South America, lies north of the equator along the Atlantic Ocean. The climate is equatorial, characterized by high temperature and humidity. Relative humidity rarely drops below 80% and average monthly temperatures remain around 27 °C throughout the year. Rainfall exhibits significant seasonal variation due to the influence of migration from the Intertropical Convergence Zone (ZCIT). The average annual rainfall is approximately 3 m, with the majority occurring during the rainy season from December to June. Conversely, the dry season, spanning from July to November, experiences lower levels of rainfall. There are considerable fluctuations in annual rainfall accumulation from year to year, influenced by El Niño. During El Niño periods, French Guiana experiences precipitation deficits, while La Niña events lead to higher annual rainfall accumulations [34,35]. Environmental conditions exhibit more heterogeneity across the territory compared to meteorological conditions. Ninety percent of the territory is covered by the Amazon rainforest, contributing to the region's rich biodiversity. Most of the population resides along the coast, where the three main municipalities are located: Cayenne and its surroundings, Kourou, and Saint-Laurent-du-Maroni municipalities.

**2.2. Data Collection**

We conducted an entomological monitoring study over an 18-month period from September 2011 to February 2013 in the urban center of Matoury, (Figure 1), an adjacent municipality to the main city of Cayenne. At the time of the study, Matoury had a population of 29,712 inhabitants, covering an area of 137 km², resulting in a population density of 217 inhabitants/km². The town's urbanization was heterogeneous, surrounded by dense forests, with the central area comprising 875 easily accessible houses, isolated from the rest of the town, representing urbanization in French Guiana. The monitoring period included 16 months of inter-epidemic phase and 2 months of dengue epidemics at the end of the study period (Figure 2). Each month, a 3-day monitoring session was conducted, during which we randomly selected twenty houses using a sampling plan with replacement. The objective was to identify and quantify the presence and types of water-holding containers, as well as to assess the presence or absence of larvae in both domestic and peridomestic containers.

During the study period, we collected daily meteorological data from the Météo France "Félix Eboué" station at Matoury airport including rainfall, minimum, average, and maximum temperatures, temperature amplitude, relative and absolute humidity, sea-level pressure, vapor pressure, global radiation, and sunshine duration, among others. For each trapping date, we averaged or cumulated 243 variables over different previous periods: 5, 10, 15, 20, 25, and 30 days, considering lag effects by integrating characteristics calculated over the last 5-10, 5-15, 5-20, 5-25, and 5-30 days. To enhance spatial rainfall information, we used satellite observations from the Tropical Rainfall Measuring Mission (TRMM). We obtained data from the nearest pixel to the area or the average of surrounding pixels.

The Pleiades satellite constellation provided high-resolution images (50 cm) with frequent updates, covering every point on the globe every three to four days. We acquired two pairs of optical Pleiades images, covering the dry season and the beginning of the rainy season, on 23 September 2012 and 10 December 2013, respectively. These images were
projected in WGS 84, UTM Zone 22 N, and geometrically corrected using the SRTM 90 m database by NASA. We applied radiometric, atmospheric, and geometric corrections to the images, processed using ENVI 5.1 and ENVI EX (Exelis Visual Information Solutions).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** French Guiana, the studied area, and the center of Matoury. Source IGN—Open Street Map.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Number of confirmed cases collected by the regional system of epidemiological surveillance and entomological collection period in Matoury, French Guiana.</figcaption></figure>
We used QGIS version 2.12.3 with Orfeo Toolbox (OTB) to calculate 14 indices grouped into three categories, including vegetation (NDVI, TNDVI, RVI, SAVI, MSAVI2, GEMI, ARVI, AVI, IPVI), land use (RI, IB, IB2), and water (NDWI2, NDTI) based on spectral characteristics related to physical and biological features. We calculated atmospheric-effect indices (e.g., GEMI) before applying atmospheric correction to the images.

For land cover maps, we classified each image into eight classes: building, pool, asphalt, water, bare ground, light vegetation, green grass, and tree. We employed supervised maximum likelihood pixel classification in ENVI 5.1, using training polygons digitized by an operator based on spectral signatures of each class. Each pixel was assigned to the class with the highest probability based on these spectral signatures. The kappa coefficient, a measure of classification accuracy, was 0.93 for each image, indicating good agreement between the resulting classes and validation areas.

We performed object-oriented classification to enhance classification accuracy and avoid class confusion. This approach considers spectral information and spatial environment. We conducted segmentation and merging of objects based on rule operations in the ENVI FX module. Using confusion trees between pixel and object classifications, we better classified objects like blue or white roofs and swimming pools. For each house, identified by its centroid (from IGN cadastral data), we extracted various information from the immediate environment, calculating averages of all indices and determining land use class zones for buffer zones at 10, 20, 30, 40, 50, 100, 200, 300, 400, and 500 m around the house. We also calculated the distance to the first pixel of each class and the area of each land use class around the experimental unit.

We obtained epidemiological data of confirmed dengue cases, spanning from November 2019 to August 2020, within the Cayenne municipality and aggregated into 100 m × 100 m pixels, from the regional epidemiology unit of Santé Publique France as part of routine surveillance and epidemic management activities.

## 2.3. Statistical and Spatial Analyses

We compiled a large and diverse database, consolidating entomological data with 747 environmental and 195 meteorological factors. We employed univariate logistic regressions to model the probability of larvae presence, testing each of the 942 variables as explanatory factors. Variables with *p*-value < 0.25 were retained in the initial set of Boosted Regression Tree (BRT) variables, resulting in 317 explanatory variables. We conducted the analysis of larvae presence/absence using a BRT model, fitting a Bernoulli distribution with a flexible probability structure to characterize the relationship between the variable of interest and numerous explanatory variables. To refine the model and improve interpretability, we prioritized variables that had the most interpretable biological or physical significance and reduced those that were highly correlated with other variables in the model. This process ensured that the retained variables contributed unique and meaningful information to the model. We set specific parameters for the BRT model, including a tree complexity (tc) value of 5, a learning rate (lr) value of 0.001, and an initial number of trees (nt) of 50. Cross-validation was performed with a bag fraction of 0.5, using a portion of observations not used to build the model.

We measured the strength of association between each explanatory variable and the outcome using the ROC (Receiver Operating Characteristic) curve. We extrapolated the resulting model to the municipalities of Cayenne, Rémire-Montjoly, and Matoury during the last dengue epidemic period (from November 2019 to August 2020) to compare the locations of reported cases with the hotspots identified by the model.

We obtained a complete meteorological dataset without missing data from meteorological stations. We used a Pleiades image taken on 20 December 2017 for analysis. We correlated the predicted entomological data and epidemiological variables (geolocated biologically confirmed cases from November 2019 to August 2020) using different grid sizes, including 500 × 500 m and 1000 × 1000 m grids, as well as IRIS (homogeneous geographical and demographic micro-neighborhoods). We visualized these correlations using
choropleth maps, representing the data by discretized ranges of values. We performed all
statistical analyses using R software version 3.0.2.

3. Results

3.1. Overview of Larval Data and Modeling Approach

We collected larval data from 261 different houses, surveyed between one and three times, resulting in 333 observations, with 27 houses declining to participate across the various sessions. Among the participating houses, 41% were positive for *Ae. aegypti* larvae. We developed the final Boosted Regression Tree (BRT) model using six explanatory variables, narrowed down from an initial set of three hundred seventeen obtained via logistic regression. The model incorporated three environmental variables accounting for 59% of the cumulative influence: surface of low vegetation within a 100 m radius (26%), Mean Redness Index (RI) within a 50 m radius (18%), and Normalized Difference Tillage Index (NDTI) within a 200 m radius (15%). Additionally, three weather variables contributed to 41% of the model's cumulative influence: mean temperature over the previous 5 days (19%), minimum temperature over the previous 25 days (12%), and rainfall over the previous 3 days (10%).

Figure 3 illustrates the directions of the associations and thresholds at which predictions change. For instance, a mean temperature of 25.5 °C over the previous 5 days decreases the probability of *Ae. aegypti* larvae presence. Conversely, for RI, beyond the cutoff point of 1.5, the probability of larvae presence increases. Cumulative rainfall over the previous 3 days constitutes a high-probability plateau for the presence of *Ae. aegypti* larvae, between 10 mm and 100 mm.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Graphs of associations between explanatory variables and the presence of larvae.</figcaption></figure>
To account for the dynamic nature of mosquito populations and their environments, we added an entomological situation variable to the six other variables. This variable, representing the average probabilities of larvae presence in houses within a 30 m radius over the previous 7 days, was created using the initial model run and then incorporated into the second modeling run to predict larvae presence/absence for each geolocated house. The area under the ROC curve was 0.72.

## 3.2. Larvae Risk Mapping

During the 10 months of the studied epidemic period (2019–2020), the extrapolated model (Figure 4) assessed the risk of *Ae. aegypti* larvae presence (red) or absence (blue) for each house. No house was at risk every day of the study period, while all buildings appeared at risk for at least 100 days. The positivity rate varied between 3.8% and 99.6% (Table 1). The highest entomological risk occurred in May and June, with larvae presence positivity rates above 99%. Conversely, at the beginning and end of the epidemic period, the positivity rates were lowest, with November at 3.8% and August at 6.6%.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Monthly entomological risk maps from the modeling experiment based on data from November 2019 until August 2020 in the Cayenne area (Cayenne, Matoury, and Rémire-Montjoly municipalities).</figcaption></figure>

<table><caption>**Table 1.** Correlation rate between entomological and epidemiological data.</caption><thead><tr><th></th><th>Number of Positive Cases</th><th>Number of Positive Houses</th><th>Proportion of Predicted Positive Houses (%)</th><th>Correlation Rate</th></tr></thead><tbody><tr><td>November 2019</td><td>3</td><td>1125</td><td>3.8</td><td>0.2</td></tr><tr><td>December 2019</td><td>0</td><td>2339</td><td>7.9</td><td>-</td></tr><tr><td>January 2020</td><td>8</td><td>15,796</td><td>58.6</td><td>0.3</td></tr><tr><td>February 2020</td><td>11</td><td>10,730</td><td>38.6</td><td>0.2</td></tr><tr><td>March 2020</td><td>29</td><td>3707</td><td>9.0</td><td>0.2</td></tr><tr><td>April 2020</td><td>79</td><td>20,874</td><td>72.8</td><td>0.3</td></tr><tr><td>May 2020</td><td>121</td><td>25,610</td><td>99.6</td><td>0.4</td></tr><tr><td>June 2020</td><td>156</td><td>25,608</td><td>99.6</td><td>0.6</td></tr><tr><td>July 2020</td><td>168</td><td>6890</td><td>24.6</td><td>0.5</td></tr><tr><td>August 2020</td><td>122</td><td>2033</td><td>6.6</td><td>0.5</td></tr></tbody></table>
The highest positive correlations between entomological and epidemiological variables were found using a 1000 m × 1000 m grid, with correlations peaking between June and August 2020 (Table 1). Choropleth maps (Figure 5) revealed spatial heterogeneity in entomological and epidemiological results. In November and December 2019, at the epidemic's onset, the number of cases and predicted positive houses was low, with no clear patterns emerging. As case numbers and predicted positive houses increased, distinct patterns appeared, with April, May, and June showing some grids as hotspot areas.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Correlation between the number of cases and the number of positive houses for the months where the model was extrapolated, Cayenne area (Cayenne, Matoury, and Rémire-Montjoly municipalities).</figcaption></figure>

**4. Discussion**

We successfully utilized very high-resolution satellite imagery and daily meteorological field data to identify environmental variables that determine the presence of *Ae. aegypti* larvae. We investigated the close and more distant environment of the houses, considering radius areas of 50, 100, and 200 m. Interestingly, no larger-scale variables were found to be associated with larvae risk, indicating that the biological and physical phenomena leading to the presence of larvae occur primarily in close proximity to houses. During the “wet” season, we observed that the probability of larvae presence increased with the surface area of lightly vegetated land within a 100 m radius. This effect has been previously demonstrated in Martinique, where “sparsely vegetated soil” was identified as a risk factor for water lodges at the same scale [25]. Additionally, our results showed that land characteristics influenced vector risks within a 50 m radius, illustrating that the characteristics of the surroundings of a house, such as the yard or garden, impact the presence or absence of water deposits and the ability of larvae to develop in those deposits. The well-known relationship between land characteristics and vegetation can influence the presence of *Ae. aegypti* by providing shade for their roost [36]. Moreover, moist environments at low temperatures were found to favor larval development, while the supply of nutrients in water breeding sites attracts females for egg laying and larvae feeding [37], resulting in better survival rates and larger females [38,39]. A higher density
of larvae in a location is evidently partly related to the presence of a high density of females. The NDTI in a 200 m buffer indicates bare rooftops and unvegetated surfaces, signifying urbanization. This index plays a positive role up to a threshold, as urbanization provides artificial lodgings and blood meal within this 200 m scale. The identification of NDTI as a significant predictor highlights the importance of urbanization in creating environments conducive to *Ae. aegypti* breeding [40,41]. By capturing the extent of bare surfaces and built-up areas, NDTI serves as a proxy for urban density, where anthropogenic activities increase the availability of breeding sites and suitable conditions for mosquito survival.

Furthermore, we identified three meteorological variables that remain constant in space and represent both the situation on previous days and the seasonal level. The average mean temperature in the preceding five days is a key factor that affects the likelihood of larvae presence, which is highest at around 25.5 °C. High air temperatures are not conducive to adult survival. The average minimum temperature in the preceding 25 days was positively associated with the presence of larvae, with two thresholds observed at approximately 22.5 and 23.5 °C, highlighting that cooler seasons are less favorable for vector presence. Additionally, the cumulative rainfall in the previous 3 days impacted the presence of larvae, quickly reaching a limit. When rainfall exceeded ten millimeters in three days, the probability of the presence of larvae reached its maximum. This finding is particularly relevant in our study area, where breeding sites such as cups, pots, and watering cans were small and quickly filled with water. Overall, the biological and physical significance of the identified predictors underscores their relevance in shaping the spatial and temporal dynamics of *Ae. aegypti* larval presence. These insights provide valuable guidance for targeted vector control strategies, especially in urbanized areas where intervention efforts can be concentrated based on predictive modeling outcomes.

The extrapolation of the model showed that March was less favorable for the presence of larvae than the other months of the rainy season. This month corresponds to the “short March summer”, when rainfall decreases significantly, between the short rainy season (January–February) and the long rainy season (April–June). As soon as the main rainy season starts, the probability of larval presence increases again.

The spatiotemporal risk maps produced in this study are a significant contribution to the development of targeted operational control systems for arboviral diseases in urban areas. While it is widely accepted that entomological risk alone cannot fully predict the occurrence of human cases of arbovirus, as other factors such as socioeconomic level, virus circulation, and protective behaviors also play crucial roles, the presence of vectors remains an essential element in assessing the risk of transmission, especially when associated with virus circulation. Our proposed approach aims to identify areas at high risk of transmission by combining the presence of vector risk and the presence of the virus, as depicted in the choropleth maps. These maps can effectively direct vector control actions to hotspots areas. The core strength of our approach lies in its flexibility. While the specific parameters of the model (such as temperature thresholds and rainfall patterns) are influenced by the local climate of the study area, the overall framework is designed to be adaptable. By recalibrating the model with region-specific data, such as different temperature ranges, rainfall patterns, and vegetation cover, it can be applied to other tropical regions, including those with a tropical monsoon climate. However, there are limitations in the collected entomological data. For instance, we did not consider the productivity of breeding sites for larvae, and it might have been interesting to model only the presence of larvae at the most productive sites. Nevertheless, this may not be a significant issue in the specific case of French Guiana, as productivity was fairly dichotomous: close to 20% at positive sites with buckets or watering can containers and plants, versus close to 0% for all other types of sites. Another limitation of the study pertains to the methodology used to explore the associations between modeled entomological risk and case occurrence. While our study relied on the surveillance system for biologically confirmed cases, it is important to note that this system may not detect all dengue cases, potentially underestimating the actual number of cases in the control houses.
Looking ahead, it will be intriguing to identify geolocalized indicators that are more representative of epidemiological risk than individual confirmed cases, allowing us to better highlight the link between entomological risk and epidemiological risk. Future steps include implementing new strategies to produce high-resolution local maps across the country and integrating the proposed approach into the master plan of information systems dedicated to the surveillance of priority diseases in French Guiana, enabling precision public health actions to be deployed by control operators.

**Author Contributions:** Conceptualization, C.F. (Claude Flamand), V.M., P.P., J.-P.L. and P.Q.; Methodology, C.F. (Claude Flamand), V.M., S.B. (Sarah Bailly) and S.B. (Samuel Beneteau); Formal analysis, V.M., S.B. (Sarah Bailly), S.B. (Samuel Beneteau) and C.F. (Claude Flamand); Investigation, R.G.; Supervision, J.-P.L., P.Q. and C.F. (Claude Flamand); Data curation, V.M., C.F. (Camille Fritzell), S.B. (Sarah Bailly) and S.B. (Samuel Beneteau); Writing—original draft preparation, S.B. (Sarah Bailly), V.M., S.B. (Samuel Beneteau) and C.F. (Claude Flamand); Writing—review and editing: S.B. (Sarah Bailly), V.M., S.B. (Samuel Beneteau), C.F. (Claude Flamand), P.Q., P.P. and R.G.; Visualization: S.B. (Sarah Bailly) and V.M.; Resources: J.-P.L. and R.G.; Funding acquisition: C.F. (Claude Flamand). All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Centre National d'Etudes Spatiales-Terre solide, Océan, Surfaces Continentales, Atmosphère Fund (grant number CNES-TOSCA-4800000720). C.F. acknowledges the European Regional Development Fund under EPI-ARBO grant agreement GY0008695.

**Institutional Review Board Statement:** The collection of epidemiological data was sourced from the national surveillance system and received approval from the Advisory Committee on Information Processing in Material Research in the Field of Health (N°07.148) and authorization by regulatory authorities (CNIL-N°1213498). No individual data or biological samples were collected in this study and used data were systematically aggregated, preventing the identification of individuals included in the study.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data used in this study are available upon request from the authors or the Institut Pasteur of French Guiana—Tel: +594(0)594292600; +594(0)594292618, Contact: https://www.pasteur-cayenne.fr/contact/, accessed on 1 August 2024).

**Acknowledgments:** We extend our deepest gratitude to Vanessa Ardillon from Santé Publique France for her invaluable support in providing the epidemiological data essential for this study. We would like to thank all the participants from the city of Matoury who graciously opened their homes for mosquito collections.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## References

1. Bhatt, S.; Gething, P.W.; Brady, O.J.; Messina, J.P.; Farlow, A.W.; Moyes, C.L.; Drake, J.M.; Brownstein, J.S.; Hoen, A.G.; Sankoh, O.; et al. The Global Distribution and Burden of Dengue. *Nature* **2013**, *496*, 504–507. [CrossRef] [PubMed]

2. Gould, E.; Pettersson, J.; Higgs, S.; Charrel, R.; de Lamballerie, X. Emerging Arboviruses: Why Today? *One Health* **2017**, *4*, 1–13. [CrossRef] [PubMed]

3. Kraemer, M.U.G.; Reiner, R.C.; Brady, O.J.; Messina, J.P.; Gilbert, M.; Pigott, D.M.; Yi, D.; Johnson, K.; Earl, L.; Marczak, L.B.; et al. Past and Future Spread of the Arbovirus Vectors Aedes Aegypti and Aedes Albopictus. *Nat. Microbiol.* **2019**, *4*, 854–863. [CrossRef]

4. Pettersson, J.H.-O.; Eldholm, V.; Seligman, S.J.; Lundkvist, Å.; Falconar, A.K.; Gaunt, M.W.; Musso, D.; Nougairède, A.; Charrel, R.; Gould, E.A.; et al. How Did Zika Virus Emerge in the Pacific Islands and Latin America? *MBio* **2016**, *7*, 10–1128. [CrossRef]

5. Adde, A.; Roucou, P.; Mangeas, M.; Ardillon, V.; Desenclos, J.-C.; Rousset, D.; Girod, R.; Briolant, S.; Quenel, P.; Flamand, C. Predicting Dengue Fever Outbreaks in French Guiana Using Climate Indicators. *PLoS Negl. Trop. Dis.* **2016**, *10*, e0004681. [CrossRef] [PubMed]

6. Flamand, C.; Fabregue, M.; Bringay, S.; Ardillon, V.; Quénel, P.; Desenclos, J.C.; Teisseire, M. Mining Local Climate Data to Assess Spatiotemporal Dengue Fever Epidemic Patterns in French Guiana. *J. Am. Med. Inform. Assoc.* **2014**, *21*, e232–e240. [CrossRef] [PubMed]

7. L'Azou, M.; Taurel, A.F.; Flamand, C.; Quénel, P. Recent Epidemiological Trends of Dengue in the French Territories of the Americas (2000–2012): A Systematic Literature Review. *PLoS Negl. Trop. Dis.* **2014**, *8*, e3235. [CrossRef]
8. Flamand, C.; Bailly, S.; Fritzell, C.; Berthelot, L.; Vanhomwegen, J.; Salje, H.; Paireau, J.; Matheus, S.; Enfissi, A.; Fernandes-Pellerin, S.; et al. Impact of Zika Virus Emergence in French Guiana: A Large General Population Seroprevalence Survey. *J. Infect. Dis.* **2019**, *220*, 1915–1925. [CrossRef]

9. Flamand, C.; Fritzell, C.; Prince, C.; Abboud, P.; Ardillon, V.; Carvalho, L.; Demar, M.; Boukhari, R.; Papaix-Puech, M.; Elenga, N.; et al. Epidemiological Assessment of the Severity of Dengue Epidemics in French Guiana. *PLoS ONE* **2017**, *12*, e0172267. [CrossRef]

10. Fritzell, C.; Raude, J.; Adde, A.; Dusfour, I.; Quenel, P.; Flamand, C. Knowledge, Attitude and Practices of Vector-Borne Disease Prevention during the Emergence of a New Arbovirus: Implications for the Control of Chikungunya Virus in French Guiana. *PLoS Negl. Trop. Dis.* **2016**, *10*, e0005081. [CrossRef]

11. Flamand, C.; Fritzell, C.; Matheus, S.; Dueymes, M.; Carles, G.; Favre, A.; Enfissi, A.; Adde, A.; Demar, M.; Kazanji, M.; et al. The Proportion of Asymptomatic Infections and Spectrum of Disease among Pregnant Women Infected by Zika Virus: Systematic Monitoring in French Guiana, 2016. *Euro Surveill.* **2017**, *22*, 17-00102. [CrossRef]

12. Bailly, S.; Rousset, D.; Fritzell, C.; Hozé, N.; Ben Achour, S.; Berthelot, L.; Enfissi, A.; Vanhomwegen, J.; Salje, H.; Fernandes-Pellerin, S.; et al. Spatial Distribution and Burden of Emerging Arboviruses in French Guiana. *Viruses* **2021**, *13*, 1299. [CrossRef]

13. Nacher, M.; Douine, M.; Gaillet, M.; Flamand, C.; Rousset, D.; Rousseau, C.; Mahdaoui, C.; Carroll, S.; Valdes, A.; Passard, N.; et al. Simultaneous Dengue and COVID-19 Epidemics: Difficult Days Ahead? *PLoS Negl. Trop. Dis.* **2020**, *14*, e0008426. [CrossRef] [PubMed]

14. Santé Publique France Situation Épidémiologique de la Dengue en Guyane. Point au 6 mai 2021. Available online: https://www.santepubliquefrance.fr/regions/guyane/documents/bulletin-regional/2021/situation-epidemiologique-de-la-dengue-en-guyane-point-au-6-mai-2021 (accessed on 22 April 2023).

15. Dusfour, I.; Zorrilla, P.; Guidez, A.; Issaly, J.; Girod, R.; Guillaumot, L.; Robello, C.; Strode, C. Deltamethrin Resistance Mechanisms in Aedes Aegypti Populations from Three French Overseas Territories Worldwide. *PLoS Neglected Trop. Dis.* **2015**, *9*, e0004226. [CrossRef] [PubMed]

16. Epelboin, Y.; Chaney, S.C.; Guidez, A.; Habchi-Hanriot, N.; Talaga, S.; Wang, L.; Dusfour, I. Successes and Failures of Sixty Years of Vector Control in French Guiana: What Is the next Step? *Mem. Inst. Oswaldo Cruz* **2018**, *113*, e170398. [CrossRef] [PubMed]

17. MacCormack-Gelles, B.; Lima Neto, A.S.; Sousa, G.S.; do Nascimento, O.J.; Castro, M.C. Evaluation of the Usefulness of Aedes Aegypti Rapid Larval Surveys to Anticipate Seasonal Dengue Transmission between 2012–2015 in Fortaleza, Brazil. *Acta Trop.* **2020**, *205*, 105391. [CrossRef]

18. Bowman, L.R.; Runge-Ranzinger, S.; McCall, P.J. Assessing the Relationship between Vector Indices and Dengue Transmission: A Systematic Review of the Evidence. *PLoS Negl. Trop. Dis.* **2014**, *8*, e2848. [CrossRef] [PubMed]

19. Pham, H.V.; Doan, H.T.; Phan, T.T.; Minh, N.N. Ecological Factors Associated with Dengue Fever in a Central Highlands Province, Vietnam. *BMC Infect. Dis.* **2011**, *11*, 172. [CrossRef]

20. Sanchez, L.; Vanlerberghe, V.; Alfonso, L.; del Marquetti, M.C.; Guzman, M.G.; Bisset, J.; van der Stuyft, P. Aedes Aegypti Larval Indices and Risk for Dengue Epidemics. *Emerg. Infect. Dis.* **2006**, *12*, 800–806. [CrossRef]

21. Liyanage, P.; Tozan, Y.; Tissera, H.A.; Overgaard, H.J.; Rocklöv, J. Assessing the Associations between Aedes Larval Indices and Dengue Risk in Kalutara District, Sri Lanka: A Hierarchical Time Series Analysis from 2010 to 2019. *Parasit Vectors* **2022**, *15*, 277. [CrossRef]

22. de Thoisy, B.; Duron, O.; Epelboin, L.; Musset, L.; Quénel, P.; Roche, B.; Binetruy, F.; Briolant, S.; Carvalho, L.; Chavy, A.; et al. Ecology, Evolution, and Epidemiology of Zoonotic and Vector-Borne Infectious Diseases in French Guiana: Transdisciplinarity Does Matter to Tackle New Emerging Threats. *Infect. Genet. Evol.* **2021**, *93*, 104916. [CrossRef] [PubMed]

23. Eisen, L.; Lozano-Fuentes, S. Use of Mapping and Spatial and Space-Time Modeling Approaches in Operational Control of Aedes Aegypti and Dengue. *PLoS Negl. Trop. Dis.* **2009**, *3*, e411. [CrossRef] [PubMed]

24. Gao, P.; Pilot, E.; Rehbock, C.; Gontariuk, M.; Doreleijers, S.; Wang, L.; Krafft, T.; Martens, P.; Liu, Q. Land Use and Land Cover Change and Its Impacts on Dengue Dynamics in China: A Systematic Review. *PLoS Negl. Trop. Dis.* **2021**, *15*, e0009879. [CrossRef] [PubMed]

25. Machault, V.; Yébakima, A.; Etienne, M.; Vignolles, C.; Palany, P.; Tourre, Y.M.; Guérécheau, M.; Lacaux, J.-P. Mapping Entomological Dengue Risk Levels in Martinique Using High-Resolution Remote-Sensing Environmental Data. *ISPRS Int. J. Geo-Inf.* **2014**, *3*, 1352–1371. [CrossRef]

26. Sarfraz, M.S.; Tripathi, N.K.; Tipdecho, T.; Thongbu, T.; Kerdthong, P.; Souris, M. Analyzing the Spatio-Temporal Relationship between Dengue Vector Larval Density and Land-Use Using Factor Analysis and Spatial Ring Mapping. *BMC Public Health* **2012**, *12*, 1–19. [CrossRef]

27. Lana, R.M.; Carneiro, T.G.S.; Honório, N.A.; Codeço, C.T. Seasonal and Nonseasonal Dynamics of Aedes Aegypti in Rio de Janeiro, Brazil: Fitting Mathematical Models to Trap Data. *Acta Trop.* **2014**, *129*, 25–32. [CrossRef]

28. de Melo, D.P.O.; Scherrer, L.R.; Eiras, Á.E. Dengue Fever Occurrence and Vector Detection by Larval Survey, Ovitrap and MosquiTRAP: A Space-Time Clusters Analysis. *PLoS ONE* **2012**, *7*, e42125. [CrossRef]

29. Flamand, C.; Quenel, P.; Ardillon, V.; Carvalho, L.; Bringay, S.; Teisseire, M. The Epidemiologic Surveillance of Dengue-Fever in French Guiana: When Achievements Trigger Higher Goals. *Stud. Health Technol. Inform.* **2011**, *169*, 629–633.
30. Hanf, M.; Friedman, E.; Basurko, C.; Roger, A.; Bruncher, P.; Dussart, P.; Flamand, C.; Carles, G.; Buekens, P.; Breart, G.; et al. Dengue Epidemics and Adverse Obstetrical Outcomes in French Guiana: A Semi-Ecological Study. *Trop. Med. Int. Health* **2014**, *19*, 153–158. [CrossRef]

31. Hallet, E.; Flamand, C.; Rousset, D.; Bonifay, T.; Fritzell, C.; Matheus, S.; Dueymes, M.; Ntab, B.; Nacher, M. ZIKA Virus Infection in Pregnant Women in French Guiana: More Precarious-More at Risk. *PLoS Negl. Trop. Dis.* **2020**, *14*, e0008193. [CrossRef]

32. Flamand, C.; Fritzell, C.; Obale, P.; Quenel, P.; Raude, J. The Role of Risk Proximity in the Beliefs and Behaviors Related to Mosquito-Borne Diseases: The Case of Chikungunya in French Guiana. *Am. J. Trop. Med. Hyg.* **2017**, *97*, 344–355. [CrossRef]

33. Bonifay, T.; Le Turnier, P.; Epelboin, Y.; Carvalho, L.; De Thoisy, B.; Djossou, F.; Duchemin, J.-B.; Dussart, P.; Enfissi, A.; Lavergne, A.; et al. Review on Main Arboviruses Circulating on French Guiana, An Ultra-Peripheric European Region in South America. *Viruses* **2023**, *15*, 1268. [CrossRef]

34. Aceituno, P. On the Functioning of the Southern Oscillation in the South American Sector. Part I: Surface Climate. *Mon. Wea. Rev.* **1988**, *116*, 505–524. [CrossRef]

35. Moron, V.; Bigot, S.; Roucou, P. Rainfall Variability in Subequatorial America and Africa and Relationships with the Main Sea-Surface Temperature Modes (1951–1990). *Oceanogr. Lit. Rev.* **1996**, *8*, 771. [CrossRef]

36. Vezzani, D.; Albicócco, A.P. The Effect of Shade on the Container Index and Pupal Productivity of the Mosquitoes Aedes Aegypti and Culex Pipiens Breeding in Artificial Containers. *Med. Vet. Entomol.* **2009**, *23*, 78–84. [CrossRef] [PubMed]

37. Reiskind, M.H.; Greene, K.L.; Lounibos, L.P. Leaf Species Identity and Combination Affect Performance and Oviposition Choice of Two Container Mosquito Species. *Ecol. Entomol.* **2009**, *34*, 447–456. [CrossRef]

38. Tun-Lin, W.; Burkot, T.R.; Kay, B.H. Effects of Temperature and Larval Diet on Development Rates and Survival of the Dengue Vector Aedes Aegypti in North Queensland, Australia. *Med. Vet. Entomol.* **2000**, *14*, 31–37. [CrossRef]

39. Yan, J.; Kibech, R.; Stone, C.M. Differential Effects of Larval and Adult Nutrition on Female Survival, Fecundity, and Size of the Yellow Fever Mosquito, Aedes Aegypti. *Front. Zool.* **2021**, *18*, 10. [CrossRef]

40. Zahouli, J.B.Z.; Koudou, B.G.; Müller, P.; Malone, D.; Tano, Y.; Utzinger, J. Urbanization Is a Main Driver for the Larval Ecology of Aedes Mosquitoes in Arbovirus-Endemic Settings in South-Eastern Côte d'Ivoire. *PLoS Neglected Trop. Dis.* **2017**, *11*, e0005751. [CrossRef]

41. Wilke, A.B.B.; Vasquez, C.; Carvajal, A.; Moreno, M.; Fuller, D.O.; Cardenas, G.; Petrie, W.D.; Beier, J.C. Urbanization Favors the Proliferation of Aedes Aegypti and Culex Quinquefasciatus in Urban Areas of Miami-Dade County, Florida. *Sci. Rep.* **2021**, *11*, 22989. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.