Skip to main content

Rift Valley fever vector diversity and impact of meteorological and environmental factors on Culex pipiens dynamics in the Okavango Delta, Botswana



In Northern Botswana, rural communities, livestock, wildlife and large numbers of mosquitoes cohabitate around permanent waters of the Okavango Delta. As in other regions of sub-Saharan Africa, Rift Valley Fever (RVF) virus is known to circulate in that area among wild and domestic animals. However, the diversity and composition of potential RVF mosquito vectors in that area are unknown as well as the climatic and ecological drivers susceptible to affect their population dynamics.


Using net traps baited with carbon dioxide, monthly mosquito catches were implemented over four sites surrounding cattle corrals at the northwestern border of the Okavango Delta between 2011 and 2012. The collected mosquito species were identified and analysed for the presence of RVF virus by molecular methods. In addition, a mechanistic model was developed to assess the qualitative influence of meteorological and environmental factors such as temperature, rainfall and flooding levels, on the population dynamics of the most abundant species detected (Culex pipiens).


More than 25,000 mosquitoes from 32 different species were captured with an overabundance of Cx. pipiens (69,39 %), followed by Mansonia uniformis (20,67 %) and a very low detection of Aedes spp. (0.51 %). No RVF virus was detected in our mosquito pooled samples. The model fitted well the Cx. pipiens catching results (ρ = 0.94, P = 0.017). The spatial distribution of its abundance was well represented when using local rainfall and flooding measures (ρ = 1, P = 0.083). The global population dynamics were mainly influenced by temperature, but both rainfall and flooding presented a significant influence. The best and worst suitable periods for mosquito abundance were around March to May and June to October, respectively.


Our study provides the first available data on the presence of potential RVF vectors that could contribute to the maintenance and dissemination of RVF virus in the Okavango Delta. Our model allowed us to understand the dynamics of Cx. pipiens, the most abundant vector identified in this area. Potential predictions of peaks in abundance of this vector could allow the identification of the most suitable periods for disease occurrence and provide recommendations for vectorial and disease surveillance and control strategies.


Mosquitoes are a major source of nuisance throughout the world, particularly because of the abundance of pathogens transmitted by their bites and their implications for animal and human health [1, 2]. Rift Valley fever virus (RVFV), an arbovirus belonging to the genus Phlebovirus, family Bunyaviridae, is one such pathogen. Transmitted to vertebrates by mosquitoes, mainly of the genera Culex and Aedes, or direct contact with viraemic animal products, it is responsible for Rift Valley fever (RVF), an acute disease considered as a significant global threat to both humans and animals [3, 4]. The virus has caused epizootics and human epidemics throughout Africa, the Arabian Peninsula and several islands of the Indian Ocean [5]. In ruminant livestock, especially sheep and cattle, the disease is characterized by high abortion and mortality rates (100 % in neonatal animals and from 10 % to 20 % among adult animals [68]). Human infections are characterized by blurred vision, retinal lesions, headache, loss of memory, lethargy, myalgia, fever and haemorrhages, although a minority (approximately 1 %) of the patients can develop complications such as encephalitis, ocular disease, retinitis or fatal haemorrhagic fever [9]. RVF, which significantly affects the health of animals and/or humans, induces a very heavy economic impact in the societies where it is present, and particularly in developing countries of tropical and sub-tropical areas [10]. In Botswana, as in other regions of sub-Saharan Africa, RVF virus is suspected to circulate among animals without reports of clinical outbreaks [11]. Indeed, RVF is often ignored and under reported in Botswana because clinical cases are not monitored and detected. To prevent and control RVF, effective knowledge about RVF risk areas and risk periods is required. As a mosquito-borne disease, RVF distribution in space and time is closely associated to the geographical distribution and dynamics of its vectors and hosts, which are themselves affected by climatic and landscape features [1215]. Indeed, the ecology of vectors is highly dependent on environmental conditions. Each mosquito species requires specific environmental conditions to develop and survive, such as water availability to lay eggs, optimal temperature for aquatic stage development, limited wind (to facilitate host-seeking or breeding-site-seeking behaviour), or specific vegetation for some species [16].

The waters of the Okavango Delta allow the establishment of various mosquito species, such as Culex spp. breeding in the margins of permanent water bodies or Aedes spp. breeding in the marginal pools created by fluctuation of waters due to flooding. To date, in this region, no outbreaks of RVF have ever been reported. However, substantial RVF outbreaks have occurred in neighbouring countries (South Africa, Namibia, Zambia) [17] and serological investigations performed close to the time of this study provided evidence for the first time of RVFV antibody circulation among wild and domestic bovids in northern Botswana [11]. Consequently, the present study was structured around two main objectives: (i) Identifying the main potential vectors of RVFV present in the Okavango Delta; and (ii) Adapting and using a climate driven abundance model [18] under different biological scenarios to appreciate the potential influence of climatic and environmental factors, such as temperature, rainfall and water availability, on Culex pipiens populations, the main RVF mosquito vector species present at the livestock-wildlife interface in the Okavango Delta. It is commonly accepted that temperature is the main driver of transition rates between eggs, pupae, nymphs, and adult stages, and of mortality in mosquitoes [19]. This work specifically assessed the importance of the water availability and the respective roles of flood and rainfall. Therefore, we tested all possible combinations of these two ecological variables: temperature and water availability.


Study area

The Okavango River basin, a unique wetland environment shared by Botswana, Angola and Namibia, supports a fragile and extremely complex ecosystem [20]. The study was carried out at the western border of the Delta located in semi-arid north-west Botswana, an interface between wetland and dry-land, precisely situated between 22°08′ and 22°18′E, and 19°07′ and 19°26′S (Fig. 1).

Fig. 1
figure 1

Map of the study area in north-western region of the Okavango Delta

The hydrological system of the Delta is very complex and dominated by annual flood from the Okavango River. With the rivers that rise in the well-watered Angolan highlands and through Namibia, ‘disappearing’ into the sands of the Kalahari Desert and the “thirsty” atmosphere above the Okavango Delta, the whole Okavango system is endoreic in nature. Most rainfall occurs as summer downpours between October and March (Fig. 2a) ranging between 1,300 mm in parts of the catchment furthest to the north-west and steadily declining to 450 mm in the lowest reaches of the Delta [21]. The dry season, between April and November, is characterized by the almost total absence of rains. Maximum temperatures get warmer from September-October, soaring up to 40 °C, and dropping considerably towards the end of November or in early December (Fig. 2b). Nevertheless, until the end of February or early March, summer temperatures stay relatively warm, ranging from 38 °C during the day to 20 °C at night. During the dry winter season (April to September), day temperatures usually stand around 25 °C and evening temperatures can be as low as 2 °C and even go below freezing. The flood events are primarily caused by a flood wave induced by the high rainfall events happening in the northern catchment which lies entirely in the Angolan part of the Okavango River basin. During these events, depending on the size of the flood, the extent of the inundated area increases from 5,000 km2 to 6,000-12,000 km2 [22]. In addition, the flooding dynamics depends on the highly seasonal rainfall to which the Delta area is exposed between October and March (Fig. 2a,c). The spatial variability in available water is presented later in this paper.

Fig. 2
figure 2

Meteorological and environmental data in the study area, Okavango Delta, 2005–2012. a Rainfall. b Mean daily temperature. c Proportion of flooded area dynamics over time

Few people live permanently in the area of the inner Delta and the remote areas of western and northern Ngamiland (district located around the Delta in northwestern Botswana). Most of the existing settlements are very small, with a high proportion of the population clustered in a handful of larger villages and in the district capital of Maun (Fig. 1). In total, the Ngamiland District population is composed of about 150,000 persons [23]. Most land is used for communal farming, private ranching, wildlife conservation and tourism. The extremely rich biodiversity of the Delta and high abundance of large mammals results in the great diversity of habitats ranging from different kinds of drylands to various wetlands and the obvious availability of a perennial supply of water. The conditions are poorly suited to agriculture and ranching because soils have low fertility, rainfall varies too much to produce good yields and pastures cannot support large numbers of livestock, all the more so because livestock diseases and crop pests limit production. Most livestock is kept on communal or tribal land, and adjacent to the buffalo veterinary fence erected to limit the spread of disease between wildlife and livestock. Cattle and goats are the most abundant stock, with smaller numbers of sheep and donkeys [20]. Cattle numbers in the western part of the Delta, where this research was implemented, are estimated to be approximately 20,000 head and their density (32 cattle heads/km2) is the highest among the wildlife/livestock interface of Northern Botswana [24].


Mosquito sampling

Sampling sites to collect mosquitoes were selected on the basis of the detection of RVFV antibodies in livestock during a preliminary serological survey undertaken in buffalo and livestock in April 2010. The sampling design is described in Jori et al. [11], see summary of methods and results in Additional file 1. This survey enabled us to identify, in the study area, cattle corrals (isolated farms) with significant seroprevalence rates for antibodies against RVFV (Table 1). In our study, those seroprevalence rates were used as a proxy of mosquito abundance. Sampling was performed by disposing the traps between 25 and 50 m from cattle corrals located in the areas named Etsha 1, Etsha 6, Kandalangondi and Gumare (Fig. 1). Despite measures of meteorological parameters on the capture sites were not taken, climate, topography and landscape in the 25 km2 area seemed quite homogeneous and we assumed that the four mosquito capture sites had very similar environmental and ecological conditions. For instance, the altitude in the four sites were very similar: Etsha 1 is located 973 m above sea level (asl), Etsha 6 is 974 m asl, Kandalangondi is 966 m asl and Gumare is 964 m asl.

Table 1 Seroprevalence of IgG antibodies detected in the five different crush pensb located on the Western boundary of the Okavango Delta

Trapping episodes took place from sunset (between 15:30 and 18:00) to sunrise (6:00 to 8:00). The traps consisted of a rectangular “tent” of fine netting spread over 4 aluminum corner poles and had dimensions of 2 m long × 1.5 m wide × 1.6 m high. The trapping method is based on the attraction of mosquitoes to carbon dioxide and RVFV hosts proximity [25]. Dry ice was used as a source of carbon dioxide to bait mosquitoes just before sunset in order to collect the mainly crepuscular/nocturnal vector species [26]. The netting was rolled up to approximately 15 cm from the ground to allow mosquito entry and left overnight. The tent traps were cleared of mosquitoes just before sunrise, by entering the trap, rolling the sides down to ground level and collecting the mosquitoes with a mechanical aspirator. Mosquitoes were transferred to small cages for transport to the base station, where they were killed by freezing for 20 min at -20 °C, sorted according to genus and stored in 1.8 ml cryotubes in lots of up to 200 per tube. The sorted mosquitoes were then frozen at -20 °C or lower for preservation and transported to the laboratory on dry ice for species identification and virological testing. The sampling extended from April 2011 to March 2012 during the summer rainfall season excluding the winter months from May to September. During the sampling period, in each location one night per month, 6 traps were set; 4 traps were disposed at the 4 cardinal points surrounding the cattle corral and the 2 remaining traps were disposed around the cardinal point on the closest side from the Delta. One trapping episode failed because of heavy rain during the night. In total, 23 trapping episodes were undertaken under good meteorological conditions (Table 2).

Table 2 Global mosquito dynamics. Trapping summary for the four sites with the date of each trapping episode, the number of tubes collected (Nb) and the number of collected Culex pipiens among all mosquitoes collected by trapping session. The date format used was 

Identification of the trapped mosquito specimens was performed according to the keys and descriptions of Jupp [27], Edwards [28], Gillies & de Meillon [29] and Gillies & Coetzee [30]. To test the presence of RVFV in vectors, a total of 18,259 mosquitoes divided in mosquito pools were processed to obtain supernatant fluid as described by Jupp et al. [31], which was inoculated in infant mice (NHLS Animal Ethics Clearance Certificate No. 124/11). Those mice were tested by RT-PCR. We concentrated on species known as vectors of RVFV [3234].


Daily temperature data were recorded by the Department of Meteorological Services of Botswana at Maun airport located at the Delta’s southern fringe, which is the closest meteorological station to the study area (Fig. 1) and has been recording daily climatic data since 1921.


Daily rainfall data were estimated from satellite Tropical Rainfall Measuring Mission (TRMM) at 0.25° spatial resolution, downloaded from NASA’s Goddard Earth Sciences Data and Information Services Center from 2005 to 2012 [35].

Mapping flooding extent

To assess the hydrological dynamics, we used inundation extent measures derived from Moderate Resolution Imaging Spectroradiometer (MODIS) satellite imagery. The daily proportion of flooded areas was computed around each trapping location (Etsha 1, Etsha 6, Gumare, Kandalangondi) within a 2 km radius and at a study area scale defined as the area around the centroid of the four sites within a 25 km buffer radius (Fig. 1). Computation was performed using Geographic Information System (GIS) functionalities following the method detailed in Additional file 2.

Modeling mosquito population dynamics

Model design and analysis

A continuous time and stage structured model designed by Cailly et al. [18] to study mosquito population dynamics and applied in European wetlands was adapted in order to simulate the population dynamics of the main RVFV vector in the north-western region of the Okavango Delta.

The model considers the entire mosquito life-cycle (Fig. 3). It uses an a priori mechanistic mathematical description of all processes of mosquito population dynamics and a deterministic representation of the average behaviour of the population. The model includes meteorological and environmental conditions, such as water presence and temperature, and specific processes of mosquito life-cycle, such as diapause.

Fig. 3
figure 3

Diagram of the generic model of mosquito population dynamics based on life-cycle inspired from Cailly et al. [18], succession of stages (not italicized black text) and events (italicized blue text). Mosquito life-cycle contains a complete metamorphosis between aquatic juvenile stages drawn on the left of the dotted line and terrestrial adult stages on the right. As in Cailly et al. [18], females were divided into nulliparous (which have never laid eggs) and parous (which have laid eggs at least once). The green dotted box indicates the females which moved to seek a host or an oviposition site. Culex species enter diapause as young adults (top-right-arrow) while Aedes species enter diapause as egg (bottom-right-arrow)

The model was run through seven consecutive years. The dynamics in year n + 1 explicitly depends on the dynamics in year n and on survival rates during the unfavourable seasons. As in Cailly et al. [18], ten different stages were considered: 3 juvenile aquatic stages (E, eggs; L, larvae; P, pupae), 1 emerging adult stage A em ), 3 adult nulliparous female stages (A 1h , A 1g , A 1o , and 3 adult parous female stages (A 2h , A 2g , A 2o ) (Fig. 3). Adults were subdivided regarding behaviours during the gonotrophic cycle, which could affect their survival (h, host-seeking; g, transition from engorged to gravid; o, oviposition site seeking). Individual transitions between stages were due to different demographic events such as egg mortality and hatching, larval mortality, pupation (moult of larvae to pupae), pupa mortality, adult emergence, adult mortality, adult engorgement, egg maturing or oviposition (Fig. 3). Density-dependent mortality is assumed at the larval and egg stages [16, 36]. The success of adult emergence is considered dependent and negatively correlated to pupal density [37]. Because of their inability to feed on blood excludes them from the pathogen transmission process, adult males were excluded from future computations.

To take on board the environmental influence and the seasonality of the mosquito life-cycle, Cailly et al. [18] based their model on two systems of ordinary differential equations (ODE): one for the favourable period, during which mosquitoes were active, and the other for the unfavourable period, during which, according to the temperature level and the length of the day, diapause occurred [19]. The diapause of Culex pipiens in southern Africa consists of a combination of winter die-off and survival of a portion of the adult females in refuges such as reed beds (natural wetland habitat widely occurring in the Okavango Delta region and characterized by a vegetation dominated by reeds). Larvae may also survive in the reed beds [3840]. For Culex and Anopheles species, for which nulliparous adults were those surviving during the diapause, the ODE systems were detailed in Additional file 3.

Adaptation of the mosquito population model to the main RVFV vector species of Okavango Delta

Once the main species present in our study area were identified, we defined the model parameters using published information and expert opinion. Systematic literature reviews were performed in order to define the values of model parameters for the most abundant species, using words related to the life-cycle of mosquitoes, the name of the species and to this specific region in southern Africa. Knowledge gaps in the literature were addressed by organizing discussions with entomologists working in the region (Botswana, South Africa) in order to discuss the missing parameter values. According to local expert opinion and literature review, we defined the favourable period beginning on the 1st September and ending on the 30th June.

The model was used with daily temperature, flooding and rainfall as inputs to highlight the main meteorological and environmental drivers of mosquito population dynamics in the study area. A daily temperature dataset collected from Maun meteorological station was the same for the four study sites, assuming a low variability in temperature among the study sites because of the very low difference in altitude between sites (maximum 10 m) and with Maun (maximum 26 m). As we assumed, differences between the rainfall measured in Maun meteorological station and local rainfall in the study area, the latter was estimated from satellite estimates. Given the proximity between the four sites and the absence of relief, rainfall was considered homogeneous on the study area. Nevertheless, spatial variations in flooding were accounted for by the use of local satellite-derived datasets - flooding was estimated independently for each site (site-scale). We tested four scenarios based on different meteorological and environmental factor influence on the population dynamics of main RVF vector in Okavango Delta: Scenario 1) Temperature alone is driving the population dynamics; Scenario 2) A combination of temperature and rainfall is sufficient to forecast the population dynamics; Scenario 3) A combination of temperature and flooding is sufficient to forecast the population dynamics; and  Scenario 4) A combination of temperature, rainfall and flooding is required to forecast the population dynamics.

Initial conditions

Each run was initiated during the unfavourable season, when mosquito numbers were supposed to be the lowest. In our study area, this date corresponds with the month of July. Therefore, the model was initialized in this month with an initial population of 107 emerging adults. The outputs of the model were the simulated daily numbers of individuals at each stage.

Model validation

To validate the model, we compared model outputs of simulated host-seeking adult abundance (A 1h + A 2h ) for the years 2011–2012 with field data. For each scenario, the degree of association between observed and simulated data at the time of trapping was assessed by calculating the Spearman correlation coefficient ρ. We used an aggregation of field data from the four sites and ran the scenarios 3 and 4 using flooding extents estimated at the study area scale. Then, in order to assess the ability of the model to predict the spatial variation of the abundance and still using the Spearman correlation coefficient ρ, we compared the simulated - from scenario 4 - and observed maximal abundance for each site independently. In order to do this, we used field data from each site independently of each other and compared them with abundances simulated independently for each site using site-scaled estimates of flooding extents (temperature and rainfall being assumed similar on the four sites).

Variance-based sensitivity analysis

We used a global sensitivity analysis, varying simultaneously all of the model parameters [41], to assess the dependence of the model to these different factors. The method tested the influence of the input variation on the normally distributed aggregated outputs. This analysis enabled us (i) to determine the key parameters involved in the biological dynamic system; and (ii) to adjust the functions together with the model.

We determined four aggregated outputs per year: (i) the maximal number of mosquito females (abundance peak); (ii) the proportion of host-seeking females per day, averaged in a time window of 20 days around the abundance peak (host-seeking rate); (iii) the parity rate, which is the proportion of parous females per day, averaged in a time window of 20 days around the abundance peak; and (iv) the date of the emergence of most mosquitoes from aquatic and immature to adult and aerial stages. The abundance peak enabled us to characterize the mosquito dynamics while the parity rate provided an estimation of the proportion of mosquitoes having already taken a blood meal, therefore, potentially infected - of course, the proportion of infected mosquitoes will usually be less than the proportion that have taken a blood meal and the proportion of infected mosquitoes that transmit will be less than the number infected.

We produced a data set where each parameter of the model varied independently by approximately 10 %, 25 % and 50 %. Simulations using this data set enabled us to estimate the contribution of each single parameter to a minimum of 90 % of the variance. Therefore, only the main factors or interactions accounting for these 90 % of the output variance were retained as key parameters for the model [42]. The contribution of the variation factors to the output variability was evaluated using a linear regression, with the aggregated output as response and the parameters as explanatory variables. Demographic simulations were performed using Scilab (version 5.4.1) [43], and statistical simulations and plotting were performed using R (version 3.1.2) [44], free and open-source software for numerical computation and graphics.


Mosquito capture results

A total of 26,289 mosquitoes from 32 different species were identified (Table 3). Our trapping records revealed a large difference in the number of trapped mosquitoes. Indeed, while an important number of mosquitoes (> 24,900) was collected in the locations of Etsha 1 and Kandalangondi (11,927 and 13,044, respectively), a much lower number of mosquitoes (<1,500) was trapped in Etsha 6 and Gumare (1,102 and 216, respectively). The three main species representing respectively 69 %, 21 % and 4 % of the total number of identified mosquitoes were Culex pipiens with 18,267 individuals, Mansonia (Mansonioides) uniformis with 5,429 specimens and Mansonia (Mansonioides) africana with 1,140 individuals. It should be noticed that a very low number of Aedes (Aedes) was detected (4 species including 134 individuals) mostly represented by Ae. (Neomelaniconion) mcintoshi and Ae. (Neomelaniconion) unidentatus. Despite all these species being able to transmit RVFV [45, 46], no evidence of RVFV or other arboviruses was found in the mosquito pools inoculated to mice.

Table 3 Summary of the identified species of mosquito

We decided to base the model on Cx. pipiens, because this species presented the strongest variation in abundance over time with 3,806 mosquitoes caught in April 2011, 4 in September 2011, 6 in November 2011, 301 in December 2011, 122 in January 2012 and 14,028 in March 2012 (Table 2) and it is also the one species in our collections most strongly associated with transmission of RVFV elsewhere [33, 47, 48].

Estimation of environmental and meteorological environment

The daily variations in rainfall, temperature and flooding on the global study area were presented in Fig. 2a, b, c, respectively. The estimated spatial variation of flooding illustrated in Fig. 4 showed more abundant water flock in Kandalangondi and Etsha 1 (dark and light blue lines) than in Gumare and Etsha 6 (red and orange lines) during the dry season.

Fig. 4
figure 4

Proportion of flooded area dynamics over time in the four trapping sites, Okavango Delta, 2005–2012

Definition of the parameters and functions of the model for Cx. pipiens

A list of 21 constant parameters related to the bio-ecology of Cx. pipiens in subtropical areas was computed from existing literature and expert opinion (Table 4). Our estimation of carrying capacity was based on the maximum density of individuals (larvae or pupae) and on the available breeding surface in a hectare (considering adult females lay their eggs only on the border of the water surface).

Table 4 Description of the different parameters used in the model

According to bibliographic review and expert opinion on the bio-ecology of Cx. pipiens, we defined the impact of meteorological and environmental variables as follows: transitions between successive stages (f X) were all temperature-driven for the aquatic stages, as well as the duration of egg maturation (transition from engorged to gravid: f Ag ) [38]. Moreover, all mortality rates were temperature-driven, except the egg mortality [37]. We integrated a relationship between adult mortality related to seeking behaviour (mr) and the environmental variables, considering that mr was inversely correlated with the water availability, indeed the farther mosquitoes need to pull away from the water to find a host, the higher is their probability of dying. We also considered that water availability impacts the environment’s carrying capacity of aquatic stages (k L and k P ), flooding increasing the number of breeding sites available. Definitions of the model functions for Cx. pipiens species were provided in Table 5.

Table 5 Functions describing Cx. pipiens life-cycle

All of these functions were modified according to the four hypotheses on the meteorological and environmental factors suspected to influence on the population dynamics of Cx. pipiens in Okavango Delta. Those functions can take into account daily temperature (T) (all scenarios) and different rates of water presence (W), including either 8 days cumulated rainfall (Scenario 2), or the flooding area proportion (Scenario 3), or the mean of both (Scenario 4).

Simulation and validation

The model was run for a seven-year period (from the beginning of the favourable period of 2005 to the one of 2012). Figure 5 shows the host-seeking adult mosquitoes abundance predicted by the four scenarios plotted with a temporal superimposition of trapping records (mean on the four sites per month), simulated abundances and climatic and environmental variations. The simulated population dynamics present a strong temporal variability (Fig. 5a) following environmental and meteorological temporal variation (Fig. 5b). According to the used scenario, the starting date of adult emergence moves back or forward in the year: October in scenario 1, January in scenario 3, February in scenarios 2 and 4; nevertheless the date of the peak of mosquitoes seemed to be subject to variation.

Fig. 5
figure 5

Confrontation between aggregated field data (red stars), simulations (coloured lines on the upper graphic) under the four scenarios (a) and environmental variations (b). Scenario 1: Temperature used as unique input to describe the Culex pipiens population dynamics. Scenario 2: Combination of temperature and rainfall used as inputs to describe the Culex pipiens population dynamics. Scenario 3: Combination of temperature and flooding used as inputs to describe the Culex pipiens population dynamics. Scenario 4: Combination of temperature, rainfall and flooding used as inputs to describe the Culex pipiens population dynamics

Scenario 1 - temperature alone has an effect on mosquito dynamics

The simulated number of host-seeking mosquitoes under the hypothesis of temperature as the only driver of mosquito dynamics (orange line in Fig. 5) increased as soon as the favourable period begun according to the increase of temperature in November. The global mosquito predicted abundance was poorly correlated with trapping collection results. The Spearman correlation coefficient (ρ) was 0.26 (P = 0.658).

In the same way, tests of comparison applied to each site separately showed low correlations: ρ = 0.06 for Etsha 1 (P = 0.913), ρ = -0.20 for Kandalangondi (P = 0.700), ρ = -0.12 for Gumare (P = 0.827), ρ = 0.15 for Etsha 6 (P = 0.805).

Scenario 2 - combined effects of temperature and rainfall have a significant influence on mosquito dynamics

The host-seeking headcount simulated using rainfall amounts as a second driver of the Culex pipiens dynamics in addition to temperature (green line in Fig. 5) increased several months after the beginning of the favourable period (April) and did not reach a peak as high as the other simulations. Simulations were significantly more correlated with the global Cx. pipiens population dynamics on the four sites. The Spearman correlation coefficient (ρ) was 0.89 (P = 0.033).

The Spearman’s tests applied to each site separately showed distinct patterns. The correlation for the sites located close to the Okavango River, Etsha 1 (ρ = 0.87, P = 0.024) and Etsha 6 (ρ = 0.97, P = 0.005), was higher than for the two further sites: Kandalangondi (ρ = 0.61, P = 0.200) and Gumare (ρ = 0.38, P = 0.461).

Scenario 3 - combined effects of temperature and flooding influence mosquito dynamics

Scenario 3 used the proportion of flooded surface estimated in the smallest area including the four sites (study area scale) as the second driver of the mosquito dynamics, in addition to temperature. The simulated number of host-seeking mosquitoes (cyan line in Fig. 5) increased a few months after the beginning of the favourable period (February-March) and quickly reached a high quantity. Simulations well fitted to the field data: the Spearman correlation coefficient was equal to 0.83 (P = 0.058).

The Spearman’s test applied to each site separately was quite close to those computed from scenario 2 simulations. The rho coefficient were very significant for the two sites with trapping locations closer to the Okavango River and more impacted by flooding: Etsha 1 (ρ = 0.93, P = 0.008) and Etsha 6 (ρ = 0.054). The correlations with southern sites: Kandalangondi (ρ = 0.46, P = 0.354) and Gumare (ρ = 0.12, P = 0.827) farther from the river were not significant.

Scenario 4 - temperature, rainfall and flooding, have a significant effect on mosquito dynamics

The number of host-seeking mosquitoes simulated according to the last scenario involving proportion of the flooded surface on the Okavango Delta, rainfall and temperature as three significant drivers of the mosquito dynamics (purple line in Fig. 5) presented an intermediate pattern between scenarios 2 and 3. Spearman’s coefficient between predicted values and global field data being the highest (ρ = 0.94, P = 0.017) testified a strong correlation.

Considering each site separately, using temperature, flooded area proportion at site-scale and rainfall as drivers of the population dynamics, the global correlation between predictions and field data produces the strongest coefficients of determination for Etsha 1: ρ = 0.87 (P = 0.024) and Etsha 6: ρ = 0.97 (P = 0.005). The correlations with trapping records from Kandalangondi: ρ = 0.67 (P = 0.148) and Gumare: ρ = 0.38 (P = 0.461) were also relatively high but still less significant than for the other sites.

Spatial dynamics

Simulations using the proportion of flooded surface estimated at the local scale (site-scale) in scenario 4 seemed able to qualitatively predict locations where mosquitoes are going to be more abundant than somewhere else. Figure 6 illustrates the correlation between the maximum number of mosquitoes collected in each of the four study sites, and the maximum number of mosquitoes as predicted by the model, using a bi-dimensionnal representation of the data. Using scenario 4 simulations, the Spearman’s correlation between maximal abundance of collected mosquitoes and the predicted number of mosquitoes over one year was extremely significant (ρ = 1, P = 0.083; Fig. 6). The model was able to accurately replicate the spatial variation of the number of mosquitoes between the different sites.

Fig. 6
figure 6

Separate confrontation between abundance of mosquitoes collected on each field mosquito capture site and the abundance expected by site-scale simulation under scenario 4 with the Spearman’s rho coefficient

Sensitivity analysis

The sensitivity analysis computed on the aggregated outputs (the abundance peak value, the host-seeking rate, the parity rate and the date of the peak of emergence of adults) of the model run under scenario 4 allowed us to identify six key parameters contributing to the aggregated output variance: (i) the adult mortality rate at the emergence; (ii) the date of the beginning of the favourable period; (iii) the development rate of emerging adults per day; (iv) the maximum adult mortality rate related to seeking behaviour; (v) the number of eggs laid by oviposition of nulliparous female; and (vi) the sex-ratio at the emergence (Fig. 7). The peak value and the host-seeking rate were affected by the same parameters. The date of the emergence was related to the beginning of the favourable period. The parity rate was mostly driven by adult mortalities related to the host seeking behaviour.

Fig. 7
figure 7

Key parameters contributing to aggregated outputs variance in scenario 4. Adult mortality rates, development rate of emerging adult and sex ratio at the emergence, number of eggs laid per nulliparous female and beginning of the favourable period significantly affect the outputs of our model


An important observation relative to the trapping results was the overabundance of Cx. pipiens (> 8,000) compared to the other species (> 60 %), and especially, compared to the small number of Aedes spp. (< 150; < 1 %) collected. The overabundance of Cx. pipiens and the very low numbers of Aedes spp. could involve different processes. A possible explanation was that the Okavango Delta, with permanent water all year round, would be unfavourable to Aedes spp. which lay its eggs on damp soil and requires an event of desiccation for the eggs to hatch. Finally, Aedes eggs hatching irregularly triggered by flooding events [4951], a single trapping session per month could have failed to capture this species of mosquitoes at the moment when these species were more abundant. The implementation of more frequent trapping sessions over the year, without a break during the dry season may enable this hypothesis to be tested. Similar observations regarding disparity between Culex spp. and Aedes spp. were made in other delta areas such as the Senegal River [52]. Another hypothesis was that the lack of field data relative to the winter months, between April and September, could have hidden a higher presence of Aedes spp. (Alan Kemp personal communication, [26, 53, 54]).

A strong spatial variability in catches was observed between the four trapping sites (Table 2), with a very high number of mosquitoes (> 24,974) collected in two sites (Etsha 1 and Kandalangondi) and a much lower number of mosquitoes (< 1,500) trapped in the other two sites (Etsha 6 and Gumare). We hypothesized that these local differences were probably due to the distance between the trapping sites and the edge of the Delta. Indeed, Kandalangondi and Etsha 1 being closer to the Delta they benefit from higher and more regular water availability (Figs. 1 and 4). It was worth noting that previous studies on the spatial distribution of Cx. pipiens in Europe showed weak correlations between land cover and Cx. pipiens abundance [55, 56], stressing the ability of this species to colonize diverse breeding sites. Our results suggest that it might not be the case in constrained environments, where water presence determines the availability of mosquito breeding sites. Indeed, predicted and observed abundance values poorly fitted under scenario 1 (temperature alone) although the fit was better under all other scenarios. Both, the presence of permanent Delta waters and seasonal regime of rainfall, were likely to influence the population dynamics of this species by affecting (i) the environmental carrying capacity of larvae and pupae, and (ii) the adult mortality rate related to seeking behaviour. The efficiency of scenario 4 (temperature, rainfall and flooding), for which the correlation coefficient between observed and predicted values was the highest, suggests a strong influence of both flooding and rainfall on Cx. pipiens dynamics in the Okavango Delta area. The impact of flooding is expected to be higher for the sites close to the Delta River and the impact of rainfall higher for the sites located further from the Okavango River. Our results confirmed these trends, showing a better fit between predicted and observed values obtained in Scenario 2 (temperature and rainfall amounts) than in Scenario 3 (temperature and flooding) for Kandalangondi and Gumare sites, both located further south from Okavango River. Moreover, taking into account local flooding seemed to explain the spatial variability in catches observed among the four trapping sites (Fig. 6) and supports the hypothesis that local water variations have an important impact on Cx. pipiens dynamics and abundance. To test this hypothesis in the field, measurements of water variations could be conducted around each study site. Additional local measures of meteorological conditions should also be conducted around each site to confirm the hypothesis that there is a low variability in temperature among the study sites.

The sensitivity analysis allowed the identification of the key parameters driving Cx. pipiens population dynamics in the Delta (Fig. 7). As for other species, an improved knowledge of those parameters through laboratory or field studies would be required to increase the precision of the model predictions [42]. Moreover, those parameters allow targeting potential control points in the biological system which can be used to discuss possible scenarios of vector control strategies [57]. Key parameters are dependent on the species but also on the environment, for instance, in Osijek, Croatia (see Loncaric et al. [57]) highlighted the pupae and adults as sensitive stages of Cx. pipiens life-cycle. In our study area, control strategies focused on adult stages are likely to have a bigger effect than those targeting aquatic stages (pupa or eggs).

The increased abundance of Cx. pipiens simulated in December was consistent with the higher seroprevalences of RVFV in human samples collected during the same month in another area of Northern Botswana observed by Tessier et al. [58]. In other words, the increase of transmission rate in humans in Nov-Dec-Jan suggested by Teissier et al. [58] could be explained by an increased Cx. pipiens abundance. The relationships between mosquito-borne disease transmission and vector dynamics have commonly been established [5962]. Nevertheless, no RVFV could be detected despite processing 65 % of the mosquito catches and using highly sensitive methods for the detection of arboviruses. Analysis concentrated on those vector species known to play a role in the epidemiology of the diseases elsewhere. This absence of virus in mosquitoes was not surprising according to the absence of active circulation of virus. Indeed, in our study area, no clinical outbreaks have ever been detected despite the reported presence of antibodies in humans, livestock and free ranging African buffalo [11, 58]. Under those circumstances, the probability of detecting RVFV in mosquitoes in the inter-epidemic period is extremely low [63]. Moreover, since the detected antibodies were only indicative of past and not recent infections, further epidemiological surveys in human and animal patients using serological tests able to detect recent infections (IgM antibodies) are necessary. In Southern Africa, RVF is more commonly controlled by vaccination of animal populations rather than through vector control. However, our results could be useful to understand patterns of Cx. pipiens dynamics and the vector control efforts should be focused on adult stages rather than aquatic stages. In the case of an emergence of a RVF outbreak, public health agencies could use the temporal and spatial predictions of mosquito abundance provided by our study based on real-time collected meteorological data to set up adapted surveillance and control plans. Agreeing to previous studies, extreme weather events might create the necessary conditions for vector borne diseases such as RVF, to expand beyond its geographical area, triggering unexpected impacts on the animal and human health of newly affected countries [6466].

The population dynamics model used in this study has been implemented so far, in three different and very diverse ecological contexts: (i) the French Camargue region to analyse Anopheles, Culex and Aedes mosquito population dynamics in a temperate wetland [18, 67]; (ii) the French coastal region of Cote d’Azur to study Aedes mosquito population dynamics in an urban area [42]; and (iii) our study, where Cx. pipiens population dynamics was described in a tropical wetland environment. This study verified versatility of the model and suggested that it could possibly be used in a wider range of applications. The model predictability under conditions of inter-annual variability had not been assessed and could be explored by setting up mosquito trapping campaigns through several years or at least each month during one complete year. To improve the model predictions and better explain the observed spatial variability in the vector abundance, information on land cover and vegetation could be accounted for in the model, provided that land cover maps are available on the study area, or derived from satellite imagery [68, 69]. Moreover, to assess the response of the model for the entire life-cycle, additional sampling efforts catching different mosquito life stages (adult/eggs/larvae/pupae) with a higher sampling frequency are required. Finally, host density (cattle, wildlife) may also be an important factor to be taken into account, because of its potential impact on vector densities.

The vector dynamics simulated by our model could be used as inputs for a compartmental (SIR) model [70], other transmission models [71, 72]. Our model could also be coupled to predictive climate models such as proposed by Guis et al. (2011) [73]. For instance, the three identified climatic factors were relevant to predict vector population increase and could provide the first step for establishing risk maps of RVF emergence or entomological risk, taking into account the seasonal variations of host and vector distributions, and allowing the establishment of assumptions about the effect of the environmental factors on RVF transmission. This kind of application could bring better knowledge but also enable testing and improving of control strategies (vaccinations or vector-control strategies [74]).


In this paper, we provide original data describing, for the first time, the diversity of mosquito populations in the region of Okavango Delta in northern Botswana where RVFV is suspected to circulate without reports of clinical outbreaks. Our work also highlights a large overabundance of Cx. pipiens compared to other mosquito species in the most populated areas of the western side of the Okavango Delta. The rare occurrence of Aedes spp. could be explained by an unfavourable environment (permanent water) avoiding the desiccation of Aedes spp. eggs required for their hatching. However, our results should be considered with caution, since potential biases could have occurred due to the mosquito capture design adopted. Further sampling during the dry season and over several years is strongly recommended.

The mosquito population model developed and validated by our field trapping data was able to reproduce Cx. pipiens abundance and dynamics using three environmental and meteorological inputs (temperature, rainfall and proportion of flooded surface). It highlights the significant role of permanent water in this particular region and the relative influence of both rainfall and flooded surface at local level depending on the distance to the banks of the Delta. To the best of our knowledge, it was also the first time that the population dynamics of Cx. pipiens in a semi-arid environment was modelled using a mechanistic approach. Our model allowed testing different assumptions on the main drivers of Cx. pipiens dynamics in the Okavango region and, despite the low number of catching episodes, it accurately predicted the spatial variability of Cx. pipiens abundance on the western side of the Delta and the major trends in the annual fluctuations of host-seeking adults in the two trapping sites closer to the Okavango River which are probably more strongly influenced by flooding.

Our study is a necessary first step in the understanding the potential dynamics of RVF in the Okavango Delta. However, increased knowledge on RVFV epidemiology requires additional field data from exposed human and animal populations, complementary vector sampling and continued development of modelling techniques for exploring plausible disease transmission scenarios, emergence mechanisms and the potential impact of intervention strategies.


asl, above sea level; MODIS, moderate resolution imaging apectroradiometer; ODE, ordinary differential equations; RVF, Rift Valley fever; RVFV, Rift Valley fever virus; TRMM, tropical rainfall measuring.


  1. Balkhy HH, Memish ZA. Rift Valley fever: an uninvited zoonosis in the Arabian peninsula. Int J Antimicrob Agents. 2003;21(2):153–7. Available from:

    Article  CAS  PubMed  Google Scholar 

  2. Linthicum KJ, Anyamba A, Britch SC, Chretien JP, Erickson RL, Small J, et al. A Rift Valley fever risk surveillance system for Africa using remotely sensed data: potential for use on other continents. Vet Ital. 2007;43(3):663–74.

    PubMed  Google Scholar 

  3. Pepin M, Bouloy M, Bird BH, Kemp A, Paweska J. Rift Valley fever virus (Bunyaviridae: Phlebovirus): an update on pathogenesis, molecular epidemiology, vectors, diagnostics and prevention. Vet Res. 2010;41(6):61.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Nicolas G, Durand B, Rakotoarimanana TT, Lacote S, Chevalier V, Marianneau P. A 3-year serological and virological cattle follow-up in Madagascar highlands suggests a non-classical transmission route of Rift Valley fever virus. Am J Trop Med Hyg. 2014;90(2):265–6.

  5. Nanyingi MO, Munyua P, Kiama SG, Muchemi GM, Thumbi SM, Bitek AO, et al. A systematic review of Rift Valley fever epidemiology 1931–2014. Infect Ecol Epidemiol. 2015;5.

  6. Coetzer J. The pathology of Rift Valley fever. I. Lesions occurring in natural cases in new-born lambs. Onderstepoort J Vet Res. 1977;44(4):205–11.

    CAS  PubMed  Google Scholar 

  7. Swanepoel R, Coetzer J. Rift valley fever. Infect Dis Livestock. 2004;2:1037–70.

    Google Scholar 

  8. Bird BH, Ksiazek TG, Nichol ST, Maclachlan NJ. Rift Valley fever virus. J Am Vet Med Assoc. 2009;234(7):883–93.

    Article  PubMed  Google Scholar 

  9. WHO. Fact sheet. World Health Organisation; 2010. 207.

  10. Peyre M, Chevalier V, Abdo-Salem S, Velthuis A, Antoine-Moussiaux N, Thiry E, et al. A systematic scoping study of the socio-economic impact of Rift Valley fever: Research gaps and needs. Zoonoses Public Health. 2014;p. 1–17. Available from:

  11. Jori F, Alexander KA, Mokopasetso M, Munstermann S, Moagabo K, Paweska JT. Serological evidence of Rift Valley fever virus circulation in domestic cattle and African buffalo in Northern Botswana (2010–2011). Front Vet Sci. 2015;2(63):1–7.

  12. Davies FG, Linthicum KJ, James AD. Rainfall and epizootic Rift Valley fever. Bull World Health Organ. 1985;63(5):941–3. PMID: 3879206 PMCID: PMC2536443. Available from:

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Linthicum KJ, Anyamba A, Tucker CJ, Kelley PW, Myers MF, Peters CJ. Climate and satellite indicators to forecast Rift Valley fever epidemics in Kenya. Sci (New York, NY). 1999;285(5426):397–400.

    Article  CAS  Google Scholar 

  14. Swanepoel R, Paweska J. Zoonoses: Biology, Clinical Practice, and Public Health Control - Chapter 36 Rift Valley fever. In: Oxford textbook of zoonoses: biology, clinical practice, and public health control. New York: Oxford University Press; 2011. p. 423-31.

  15. Mweya CN, Kimera SI, Mellau LS, Mboera LE. Inter-epidemic abundance and distribution of potential mosquito vectors for Rift Valley fever virus in Ngorongoro district, Tanzania. Global Health Action. 2015;8

  16. Clements AN. The Biology of Mosquitoes: Development, Nutrition and Reproduction. London, UK: Chapman & Hall; 1992.

  17. Manore CA, Beechler BR. Inter-epidemic and between-season persistence of rift valley fever: vertical transmission or cryptic cycling? Transbound Emerg Dis. 2015;62(1):13-23.

  18. Cailly P, Tran A, Balenghien T, L’Ambert G, Toty C, Ezanno P. A climate-driven abundance model to assess mosquito control strategies. Ecol Modell. 2012;227:7–17. Available from:

    Article  Google Scholar 

  19. Loetti V, Schweigmann N, Burroni N. Development rates, larval survivorship and wing length of Culex pipiens (Diptera: Culicidae) at constant temperatures. J Nat Hist. 2011;45(35–36):2203–13.

  20. Mendelson J, Vanderpost C, Ramberg L, Hudson MM, Mosepele K. Okavango Delta: Floods of Life. World Conservation Union. RAISON for IUCN and HOORC: Windhoek, Namibia; 2010.

    Google Scholar 

  21. Ioris AAR. Tropical wetland management: The South-American Pantanal and the international experience. Ioris AAR, editor. Routledge: University of Edinburgh, UK; 2017.

    Google Scholar 

  22. Wolski P, Murray-Hudson M, Gumbricht T. Modelling of the hydrology of the Okavango Delta - Water and Ecosystem Resources for Regional Development. Harry Oppenheimer Okavango Research Centre - University of Botswana; 2005

  23. Botswana S. Botswana population and housing census. Dep Printing Pub Serv: Gaborone. 2011;211:4.

    Google Scholar 

  24. Chase M. Dry season fixed-wing aerial survey of elephants and wildlife in northern Botswana. Elephants without Borders, San Diego Zoological Society, Department of Wildlife and National Parks Report. 2011;p. 138.

  25. Tchouassi DP, Sang R, Sole CL, Bastos AD, Teal PE, Borgemeister C, et al. Common host-derived chemicals increase catches of disease-transmitting mosquitoes and can improve early warning systems for Rift Valley fever virus. PLoS NeglTrop Dis. 2013;7(1), e2007.

    Article  Google Scholar 

  26. Magnarelli L. Relative abundance and parity of mosquitoes collected in dry-ice baited and unbaited CDC Miniature Light Traps. Mosq News. 1975;35(3):350–3.

    Google Scholar 

  27. Jupp PG. Mosquitoes of southern Africa: Culicinae and Toxorhynchitinae. Hartebeespoort: Ekogilde Publishers; 1996.

    Google Scholar 

  28. Edwards FW, Hopkins GHE, Evans AM. Mosquitoes of the Ethiopian Region III: Culicine adults and pupae. London, UK: Trustees, British Museum (Natural History), B. Quaritch, Limited; 1941.

    Google Scholar 

  29. Gillies M, De Meillon B. The Anophelinae of Africa South or the Sahara (Ethiopian Zoogeographical Region). S Afr Inst Med Res. 1968;54:1–343.

  30. Gillies M, Coetzee M. A Supplement to the Anophelinae of Africa South of the Sahara. S Afr Inst Med Res. 1987;55:1–143.

    Google Scholar 

  31. Jupp PG, Kemp A, Grobbelaar A, Lema P, Burt FJ, Alahmed AM, et al. The 2000 epidemic of Rift Valley fever in Saudi Arabia: mosquito vector studies. Med Vet Entomol. 2002;16(3):245–52.

    Article  CAS  PubMed  Google Scholar 

  32. McIntosh B, Jupp P. Epidemiological aspects of Rift Valley fever in South Africa with reference to vectors. Proceed Rift Valley Fever Workshop. 1981;9:92.

    Google Scholar 

  33. Turell MJ, Linthicum KJ, Patrican LA, Davies FG, Kairo A, Bailey CL. Vector competence of selected African mosquito (Diptera: Culicidae) species for Rift Valley fever virus. J Med Entomol. 2008;45(1):102–8.

    Article  PubMed  Google Scholar 

  34. Arum SO, Weldon CW, Orindi B, Landmann T, Tchouassi DP, Affognon HD, et al. Distribution and diversity of the vectors of Rift Valley fever along the livestock movement routes in the northeastern and coastal regions of Kenya. Parasit Vectors. 2015;8(1):294.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Huffman G, Bolvin D. TRMM and Other Data Precipitation Data Set Documentation. 2014.

  36. Alto BW, Muturi EJ, Lampman RL. Effects of nutrition and density in Culex pipiens. Med Vet Entomol. 2012;26(4):396–406.

  37. Jetten TH, Takken W. Anophelism Without Malaria in Europe: A review of the ecology and distribution of the genus Anopheles in Europe. Wageningen Agricultural University papers. 6708 PB Wageningen: Wageningen Agricultural University; 1994.

  38. Jupp P, McIntosh B. Ecological studies on Sindbis and West Nile viruses in South Africa. II. Mosquito bionomics. S Afr J Med Sci. 1967;32(1):15.

    CAS  PubMed  Google Scholar 

  39. Jupp P. Preliminary studies on the overwintering stages of Culex mosquitoes (Diptera: Culicidae) in the highveld region of South Africa. Entomol Soc Southern Afr J. 1969;32:91–8.

  40. Jupp P. Further studies on the overwintering stages of Culex mosquitoes (Diptera: Culicidae) in the highveld region of South Africa. J Entomol Soc South Afr. 1975;38(1).

  41. Saltelli A, Chan K, Scott EM, et al. Sensitivity analysis, vol. 134. Wiley New York: New York, USA; 2000.

    Google Scholar 

  42. Tran A, L’Ambert G, Lacour G, Benoît R, Demarchi M, Cros M, et al. A rainfall- and temperature-driven abundance model for Aedes albopictus populations. Int J Environ Res Public Health. 2013;10(5):1698–719. PMID: 23624579 PMCID: PMC3709343.

  43. Scilab Entreprises. Scilab: Free and Open Source software for numerical computation (OS Version 5.4.1) [Software]. Available from: 2012.

  44. R Development Core Team. R: A language and environment for statistical computing. Vienna, Austria: 1999–2012 R Foundation for Statistical Computing; 2011.

  45. Prehaud C, Bouloy M. La fièvre de la vallée du Rift : Un modèle d’étude des fièvres hémorragiques virales. Annales de l’Institut Pasteur Actualités. 2007;8(3):233–44. Available from:

    Article  Google Scholar 

  46. Modou M. La fièvre de la Vallée du Rift dans l’Afrique Sub-saharienne. 2012. Institut Sénégalais de Recherches Agricoles.

    Google Scholar 

  47. Turell MJ. Gargan 2nd T, Bailey CL. Replication and dissemination of Rift Valley fever virus in Culex pipiens. Am J Trop Med Hyg. 1984;33(1):176–81.

  48. Amraoui F, Krida G, Bouattour A, Rhim A, Daaboub J, Harrat Z, et al. Culex pipiens, an experimental efficient vector of West Nile and Rift Valley fever viruses in the Maghreb region. PLoS One. 2012;7(5), e36757.

  49. Gillett J. Contributions to the oviposition-cycle by the individual mosquitoes in a population. J Insect Physiol. 1962;8(6):665–81. Available from:

    Article  Google Scholar 

  50. In Ecology of Mosquitoes: proceedings of a workshop (L P Lounibos JRR, H Frank E. The nature of egg hatching in Aedes triseriatus: ecological implications and evolutionary consequences. 1985.

  51. Ahumada JA, Lapointe D, Samuel MD. Modeling the population dynamics of Culex quinquefasciatus (Diptera: Culicidae), along an elevational gradient in Hawaii. J Med Entomol. 2004;41(6):1157–70.

  52. Fall AG, Diaté A, Lancelot R, Tran A, Soti V, Etter E, et al. Feeding behaviour of potential vectors of West Nile virus in Senegal. Parasit Vectors. 2011;4(1).

  53. Jones R, Barker-Hudson P, Kay B. Comparison of dry ice baited light traps with human bait collections for surveillance of mosquitoes in northern Queensland. Australia J Am Mosq Control Assoc. 1991;7(3):387–94. Available from:

    CAS  PubMed  Google Scholar 

  54. Morris JA, Lampman RL, Ballmes G, Funes J, Halvorsen J, Novak RJ. First record of Aedes japonicus japonicus in Illinois: defining its spatial distribution and associated mosquito species. J Am Mosq Control Assoc. 2007;23(3):243–51. Available from:

  55. Cailly P, Balenghien T, Ezanno P, Fontenille D, Toty C, Tran A. Role of the repartition of wetland breeding sites on the spatial distribution of Anopheles and Culex, human disease vectors in Southern France. Parasit Vectors. 2011;6:65.

  56. Rosà R, Marini G, Bolzoni L, Neteler M, Metz M, Delucchi L, et al. Early warning of West Nile virus mosquito vector: climate and land use models successfully explain phenology and abundance of Culex pipiens mosquitoes in north-western Italy. Parasit Vectors. 2014;7(1):269.

  57. Loncaric Ž, Hackenberger BK. Stage and age structured Aedes vexans and Culex pipiens (Diptera: Culicidae) climate-dependent matrix population model. Theor Popul Biol. 2013;83:82–94.

  58. Tessier SF, Rollin PE, Sureau P. Viral haemorrhagic fever survey in Chobe (Botswana). Trans R Soc Trop Med Hyg. 1987;81(4):699–700. Available from:

  59. Barker CM. Spatial and temporal patterns in mosquito abundance and virus transmission in California. ProQuest: University of California, USA; 2008.

    Google Scholar 

  60. Sang R, Kioko E, Lutomiah J, Warigia M, Ochieng C, O’Guinn M, et al. Rift Valley fever virus epidemic in Kenya, 2006/2007: the entomologic investigations. Am J Trop Med Hyg. 2010;83(2 Suppl):28–37.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Rolin AI, Berrang-Ford L, Kulkarni MA. The risk of Rift Valley fever virus introduction and establishment in the United States and European Union. Emerg Microbes Infect. 2013;2(12), e81.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Lutomiah J, Bast J, Clark J, Richardson J, Yalwala S, Oullo D, et al. Abundance, diversity, and distribution of mosquito vectors in selected ecological regions of Kenya: public health implications. J Vector Ecol. 2013;38(1):134–42.

    Article  PubMed  Google Scholar 

  63. McIntosh B, Jupp P, Dos Santos I, Barnard B, et al. Vector studies on Rift Valley fever virus in South Africa. S Afr Med J. 1980;58(3):127–32.

    CAS  PubMed  Google Scholar 

  64. Martin V, Chevalier V, Ceccato P, Anyamba A, De Simone L, Lubroth J, et al. The impact of climate change on the epidemiology and control of Rift Valley fever. Rev Sci Tech OIE. 2008;27(2):413–26.

    Article  CAS  Google Scholar 

  65. Tourre YM, Lacaux JP, Vignolles C, Lafaye M. Climate impacts on environmental risks evaluated from space: a conceptual approach to the case of Rift Valley Fever in Senegal. Global health action. 2009;2. PMID: 20052381 PMCID: PMC2799327.

  66. Konrad SK, Miller SN. A temperature-limited assessment of the risk of Rift Valley fever transmission and establishment in the continental United States of America. Geospat Health. 2012;6(2):161–70.

    Article  PubMed  Google Scholar 

  67. Ezanno P, Aubry-Kientz M, Arnoux S, Cailly P, L'Ambert G, Toty C, Balenghien T, Tran A. A generic weather-driven model to predict mosquito population dynamics applied to species of Anopheles, Culex and Aedes genera of southern France. Prev Vet Med. 2015;120(1):39-50.

  68. Vanwambeke SO, Somboon P, Harbach RE, Isenstadt M, Lambin EF, Walton C, et al. Landscape and land cover factors influence the presence of Aedes and Anopheles larvae. J Med Entomol. 2007;44(1):133–44.

  69. Tran A, Ponçon N, Toty C, Linard C, Guis H, Ferré 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 Geographics. 2008;7(1):9.

  70. McMahon B, Manore C, Hyman J, LaBute M, Fair JM. Coupling vector-host dynamics with weather geography and mitigation measures to model Rift Valley fever in Africa. Math Model Nat Phenom. 2014;9(02):161–77.

  71. Mweya CN, Holst N, Mboera LEG, Kimera SI. Simulation modelling of population dynamics of mosquito vectors for Rift Valley fever virus in a disease epidemic setting. PLoS One. 2014;9(9), e108430. Available from:

  72. Nicolas G, Chevalier V, Tantely LM, Fontenille D, Durand B. A Spatially explicit metapopulation model and cattle trade analysis suggests key determinants for the recurrent circulation of Rift Valley fever virus in a pilot area of Madagascar highlands. PLoS NeglTrop Dis. 2014;8(12):e3346.

  73. Guis H, Caminade C, Calvete C, Morse AP, Tran A, Baylis M. Modelling the effects of past and future climate on the risk of bluetongue emergence in Europe. J R Soc Interface. 2011;9(67):339–50. PMID: 21697167 PMCID: PMC3243388.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Dumont Y, Chiroleu F. Vector control for the Chikungunya disease. Math Biosci Eng. 2010;7(2):313–45.

    Article  PubMed  Google Scholar 

  75. Vinogradova EB. Culex pipiens pipiens mosquitoes: Taxonomy, distribution, ecology, physiology, genetics, applied importance and control. Sofia, Moscow: Pensoft Publishers; 2000.

  76. Vasuki V. Effect of insect growth regulators on hatching of eggs of three vector mosquito species. Proc Anim Sci. 1990;99(6):477–82. Available from:

    Article  CAS  Google Scholar 

  77. Zequi JAC, Lopes J. Development of the immature stages of Culex (Culex) saltanensis Dyar (Diptera, Culicidae) under laboratory conditions. Rev Bras Entomol. 2012;56(1):106–10. Available from:

  78. Lutambi AM, Penny MA, Smith T, Chitnis N. Mathematical modelling of mosquito dispersal in a heterogeneous environment. Math Biosci. 2013;241(2):198–216.

    Article  PubMed  Google Scholar 

  79. Craig M, Snow R, Le Sueur D. A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitol Today. 1999;15(3):105–11. Available from:

    Article  CAS  PubMed  Google Scholar 

  80. Shaman J, Spiegelman M, Cane M, Stieglitz M. A hydrologically driven model of swamp water mosquito population dynamics. Ecol Modell. 2006;194(4):395–404. Available from:

    Article  Google Scholar 

Download references


We are very grateful to the Entomology Unit from the National Veterinary Laboratory of Botswana and the University of Montpellier (France) for supporting this project. We are also so grateful to M. Poitr Wolski from the Okavango Research Institute and to the Maun Weather station for their help on accessing climatic data and inputs on the complex hydrological system of the Delta and to Mrs. Keitumetse Gladys Mangate from the Botswana National Veterinary Laboratory for her instrumental assistance in the field. Finally, we would like to acknowledge Dr. Catherine Sole from the Department of Zoology and Entomology for her help in organizing the trap design and the provision of entomological references.


Funding was provided by FAO. Award Number: LoA OSRO/INT/602/USA/B1.

Availability of data and material

Data and materials available on request.

Authors’ contributions

HP performed the model adaptation and drafted the manuscript. TA participated in the design of the study and the model adaptation and performed the flooding data acquisition and drafted the manuscript. KA performed the virological diagnostics, participated in the mosquito identification and provided entomological interpretation of the results. TP carried out the mosquitoe identifications. KP helped designing the sampling design and carried out the mosquitoes catches. CV helped to draft the manuscript and provided RVF knowledge. PJ coordinated the mosquito identification and viral detection activities. JF designed and coordinated the study, he also coordinated and participated in field activities and drafted the manuscript. All authors read and approved the final manuscript.

Authors’ information

CV and PJ are Rift Valley Fever specialists holding several publications about this disease in sub-Saharan Africa.

KA is a recognized clinical entomologist in Southern Africa and has monitored several outbreaks of RVF in Southern Africa.

TA is an entomological modeler with several publications in this field in tropical and temperate environments. She participated in building the original model with Cailly Priscilla and in further adaptations.

Competing interests

The authors declare that they have no competing interests (financial and non-financial).

Consent for publication

Not applicable.

Ethics approval and consent to participate

The use of laboratory mice (strain NIH) for inoculation of mosquito homogenates for virus isolation was approved by the Animal Ethics Committee of the South African National Health Laboratory Service (Clearance Certificate No. 124/11), according to the South African National Standard for The Care and Use of Animals for Scientific Purposes, SANS 10386:2008 (ISBN 978-0-626-22296-3).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Hammami Pachka.

Additional files

Additional file 1:

Serological analysis (Sampling strategy, laboratory analysis and results). (PDF 71 kb)

Additional file 2:

Mapping flooding extent method. Figure in Additional file 2. Maps of Modified Normalized Difference Water Index (MNDWI) derived from MODIS imagery at different dates corresponding to the study period. (PDF 3185 kb)

Additional file 3

Detail of the ordinary differential equation system. (PDF 139 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Pachka, H., Annelise, T., Alan, K. et al. Rift Valley fever vector diversity and impact of meteorological and environmental factors on Culex pipiens dynamics in the Okavango Delta, Botswana. Parasites Vectors 9, 434 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: