Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Dynamical Mapping of Anopheles darlingi Densities in a Residual Malaria Transmission Area of French Guiana by Using Remote Sensing and Meteorological Data

  • Antoine Adde ,

    adde.antoine@gmail.com

    Affiliation Unité d’Entomologie Médicale, Institut Pasteur de la Guyane, Cayenne, French Guiana

  • Emmanuel Roux,

    Affiliation UMR ESPACE-DEV, Institut de Recherche pour le Développement, Montpellier, France

  • Morgan Mangeas,

    Affiliation UMR ESPACE-DEV, Institut de Recherche pour le Développement, Montpellier, France

  • Nadine Dessay,

    Affiliation UMR ESPACE-DEV, Institut de Recherche pour le Développement, Montpellier, France

  • Mathieu Nacher,

    Affiliation Centre d’Investigation Clinique et Epidémiologie Clinique Antilles Guyane, Centre hospitalier Andrée-Rosemon, Cayenne, French Guiana

  • Isabelle Dusfour,

    Affiliation Unité d’Entomologie Médicale, Institut Pasteur de la Guyane, Cayenne, French Guiana

  • Romain Girod,

    Affiliation Unité d’Entomologie Médicale, Institut Pasteur de la Guyane, Cayenne, French Guiana

  • Sébastien Briolant

    Affiliations Unité d’Entomologie Médicale, Institut Pasteur de la Guyane, Cayenne, French Guiana, Direction Interarmées du Service de Santé en Guyane, Cayenne, French Guiana, Unité de Parasitologie et d’Entomologie Médicale, Institut de Recherche Biomédicale des Armées, Marseille, France, Unité de Recherche en Maladies Infectieuses Tropicales Emergentes, Faculté de Médecine La Timone, Marseille, France

Abstract

Local variation in the density of Anopheles mosquitoes and the risk of exposure to bites are essential to explain the spatial and temporal heterogeneities in the transmission of malaria. Vector distribution is driven by environmental factors. Based on variables derived from satellite imagery and meteorological observations, this study aimed to dynamically model and map the densities of Anopheles darlingi in the municipality of Saint-Georges de l’Oyapock (French Guiana). Longitudinal sampling sessions of An. darlingi densities were conducted between September 2012 and October 2014. Landscape and meteorological data were collected and processed to extract a panel of variables that were potentially related to An. darlingi ecology. Based on these data, a robust methodology was formed to estimate a statistical predictive model of the spatial-temporal variations in the densities of An. darlingi in Saint-Georges de l’Oyapock. The final cross-validated model integrated two landscape variables—dense forest surface and built surface—together with four meteorological variables related to rainfall, evapotranspiration, and the minimal and maximal temperatures. Extrapolation of the model allowed the generation of predictive weekly maps of An. darlingi densities at a resolution of 10-m. Our results supported the use of satellite imagery and meteorological data to predict malaria vector densities. Such fine-scale modeling approach might be a useful tool for health authorities to plan control strategies and social communication in a cost-effective, targeted, and timely manner.

Introduction

Nearly half of the world’s population, i.e., 3.2 billion people, is exposed to the risk of malaria [1]. Malaria is caused by Plasmodium parasites transmitted to humans through the bite of infected Anopheles mosquitoes. At every level of endemicity, the spatial and temporal distributions of malaria transmission are heterogeneous [2]. Local variation in abundance and exposure to the vectors is the key to explaining this heterogeneity [3, 4]; the more individuals that are bitten, the more likely are they to become infected and enhance transmission by infecting new mosquitoes. Accurate information on where and when malaria vectors proliferate is essential for malaria surveillance and elimination since it allows targeted interventions that remarkably increase the efficiency of control measures [5].

Malaria vector distribution is strongly influenced by environmental factors that determine the availability and productivity of Anopheles habitats. Environmental factors can be accurately detected and spatially analyzed using remote sensing techniques to characterize Anopheles ecological preferences, model population densities, and produce hazard maps [69]. The visual nature of maps makes them helpful tools to identify locations where interventions can be targeted [10]. Demand for maps of worldwide malaria-vulnerable areas, e.g., to steer World Health Organization (WHO) international actions, has led to an increasing number of studies aiming to circumscribe the geographical distribution of vectors [1116]. Global maps from these studies constitute intelligent support, depicting the schematic large-scale distribution of the species. However, the spatial resolution of current maps is too coarse to capture the local spatial and temporal dynamics of Anopheles, which are very heterogeneous, even at fine scales [1719]. The high heterogeneity of the Anopheles distribution, combined with major knowledge gaps in anopheline bio-ecology and experience from previous control efforts, suggest that ecological studies and control policies should be tailored to individual areas [20].

French Guiana is an overseas French territory of 230,000 inhabitants located in northern South America. Despite a continual decrease in the annual cases over the past decade (1.8 cases per 1,000 inhabitants in 2015), malaria remains a public health issue [21, 22]. Plasmodium vivax is the predominant malaria parasite species and was responsible for 80% of the diagnosed cases in 2015, with the remainder mainly due to P. falciparum [2124]. Most malaria cases are reported in villages located along the main rivers flowing through the territory and in illegal gold mining areas, which are propitious places for malaria transmission [2325]. As in a large part of the Americas, the main malaria vector in French Guiana has historically been considered to be An. darlingi [2632], from its natural infectivity, anthropophilic behavior, high density during malaria transmission periods, and wide distribution [3337]. Although An. darlingi has been the target of numerous studies, many aspects of its ecology and biology are still unknown. This species is highly adaptable, behaviorally variable, and has been observed in heterogeneous breeding sites in a wide range of ecosystems, including forest, savannah, swamps, and human-disturbed environments [12, 26, 34, 36, 3842].

This study was conducted in Saint-Georges de l’Oyapock, a persistent malaria-endemic municipality of French Guiana where An. darlingi have been found infected by P. vivax [43]. Anopheles darlingi densities were longitudinally monitored during the malaria transmission period (i.e., the September–November dry season) in 2012, 2013, and 2014 at eight different sites. In parallel, appropriate satellite imagery and meteorological data were selected and processed in order to extract a panel of variables potentially related to An. darlingi ecology. Based on these data, this study aimed to dynamically model and map densities of An. darlingi in the municipality of Saint-Georges de l’Oyapock at weekly intervals during the malaria transmission period. Models and density maps of An. darlingi contribute toward controlling malaria transmission by bridging the knowledge gaps in An. darlingi bio-ecology and providing tools for health authorities to effectively apply vector controls.

Materials and Methods

Geography of the study area

Saint-Georges de l’Oyapock (3.89° N, 51.8° W) is a municipality of French Guiana located in the Amazon forest at the eastern border with Brazil along the Oyapock River. The 3,855 inhabitants (data from the French Institute of Statistics and Economic Studies, 2012) are mainly concentrated in the city center and more sporadically in the peripheral Amerindian villages. Contrasts between neighborhoods can be important because the population stems from a great diversity of ethnic groups. Therefore, local habits such as agricultural practices or housing types can vary widely from one area to another. In the “traditional” Amerindian places, wooden houses are mostly built on stilts, and fishing or slash-and-burn agriculture is the main livelihood activities of the populations. In the “modern” places, populations have adopted a western lifestyle. Houses are built with solid materials (concrete, bricks, etc.) and are surrounded by planted and tended gardens. Landscape of the study area is characterized by four main elements—the forest, savanna, urban, and river. The eastern side of the city is bordered by the Oyapock River, and the western side is open to the savanna. The northern and southern boundaries of the town are delimited by dense primary forest. The climate is equatorial: hot, wet, and rainy. The average temperatures range between 26.05°C in February and 27.85°C in October. The mean annual cumulative rainfall is 3,345 mm, with four alternating seasons: a long, rainy season from April to June; a long, dry season from July to December; a short, rainy season from January to February; and a short, dry season in March.

Mosquito collection

Mosquitoes were collected during the dry seasons in 2012 (September 3–November 25; 12 weeks), 2013 (September 2–November 24; 12 weeks), and 2014 (September 1–October 26; 8 weeks) from eight different sites that varied by year (Fig 1). Sites were selected from field expertise and followed a longitudinal transect that maximized heterogeneity of the landscape and lifestyle characteristics of the population (Table 1). Landscape heterogeneities were decisive elements for sampling site selection. Photo interpretation and field expertise were used to select sampling sites such that the environmental richness of the study area was covered as best as possible. Since the lifestyle of the populations impacts the local environment (housing type or livelihood practices described in “Geography of the study area” section), this criterion was also considered. The eight sampling sites were finally decided by considering the best possible environmental criteria and field constraints (field hostility, access to the site, low risk of trapping material theft, etc.). Importantly, the choice of sampling sites was not oriented by the identification of places where we were confident to find numerous An. darlingi. Trapping sessions were conducted on two non-consecutive nights per week, between 6:00 pm and 8:00 am, by using Mosquito Magnet® traps (Woodstream Corporation, Lititz, PA) baited with octenol (MMoct). Such traps have previously been proven useful to monitor the spatial and temporal abundance of malaria vectors in French Guiana [38, 44]. Mosquitoes were stored at -20°C until counting and morphological identification [31, 45, 46]. To achieve a weekly cartographic output and to smooth variation between collection nights, the two bi-weekly observations were averaged to approximate a daily number of specimens and then multiplied by seven to approximate a weekly number of specimens. In all, 165 weekly approximated An. darlingi density records across the eight study sites was finally available for analyses, which was slightly less than our expected maximum of 184, owing to the mechanical failure of traps. Since our activities were not conducted in a protected area (not a national park or a regional nature reserve), no specific permission was required, and field studies did not involve endangered or protected species.

thumbnail
Fig 1. Study area and sampling sites.

Location of collection sites with information about their sampling period, Saint-Georges de l’Oyapock, French Guiana. Aerial photograph acquired in 2006 (BD Ortho® product from IGN, the French National Institute of Geographic and Forest Information).

https://doi.org/10.1371/journal.pone.0164685.g001

Landscape characterization

A SPOT-5 image acquired on October 14, 2012, with four color channels (red, green, near-infra-red, and middle infrared) at 10-m spatial resolution was selected to characterize the landscape of the study area. One image was sufficient to cover the eight trapping sites. However, the presence of clouds required the posterior use of a second SPOT-5 image (July 22, 2013) in order to fill the missing data (5% of the total study area) and obtain a spatially complete product. A land cover map of the study area was produced based on field observations and a supervised training approach with maximum likelihood classification. The classification included five land cover types identified as “built, roads, and bare soils”; “low vegetation”; “forest”; “very dense forest”; and “water.” A 7 × 7 pixel mode filter (i.e., each pixel value being replaced by its most common neighbor in a 7 × 7 cell moving window) was applied to the classification to reduce noise. BD-Topo® 2012 (IGN, the French National Institute of Geographic and Forest Information) was used to separate the “built” surfaces and “roads and bare soil” surfaces, resulting in a six-class land cover map (Fig 2). Land cover map quality was assessed by identifying the actual land cover of the training pixels by using photo-interpretation and local field expertise and applying a 5-fold cross validation procedure. The resulting mean Kappa coefficient was 0.84. A simplified binary classification, with two classes “forest” and “non-forest,” was also produced to provide complementary information for subsequent analyses. Biases from incomplete detection of water bodies, owing to forest canopy and thick vegetation covering the surface of water points were avoided by excluding the “water” land cover class from further statistical analyses. The remaining classes refer to perennial environmental facies. Therefore, changes in land cover during the study period were considered negligible, supporting the use of two satellite images acquired during different months to classify the study area. The Oyapock River and its main tributaries were excluded from the analysis since they are not a suitable habitat for An. darlingi breeding sites. Indeed, the estuarine situation of Saint-Georges de l’Oyapock implied a daily tidal influence and recurrent leaching of the riverbanks.

thumbnail
Fig 2. Land cover map of the study area, Saint-Georges de l’Oyapock, French Guiana.

https://doi.org/10.1371/journal.pone.0164685.g002

Given the lack of knowledge on the specific bio-ecology and behavior of An. darlingi in Saint-Georges de l’Oyapock, an exploratory approach was adopted by extracting a wide range of variables to characterize the landscape. These variables were extracted in a 200-m radius buffer around each trap by using FRAGSTATS software [47]. The radius represented a compromise between a relevant landscape characterization according to the satellite image spatial resolution, the overlap between neighboring buffers (to avoid information redundancy and artificial spatial auto-correlation), and the MMoct constructor attraction radius (up to 50 m). In total, 22 landscape variables were computed for each sampling site (Table 2). Four variables (PAFRAC, SPLIT, TE, and PRD) were computed for each of the two classifications (i.e., the final and simplified forest/non-forest binary classifications) by considering all classes together. Two variables (AREA and EDGE) were computed for each of the land cover classes (i.e., “built,” “roads and bare soils,” “low vegetation,” “forest,” “very dense forest,” “binary forest,” and “binary non-forest”).

Meteorological variables

Daily meteorological records, including rainfall, temperature, relative humidity, solar radiation, and evapotranspiration (ETP) for the years 2012–2014 were obtained from the Meteo-France weather station located in the city center of Saint-Georges de l’Oyapock. ETP was approximated using the Penman–Monteith equation [48]. For each meteorological parameter, a panel of variables (Table 3) was extracted using two different time interval schemes: a cumulative 7-day time interval scheme (i.e., past weather from days 0 to 6, days 0 to 13, days 0 to 20, etc.) and a non-cumulative 7-day time interval scheme (i.e., past weather from days 0 to 6, days 7 to 13, days 14 to 20, etc.). Extraction of meteorological variables ended at the ninth week (day 63) before the beginning of the trapping session to remain within the dry season. A total of 2,214 meteorological variables were extracted for each of the 32 weeks studied.

Statistical model

The spatial and temporal dynamics of An. darlingi densities in response to landscape and meteorology in Saint-Georges de l’Oyapock were assessed using a cumulative link mixed model (CLMM) [49]. The CLMM was estimated using maximum likelihood by using Laplace approximation with an ordinal response. Weekly An. darlingi densities were allocated to three classes following the tercile method: “Low densities” (first tercile), “Medium densities” (second tercile), and “High densities” (third tercile). The sampling scheme used potentially resulted in grouping structures as each trap and week providing several observations (replicates). In other words, An. darlingi densities observed in the same trap or in the same week were more likely to be similar than observations from different traps or weeks. CLMM allows multiple random effects to account for grouping variables and handling replications [50]. The ordinal response of An. darlingi density was analyzed using CLMM with logit link, two crossed random effects (trap position and catching week) and C = 3 categories (“Low,” “Medium,” and “High” densities): where k = 1, 2, …, nk is the week index and j = 1, 2, …, nj is the trap, whereas γcjk is the cumulative probability up to cth category for observation in trap j at week k. The covariate vector xjkincludes the observation characteristic related to the regression coefficient (β), whereas the terms uj and uk are the random effects representing unobserved factors at the trap and week levels, respectively. The parameter αc is the threshold (cutpoint). The probability of the category (πc) is obtained by the difference as follows:

Model selection

Three models were constructed. The spatial and temporal dynamics of An. darlingi densities were better and more easily evaluated by separately investigating the relative effects of the landscape and meteorological variables in the first and second models, respectively. Targeted landscape and meteorological variables were then merged into a third dynamic spatial-temporal model. Univariate CLMMs of An. darlingi densities were fitted using the 22 landscape and 2,214 meteorological features as explanatory variables. From the univariate analysis, landscape variables with p-values below 0.20 were retained for multivariate analyses. Because of the numerous meteorological variables, a different strategy was used: only variables listed in both the hundred lowest p-values and hundred highest log-likelihoods were retained for multivariate analyses. For collinear explanatory variables, we selected those that maximized the log-likelihood and were the most easily interpretable from a bio-ecological point of view. All possible multivariate combinations with the remaining variables were tested separately for landscape and meteorological features. Selection of the best spatial and temporal An. darlingi density models was based on statistical indicators, including minimization of the Akaike information criterion (AIC) [51], minimization of the random effects variance (RE), and maximization of the area under the curve (AUC) computed from a receiver operating characteristic (ROC) analysis [52]. Once defined separately, the landscape (spatial) and meteorological (temporal) variables were used together to fit a single spatial-temporal model. Weekly maps of An. darlingi densities in Saint-Georges de l’Oyapock were produced by calculating the probabilities of each pixel in the study area belonging to “Low,” “Medium,” or “High” classes according to the landscape and meteorological conditions for each of the 32 weeks of the study. For a given pixel, the class with the highest probability was assigned.

Model validation

Model validity was assessed using a 10-fold cross-validation [53], which measures the stability of a model. The dataset was randomly partitioned into 10 parts, and the model was refitted 10 times with the partitions temporarily removed in rotation. Quality was evaluated by computing a cross-validated AUC (CVAUC), which was calculated as the average of the 10 AUC of the refitted models.

Results

Anopheles darlingi collections

Observed densities of An. darlingi in Saint-Georges de l’Oyapock showed high spatial and temporal heterogeneity (Table 4). The highest weekly mean densities were found in the peripheral Amerindian villages of Blondin and, to a lesser extent, in Martin, with 260.5 and 57.7 specimens, respectively. In the traps of Savane, Onozo, Adimo, Esperance, Maripa, and Bernet, near the city center, the weekly mean densities were lower, below 15 specimens. Temporal dynamics were characterized by maximal weekly densities observed in September (43.7) or October (69.0), depending on the study site. Minimal weekly densities were recorded at the end of the dry season in November (1.7). The spatial-temporal variation in An. darlingi densities is shown in detail in S1 Fig. The distribution of densities was characterized by a high number of null observations (Table 4). For example, the upper limit of the first tercile was 0 for the five traps in Adimo, Esperance, Maripa, Onozo, and Savane. Therefore, further modeling procedures used observations classified into three densities following the tercile method: “Low densities” ([0.0]), “Medium densities” ([0.0–12.2]), and “High densities” ([12.2–1,386.0]).

Univariate analyses

Univariate predictive CLMMs of An. darlingi densities were fitted using each of the 22 landscape and 2,214 meteorological features as explanatory variables. Seven noncollinear landscape variables showed p-values lower than 0.20, and seven noncollinear meteorological variables were listed in both the hundred lowest p-values and the hundred highest log-likelihood values (Table 5). Following the model selection procedures, these variables were retained for multivariate analyses. For the landscape features, forest-related variables were positively associated, whereas all the other variables were negatively associated with An. darlingi densities. The selected meteorological variables included four related to temperature (T) and one each related to evapotranspiration (ETP), humidity (H), and rainfall (Rain). With the exception of “TN_MaxNbConsecutiveDays_56–0_<22.5” all the other variables were negatively associated with An. darlingi densities.

thumbnail
Table 5. Landscape and meteorological variables selected for multivariate CLMMs of An. darlingi densities during the malaria transmission period (i.e., September–November dry season) in Saint-Georges de l’Oyapock, French Guiana.

https://doi.org/10.1371/journal.pone.0164685.t005

Multivariate analyses

All possible multivariate combinations, separately for the landscape and meteorological variables listed in Table 5, were tested. Statistical performances of the significant predictive multivariate landscape- and meteorology-based models are listed in Table 6. It resulted in two landscape-based models and 14 meteorology-based models (only the five best models, i.e., with the highest explanatory value, are reported).

thumbnail
Table 6. Statistical performances of the best predictive multivariate landscape- and meteorology-based cumulative link mixed models of An. darlingi densities during the malaria transmission period (i.e., the September–November dry season) in Saint-Georges de l’Oyapock, French Guiana.

https://doi.org/10.1371/journal.pone.0164685.t006

Selection of the final multivariate landscape- and meteorology-based CLMMs is detailed in S1 Text. For landscape, the LAND_2 model was selected. In this model, an increase in the built surface and dense forest surface in a 200-m radius around a trap was associated with a decrease and increase in An. darlingi densities, respectively. For meteorology, the METEO_1 model was selected showing the following: the maximal value of ETP during the 28 days preceding the collection (“ETP_max_28–0”) was negatively associated with An. darlingi densities; the maximum number of consecutive days without rainfall during the 49 days preceding the collection (“MaxNbConsecutiveDaysNoRain_49–0”) was negatively associated with An. darlingi densities; the maximum number of consecutive days with a minimal temperature below 22.5°C (25th percentile) during the 56 days preceding the collection (“TN_MaxNbConsecutiveDays_56–0_<22.5”) was positively associated with An. darlingi densities; and the maximum number of consecutive days with a maximal temperature above 33.2°C (75th percentile) during the ninth week preceding the collection (“TX_MaxNbConsecutiveDays_63–57_>33.2”) was negatively associated with An. darlingi densities.

Spatial-temporal model

The targeted landscape and meteorological variables were used together to fit a single final spatial-temporal predictive model (Table 7). The AIC of the model was 287.05, and the AUC of the “Low,” “Medium,” and “High” An. darlingi density classes were 0.78, 0.64, and 0.80, respectively.

thumbnail
Table 7. Parameters of the predictive spatial-temporal cumulative link mixed model of An. darlingi densities during the malaria transmission period (i.e., the September–November dry season) in Saint-Georges de l’Oyapock, French Guiana.

https://doi.org/10.1371/journal.pone.0164685.t007

Mapping An. darlingi densities

Extrapolation of the landscape-based model produced a static map of An. darlingi densities during the dry season (September–November) for the entire municipality of Saint-Georges de l’Oyapock (Fig 3), depicting the background influence of the landscape on An. darlingi densities, irrespective of intra-seasonal variation. In agreement with the model parameters, the densely forested surroundings of Saint-Georges de l’Oyapock were predicted as a high An. darlingi density area. This area encompassed the Amerindian villages of Blondin and Martin and included locations closer to the city center, such as the northern parts of the Onozo and Savane neighborhoods and a major part of the Bernet neighborhood, which is a military camp. The savannah plain in the northeast part of the city, mainly characterized by low vegetation (Fig 2), as well as the aerodrome neighborhood dominated by herbaceous vegetation, was predicted as medium An. darlingi density area. The most urbanized areas, i.e., city center neighborhoods, were predicted as low An. darlingi density areas.

thumbnail
Fig 3. Landscape-based model map of predicted An. darlingi densities during the malaria transmission period (i.e., the September–November dry season) in Saint-Georges de l’Oyapock, French Guiana.

https://doi.org/10.1371/journal.pone.0164685.g003

Extrapolation of the spatial-temporal model (Table 7), based on both the land cover map of the study area and meteorological observations, allowed to dynamically map An. darlingi densities at weekly time intervals (S1 Movie). For the three sampling periods of 2012, 2013, and 2014, maps showed that the model predicted a global dynamic characterized by high An. darlingi densities in September across almost the entire study area, except most of the urbanized neighborhoods. Densities of An. darlingi then progressively decreased until the end of the dry season and maintained a spatial gradient from the city center neighborhoods (lowest densities) to the densely forested areas (highest densities).

Model validation

AUC values were previously computed using the entire data set to guide model selection. The quality of the model was assessed by calculating a CVAUC from a 10-fold cross validation. The CVAUC means (standard deviations) were 0.78 (< 0.01), 0.64 (0.01), and 0.79 (0.02) for the “Low,” “Medium,” and “High” An. darlingi density classes, respectively. These values were very close to those computed using the entire data set (0.78, 0.64, and 0.80) and had low standard deviations. These results confirmed the stability of the model and its predictive value for the “Low” and “High” An. darlingi density classes.

Discussion

Although An. darlingi is one of the most important malaria vectors in the Americas, in many areas, the influence of ecological factors on local population dynamics is not clearly understood. To our knowledge, this is the first study to implement a dynamic model of An. darlingi densities in French Guiana to assist the development of vector control strategies. An original and powerful modeling method was proposed based on a CLMM that exploits the ordered nature of observations and offers a flexible regression framework. The model selection procedures were guided by objective statistical criteria while retaining the central role of entomological expertise. The final model was based on two landscape variables (dense forest surface and built surface) and four meteorological variables (absence of rainfall, maximal ETP, and minimal and maximal temperatures). The accuracy of predicting the spatial-temporal densities of An. darlingi was evaluated by using a CVAUC, which was 0.74 in average for the three density classes. The average performance was reduced by the “Medium” density class, which had an AUC of 0.64; therefore, the predictions for intermediate density areas remain uncertain. More importantly, the best AUC value (0.80) was observed for the “High” density class. In other words, the model was the most accurate for the areas needing the highest malaria vector control.

Dense forest surface was an efficient predictor of high An. darlingi densities. However, the classification of the unpopulated dense forest surrounding Saint-Georges de l’Oyapock as a high An. darlingi density area should be discussed since this species is known for its anthropophilic behavior [36, 37, 54, 55]. Several specific characteristics of the study, combined with the ecological preferences of An. darlingi, might explain this model prediction. First, the hostility of the environment, i.e., the equatorial primary forest, prevented the placing of traps far from the edges of the dense forest; hence, all the sampling sites were effectively located close to human settlements. Thus, model predictions depicted that villages encompassing large patches of dense forest at their borders are highly exposed to malaria vectors. Second, the extrapolation of the model to the entire study area was not random. Indeed, the study was focused on a circumscribed area so that any pixel of the predictive map was not further than 2 km from human settlements, i.e., from An. darlingi blood meal sources. Although little is known about the flight range of An. darlingi, densely forested places in the study area could very less likely be the hotspots for An. darlingi. Forest ecosystems with human settlements at their borders are reported to provide optimal breeding, resting, and feeding sites that favor a high density and survival of human malaria vectors [9, 56].

Built surfaces were negatively associated with An. darlingi densities. Densely built-up areas favor runoff and are not likely to offer natural water bodies suitable for An. darlingi breeding or resting [26]. In addition, the increasing density of human settlement from the periphery of the municipality to the city center constitutes a behavioral barrier for mosquitoes. Females do not need to fly far into the city area to secure a blood meal. Similar observations and hypotheses have been made for An. gambiae in Africa with the number of specimens being inversely proportional to the level of urbanization in the area [57, 58].

A main limitation of this study is the absence of landscape variables directly related to the water bodies that are essential for breeding mosquitoes. This is an inherent flaw in the Amazonian region, which prevents the effective use of optical satellite imagery to exhaustively detect water bodies under the thick forest cover. The dense forest surface variable used to model An. darlingi densities might be considered as proxy for the presence of water. Indeed, the majority of the hydrographic network could likely be hidden by the canopy. Thus, although this research supports the use of optical satellite imagery to predict An. darlingi densities, it is limited in comparison with similar studies conducted in open landscapes, such as in Sahelian Africa [57] or Europe [59]. Further studies that evaluate the potential of radar imagery to detect water bodies under the canopy might help rectify this limitation.

Another limitation implied by the use of optic satellite imagery in the equatorial area was the scarcity of cloud-free image. Thus, the landscape of the study area was characterized based on a single time point land cover map. From the very regular presence of the authors in Saint-Georges de l’Oyapock for field collections, environmental changes during the entire study period were considered negligible. Despite this, fine-scale landscape modifications might have occurred between September 2012 and November 2014. Further, the landscape changes within a season (September–November) were not objectively evaluated. These two issues might result in incorrect predictions of local An. darlingi densities. However, since the land cover classes targeted in the present study refer to coarse environmental facies, the potential impacts of environmental modifications on the results were considered as minor.

The maximal value of ETP during the 28 days preceding the collection period was negatively associated with An. darlingi densities. Mosquitoes exposed to high ETP are potentially subjected to desiccation. Further, ETP might be considered a proxy for the loss of breeding sites. The maximum number of consecutive days without rainfall during the 49 days preceding the collection period was also negatively associated with An. darlingi densities. Long, dry spells during the proliferation season of An. darlingi (i.e., the dry season) most likely caused the reduction of aquatic breeding sites and interrupted population dynamics, leading to low densities in the following weeks. The maximum number of consecutive days with a minimal temperature below 22.5°C and a maximal temperature above 33.2°C were positively and negatively associated with An. darlingi densities, respectively. With regard to ETP, very high temperatures might cause mosquito desiccation and accelerate the loss of aquatic breeding sites. In contrast, lower temperatures appeared to benefit mosquito breeding. In the study area, minimal temperatures never fell below 20°C and were always suitable for An. darlingi. In addition, during the dry season, the lowest daily temperatures were most often recorded on rainy days. Therefore, the relationship between the maximum number of consecutive days with a minimal temperature below 22.5°C and An. darlingi densities might indicate the impact of rainy episodes during the dry season on mosquito density. Meteorological variables selected in the final model were related to the weather conditions observed several weeks before the sampling sessions (from four to nine weeks), which might exceed the lifespan of An. darlingi females. Therefore, the impact of meteorological factors on An. darlingi should be interpreted over several generations, i.e., adverse weather conditions during a given generation might impact subsequent ones.

By using entomological data sampled both in the dry and wet seasons, Adde et al. (unpublished data) observed a negative correlation between monthly rainfall and An. darlingi densities. Further, other studies [60, 61] have suggested that this negative relationship highlights the local seasonal framework conducive to high An. darlingi densities: the driest months of September–November. However, when focusing on the dry season and downsizing to a weekly scale, the relationship appears to be reversed: rainy events during the dry season are beneficial for the proliferation of An. darlingi. One hypothesis is that, during the wet season (April–June), the breeding sites are flooded and are regularly flushed, resulting in a drastic reduction in mosquito densities [36]. Dry conditions during the offseason (July–August) and at the onset of the dry season (September–mid October) might decrease the level of water bodies suitable for breeding, and flushing events become rarer, resulting in high densities of mosquitoes. Drying of the water bodies continues through the following weeks, and thus the availability of suitable breeding places becomes limited again toward the end of the dry season (end-October–November), resulting in a substantial reduction in An. darlingi densities.

Significant technical, human, and financial resources were required to achieve the considerable entomological collection in this study. The collection period was restricted to extend spatial sampling, which was not possible to continue throughout the year. This study focused on the middle of the dry season (September–November), which is the major period of malaria transmission and the season when An. darlingi is abundant. Unfortunately, this collection period did not allow the measurement of the increase in mosquito densities after the wet season: it begins in September when the densities are almost maximal. Given the hypothetical differentiated impact of rainfall on An. darlingi densities across seasons outlined above, the four different models (Dry season; Offseason 1; Wet season; Offseason 2) would be necessary to capture the dynamics over the entire year.

Mapping vector-borne disease determinants and entomological hazards is highlighted by the WHO as a central feature for efficient and integrated vector control management [62]. Hazard maps should be considered as predictive support that provide relevant information about where and when vector control interventions could be focused. This information is essential for the success of malaria elimination programs [6367]. This study was conducted at a local level to provide operational results consistent with the intended scale of spatial intervention. Inevitably, this limits the generality and reproducibility of the results in other contexts. However, considering the high heterogeneity of malaria vector dynamics, focusing on local specificities is necessary to devise sustainable and effective control tools, which are not possible to implement from synoptic or large-scale approaches.

Conclusions

The analysis of the relationships between landscape parameters derived from SPOT-5 satellite imagery, meteorological observations, and An. darlingi densities allowed developing a robust and operational methodology to dynamically map malaria vector density in Saint-Georges de l’Oyapock, an area of French Guiana where residual malaria transmission still exists. Such a fine-scale modeling approach fills the knowledge gap in local An. darlingi bio-ecology and provides a tool for health authorities to establish effective vector control. Production of near real-time An. darlingi density maps can be implemented by automating satellite and meteorological data processing to aid public health authorities in planning control strategies and preparing social communication in a cost-effective, targeted, and timely manner.

Supporting Information

S1 Fig. Anopheles darlingi weekly densities. Weekly densities of An. darlingi in the eight study sites according to their respective sampling periods.

Black crosses represent the non-sampled weeks owing to the technical failure of the traps.

https://doi.org/10.1371/journal.pone.0164685.s001

(TIF)

S1 Text. Selection of the multivariate landscape and meteorological models of Anopheles darlingi densities during the malaria transmission period (i.e., the September–November dry season) in Saint-Georges de l’Oyapock, French Guiana.

https://doi.org/10.1371/journal.pone.0164685.s002

(DOCX)

S1 Movie. Dynamic mapping of the predicted Anopheles darlingi densities during the malaria transmission period (i.e., the September–November dry season) of 2012, 2013, and 2014 in Saint-Georges de l’Oyapock, French Guiana.

https://doi.org/10.1371/journal.pone.0164685.s003

(WMV)

Acknowledgments

The authors thank the technical team of the Medical Entomology Unit of the Pasteur Institute of French Guiana—Romuald Carinci, Jean Issaly, Pascal Gaborit, and Christophe Nguyen—for their help in field collections and Anopheles morphological identification. We are grateful to Cécile Vignolles (CNES, Toulouse, France), Vanessa Machault (OMP, Toulouse, France), Yves Tourre (OMP, Toulouse, France), and Jean-Pierre Lacaux (OMP, Toulouse, France) for their support in this study. We also thank the DAAC association and the Camp Bernet staff of Saint-Georges de l’Oyapock for their assistance.

Author Contributions

  1. Conceptualization: SB AA.
  2. Data curation: SB AA.
  3. Formal analysis: AA.
  4. Funding acquisition: SB.
  5. Investigation: SB ID RG AA.
  6. Methodology: SB ER MM ND ID RG AA.
  7. Project administration: SB.
  8. Resources: SB RG.
  9. Software: SB ER MM AA.
  10. Supervision: SB.
  11. Validation: SB ER MM ND MN ID RG AA.
  12. Visualization: SB AA.
  13. Writing – original draft: AA.
  14. Writing – review & editing: SB ER MM ND MN ID RG AA.

References

  1. 1. WHO. World Malaria Report 2015. Geneva: World Health Organization. 2015. Available: http://www.who.int/malaria/publications/world-malaria-report-2015/report/en/.
  2. 2. Bousema T, Griffin JT, Sauerwein RW, Smith DL, Churcher TS, Takken W, et al. Hitting hotspots: spatial targeting of malaria for control and elimination. PLoS Med. 2012;9:e1001165. pmid:22303287
  3. 3. Bousema T, Drakeley C, Gesase S, Hashim R, Magesa S, Mosha F, et al. Identification of hot spots of malaria transmission for targeted malaria control. J Infect Dis. 2010;201:1764–1774. pmid:20415536
  4. 4. Tanner M, Takken W, Lyimo E, Charlwood J. The level of anti-sporozoite antibodies in a highly endemic malaria area and its relationship with exposure to mosquitoes. Trans R Soc Trop Med Hyg. 1992;86:499–504. pmid:1362013
  5. 5. Carter R, Mendis KN, Roberts D. Spatial targeting of interventions against malaria. Bull World Health Organ. 2000;78:1401–1411. pmid:11196487
  6. 6. Machault V, Vignolles C, Borchi F, Vounatsou P, Pages F, Briolant S, et al. The use of remotely sensed environmental data in the study of malaria. Geospat Health. 2011;5:151–168. pmid:21590665
  7. 7. Rogers DJ, Randolph SE, Snow RW, Hay SI. Satellite imagery in the study and forecast of malaria. Nature. 2002;415:710–715. pmid:11832960
  8. 8. Ceccato P, Connor SJ, Jeanne I, Thomson MC. Application of geographical information systems and remote sensing technologies for assessing and monitoring malaria risk. Parassitologia. 2005;47:81–96. pmid:16044677
  9. 9. Stefani A, Dusfour I, Corrêa AP, Cruz MC, Dessay N, Galardo AK, et al. Land cover, land use and malaria in the Amazon: a systematic literature review of studies using remotely sensed data. Malar J. 2013;12:192. pmid:23758827
  10. 10. Kraemer MU, Hay SI, Pigott DM, Smith DL, Wint GR, Golding N. Progress and challenges in infectious disease cartography. Trends Parasitol. 2016;32:19–29. pmid:26604163
  11. 11. Mouchet J, Carnevale P, Coosemans M, Julvez J, Manguin S, Richard-Lenoble D, et al. Biodiversité du paludisme dans le monde. Montrouge: John Libbey Eurotext. 2004.
  12. 12. Sinka ME, Bangs MJ, Manguin S, Rubio-Palis Y, Chareonviriyaphap T, Coetzee M, et al. A global map of dominant malaria vectors. Parasit Vectors. 2012;5:69. pmid:22475528
  13. 13. Kiszewski A, Mellinger A, Spielman A, Malaney P, Sachs SE, Sachs J. A global index representing the stability of malaria transmission. Am J Trop Med Hyg. 2004;70:486–498. pmid:15155980
  14. 14. May JM. Map of the world distribution of malaria vectors. Geogr Rev. 1951;41:638–639.
  15. 15. White G. Malaria. Geographical distribution of arthropod-borne diseases and their principal vectors. Geneva: World Health Organization, Division of Vector Biology and Control. 1989. pp. 7–22.
  16. 16. Hay SI, Sinka ME, Okara RM, Kabaria CW, Mbithi PM, Tago CC, et al. Developing global maps of the dominant Anopheles vectors of human malaria. PLoS Med. 2010;7:e1000209. pmid:20161718
  17. 17. Mbogo CM, Mwangangi JM, Nzovu J, Gu W, Yan G, Gunter JT, et al. Spatial and temporal heterogeneity of Anopheles mosquitoes and Plasmodium falciparum transmission along the Kenyan coast. Am J Trop Med Hyg. 2003;68:734–742. pmid:12887036
  18. 18. Trung HD, Bortel WV, Sochantha T, Keokenchanh K, Briët OJ, Coosemans M. Behavioural heterogeneity of Anopheles species in ecologically different localities in Southeast Asia: a challenge for vector control. Trop Med Int Health. 2005;10:251–262. pmid:15730510
  19. 19. Lounibos LP, Conn JE. Malaria vector heterogeneity in South America. Am Entomol. 2000;46:238–249.
  20. 20. Alimi TO, Fuller DO, Quinones ML, Xue R-D, Herrera SV, Arevalo-Herrera M, et al. Prospects and recommendations for risk mapping to improve strategies for effective malaria vector control interventions in Latin America. Malar J. 2015;14:519. pmid:26694047
  21. 21. Musset L, Pelleau S, Girod R, Ardillon V, Carvalho L, Dusfour I, et al. Malaria on the Guiana Shield: a review of the situation in French Guiana. Mem Inst Oswaldo Cruz. 2014;109:525–533. pmid:25184998
  22. 22. Petit-Sinturel M, Carvalho L, Andrieu A, Prince C, Abboud P, Djossou F et al. Situation du paludisme en Guyane française en 2015. Bulletin de veille sanitaire Antilles Guyane. 2016;2:6–10.
  23. 23. Pommier de Santi V, Dia A, Adde A, Hyvert G, Galant J, Mazevet M, et al. Malaria in French Guiana linked to illegal gold mining. Emerg Infect Dis. 2016;22:344–346. pmid:26811971
  24. 24. Pommier de Santi V, Girod R, Mura M, Dia A, Briolant S, Djossou F, et al. Epidemiological and entomological studies of a malaria outbreak among French armed forces deployed at illegal gold mining sites reveal new aspects of the disease’s transmission in French Guiana. Malar J. 2016;15:35. pmid:26801629
  25. 25. Pommier de Santi V, Djossou F, Barthes N, Bogreau H, Hyvert G, Nguyen C, et al. Malaria hyperendemicity and risk for artemisinin resistance among illegal gold miners, French Guiana. Emerg Infect Dis. 2016;22:903–906. pmid:27089004
  26. 26. Hiwat H, Bretas G. Ecology of Anopheles darlingi Root with respect to vector importance: a review. Parasit Vectors. 2011;4:177. pmid:21923902
  27. 27. Sinka ME, Rubio-Palis Y, Manguin S, Patil AP, Temperley WH, Gething PW, et al. The dominant Anopheles vectors of human malaria in the Americas: occurrence data, distribution maps and bionomic précis. Parasit Vectors. 2010;3:72. pmid:20712879
  28. 28. Laporta GZ, Linton Y-M, Wilkerson RC, Bergo ES, Nagaki SS, Sant’Ana DC, et al. Malaria vectors in South America: current and future scenarios. Parasit Vectors. 2015;8:426. pmid:26283539
  29. 29. Floch H. Antimalarial campaign in French Guiana. I. Anophelism. Riv Malariol. 1955;34:57–65. pmid:13255704
  30. 30. Pajot F-X, Le Pont F, Molez J-F, Degallier N. Agressivité d'Anopheles (Nyssorhynchus) darlingi Root, 1926 (Diptera, Culicidae) en Guyane Française. Cah ORSTOM Ser Ent Med Parasitol. 1977;15:15–22.
  31. 31. Floch H, Abonnenc E. Anophelines in French Guiana. Arch Inst Pasteur Guyane. 1951;236:1–92.
  32. 32. Floch H, Abonnenc E. Sur le role de A. darlingi Root 1926 dans la transmission du paludisme en Guyane Francaise. Arch Inst Pasteur Guyane. 1943;7:1–10.
  33. 33. Fouque F, Gaborit P, Carinci R, Issaly J, Girod R. Annual variations in the number of malaria cases related to two different patterns of Anopheles darlingi transmission potential in the Maroni area of French Guiana. Malar J. 2010;9:80. pmid:20307300
  34. 34. Girod R, Roux E, Berger F, Stefani A, Gaborit P, Carinci R, et al. Unravelling the relationships between Anopheles darlingi (Diptera: Culicidae) densities, environmental factors and malaria incidence: understanding the variable patterns of malarial transmission in French Guiana (South America). Ann Trop Med Parasitol. 2011;105:107–122. pmid:21396247
  35. 35. Claustre J, Venturin C, Nadiré M, Fauran P. Vecteurs de paludisme en Guyane Française: étude dans un foyer épidémique proche de Cayenne (1989–1998). Bull Soc Path Exot. 2001;94:353–357.
  36. 36. Hiwat H, Issaly J, Gaborit P, Somai A, Samjhawan A, Sardjoe P, et al. Behavioral heterogeneity of Anopheles darlingi (Diptera: Culicidae) and malaria transmission dynamics along the Maroni River, Suriname, French Guiana. Trans R Soc Trop Med Hyg. 2010;104:207–213. pmid:19732925
  37. 37. Girod R, Gaborit P, Carinci R, Issaly J, Fouque F. Anopheles darlingi bionomics and transmission of Plasmodium falciparum, Plasmodium vivax and Plasmodium malariae in Amerindian villages of the Upper-Maroni Amazonian forest, French Guiana. Mem Inst Oswaldo Cruz. 2008;103:702–710. pmid:19057822
  38. 38. Vezenegho SB, Carinci R, Gaborit P, Issaly J, Dusfour I, Briolant S, et al. Anopheles darlingi (Diptera: Culicidae) dynamics in relation to meteorological data in a cattle farm located in the coastal region of French Guiana: advantage of Mosquito Magnet trap. Environ Entomol. 2015;44:454–462. pmid:26313950
  39. 39. Vittor AY, Pan W, Gilman RH, Tielsch J, Glass G, Shields T, et al. Linking deforestation to malaria in the Amazon: characterization of the breeding habitat of the principal malaria vector, Anopheles darlingi. Am J Trop Med Hyg. 2009;81:5–12. pmid:19556558
  40. 40. Tadei WP, Thatcher BD, Santos JM, Scarpassa VM, Rodrigues IB, Rafael MS. Ecologic observations on anopheline vectors of malaria in the Brazilian Amazon. Am J Trop Med Hyg. 1998;59:325–335. pmid:9715956
  41. 41. Zeilhofer P, dos Santos ES, Ribeiro AL, Miyazaki RD, dos Santos MA. Habitat suitability mapping of Anopheles darlingi in the surroundings of the Manso hydropower plant reservoir, Mato Grosso, Central Brazil. Int J Health Geogr. 2007;6:7. pmid:17343728
  42. 42. Achee NL, Grieco JP, Masuoka P, Andre RG, Roberts DR, Thomas J, et al. Use of remote sensing and geographic information systems to predict locations of Anopheles darlingi-positive breeding sites within the Sibun River in Belize, Central America. J Med Entomol. 2006;43:382–392. pmid:16619625
  43. 43. Dusfour I, Issaly J, Carinci R, Gaborit P, Girod R. Incrimination of Anopheles (Anopheles) intermedius Peryassú, An.(Nyssorhynchus) nuneztovari Gabaldón, An. (Nys.) oswaldoi Peryassú as natural vectors of Plasmodium falciparum in French Guiana. Mem Inst Oswaldo Cruz. 2012;107:429–432. pmid:22510842
  44. 44. Vezenegho SB, Adde A, Gaborit P, Carinci R, Issaly J, Pommier de Santi V, et al. Mosquito magnet® liberty plus trap baited with octenol confirmed best candidate for Anopheles surveillance and proved promising in predicting risk of malaria transmission in French Guiana. Malar J. 2014;13:384. pmid:25260354
  45. 45. Forattini OP. Entomologia médica. São Paulo: Univ. de São Paulo, Faculdade de Higiene e Saúde Pública, Departamento de Parasitololgia. 1962.
  46. 46. Faran ME, Linthicum KJ. A handbook of the Amazonian species of Anopheles (Nyssorhynchus) (Diptera: Culcidae). Mosq Syst. 1981;13:1–81.
  47. 47. McGarigal K, Cushman SA, Ene E. FRAGSTATS v4: spatial pattern analysis program for categorical and continuous maps. Computer software program produced by the authors at the University of Massachusetts, Amherst. Available: http://www.umass.edu/landeco/research/fragstats/fragstats.html
  48. 48. Monteith JL. Evaporation and environment. In: The state and movement of water in living organisms. Cambridge, U. K: B. D. Fogg Cambridge Univ. Press. 1964. pp. 205–234.
  49. 49. Christensen R. Package ordinal—Regression Models for Ordinal Data.R package version 2015. 2015:1–21. Available: https://cran.r-project.org/web/packages/ordinal/ordinal.pdf.
  50. 50. Agresti A, Natarajan R. Modeling clustered ordered categorical data: A survey. Int Stat Rev. 2001;69:345–371.
  51. 51. Akaike H. Information theory and an extension of the maximum likelihood principle. Second International Symposium on Information Theory. pp. 267–281. Budapest: Akadémiai Kiado. 1973.
  52. 52. Metz CE. Basic principles of ROC analysis. Semin Nucl Med. 1978;8:283–298. pmid:112681
  53. 53. Arlot S, Celisse A. A survey of cross-validation procedures for model selection. Statist Surv. 2010;4:40–79.
  54. 54. Zimmerman RH, Galardo AK, Lounibos LP, Arruda M, Wirtz R. Bloodmeal hosts of Anopheles species (Diptera: Culicidae) in a malaria-endemic area of the Brazilian Amazon. J Med Entomol. 2006;43:947–956. pmid:17017232
  55. 55. Charlwood JD, Alecrim WA. Capture-recapture studies with the South American malaria vector Anopheles darlingi, Root. Ann Trop Med Parasitol. 1989;83:569–576. pmid:2619371
  56. 56. Kar NP, Kumar A, Singh OP, Carlton JM, Nanda N. A review of malaria transmission dynamics in forest ecosystems. Parasit Vectors. 2014;7:265. pmid:24912923
  57. 57. Machault V, Vignolles C, Pagès F, Gadiaga L, Tourre YM, Gaye A, et al. Risk mapping of Anopheles gambiae s.l. densities using remotely-sensed environmental and meteorological data in an urban area: Dakar, Senegal. PLoS ONE. 2012;7:e50674. pmid:23226351
  58. 58. Robert V, Gazin P, Benasseni R, Carnevale P. Le paludisme urbain à Bobo-Dioulasso (Burkina Faso). In: Urbanisation et santé dans le Tiers Monde: transition épidémiologique, changement social et soins de santé primaires. Paris: ORSTROM. 1989. pp. 181–185.
  59. 59. Tran A, Poncon N, Toty C, Linard C, Guis H, Ferre JB, et al. Using remote sensing to map larval and adult populations of Anopheles hyrcanus (Diptera: Culicidae) a potential malaria vector in Southern France. Int J Health Geogr. 2008;7:9. pmid:18302749
  60. 60. Forattini OP. Exophilic behavior of Anopheles darlingi Root in a southern region of Brazil. Rev Saude Publica. 1987;21:291–304. pmid:3445112
  61. 61. Hudson JE. Anopheles darlingi Root (Diptera: Culicidae) in the Suriname rain forest. Bull Entomol Res. 1984;74:129–142. http://dx.doi.org/10.1017/S0007485300010002.
  62. 62. WHO. Handbook for integrated vector management. Geneva: World Health Organization. 2012. Available: http://apps.who.int/iris/bitstream/10665/44768/1/9789241502801_eng.pdf.
  63. 63. Takken W, Knols BG. Malaria vector control: current and future strategies. Trends Parasitol. 2009;25:101–104. pmid:19168392
  64. 64. Karunamoorthi K. Vector control: a cornerstone in the malaria elimination campaign. Clin Microbiol Infect. 2011;17:1608–1616. pmid:21996100
  65. 65. Enayati A, Hemingway J. Malaria management: past, present, and future. Annu Rev Entomol. 2010;55:569–591. pmid:19754246
  66. 66. Hemingway J. The role of vector control in stopping the transmission of malaria: threats and opportunities. Philos Trans R Soc Lond B Biol Sci. 2014;369:20130431. pmid:24821917
  67. 67. Mnzava AP, Macdonald MB, Knox TB, Temu EA, Shiff CJ. Malaria vector control at a crossroads: public health entomology and the drive to elimination. Trans R Soc Trop Med Hyg. 2014;108:550–554. pmid:25009173