Next Article in Journal
Current Knowledge on Chemosensory-Related Candidate Molecules Potentially Involved in Tick Olfaction via Haller’s Organ
Next Article in Special Issue
Machine Learning Modeling of Aedes albopictus Habitat Suitability in the 21st Century
Previous Article in Journal
A Review of the Host Plant Location and Recognition Mechanisms of Asian Longhorn Beetle
Previous Article in Special Issue
The Effect of Environmental Factors on the Diversity of Crane Flies (Tipulidae) in Mountainous and Non-Mountainous Regions of the Qinghai-Tibet Plateau and Surrounding Areas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predicting Culex pipiens/restuans Population Dynamics Using a Weather-Driven Dynamic Compartmental Population Model

1
Institute for Medical Microbiology and Hygiene, AGES—Austrian Agency for Health and Food Safety, 1090 Vienna, Austria
2
Department of Veterinary and Animal Sciences, University of Copenhagen, Frederiksberg Campus, 1870 Copenhagen, Denmark
3
Complexity Science Hub Vienna, 1080 Vienna, Austria
*
Author to whom correspondence should be addressed.
Insects 2023, 14(3), 293; https://doi.org/10.3390/insects14030293
Submission received: 28 February 2023 / Revised: 14 March 2023 / Accepted: 15 March 2023 / Published: 17 March 2023
(This article belongs to the Special Issue Climate Sensitive Ecological and Dynamical Models of Insects)

Abstract

:

Simple Summary

In this study, we present a compartmental population model for Cx. pipiens/restuans, incorporating mosquito life cycle parameters, as well as temperature, precipitation, and geographic latitude. The model was validated against mosquito count data from Cook County (IL, USA). The model fitted the observation data and was able to reproduce between-year differences in the abundance of the Cx. pipiens/restuans mosquitoes, as well as the different seasonal trends.

Abstract

Mosquitoes of the genus Culex are important vectors of a variety of arthropod-borne viral infections. In most of the northern parts of the USA, Cx. pipiens/restuans is the predominant representative of this genus. As vectors, they play a key role in the spreading of arboviruses and thus, knowledge of the population dynamic of mosquitoes is important to understand the disease ecology of these viruses. As poikilotherm animals, the vital rates of mosquitoes are highly dependent on ambient temperature, and also on precipitation. We present a compartmental model for the population dynamics of Cx. pipiens/restuans. The model is driven by temperature, precipitation, and daytime length (which can be calculated from the geographic latitude). For model evaluation, we used long-term mosquito capture data, which were averaged from multiple sites in Cook County, Illinois. The model fitted the observation data and was able to reproduce between-year differences in the abundance of the Cx. pipiens/restuans mosquitoes, as well as the different seasonal trends. Using this model, we evaluated the effectiveness of targeting different vital rates for mosquito control strategies. The final model is able to reproduce the weekly mean Cx. pipiens/restuans abundance for Cook County with a high accuracy, and over a long time period of 20 years.

1. Introduction

Mosquitoes of the genus Culex (family Culicidae) are important vectors of a variety of infectious diseases, e.g., Rift Valley fever, St. Louis encephalitis, Japanese encephalitis, Western equine encephalitis, and West Nile fever [1,2,3,4]. As vectors, they play a key role in virus transmission cycles, and thus in spreading viral infections. Therefore, in order to understand the epidemiology of mosquito borne arboviruses, knowledge of vector population dynamics is necessary. The main vectors for West Nile virus (WNV) in North America are mosquitoes of the genus Culex—especially Cx. pipiens, Cx. tarsalis, Cx. quinquefasciatus, and Cx. restuans—but other Culicidae have also been shown to transmit WNV [4,5]. According to their distribution, Cx. pipiens is the main vector in the north, Cx. tarsalis is the main vector in the west, and Cx. quinquefasciatus is the main vector in the south of North America [5,6,7].
Being poikilotherm animals, the population dynamics of Culex mosquitoes are highly influenced by ambient temperature [8,9,10]. Laboratory studies have shown the influence of temperature on the development time and survival rates of eggs, larvae, and pupae [11,12,13,14,15]. In adult mosquitoes, temperature affects their lifespan, length of the gonotrophic cycle, and virus transmission rate [11,12,14]. Some Culex mosquitoes undergo reproductive diapause, which is also, along with daytime length, affected by temperature [16,17]. Precipitation has also been shown to have a strong impact on Culex abundance: a high abundance of Culex mosquitoes has been found to be associated with high precipitation occurring several weeks before a capture event [9,10]. The authors of these studies suggested that rainfall increases the water’s surface area and, therefore, possible oviposition sites. Other factors which have been shown to influence the abundance of Culex mosquitoes are landscape patterns and land use [9,18,19].
In general, there are two types of models which have been used to describe Culex population dynamics: regression models and compartmental models. Regression models have proven to be useful tools in discovering connections between population densities and environmental factors [9,20,21,22]. However, they cannot reveal the biological mechanisms behind those relationships. Regression models are, therefore, of limited use if we want to discover the mechanisms of population dynamics, which is necessary knowledge if we want to manipulate a species population, such as for pest control. Compartmental models, on the other hand, are aimed to represent the mechanisms behind population dynamics. This makes compartmental models ideal for testing as to how the population dynamics of a particular species could be influenced. However, the downside of compartmental models is their need for detailed information on the vital rates of the investigated species. Examples for compartmental Culex models are rather rare, requiring detailed information on the location the model is applied to (e.g., the water volume of the breeding site), and/or have only been validated for short time periods [23,24,25]. Therefore, it would be difficult to apply these models to other locations.
In epidemiology, simplified compartmental Culex models have been integrated in susceptible, infectious, and/or recovered (SIR) models of mosquito-borne infectious diseases [26,27,28]. In those studies, however, abundance predictions were not compared to actual mosquito capture data, and, thus, these models were not validated. Yu et al. [29] developed a compartmental abundance model for Cx. pipiens/restuans that did validate model results against actual mosquito abundance data. Besides mosquito life cycle parameters obtained from field and laboratory studies, this model was, however, solely temperature-driven. The researchers suggested including other factors impacting mosquito abundance, such as the precipitation, landscape, and wind, in order to help improve the predictive power of future models [29].
In this study, we present a compartmental population model for Cx. pipiens/restuans, incorporating mosquito life cycle parameters, as well as temperature, precipitation, and geographic latitude. The model was validated against mosquito count data from Cook County (IL, USA). The aim was not to generate a model for a specific area, but rather to produce a model which generates mean values for mosquito abundance that are valid for larger areas (here, we use the county level), which could be used for other areas as well. Thus, we use only easily accessible input parameters (temperature, precipitation, and geographic latitude). This model could therefore be simply integrated in SIR models of Culex-borne infectious diseases. Further, we test how manipulating the vital rates influences population size, which is important for mosquito control strategies. Control strategies could be aimed at increasing the mortality rate of certain stages, by using insecticides or biological pest control measurements [30], or, they could be aimed at decreasing the birth rate of mosquitoes, e.g., by releasing sterile individuals (Sterile Insect Technique) [31].

2. Materials and Methods

2.1. Cx. pipiens/restuans Population Model

To simulate Cx. pipiens/restuans population dynamics, we developed a dynamic compartmental population model in R 4.1.2 [32], consisting of a system of ordinary differential equations (OED) with daily time steps. The model is based on the models presented by Laperriere et al. [26] and Rubel et al. [27].
Our population model consists of four stage variables: E—egg stage, L—larval stage, including pupal stage, A—active adult mosquitoes, H—adult mosquitoes in diapause (Figure 1). Furthermore, we only allowed the proportion of surviving individuals of one stage to have the possibility to change into the next stage.
The daily changes of the number of animals within each stage are described by Equations (1)–(4), where definitions of the parameters used in the equations are described in Table 1.
d   E d   t = b E A m E E b L E 1 m E
d   L d   t = b L E 1 m E m L K L / L L b A L 1 m L K L / L
d   A d   t = b A L 1 m L K L / L 1 m A h A + 1 m H 1 h H
d   H d   t = m H H + 1 m A h A 1 m H 1 h H
The transition parameters used (Table 1) were estimated from literature values and were influenced only by temperature (T), precipitation (P), and daytime length (D). We combined data from different sources, giving values at different temperatures, and we also used data from other Culex species if there was not sufficient data for Cx. pipiens/restuans available. By visual inspection of the distribution of the data points, we appraised the form of the function and fitted functions through the data points, using the non-linear least square approach (nls package in R 4.1.2, [32], see Equations (5)–(14)). Estimates were rounded to 2 significant digits. Some functions (e.g., mortality rates) would allow for values above 1 at certain temperatures, which are, from a biological point of view, not possible. Thus, if necessary for a function, it was capped at the value of 1. In this case, we calculated the intersection with a horizontal line at 1. Those interception points represent the lower and upper temperature threshold for the validity of a specific function. Below and above those thresholds, the value was set to 1. We arbitrarily set initial values for E, L, and A to 0.0 and 0.1, for H to have a starting point. To minimize the effects of the chosen initial values, we started the model one year before the first mosquito capture data was available, on the date 1 January 1990. For the next sections, the parameters a, b, and c in the function descriptions do not have a biological meaning, but are necessary to shape the parameter functions.

2.1.1. Egg Stage—E

The birth rate of eggs (bE) represents the daily egg-laying rate. It can be calculated as the reciprocal value of the length of the gonotrophic cycle. Madder et al. [11], Vinogradova [15] and Elizondo-Quiroga et al. [33] give temperature (T) dependent values for the gonotrophic cycle in the Cx. pipiens complex, which we used to fit a logistic function:
b E T = a 1 + b exp c   T
Various authors have suggested the influence of precipitation on the birth rate [9,10], however, no quantitative data describing this relationship were available. As the availability of suitable Cx. pipiens/restuans breeding sites increases with preceding rain events, we included an additional term in order to increase bE: the mean sum of precipitation (P) of the k days before egg laying. Further, to obtain the “actual” birth rate, i.e., the number of female offspring per female per day, we multiplied bE by a value equaling half of the mean number of eggs found in an egg raft, as it has been shown that there is an equal number of male and female eggs [34]. For Cx. pipiens, this value was estimated to be ≈100 (calculated from Madder et al. [11], Oda and Ueda [35], and Oda et al. [36]).
b E T ,   P = 100 a + d   i = 1 k P t i 1 + b exp c   T
Data on the temperature dependence of the egg mortality rate (mE) were given by Madder et al. [11], Shriver and Bickley [37], and Suman et al. [14]. We assumed a polynomic function:
m E T = a + b T + c T 2 ,   if   t r l m E < T < t r u m E 1 ,   otherwise
where trl and tru are lower and upper temperature thresholds for the validity of the function, respectively, as the mortality rate cannot exceed 1.

2.1.2. Larval Stage—L

In our model, we combined the larval and pupal stages of the mosquitos into one compartment (hereafter only referred to as the “larval stage”). The larval birth rate (bL) is the daily hatching rate of eggs into larvae, and is calculated as the reciprocal value of the duration of the egg stage. The temperature-dependent larval birth rate has been well documented, and can be described by a logistic function [11,14,15,37,38,39,40]:
b L T = 1 1 + a exp b   T
The temperature-dependent mortality rate of the larval stage (mL) has been described for Cx. pipiens by Madder et al. [11], Olejniček and Gelbič [41], and Tamarina [42]). The data fitted a rational function:
m L T = a + b   T T T m i n . l a r v T T m a x . l a r v ,   if   t r l m L < T < t r u m L 1 ,   otherwise
where trl and tru are lower and upper temperature thresholds for the validity of the function, respectively, as the mortality rate cannot exceed 1.
This function required information on the minimum and maximum temperature for larval survival (Tmin.larv and Tmax.larv). Mean values found in the literature were Tmin.larv = 8.7 °C [11,43,44,45] and Tmax.larv = 36.4 °C [43,46]. Mortality rates are likely to be biased and highly underestimated under lab conditions. For Cx. tarsalis, it has been shown that only 10% of larval mortality in the field is due to abiotic factors [34]. Thus, we increased obtained mortality rates under lab conditions for larvae by a factor of 10.
The mortality of larvae has been shown to increase with larval density [11,47,48], which we incorporated in our model as well. Suitable densities for larval development are likely to depend on habitat structure, and, thus, the larval carrying capacity KL. Unfortunately, this parameter, as well as its relation to the number of adult mosquitoes, is unknown. Assuming a correlation between the number of larval and adult mosquitoes, and, thus, between the carrying capacity of the two stages, we estimated a scaling factor, which gives KL in relation to the capture rate of adult females. Using a scaling factor instead of estimating KL directly makes the results of this study applicable to other areas. We used the mean number of captured Cx. pipiens/restuans during their peak abundance (i.e., in July), which was found to be 5.75 individuals per day for our study area. As the volume of larval habitat increases with preceding rain events, we included an additional term to increase KL: the mean sum of precipitation of the k preceding days. Please note that k is the same as for Equation (6).
K L P ; t = a   5.75 + b   1 k i = 1 k P t i

2.1.3. Adult Stages—A and H

The birth rate of adult mosquitoes (bA) represents the hatching rate of mosquitoes, from larvae to imagoes. Analogous to the larval birth rate, the adult birth rate can be calculated as the reciprocal value of the duration of the larval period. Temperature-dependent durations of the larval period for the Cx. pipiens complex are given in [13,38,40,46,49]. Again, the data fitted a logistic function:
b A T = 1 1 + a exp b   T
Almir’on and Brewer [38], Gunay et al. [50], Oda et al. [51], Reisen et al. [12] and Suman et al. [14] give data on the temperature dependence of adult mortality rates (mA). We assumed a polynomic function:
m A T = a + b T + c T 2 ,   if   t r l m A < T < t r u m A 1 ,   otherwise
where trl and tru are lower and upper temperature thresholds for the validity of the function, respectively, as the mortality rate cannot exceed 1.
In winter, Cx. pipiens and Cx. restuans undergo reproductive diapause [16,17,48]. This diapause is mainly characterized by the arrest of ovarian development, but also accompanied by other metabolic and behavioral changes [15]. Ovarian development, which we will use as a proxy for the occurrence of diapause, is dependent on the temperature (T) and daytime length (D), as investigated for Cx. pipiens by Eldridge [52], and a logistic function can be fitted to the data:
  h T , D = 1 1 + a exp b   T + c   D
Therefore, h is the rate at which an adult mosquito changes from the active stage, A, into the diapausing (hibernating) stage, H, and 1 h is the rate at which adult mosquitoes terminate diapause. The mosquitoes spend hibernation in natural and artificial shelters, where they are protected from the harsh climatic conditions during winter [15]. Bailey et al. [53] reported that for diapausing Cx. pipiens, the proportions of individuals surviving the winter (14 December—8 March; equates to 84 days) were 0.157, 0.224 and 0.247. Thus, we can calculate the daily survival rate during diapause in winter, and hence the daily winter mortality rate:
m H = 1 e x p log a 84

2.2. Model Evaluation

2.2.1. Mosquito Count Data

For the evaluation of our model, we used mosquito count data from Cook County (IL, USA), provided by the Desplaines Valley Mosquito Abatement District (P. Geery, 2011, pers. comm.). This agency had long-term capture data available, as it has been capturing mosquitoes with New Jersey light traps at 25 trap sites since 1980.
However, not all trap sites were continuously operated during all years, and, thus, only those trap sites providing continuous time series for at least 20 years were selected, resulting in six capture sites (details on the trap location, data preparation and mosquito data selection are published by Lebl et al. [22]). For this study, the numbers of female mosquitoes were averaged over these six capture sites. Although the sites are not uniformly distributed within the county, the resulting average time series was considered to be representative for Cook County. Of the 81,022 Culex females caught at the trap sites between 1991 and 2010, 97% were Cx. pipiens/restuans, and the rest were Cx. territans, Cx. tarsalis, Cx. erraticus, and Cx. salinarius. The catches of the morphologically rather similar species Cx. pipiens and Cx. restuans were combined, as distinguishing them based on morphological characteristics is rather imprecise.

2.2.2. Environmental Data

The mean daily temperature (T) and daily sum of precipitation (P) data were taken from the weather station of the Chicago O’Hare International Airport (WMO No. 72530 [54]; 87.900° W, 41.983° N and altitude 205 m). The daytime length, D (in hours), is a function of the sun’s declination, E, and the geographic latitude ϕ, and was calculated using the CBM model described by Forsythe et al. [55]:
D = 24 24 π acos sin 0.8333 + sin φ sin ϵ cos φ cos ϵ
with the declination of the sun ε calculated as a function of the calendar day d:
ϵ = asin 0.39795   cos 0.2163108 + 2   atan 0.9671396   tan 0.00860   d 186  

2.2.3. Model Calibration and Validation

We used function estimates as initial values for our model. For the parameter associated with precipitation, our initial value was chosen as 0, i.e., precipitation had no effect. The model results had a temporal resolution of 1 day. Due to the high inter-day variation in the capture rates [22], we calculated the weekly mean of the capture rates and model estimates to be used for further analysis. To evaluate model fit, we calculated the root mean square error (RMSE) between the mosquito capture data and the model estimates for the active adult mosquitoes, A, as well as the Spearman’s correlation coefficient, rS (as the count data followed a non-Gaussian distribution).
The individual function parameters (a, b, and c) for the transition parameters were estimated based on literature values obtained at laboratory conditions, thus, it is possible that the parameters of these functions would be different at field conditions. Consequently, we allowed for the function parameters to vary within the standard error (SE) of their estimated value (Table 2). To divide the data set into a training and a test set, we distinguished between even and odd years. In a previous study, we used a random process to choose odd years as the test data, and we retained this division in order to make the results comparable [22].
With the training data set, a genetic algorithm was used, in order to identify the parameter combination that fitted the data best (i.e., minimizing RMSE). This algorithm is a search heuristic inspired by natural evolution that is used for optimization problems [56] where a given number of possible solutions (populations) to a given problem exists. The values generating the best solutions for the optimization problem (this study: lowest RMSE) undergo “mutations” and “recombinations” to produce a new population to be tested. This procedure is continued for a chosen number of iterations. For the calculation of the genetic algorithm, we used program R 4.1.2 [32], with the package genalg [57].
The genetic algorithm was run 30 times, each with a population size of 500 and 2000 iterations. Validation of the stability of the RMSE were performed depending on the number of iterations and runs (see Supplementary Material Figure S1). The model calibration with the lowest RMSE was chosen as the final model, and the model’s fit was assessed using the test data. To compare variations in estimates between the 30 runs, we calculated the relative standard deviation (RSD).

2.2.4. Estimating the Effect of Mosquito Control Strategies

Mosquito control strategies could either be aimed to reduce the birth rate (bE) or increase the mortality of the different stages (mE, mL, mA, mH). After calibrating the vital rate functions, we decreased the birth rate and increased the mortality rates by 5%, 10%, 25%, 50%, and 75%, both individually and with some of these changes combined. Mortality rates greater than one were set to one. To test the effect of the modified vital rates on the population size, we compared the total number of adult mosquitoes (A) during the study period 1991–2010, as predicted by our model, with the models in which the altered rates were used. The total number of mosquitoes in the model after calibration was set to be 100%, and the population size obtained with the altered vital rates was calculated as a proportion of this size.

3. Results

The initial function parameter estimates using non-linear least-squares can be seen in Table 2. Using these initial function parameters, without any calibration, produced a seasonal pattern of mosquito abundance, which fitted the capture data passably. We found (with the initial estimates) an RMSE = 5.1814 and an rS = 0.832 for the training data, and an RMSE = 5.1840 and an rS = 0.826 for the test data.
With parameter calibration (parameters are shown in Table 2) the estimated abundance of active mosquitoes (A), as calculated by the final model, fitted the observations from the training data well (RMSE = 1.3457, rS = 0.905; Figure 2).
The results from the test data set showed an equally good fit (RMSE = 1.437, rS = 0.900). For the whole observation period (1991–2010, training and test data combined), we calculated the RSME = 1.393 and rS = 0.902.
This final model was able to reproduce inter-year variability in terms of the abundance of the Cx. pipiens/restuans mosquitoes, as well as in terms of the different yearly trends. However, it overestimated mosquito abundance in years with very low capture rates, and, thus, with potentially low abundance (1994, 2005, 2007), while underestimating the abundance for some years with very high mosquito capture rates (2000, 2009).
To investigate the importance of including precipitation into the model, we set the parameters associated with precipitation to zero (KL(b) and bE(d)), while using all other parameters from the final model. This model without precipitation still resulted in a high correlation with the capture data from the test data set (rS = 0.899), but with a lower RMSE = 1.898 (Figure 3).
The estimates for the calibrated model parameters, from the 30 runs conducted using the genetic algorithm, were mostly located within a small section of their possible range (Table 2). The highest relative variation was found for parameter a of the function for h, which modulates the horizontal position of the inflection point. Parameters a and b of the KL function, which we had to choose to our best knowledge, also showed a higher variation. Surprisingly, parameter k, which is part of the bE and KL function and gives the sum of the days with precipitation, showed only a small variation, despite the large possible range we allowed for with this parameter, due to missing preliminary information.
The resulting functions for the birth rates used in the final model showed that temperature developmental time for adults (bA) showed little deviation from the literature (Figure 4), whereas temperature developmental time was lower for larvae (bL) and higher for eggs (bE) compared to the initial functions obtained from literature. For the mortality rates, the estimated function parameters for mL were nearly identical to the non-calibrated function parameters, while the mortality rates for the eggs were lower, allowing the eggs to survive at temperatures below 10°; therefore, the optimal temperature for egg survival was shifted to lower temperatures.
The highest modification was found for h, the transition rate to and from the diapausing stage. In the final model, the function for h implied that under short day conditions, most mosquitoes were in the diapausing stage, which was nearly independent from the temperature. At long daytime lengths, on the other hand, temperature had a high impact on the transition between the stages A and H.
When compared to values from literature, the function of the diapausing rate (h) was skewed to the right, suggesting that, even at higher temperatures, some mosquitoes are inactive. To reduce the population of Cx. pipiens/restuans, the vital rates would have to be influenced considerably (Table 3). The most effective modification is a simultaneous increase of mL and mA, followed by other combinations, including an increase of mL. The model simulation showed that the population size is rather robust against decreasing bE or an increasing mE.

4. Discussion

The dynamic compartmental model presented in this study is able to reproduce the weekly mean Cx. pipiens/restuans abundance for Cook County at a high accuracy and over a long time period of 20 years. It only needs a few, easily accessible input parameters, namely, temperature, precipitation, and daytime length (as a function of the geographic latitude). Due to the averaging of values over several trap locations, we were able to build a model that was independent from trap-site-specific effects on mosquito abundance, such as small water bodies. Therefore, our model can be an important tool to predict Cx. pipiens/restuans abundance, and, thus, the spread of diseases transmitted by this vector species. In a previous study, we used lagged weather data and a Poisson regression model to generate a predictive model for Cx. pipiens/restuans abundance [22]. The performance of the Poisson model is directly comparable to the compartmental model presented in this study, as the same mosquito capture data were used to generate and to evaluate the models. The fit of both model types is similar (the Poisson model had a slightly poorer performance, with an RMSE = 1.752 for the test data set), but the Poisson model needed two additional input parameters, namely, relative humidity and wind speed. A further benefit is that this compartmental model can easily be included in SIR models for pathogens transmitted by Cx. pipiens/restuans.
The final model in this study showed that the highest difference in vital rates from literature values was found for the transition between the active and diapausing adult stage. This was not surprising, as we had only scarce data available to generate our initial function, which also did not give the proportion of diapausing females per se, but instead gave proportion of the females with arrested ovarian development. Further, studies conducted under laboratory conditions usually work with constant parameters (e.g., constant temperature, light or dark), while field conditions are characterized by changing conditions, such as temperature variations during the day or dim light conditions at dawn and dusk. This variation in the conditions is likely to modify the temperature and photoperiod dependence observed in the laboratory. This does not only apply for the diapausing rate, but likely for all included vital rates [58], and could explain other, not so strongly pronounced, deviations in the vital rates when compared to literature values.
Including precipitation in our model (for estimating the egg-laying rate bE and the carrying capacity KL) considerably contributed to the fit of the final model to the observed mosquito abundances. Interestingly, although our settings would have allowed us to include accumulated rainfall from up to 20 days prior, only rainfall during the previous day was selected through the genetic algorithm. Both Cx. pipiens and Cx. restuans use a wide variety of water bodies as breeding sites, ranging from larger, semi-permanent pools to small (artificial) containers, such as tin cans [59]. As small, temporal water bodies are most affected by the most recent rainfall events, our study emphasizes the importance of those temporal water bodies as breeding sites for Cx. pipiens/restuans. As with the other environmental factors, the variation in the distribution of rainfall likely modulates the effect of the total sum of precipitation during a certain period. For example, as shown by Valdez et al. [60], a high variability increases the production of mosquitoes at a low rainfall regime, while the opposite behavior takes place in the case of high precipitation.
Importantly, the investigated mosquito population has been subject to mosquito control measures. However, during the whole study period, the control methods have been constantly changing, and so have their efficiency. Control measurements have been mainly targeting the larval stage of mosquito development, but since 2002, there has been also occasional spraying of an adulticide (P. Geery, pers. comm.). Aberrations in our model’s results from the observation data could partly be caused by yearly differences in the use or efficiency of the control strategies. We also need to emphasize that our model represents a “managed” Cx. pipiens/restuans population. This needs to be taken into account when applying the model to other areas and regions. However, in many areas (at least in the USA) where mosquito-borne diseases such as WNV occur, there is some kind of mosquito control, thus, the model might be highly applicable for such areas. By manipulating the vital rates in our model, we could show that additional control measures could result in a further reduction in mosquito abundance when compared to the use of the control measures that have been already applied in Cook County during 1991–2010. As our results showed, the resulting effects on the population size are nearly linear, thus, our results would apply to unmanaged Cx. pipiens/restuans populations as well. Mosquito management strategies typically consist of a combination of applying larvicide and adulticide, as well as efforts to reduce breeding sites, but the efficiency of these methods varies [30,61]. The results of our study suggest that focus should be on increasing the mortality rates of larvae and adults. However, to obtain a sustainable effect, the mortality rates would have to be strongly increased.
The Cx. pipiens/restuans abundance model by Yu et al. [29] resembles the model created in this study, however, it did not include precipitation, which in our model increased model fit. As with our model, their abundance model showed some restraints regarding predictive power in some years, but no RMSE or correlation coefficients were given, so we cannot directly compare our model to their work. However, as Yu et al. [29] suggest, discrepancies in observed and estimated data could be due to skewness in the surveillance data, caused by some of the traps capturing a disproportionately large number of mosquitoes relative to other traps included in the surveillance program. In our study, we aimed to reduce the bias caused by inhomogeneous traps, by altering the selection process of the data used [22], however, we cannot avoid this bias completely. Additionally, poor predictive power in our model for some years could be due to differences in mosquito control measures as described above, or a combination of unusual weather phenomena, resulting in unusually low or high mosquito numbers. Whereas control measures can be accounted for in models, given detailed information, unusual weather phenomena can be very hard to predict and capture in a model.
Studies evaluating mosquito population models with capturing face the problem that the capture data does not represent the whole (female) mosquito population size, as mosquito traps usually target host-seeking females. Further, trap types differ in their efficiency, which has to be considered when comparing model output with capture data obtained from other trap types. A further challenge will be to test this model in different areas and regions [62,63,64]. As the geographic latitude is already included in this model, it should be rather straight forward to apply this model to other areas. Further, we intend to adapt this model for other Culex species, especially for the other important WNV vectors, Cx. tarsalis in the west and Cx. quinquefasciatus in the south of the United States. Further development of this model by implementing the effect of varying environmental conditions would likely increase the predictive power.

5. Conclusions

The presented dynamic compartmental model was able to reproduce the weekly mean Cx. pipiens/restuans abundance for Cook County at high precision. As it only needs a few input parameters, namely, temperature, precipitation, and latitude, it can be adapted to other regions and easily implemented into SIR models for pathogens transmitted by Cx. pipiens/restuans mosquitoes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/insects14030293/s1, Figure S1: Validation of the stability of root-mean-square error (RMSE) depending on the number of iterations and runs.

Author Contributions

Conceptualization: K.B.-L.; methodology: K.B.-L. and B.C.; software: K.B.-L. and B.C.; validation: L.J.K.; formal analysis: K.B.-L.; investigation: K.B.-L., L.J.K. and B.C.; data curation: K.B.-L.; writing—original draft preparation: K.B.-L. and B.C.; writing—review and editing: K.B.-L., L.J.K. and B.C.; visualization: K.B.-L. and B.C.; project administration: K.B.-L. and B.C. All authors have read and agreed to the published version of the manuscript.

Funding

The open access fees were covered by the Section Animal Welfare and Disease Control, Department of Veterinary and Animal Sciences, University of Copenhagen.

Data Availability Statement

The data presented in this study are available on reasonable request from the corresponding author.

Acknowledgments

We are grateful to P. Geery from the Desplaines Valley Mosquito Abatement District for providing the mosquito capture data. We thank Dániel Kondor from the Complexity Science Hub Vienna for his support running the code on a high-performance computing server.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Calisher, C.H. Medically important arboviruses of the United States and Canada. Clin. Microbiol. Rev. 1994, 7, 89–116. [Google Scholar] [CrossRef] [PubMed]
  2. Dobler, G.; Aspöck, H. Krank Durch Arthropoden. Denisia 30., Ch. Durch Stechmücken übertragene Arboviren als Erreger von Infektionskrakheiten des Menschen; Land Oberosterreich, Biologiezentrum/Oberösterreichische Landesmuseen: Linz, Austria, 2010; pp. 501–553. [Google Scholar]
  3. Farajollahi, A.; Fonseca, D.M.; Kramer, L.D.; Marm Kilpatrick, A. ”Bird biting” mosquitoes and human disease: A review of the role of Culex pipiens complex mosquitoes in epidemiology. Infect. Genet. Evol. 2011, 11, 1577–1585. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Kramer, L.; Ebel, G. Dynamics of flavivirus infection in mosquitoes. Adv. Virus Res. 2003, 60, 187–232. [Google Scholar] [PubMed]
  5. Hayes, E.B.; Komar, N.; Nasci, R.S.; Montgomery, S.P.; O’Leary, D.R.; Campbell, G.L. Epidemiology and transmission dynamics of West Nile virus disease. Emerg. Infect. Dis. 2005, 11, 1167–1173. [Google Scholar] [CrossRef]
  6. Barr, A.R. The distribution of Culex p. pipiens and C. p. quinquefasciatus in North America. Am. J. Trop. Med. 1957, 6, 153–165. [Google Scholar] [CrossRef]
  7. Darsie, R.F.; Ward, R.A. Identification and Geographical Distribution of the Mosquitoes of North America, North of Mexico; University Press of Florida: Gainesville, FL, USA, 2004. [Google Scholar]
  8. Bolling, B.G.; Kennedy, J.H.; Zimmerman, E.G. Seasonal dynamics of four potential West Nile vector species in north-central texas. J. Vector Ecol. 2005, 30, 186–194. [Google Scholar]
  9. Chuang, T.W.; Hildreth, M.B.; Vanroekel, D.L.; Wimberly, M.C. Weather and land cover influences on mosquito populations in Sioux Falls, South Dakota. J. Med. Entomol. 2011, 48, 669–679. [Google Scholar] [CrossRef] [Green Version]
  10. Wang, J.; Ogden, N.H.; Zhu, H. The impact of weather conditions on Culex pipiens and Culex restuans (Diptera: Culicidae) Abundance: A case study in Peel Region. J. Med. Entomol. 2011, 48, 468–475. [Google Scholar] [CrossRef]
  11. Madder, D.J.; Surgeoner, G.A.; Helson, B.V. Number of generations, egg production, and developmental time of Culex pipiens and Culex restauns (Diptera: Culicidae) in southern Ontario. J. Med. Entomol. 1983, 20, 275–287. [Google Scholar] [CrossRef]
  12. Reisen, W.K. Effect of Temperature on Culex tarsalis (Diptera: Culicidae) from the Coachella and San Joaquin Valleys of California. J. Med. Entomol. 1995, 32, 636–645. [Google Scholar] [CrossRef] [PubMed]
  13. Rueda, L.M.; Patel, K.J.; Axtell, R.C.; Stinner, R.E. Temperature-dependent development and survival rates of Culex quinquefasciatus and Aedes aegypti (Diptera: Culicidae). J. Med. Entomol. 1990, 27, 892–898. [Google Scholar] [CrossRef]
  14. Suman, D.S.; Tikar, S.N.; Mendki, M.J.; Sukumaran, D.; Agrawal, O.P.; Parashar, B.D.; Prakash, S. Variations in life tables of geographically isolated strains of the mosquito Culex quinquefasciatus. Med. Vet. Entomol. 2011, 525, 276–288. [Google Scholar] [CrossRef] [PubMed]
  15. Vinogradova, E.B. Mosquitoes Culex pipiens Pipiens: Taxonomy, Distribution, Ecology, Physiology, Genetics and Control; PenSoft: Sofia, Bulgaria, 2000; 280p. [Google Scholar]
  16. Spielman, A.; Wong, J. Environmental control of ovarian diapause in Culex pipiens. Ann. Entomol. Soc. Am. 1973, 66, 905–907. [Google Scholar] [CrossRef]
  17. Wilton, D.P.; Smith, G.C. Ovarian diapause in three geographic strains of Culex pipiens (Diptera: Culicidae). J. Med. Entomol. 1985, 22, 524–528. [Google Scholar] [CrossRef]
  18. Pecoraro, H.L.; Day, H.L.; Reineke, R.; Stevens, N.; Withey, J.C.; Marzluff, J.M.; Meschke, J.S. Climatic and landscape correlates for potential West Nile virus mosquito vectors in the Seattle region. J. Vector Ecol. 2007, 32, 22–28. [Google Scholar] [CrossRef]
  19. Trawinski, P.R.; MacKay, D.S. Identification of environmental covariates of West Nile virus vector mosquito population abundance. Vector Borne Zoonotic Dis. 2010, 10, 515–526. [Google Scholar] [CrossRef]
  20. Diuk-Wasser, M.A.; Brown, H.D.; Andreadis, T.G.; Fish, D. Modeling the spatial distribution of mosquito vectors West Nile virus in Connecticut, USA. Vector Borne Zoonotic Dis. 2006, 6, 283–295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Trawinski, P.R.; MacKay, D.S. Meteorologically conditioned time-series predictions of West Nile virus vector mosquitoes. Vector Borne Zoonotic Dis. 2008, 8, 505–521. [Google Scholar] [CrossRef] [Green Version]
  22. Lebl, K.; Brugger, K.; Rubel, F. Predicting Culex pipiens/restuans population dynamics by interval lagged weather data. Parasites Vectors 2013, 6, 129. [Google Scholar] [CrossRef] [Green Version]
  23. Ahumada, J.A.; Lapointe, D.; Samuel, M.D. Modeling the population dynamics of Culex quinquefasciatus (Diptera: Culicidae), along an elevational gradient in Hawaii. J. Med. Entomol. 2004, 41, 1157–1170. [Google Scholar] [CrossRef] [Green Version]
  24. Eisenberg, J.N.; Reisen, W.; Spear, R.C. Dynamic Model Comparing the Binoimics of Two Isolated Culex tarsalis Populations: Model Development. J. Med. Entomol. 1995, 32, 83–97. [Google Scholar] [CrossRef] [PubMed]
  25. Morin, C.W.; Comrie, A.C. Modeled response of the West Nile virus vector Culex quinquefasciatus to changing climate using the dynamic mosquito simulation model. Int. J. Biometeorol. 2010, 54, 517–529. [Google Scholar] [CrossRef] [PubMed]
  26. Laperriere, V.; Brugger, K.; Rubel, F. Simulation of the seasonal cycles of bird, equine and human West Nile virus cases. Prev. Vet. Med. 2011, 98, 99–110. [Google Scholar] [CrossRef] [PubMed]
  27. Rubel, F.; Brugger, K.; Hantel, M.; Chvala-Mannsberger, S.; Bakonyi, T.; Weissenböck, H.; Nowotny, N. Explaining Usutu virus dynamics in Austria: Model development and calibration. Prev. Vet. Med. 2008, 85, 166–186. [Google Scholar] [CrossRef]
  28. Wonham, M.J.; de Camino-Beck, T.; Lewis, M.A. An epidemiological model for West Nile virus: Invasion analysis and control applications. Proc. Royal Soc. B 2004, 271, 501–507. [Google Scholar] [CrossRef] [Green Version]
  29. Yu, D.; Madras, N.; Zhu, H. Temperature-driven population abundance model for Culex pipiens and Culex restuans (Diptera: Culicidae). J. Theor. Biol. 2018, 443, 28–38. [Google Scholar] [CrossRef] [PubMed]
  30. Ballenger-Browning, K.; Elder, J. Multi-modal aedes aegypti mosquito reduction interventions and dengue fever prevention. Trop. Med. Int. Health. 2009, 14, 1542–1551. [Google Scholar] [CrossRef] [PubMed]
  31. Alphey, L.; Benedict, M.; Bellini, R.; Clark, G.; Dame, D.; Service, M.; Dobson, S. Sterile-insect methods for control of mosquito-borne diseases: An analysis. Vector Borne Zoonotic Dis. 2010, 10, 295–311. [Google Scholar] [CrossRef]
  32. R Development Core Team. R: A Language and Environment for Statistical Computing, Version 4.1.2; R Foundation for Statistical Computing: Vienna, Austria, 2022.
  33. Elizondo-Quiroga, A.; Flores-Suarez, A.; Elizondo-Quiroga, D.; Ponce-Garcia, G.; Blitvich, B.J.; Contreras-Cordero, J.F.; Gonzalez-Rojas, J.I.; Mercado-Hernandez, R.; Beaty, B.J.; Fernandez-Salas, I. Gonotrophic cycle and survivorship of Culex quinquefasciatus (Diptera: Culicidae) using sticky ovitraps in Monterrey, northeastern Mexico. J. Am. Mosq. Control Assoc. 2006, 22, 10–14. [Google Scholar]
  34. Reisen, W.K.; Milby, M.M. Studies on autogeny in Culex tarsalis: 3. Life table attributes of autogenous and anautogenous strains under laboratory conditions. J. Am. Mosq. Control Assoc. 1987, 3, 619–625. [Google Scholar]
  35. Oda, T.; Ueda, M. Seasonal changes in the oviposition activity of Culex pipiens molestus Forskal in the above-ground water. Trop. Med. 1979, 21, 139–144. [Google Scholar]
  36. Oda, T.; Mori, A.; Ueda, M.; Kurokawa, K. Effects of temperatures on the oviposition and hatching of eggs in Culex pipiens molestus and Culex pipiens quinquefasciatus. Trop. Med. 1980, 22, 167–172. [Google Scholar]
  37. Shriver, D.; Bickley, W.E. The effect of temperature on hatching eggs of the mosquito, culex pipiens quinquefasciatus say. Mosquito News 1964, 24, 137–140. [Google Scholar]
  38. Almir’on, W.R.; Brewer, M.E. Winter biology of Culex pipiens quinquefasciatus Say, (Diptera: Culicidae) from Cord´oba, Argentina. Mem. Do Inst. Oswaldo Cruz 1996, 91, 649–654. [Google Scholar] [CrossRef] [Green Version]
  39. Kirkpatrick, T.W. The Mosquitoes of Egypt; Government Press: Cairo, Egypt, 1925; p. 224.
  40. Kramer, S.D. The effect of temperature on the life cycle of Musca domestica and Culex pipiens. Science 1915, 41, 874–877. [Google Scholar] [CrossRef] [PubMed]
  41. Olejníček, J.; Gelbič, I. Differences in response to temperature and density between two strains of the mosquito, Culex pipiens molestus Forskal. J. Vector Ecol. 2000, 25, 136–145. [Google Scholar]
  42. Tamarina, N.A. The developmental cycle of culex pipiens molestus forsk. Vestn. MGU (Biol.) 1966, 2, 65–75. [Google Scholar]
  43. Bailey, S.F.; Gieke, P.A. A study of the effect of water temperatures on rice field mosquito development. In Proceedings and Papers of the Annual Conference of the California Mosquito Control Association; CMCA Press: Rockland, ME, USA, 1968; pp. 53–61. [Google Scholar]
  44. Mogi, M. Temperature and photoperiod effects on larval and ovarian development of New Zealand strains of Culex quinquefasciatus (Diptera: Culicidae). Ann. Entomol. Soc. Am. 1992, 85, 58–66. [Google Scholar] [CrossRef]
  45. Mori, A.; Oda, T.; Zaitsu, M.; Ueda, M.; Kurokawa, K. Studies on the developing period of larval stage of the Culex pipiens complex in Japan. Trop. Med. 1988, 30, 155–161. [Google Scholar]
  46. Shelton, R.M. The effect of temperatures on development of eight mosquito species. Mosq. News 1973, 33, 1–12. [Google Scholar]
  47. Reisen, W.K.; Meyer, R.P.; Shields, J.; Arbolante, C. Population ecology of preimaginal Culex tarsalis (Diptera: Culicidae) in Kern County, California. J. Med. Entomol. 1989, 26, 10–22. [Google Scholar] [CrossRef] [PubMed]
  48. Buth, J.L.; Brust RAEllis, R.A. Development time, oviposition activity and onset of diapause in Culex tarsalis, Culex restuans and Culiseta inornata in southern Manitoba. J. Am. Mosq. Control Assoc. 1990, 6, 55–63. [Google Scholar] [PubMed]
  49. Kalchenko, E.I. The biology of Culex pipiens molestus Forsk. (Diptera, Culicidae). Entomol. Obozr. 1962, 41, 94–98. [Google Scholar]
  50. Gunay, F.; Alten, B.; Ozsoy, E.D. Estimating reaction norms for predictive population parameters, age specific mortality, and mean longevity in temperature-dependent cohorts of Culex quinquefasciatus Say (Diptera: Culicidae). J. Vector Ecol. 2010, 35, 354–362. [Google Scholar] [CrossRef] [PubMed]
  51. Oda, T.; Eshita, Y.; Uchida, K.; Mine, M.; Kurokawa, K.; Ogawa, Y.; Tahara, H.K. Reproductive activity and survival of Culex pipiens pallens and Culex quinquefasciatus (Diptera: Culicidae) in Japan at high temperature. J. Med. Entomol. 2002, 39, 185–190. [Google Scholar] [CrossRef]
  52. Eldridge, B.F. Environmental Control of Ovarian Development in Mosquitoes of the Culex pipiens complex. Science 1966, 151, 826–828. [Google Scholar] [CrossRef] [PubMed]
  53. Bailey, C.M.; Faran, T.; Gargan, I.I.; Hayes, D. Winter survival of blood-fed and nonblood-fed Culex pipiens L. Am. J. Trop. Med. 1982, 31, 1054–1061. [Google Scholar] [CrossRef] [PubMed]
  54. NCDC, 2012: NNCD Climate Data Online. Available online: https://www.ncei.noaa.gov/cdo-web/ (accessed on 1 December 2022).
  55. Forsythe, W.C.; Rykiel, E.J.; Stahl, R.S.; Wu, H.; Schoolfield, R.M. A model comparison for daylength as a function of latitude and day of year. Ecol Modell. 1995, 80, 87–95. [Google Scholar] [CrossRef]
  56. Katoch, S.; Chauhan, S.S.; Kumar, V. A review on genetic algorithm: Past, present, and future. Multimed. Tools Appl. 2021, 80, 8091–8126. [Google Scholar] [CrossRef]
  57. Willighagen, E. Genalg: R Based Genetic Algorithm; version 0.1.1; 2005. Available online: http://cran.rproject.org/web/packages/genalg/index.html (accessed on 27 February 2023).
  58. Erguler, K.; Mendel, J.; Petrić, D.V.; Petrić, M.; Kavran, M.; Demirok, M.C.; Gunay, F.; Georgiades, P.; Alten, B.; Lelieveld, J. A dynamically structured matrix population model for insect life histories observed under variable environmental conditions. Sci. Rep. 2022, 12, 11587. [Google Scholar] [CrossRef]
  59. Becker, N.; Petrić, D.; Zgomba, M. Mosquitoes-Identification, Ecology and Control, 3rd ed.; Springer: Cham, Switzerland, 2020; p. 570. [Google Scholar]
  60. Valdez, L.; Sibona, G.; Díaz LContigiani, M.; Condat, C. Effects of rainfall on Culex mosquito population dynamics. J. Theor. Biol. 2017, 421, 421. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Kramer, L.; Styer, L.M.; Ebel, G.A. Global perspective on the epidemiology of the west nile virus. Annu. Rev. Entomol. 2008, 53, 61–81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Pinior, B.; Lebl, K.; Firth, C.; Rubel, F.; Fuchs, R.; Stockreiter, S.; Loitsch, A.; Köfer, J. Cost analysis of bluetongue virus serotype 8 surveillance and vaccination programmes in Austria from 2005 to 2013. Vet. J. 2015, 206, 154–160. [Google Scholar] [CrossRef] [Green Version]
  63. Pinior, B.; Brugger, K.; Köfer, J.; Schwermer, H.; Stockreiter, S.; Loitsch, A.; Rubel, F. Economic comparison of the monitoring programmes for bluetongue vectors in Austria and Switzerland. Vet. Rec. 2015, 176, 464. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Zittra, C.; Flechl, E.; Kothmayer, M.; Vitecek, S.; Rossiter, H.; Zechmeister, T.; Fuehrer, H.-P. Ecological characterization and molecular differentiation of Culex pipiens complex taxa and Culex torrentium in eastern Austria. Parasites Vectors 2016, 9, 197. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Block diagram of the Culex population model, incorporating 4 stages: E—eggs, L—larvae (inclusive pupal sage), A—active adults; H—hibernating adults. For definition of the parameters, see Table 1.
Figure 1. Block diagram of the Culex population model, incorporating 4 stages: E—eggs, L—larvae (inclusive pupal sage), A—active adults; H—hibernating adults. For definition of the parameters, see Table 1.
Insects 14 00293 g001
Figure 2. Standardized weekly averaged observation data and final model results. Years written in regular type indicate the years for the training data set; years in bold type indicate the test data set.
Figure 2. Standardized weekly averaged observation data and final model results. Years written in regular type indicate the years for the training data set; years in bold type indicate the test data set.
Insects 14 00293 g002
Figure 3. The estimated abundances of active mosquitoes, using the capture data from the test data set, without including precipitation.
Figure 3. The estimated abundances of active mosquitoes, using the capture data from the test data set, without including precipitation.
Insects 14 00293 g003
Figure 4. Functions used to describe the transition parameters (see Table 1). Dashed lines represent the estimated functions fitted to literature values (mL multiplied by a factor of 10), and solid lines represent the functions after calibration with the training data set.
Figure 4. Functions used to describe the transition parameters (see Table 1). Dashed lines represent the estimated functions fitted to literature values (mL multiplied by a factor of 10), and solid lines represent the functions after calibration with the training data set.
Insects 14 00293 g004
Table 1. Definition of the transition parameters used in the Culex population model. Model parameters are per capita rates in units days−1.
Table 1. Definition of the transition parameters used in the Culex population model. Model parameters are per capita rates in units days−1.
Parameter 1Description
bE (T, P)birth rate eggs—reciprocal value of the length of the gonotrophic cycle multiplied with the number of female eggs laid per egg raft
mE (T)mortality rate of eggs
bL (T)birth rate larvae—reciprocal value of the length of the egg stag
mL (T, KL)mortality rate of larvae and pupae combined
KL (P)carrying capacity of the larvae
bA (T)birth rate adults—reciprocal value of the length of the larval and pupal stage combined
mA (T)mortality rate of active adult mosquitoes
h (T, D)diapausing rate—rate of adult mosquito leaving the active stage and entering the diapause stage
mHmortality rate of adult mosquitoes during winter diapause
1 Information on parameters, depending on environmental factors, is given in parentheses: T—daily mean ambient temperature (°C), P—daily sum of precipitation (mm) and D—daytime length (h).
Table 2. Parameter estimates for Equations (6)–(14) for the Cx. pipiens/restuans mosquito model. The estimates are the values for a specific function fitted to literature data, and the according lower and upper margin (estimate ± SE). Values estimated to our best knowledge, i.e., not based on literature data, are indicated by an asterisk (*). Furthermore, the model calibration with the lowest RMSE (in bold), as well as the minimum, maximum, mean, standard deviation (SD) and relative standard deviation (RSD), obtained by the 30 runs of the genetic algorithm (training data set), are shown.
Table 2. Parameter estimates for Equations (6)–(14) for the Cx. pipiens/restuans mosquito model. The estimates are the values for a specific function fitted to literature data, and the according lower and upper margin (estimate ± SE). Values estimated to our best knowledge, i.e., not based on literature data, are indicated by an asterisk (*). Furthermore, the model calibration with the lowest RMSE (in bold), as well as the minimum, maximum, mean, standard deviation (SD) and relative standard deviation (RSD), obtained by the 30 runs of the genetic algorithm (training data set), are shown.
Func.Par.EstimateLowerUpperCalibMinMaxMeanSDRSD
bEa3.9 × 10−13.4 × 10−14.5 × 10−14.5 × 10−14.2 × 10−14.5 × 10−14.4 × 10−17.9 × 10−31.8 × 10−2
b9.5 × 1012.9 × 1011.6 × 1023.0 × 1012.9 × 1013.2 × 1013.0 × 1018.6 × 10−12.9 × 10−2
c−2.1 × 10−1−2.6 × 10−1−1.7 × 10−1−2.6 × 10−1−2.6 × 10−1−2.4 × 10−1−2.5 × 10−15.3 × 10−3−2.1 × 10−2
d *0.0 × 1000.0 × 1001.0 × 10−11.0 × 10−15.0 × 10−21.0 × 10−19.7 × 10−29.1 × 10−39.4 × 10−2
bE, KLk *0.0 × 1000.0 × 1002.0 × 1011.1 × 1008.4 × 10−11.8 × 1001.2 × 1002.5 × 10−12.1 × 10−1
mEa2.5 × 1002.0 × 1002.9 × 1002.0 × 1002.0 × 1002.5 × 1002.1 × 1001.1 × 10−15.4 × 10−2
b−2.0 × 10−1−2.4 × 10−1−1.6 × 10−1−2.0 × 10−1−2.2 × 10−1−1.8 × 10−1−2.0 × 10−11.1 × 10−2−5.6 × 10−2
c4.2 × 10−33.3 × 10−35.0 × 10−35.0 × 10−34.1 × 10−35.0 × 10−34.9 × 10−32.4 × 10−45.0 × 10−2
bLa2.7 × 1011.4 × 1014.0 × 1011.4 × 1011.3 × 1013.3 × 1011.6 × 1014.2 × 1002.7 × 10−1
b−1.5 × 10−1−1.7 × 10−1−1.3 × 10−1−7.0 × 10−2−1.7 × 10−1−5.9 × 10−2−1.2 × 10−13.2 × 10−2−2.7 × 10−1
mLa1.9 × 101−3.6 × 1004.1 × 1013.1 × 1012.5 × 1013.3 × 1012.8 × 1012.3 × 1008.1 × 10−2
b−2.9 × 100−3.8 × 100−2.0 × 100−3.7 × 100−3.8 × 100−3.7 × 100−3.7 × 1004.1 × 10−2−1.1 × 10−2
KLa *1.0 × 1001.0 × 1005.0 × 1011.0 × 1001.0 × 1001.8 × 1001.2 × 1002.2 × 10−11.8 × 10−1
b *1.0 × 1001.0 × 1007.0 × 1013.6 × 1013.6 × 1013.7 × 1013.6 × 1015.7 × 10−11.6 × 10−2
bAa1.1 × 1027.3 × 1011.5 × 1021.0 × 1027.5 × 1011.5 × 1021.1 × 1022.2 × 1012.1 × 10−1
b−8.5 × 10−2−9.8 × 10−2−7.3 × 10−2−7.3 × 10−2−7.4 × 10−2−7.3 × 10−2−7.3 × 10−23.9 × 10−4−5.2 × 10−3
mAa9.2 × 10−26.4 × 10−21.2 × 10−16.4 × 10−26.4 × 10−28.5 × 10−26.8 × 10−25.0 × 10−37.3 × 10−2
b−7.5 × 10−3−1.0 × 10−2−5.0 × 10−3−5.1 × 10−3−5.3 × 10−3−4.9 × 10−3−5.0 × 10−31.0 × 10−4−2.1 × 10−2
c2.0 × 10−41.5 × 10−42.6 × 10−42.6 × 10−42.6 × 10−42.6 × 10−42.6 × 10−48.7 × 10−73.4 × 10−3
mHa1.8 × 10−21.7 × 10−22.0 × 10−21.7 × 10−21.7 × 10−22.0 × 10−21.8 × 10−28.1 × 10−44.5 × 10−2
ha8.4 × 10−6−6.7 × 10−62.4 × 10−58.8 × 10−88.0 × 10−83.2 × 10−71.2 × 10−75.6 × 10−84.6 × 10−1
hb2.8 × 10−12.3 × 10−13.3 × 10−12.3 × 10−12.3 × 10−12.4 × 10−12.4 × 10−11.9 × 10−38.0 × 10−3
hc7.5 × 10−16.2 × 10−18.7 × 10−18.7 × 10−17.8 × 10−18.7 × 10−18.5 × 10−12.4 × 10−22.9 × 10−2
Table 3. Effect of mosquito control strategies, increasing/decreasing different vital rates by a certain percentage, on the relative population size. The percentage of active mosquitoes during the study period (1991–2010) was obtained by using the altered rate in relation to the number of mosquitoes predicted by the model after calibration.
Table 3. Effect of mosquito control strategies, increasing/decreasing different vital rates by a certain percentage, on the relative population size. The percentage of active mosquitoes during the study period (1991–2010) was obtained by using the altered rate in relation to the number of mosquitoes predicted by the model after calibration.
5%10%25%50%75%
decrease bE99.699.297.894.286.5
increase mE99.799.598.697.496.3
increase mL96.292.682.668.456.5
increase mA96.292.783.671.963.1
increase mH99.198.395.992.389.1
decrease bE + increase mE99.398.796.491.482.2
decrease bE + increase mL95.891.880.764.248.0
decrease bE + increase mA95.892.081.667.052.6
decrease bE + increase mH98.797.593.685.973.8
increase mE + increase mL95.992.181.566.554.2
increase mE + increase mA95.992.282.469.860.4
increase mE + increase mH98.897.794.689.885.6
increase mL + increase mA92.685.869.149.335.8
increase mL + increase mH95.391.079.162.949.7
increase mA + increase mH95.491.180.466.856.8
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bakran-Lebl, K.; Kjær, L.J.; Conrady, B. Predicting Culex pipiens/restuans Population Dynamics Using a Weather-Driven Dynamic Compartmental Population Model. Insects 2023, 14, 293. https://doi.org/10.3390/insects14030293

AMA Style

Bakran-Lebl K, Kjær LJ, Conrady B. Predicting Culex pipiens/restuans Population Dynamics Using a Weather-Driven Dynamic Compartmental Population Model. Insects. 2023; 14(3):293. https://doi.org/10.3390/insects14030293

Chicago/Turabian Style

Bakran-Lebl, Karin, Lene Jung Kjær, and Beate Conrady. 2023. "Predicting Culex pipiens/restuans Population Dynamics Using a Weather-Driven Dynamic Compartmental Population Model" Insects 14, no. 3: 293. https://doi.org/10.3390/insects14030293

APA Style

Bakran-Lebl, K., Kjær, L. J., & Conrady, B. (2023). Predicting Culex pipiens/restuans Population Dynamics Using a Weather-Driven Dynamic Compartmental Population Model. Insects, 14(3), 293. https://doi.org/10.3390/insects14030293

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop