 Research
 Open Access
 Published:
A climatedriven mechanistic population model of Aedes albopictus with diapause
Parasites & Vectors volume 9, Article number: 175 (2016)
Abstract
Background
The mosquito Aedes albopitus is a competent vector for the transmission of many bloodborne pathogens. An important factor that affects the mosquitoes’ development and spreading is climate, such as temperature, precipitation and photoperiod. Existing climatedriven mechanistic models overlook the seasonal pattern of diapause, referred to as the survival strategy of mosquito eggs being dormant and unable to hatch under extreme weather. With respect to diapause, several issues remain unaddressed, including identifying the time when diapause eggs are laid and hatched under different climatic conditions, demarcating the thresholds of diapause and nondiapause periods, and considering the mortality rate of diapause eggs.
Methods
Here we propose a generic climatedriven mechanistic population model of Ae. albopitus applicable to most Ae. albopictuscolonized areas. The new model is an improvement over the previous work by incorporating the diapause behaviors with many modifications to the stagespecific mechanism of the mosquitoes’ lifecycle. monthly Container Index (CI) of Ae. albopitus collected in two Chinese cities, Guangzhou and Shanghai is used for model validation.
Results
The simulation results by the proposed model is validated with entomological field data by the Pearson correlation coefficient r ^{2} in Guangzhou (r ^{2} = 0.84) and in Shanghai (r ^{2} = 0.90). In addition, by consolidating the effect of diapauserelated adjustments and temperaturerelated parameters in the model, the improvement is significant over the basic model.
Conclusions
The model highlights the importance of considering diapause in simulating Ae. albopitus population. It also corroborates that temperature and photoperiod are significant in affecting the population dynamics of the mosquito. By refining the relationship between Ae. albopitus population and climatic factors, the model serves to establish a mechanistic relation to the growth and decline of the species. Understanding this relationship in a better way will benefit studying the transmission and the spatiotemporal distribution of mosquitoborne epidemics and eventually facilitating the early warning and control of the diseases.
Background
Aedes albopitus (Skuse), also known as the Asian tiger mosquito, is a competent vector for the transmission of many bloodborne epidemics such as dengue fever, West Nile virus infections and Chikungunya fever [1–5]. Ae. albopictus is a species native to the tropical areas of Southeast Asia; and during the last century, it has rapidly invaded countries throughout the world. Now the species is pervasively found in the subtropical and temperate climate areas of East Asia, Europe, Africa, the Middle East and the Americas [1, 2]. The extensive spreading of Ae. albopictus is not only prompted by the increasing trend of international trade and travel [3–5] but is also mediated by climate change in a global context [6]. It has been corroborated that the growth of Ae. albopictus is constrained by the changing nature of physical environment; as a result, their population density is extrinsically impinged by a series of climatic factors including temperature, precipitation and photoperiod [7]. Because of this link to climate change, modeling the population dynamics of Ae. albopictus based on climatic factors serves a critical role for further identifying the causal relation to the transmission and control of mosquitoborne pathogens [8].
Within the realm of quantitative modeling, the climatedriven dynamics of the mosquito has been explored by two types of models: the statistical population model and the mechanistic population model [9]. The statistical population model aims to establish a mathematical correlation between the population abundance and climatic factors using data solicited from controlled experiments or field observations [10–12]. Although these statistical relationships are relatively straightforward and can be easily understood, they are flawed in describing the intrinsic biological mechanism of how the morphology of pathogen carriers is mediated by the environment. Another overlooked facet in the statistical population model pertains to the limited coverage of the species’ lifecycle stages; for example, the majority of mosquito models have only differentiated between the aquatic period and the aerial period, while other substages (e.g. eggs, larvae, pupae) have been less emphasized [13]. When climatic factors are involved in promoting or inhibiting the development, the influence on each substage of the lifecycle is a critical aspect that must be further scrutinized [14, 15]. To overcome these limitations, the mechanistic population model explores the fluctuation of population by factoring in a priori established development process in a context constrained by environment [8]. For the study of Ae. albopictus, this type of model and its variations have been applied to specific geographic regions and have been modified for extensions to other Aedes spp. [16–22]. One important work in this area is attributed to Erickson et al. [23], who developed a stagestructured population model consisting of six ordinary differential equations to correspond to different stages of Ae. albopictus’ lifecycle, with each equation measuring the mortality rate and growth rate dependent on a stagespecific temperature variable. Thereafter, Cailly et al. [24] established a generic temperaturedriven model that extended the application to different mosquito species. This model served a more general purpose to represent the complete lifecycle with ten model compartments and temperaturedriven mortality rate and growth rate [25, 26]. In a followup study, Tran et al. [27] optimized the parameters and transition functions in Cailly’s model to estimate the population dynamics of Ae. albopictus with improved accuracy and a better fit to field observations.
An overlooked area in the formulation of the mechanistic population model is the phenomenon of diapause. Diapause refers to the physiological mechanism within certain species that inhibits the development of organism as a strategy to survive unfavorable environmental conditions, such as extreme weather [28, 29]. Diapause is observed among Ae. albopictus in selected subtropical areas and temperate climate areas where the relatively cold season of a year makes the mosquito eggs become dormant and unable to hatch [7, 30]. The delay in development can only be remedied when the eggs are exposed to enough warmth or a long photoperiod [1, 7]. Many entomologists previously associated Ae. albopictus diapause with environmental conditions using controlled experiments. For example, Wang [31] and Imai & Maeda [32] discovered that the intervening factors of diapause included low temperatures and short photoperiods; and this conclusion was later confirmed by a series of controlled experiments [1, 7, 30, 33]. In addition, comparative studies identified that the unique mechanism of diapause greatly improved the survival rate of Ae. albopictus under extreme desiccation and coldstress conditions [34, 35]. These findings, however, have been rarely incorporated in the mechanistic population model [27]. Two exceptions are the models proposed by Cailly et al. [24] and Tran et al. [27], where diapause was quantified as a piecewise function dichotomizing the development into diapause period, say last September through early March, and that of nondiapause for the rest of year. Compared with Erickson et al. [23], accounting for diapause greatly improved the model performance and generated a close approximation of the field observation. However, the adjustment was based on a relative empirical assumption that diapause occurs in the same time period over different years as well as does not manifest regional differences. As diapause is a dynamic process primarily dictated by degrees of temperature and lengths of photoperiod [31, 32], it is more reasonable to assume diapause as a temperature and photoperioddriven phenomenon in a general way to consolidate the effect of temporal as well as regional differences.
Along the line of existing mechanistic population models [23, 24, 27], this paper, proposes a generic model to study Ae. albopictus population by considering diapause as a dynamic process conditional on the change of temperature and photoperiod. By restructuring the model and tuning the parameters, the study strictly captures the stagespecific mechanism of the mosquitoes’ lifecycle. In addition, the model has been validated by field data collected in two Chinese cities over a fiveyear span where the phenomenon of diapause was introduced by seasonality. The proposed model aims to: (i) define a finetuned relationship between major climatic variables (temperature and photoperiod) and diapauserelated events in Ae. albopictus’s development cycle; (ii) explicitly quantify the effect of temperature on the population dynamics during the aquatic period based on controlled laboratory experiments and related literature.
The paper is organized as follows. The Methods section first presents the two sets of field observations in Chinese cities, which are employed as evidence for model validation. Then the new model is proposed based on two existing mechanistic population models. The Results section demonstrates simulation results, which are further compared and validated with field observations. The Discussion section discusses the improvement of the model over the original model as well as the sensitivity of nonclimatic variables in the model. Lastly, the Conclusions section summarizes the contribution of the study and proposes directions for future research.
Methods
Study areas and data
To validate the model, we firstly conducted rigorous field experiments to collect Ae. albopictus larval samples. Field observations on the Ae. albopictus population were conducted in two southern Chinese cities, Guangzhou and Shanghai, over a respective period of five years. The two study areas both belong to the subtropical climate zone with a humid and hot summer and an arid and cold winter. Compared to Guangzhou, the weather in Shanghai exhibits a higher degree of seasonality with a colder spring, autumn and winter.
The field observations were quantified by the monthly Container Index (CI) of Ae. albopictus larvae in Guangzhou (2007–2011) and Shanghai (2009–2013) collected by the China Centre for Disease Control and Prevention (CDC), as shown in Table 1. To acquire the data, we distributed 50–100 sampling containers with clean water in outdoor conditions and then checked the larval density at the end of each sampling period. Generally, the sampling was conducted every 3–5 days in areas where cases of dengue fever were observed and every 15 days in other areas. The corresponding CI was calculated by Equation 1, where N ^{+}denotes number of containers with at least one Ae. albopictus larva and N denotes the number of total sampling containers [36]. Then we consolidated the CI values collected in different sampling areas into a monthly index. This monthly CI represents the monthly Ae. albopictus population abundance at the larval stage and serves as the field evidence for model validation.
To correspond to the CI data for each city, we collected daily mean temperature and precipitation data from China Meteorological Data Sharing Service System. We also derived photoperiod data using an existing source [37]. These datasets included Guangzhou over 2007–2011 and Shanghai over 2009–2013.
Basic mechanistic population model
Our proposed model is a natural extension of two basic mechanistic population models proposed by Cailly et al. [24] and Tran et al. [27]. These two models were formulated by considering the lifecycle of Ae. albopictus that consists of two principal periods by distinct habitats: the aquatic period and the aerial period, as shown in Fig. 1a. The two periods can be further divided into the following eight substages: four aquatic stages for growth (eggs, larvae, pupae and emergence) and four aerial stages for breeding among adults (mating, blood feeding, gestating and ovipositing) [25]. To simplify the entire lifecycle, our model merges the emergence stage and the mating stage into the stage of emerging adults, while other seven substages are retained, as shown in Fig. 1b. In addition, to further account for the effect of diapause, the abundance of Ae. albopictus eggs are divided into two groups: nondiapause eggs and diapause eggs to represent two distinct hatching behaviors.
Temperature plays an indispensable role in mediating the development and mortality rate of the Ae. albopictus throughout the lifecycle stages; and the aquatic period is significantly affected [38]. Erickson et al. [23] firstly developed a basic sixstaged population model to describe this temperaturedriven mechanism. Meanwhile, Cailly et al. [24] introduced diapause and also accounted for temperature as the major climatic trigger using a tenstaged population model. In addition, the second climatic mediator refers to precipitation, which contributes to the environmental carrying capacity (defined as the maximum population abundance the environment can sustain) at both the larval stage and the pupal stage [12, 39–41]. Based on Cailly et al. [24], Tran et al. [27] further considered both temperature and precipitation in the development of Ae. albopictus by a tenstage mechanistic model. Based on these three models [23, 24, 27], an equivalent sevenstage model is proposed as the basis of the paper, as given in Equation 2. In this equation, \( \dot{X} \) denotes a variable X taking a derivative with respect to time t in day of year. This equation denotes the daily variation of population abundance at each of the seven substages including eggs (E), larvae (L) and pupae (P) in the aquatic period as well as emerging adults (A _{em}), blood feeding adults (A _{b}), gestating adults (A _{g}) and ovipositing adults (A _{o}) in the aerial period. In general, this model specifies that the population at each substage is composed of (1) change of the population at the current stage (determined by the mortality rate and development rate) and (2) the population accumulated from the previous stage. In Equation 2, some parameters are dependent on climatic variables including temperature (T) and precipitation (P), as described in Table 2; other parameters independent of T and P are given by experiments from existing literature, as given in Table 3.
This proposed basic model is still in need of further scrutiny, in that it overlooks two important mechanistic facets. First, as noted by Cailly et al. [24, 27] and Tran et al. [27], diapause characterized only by a single binary variable (z _{dia}) appears to be flawed. When applied for a different region, this model is less effective by assuming a static diapause period, say late September through early March (z _{dia} = 0), and the rest of the year as the nondiapause period (z _{dia} = 1) [27]. As diapause is a climatedriven phenomenon that differs across geographic regions and arises on different days [7], a model that considers the temporal variation of diapause needs to be considered and validated for all lifecycle stages of Ae. albopictus. Secondly, when the effect of diapause is incorporated, the role of climatedependent parameters (Table 2) needs to be further scrutinized toward a better model performance. To address these two concerns, we propose an improved mechanistic population model that aims to adjust asterisked parameters in Table 2.
Improved mechanistic population model
In this section, we focus on improving the basic model (Equation 2) from two perspectives: model structure and model parameters. First, we integrate the confirmed relationship between diapauserelated effects and climatic variables (temperature and photoperiod); and then we restructured the model by quantifying the conditions when the diapause arises. Secondly, we adjust several model parameters, such as development rates, mortality rates and oviposition rate (see asterisked parameters in Table 2).
Model structure: relationship between diapause and climate
In our field work, the absence of larvae throughout most winters indicates the likelihood of diapause in the two Chinese cities (Table 1). For the study of Ae. albopictus, the pressing need is to quantify the climatic thresholds under which diapause occurs and how diapause contingently influences subsequent development. Based on former experimental findings [31, 42], a female chooses to lay diapause eggs when two following conditions are simultaneously fulfilled: (1) the temperature is below 21 °C and (2) the length of daylight hours is around 13h/14h. Based on this conclusion, we denote a binary variable z _{1}(t) as the state of diapause egg oviposition at day t of the year (Equation 3). This variable is dependent on the average temperature (T _{aver}) and the average daylight hours (D _{aver}) observed in the week before t. As the oviposition takes place normally in early autumn [7], here we consider that the oviposition process starts in the beginning of autumn (t _{1}, August 31 or the 243th day of year [DOY]) and ends in the day when the diapause period begins (t _{begin}).
For simplicity, here we denote the percentage of diapause eggs (r _{dia}) as E _{dia}/E, where E = E _{0} + E _{dia} (E _{0} is the number of nondiapause eggs and E _{dia} is the number of diapause eggs). According to Tran et al. [27], when a significant portion of diapause eggs are oviposited (r _{dia} > 0.9), the species steps into its diapause period at t _{begin}. The diapause continues and ends when these eggs begin to hatch. As observed in Toma et al. [43], the diapause eggs did not hatch until weekly mean temperature rose to 10 °C/11 °C and the daylight hours reached 11h/11.5h. Based on the evidence, we denote another binary z _{2}(t) to control for the climatic threshold for the hatching event, indicating the end of egg diapause at day t of year (Equation 4), where t _{end} is the day when the diapause period ends and t _{2} is the first day of year when r _{dia} (t) < 0.1.
During the diapause period (t _{begin} < t < t _{end}), on average Ae. albopictus adults succumb to a temperature of under 9.5 °C (γ _{ Aem } = γ _{ Ag } = γ _{ Ar } = 0 when T _{aver} < 9.5 °C) [43].To survive such extreme conditions, egg diapause occurs as an adaptive strategy to improve survival rate compared with nondiapause egg [7, 34, 35]. However, the mortality rate and development rate of diapause eggs are less explored by existing literature; and therefore we denote these two parameters as m _{dia} and f _{dia} respectively in Equation 5 and conduct their sensitivity analysis in Results section. In addition to Equation 4, the hatching of diapause eggs is subject to other environmental factors, such as water and food availability [32, 44, 45]. In accordance with Tran et al. [27], we consider that the hatching is also dependent on the first precipitation event in spring (P _{week} > P _{0}, where P _{0} = 0 mm).
With these improvements applied in Equation 2, the adjusted model is presented in Equation 5 to characterize the population dynamics of Ae. albopictus in a sevenstage development process, where the notations follow Tables 2 and 3. This model features the temperature and photoperiod driven mechanism in both the oviposition z _{1}(t) and the egg diapause z _{2}(t), while considering diapause egg survival and hatching during the diapause period.
Model parameters: temperaturedriven mechanism
To derive the asterisked parameters in Table 2, we have conducted a separate set of controlled experiments to investigate the response of Ae. albopictus to different temperatures in the aquatic period [46]. In these experiments, the development lengths of eggs, larvae and pupae, and the survival rate of larvae under different temperatures (16 °C, 21 °C, 26 °C, 31 °C and 36 °C) have been identified (Table 4). Unfortunately, due to limited laboratory conditions, several other records, such as the survival rate of pupae, could not be obtained. Therefore, we have derived these missing variables from existing literature as a compromise solution [7, 47, 48], as shown in Table 4.
We then describe these data in a scatter diagrams: the hatching rate (f _{E}), the larval development rate (f _{L}) and the pupal development rate (f _{P}) as the reciprocal of the development length under different temperatures (T), respectively (Figs. 2ac). Previous studies have explored a linear or a quadratic equation [27] to describe the relationship between the development rates (f _{X,} X = E, L, P) and temperature T. Our preliminary analysis has identified a better fit with the Gaussian function. For each stage X, f _{X} is formulated by the Gaussian curve as a dependent of T, as shown in Equation 6. Coefficients and fitting rsquares for each stage X are given in Table 5. This relationship describes a biological pattern that f _{X} gradually increases at a low temperature, reaches a maximum in an optimal temperature range (28–34 °C), and then drops to zero at a lethal high temperature.
To explore the relationship between the larva/pupa mortality rate (m _{X}, X = L, P) and temperature, we explored functions established in other literature through trial and error, such as the monotonous exponential function (m _{X} = e ^{T/2} _{+} constant) [27] and the linear function (m _{X} = aT _{+} constant, X = A) [13]. However, none of these functions fit well with our observations. Our collected data showed that the mosquito had relatively higher mortality rates under high and low temperatures (Fig. 2df) [47, 48]. Based on this finding, a new regression function between the m _{X} and T is established, where m _{X} is adjusted to the range of [0, 1], as shown in Equation 7 and Table 5. This equation provides a better fit to our observed population abundance and is relatively consistent with the regression pattern of the mosquito’s survival rate [49–52]. Appendix A describes the process of establishing this equation in greater detail.
In addition to the temperaturedriven mechanism of development rates and mortality rates, we have adopted an established relationship between the egg ovipostion rate and temperature proposed by Yang et al. [53], as given in Equation 8. This equation is employed to substitute the oviposition rate in Tran et al.’s model [27], whereas in their model the rate was formulated as a constant.
Results
The scientific significance of the proposed model can only be corroborated by rigorous validation with respect to field data. To prepare for validation, we aimed to derive the simulation results first with an initial population of 10^{6} eggs and a starting time t _{0} of January 1^{st}. Then the model was discretized using the Euler Method in MatLab 2015 [54]. Specifically, multiple rounds of simulations were conducted on a daily basis over six years in Guangzhou (2006–2011) and Shanghai (2008–2013). The results derived for the first year (i.e. 2006 for Guangzhou and 2008 for Shanghai) were not used for comparison, as they were strongly dependent on the initial setting [24, 27].
We then compared larva abundance between the monthly observed CI (Table 1) and the monthly simulated results (L _{R}, which is relative to the maximum simulated value over the entire study period) using the Pearson’s correlation coefficient (r). In this case, the study period is 2007–2011 for Guangzhou and 2009–2013 for Shanghai.
Simulated diapause periods
In the model, two parameters m _{dia} and f _{dia} related to diapause egg survival remained undetermined. To derive this set of parameters, we compared them with m _{E} and f _{E} in the original model (Equation 2): m _{dia} = a _{1} m _{E} and f _{dia} = a _{2} f _{E} where 0 < a _{1} < 1, 0 < a _{2} < 1; and then we traversed a _{1} and a _{2} in a given range (a _{1}: 0.01–1, a _{2}: 0.01–1) and derived the optimal combination with the highest r for each city, as shown in Fig. 3. The highest r was chosen as it represented the best fit between the field data and the simulation results, eventually yielding the best model performance. In Guangzhou, the r is comparatively stable within a narrow range of 0.80–0.85 (Fig. 3a), whereas in Shanghai, r fluctuates over a wider range of 0.5–0.9 (Fig. 3b). The highest r was derived at a _{1} = 0.1 and a _{2} = 0.1 for both cities (at r = 0.84 for Guangzhou and r = 0.90 for Shanghai). Based on the optimal combination, we estimated the starting day and ending day of the diapause with the new model.
Figure 4 summarizes the finding about the diapause period and the nondiapause period with starting time (t _{b}) and ending time (t _{e}) in both cities. It shows that on average, the diapause in Guangzhou began in late November on the 330st day of the year (DOY 330) and ended in early March the next year (DOY 46); and in Shanghai, it began in midOctober (DOY 294) and ended in midMarch the next year (DOY 69). This result indicates an average diapause period of 81 days in Guangzhou and 140 days in Shanghai.
Simulated population abundance
Based on the highest r derived at a _{1} = 0.1 and a _{2} = 0.1, we simulated the population abundance for each of the seven stages on a daily basis. Figure 5 shows part of the results in the aquatic period (E, L, P) and the aerial period (A _{b}, A _{g}, A _{o}). A general observation is that the simulated population at each stage is nearly two times greater in Guangzhou (2007–2011) than that in Shanghai (2009–2013). This result might be attributed to the fact that the duration of the nondiapause period is longer in Guangzhou, which creates a more favorable environment for the mosquito to survive and develop. We then estimated the favorable development period at the larvalpupal stage (L + P) and at the adult stage (A _{b} + A _{o}) for the two cities on an annual basis, as given by Table 6.
Model validation
The validation of the model is shown in Figs. 6 and 7, using the correlation coefficient (r) and zerointercept rsquare (r _{ 0 } ^{2}), respectively.
Figure 6 is the comparison between the observed CI and the simulated L _{R} over a fiveyear period in Guangzhou and that in Shanghai (Fig. 6). The simulated population is highly consistent with the field observations (CI), with r equals to 0.84 for Guangzhou (Fig. 6a) and 0.90 for Shanghai (Fig. 6b). However, a good match with field data does not fully capture the fluctuation of population, leading to underestimations (green dots) and overestimations (red dots) for certain months. Four apparent underestimations are Oct 2008, May 2010 for Guangzhou and Jun 2011, Jun 2013 for Shanghai.
Figure 7 shows the r ^{2} for the comparison within each single year. Generally, over a fiveyear period, Shanghai (Fig. 7l) has a better fit than Guangzhou (Fig. 7f) in terms of r _{0} ^{2}. Specifically, in Shanghai, the simulation result is overall very satisfactory with the best fit appearing in 2012 (r _{0} ^{2} = 0.91, Fig. 7j); in Guangzhou, the simulation result in 2007 is considerably underestimated (r _{0} ^{2} = 0.67 compared to the average r _{0} ^{2} = 0.71 in Fig. 7f).
These results share similar findings with other case studies in Shanghai [55, 56] and Guangzhou [57, 58]. Specifically, the predicted population growth periods in Shanghai (i.e. F _{LP} > 1 %, MarchNovember; F _{A} > 1 %, AprilNovember, Table 6) parallel with existing studies [55, 56].
Sensitivity analysis
In addition to model validation, another pressing need is to examine the sensitivity of the new model to the change of parameters. The nine parameters independent of climatic variables (given in Table 3) were examined by using the fraction factorial design [59]. The initial values of these parameters were assigned according to existing literature [23, 24, 27, 47]. Then we reassigned a value to each of the parameters (e.g. m _{E}) within a ±10 % range (e.g. 0.9m _{E}1.1m _{E}) and derived the maximum r in the new model (Equation 5). We followed this method by changing each of the nine parameters, one at a time, and derived nine maximum r. Fig. 8 shows the results.
It can be seen from Fig. 8 that the change of r introduced by the singlefactor sensitivity analysis is generally minor. The overall improvement of r ranges from 0.004 to 0.017 (or 0.4–2.0 %) in Guangzhou and from 0 to 0.016 (0–1.7 %) in Shanghai. The greatest change is observed on parameter σ in the case of Guangzhou and m _{E} in Shanghai. This analysis provides compelling evidence that the proposed model is relatively robust and is not sensitive to nonclimatic variables, corroborating the applicability of the model.
Discussion
This section evaluates the improvement of the new model based on several adjustments to the original model and discusses the difference in model performance between the two study areas.
Model comparison
First, we would like to examine the improvement of the new model (designated as Model AB) by comparing with the original model (Model O) and two partially adjusted models (Model A and Model B) according to the Methods section.

Model O: Original model (Equation 2)

Model A: Original model with diapauserelated structural adjustment (Equation 5)

Model B: Original model with temperaturerelated parameter adjustment (Equation 2 + Equations 6–8)

Model AB: New model with all adjustments applied (Equation 5 + Equations 6–8)
Figure 9 shows the correlation coefficient r with the two sets of field data applied into these four models. It can be seen from the comparison that although Model O is well suited for Guangzhou (r = 0.77), it does not perform well for Shanghai (r = 0.42). This result is possibly introduced by the static diapause parameters in Model O that could not fully capture the climaticdriven mechanism of diapause, degrading the performance of the model. When the diapause effect is included, Model A demonstrates a huge improvement in cases of Shanghai (r = 0.75 or an improvement of 78.6 %), where diapause shows a confirmed seasonal pattern. Adjusting the model with only temperaturerelated parameters (Model B) yields moderate improvement in both study areas; and the best performance is achieved with two facets of adjustments applied.
Regional comparison
The results above corroborate the model’s potential in simulating Ae. albopictus population. However, compared to Shanghai (31.2°N), the performance in Guangzhou (23.1°N) seems to be less satisfactory (i.e. curve fitting rate is 90 % for Shanghai and 84 % for Guangzhou, Fig. 6). This regional difference could be explained in two very different respects.
First, the effect of the diapause in the lowlatitude area is of arguable existence [7]. Guangzhou is located further south and has comparatively warm and humid winters. The mild climatic conditions could be favorable for certain Ae. albopictus to sustain and develop without taking the strategy of diapause. For example, our field experiments showed evidence that during winter times there were still a very low number of larvae, pupae and adults (Fig. 5) in Guangzhou that could circumvent the process of diapause. This observation was supported by Liu et al. [58] that found in the same region a small proportion of Ae. albopictus larvae hatched after midNovember could still survive the extreme weather and grow into nonovipositing adults. A separate evidence could be found in Fig. 9 that applying the diapauserelated structural adjustment is less effective in Guangzhou (r: 0.72 → 0.80) than in Shanghai (r: 0.42 → 0.75).
Secondly, the seasonality of Ae. albopictus is relatively complex in tropical climate areas. In an annual development cycle, Ae. albopictus population showed one peak in subtropical climate areas and two peaks with different magnitudes in tropical climate areas [60]. As Guangzhou is located in close vicinity to tropical climate areas, our field observations of CI show interchangeable patterns (i.e. 2007, 2010 and 2011 have two peaks, while 2008 and 2009 have one peak, Fig. 6a). This dynamic pattern created a certain degree of uncertainty for validation and eventually degraded model performance for cases in Guangzhou. Comparatively, Shanghai with an apparent singlepeak pattern in the years of study yielded better simulation results (Fig. 6b).
Conclusions
The phenomenon of diapause among Ae. albopictus has been confirmed in most temperate climate areas and is of arguable existence in subtropical climate areas [7]. Generally, low temperatures and short photoperiods create an unfavorable condition for Ae. albopcitus to develop. To survive the extreme conditions, especially during winter times, egg diapause occurs as an adaptive strategy to lower mortality rates [1, 7]. This climatedriven mechanism of diapause has only been explored in a limited manner [24, 27].
This paper proposes a climatedriven mechanistic population model of Ae. albopictus that accounts for the biological phenomenon of diapause. The model is a natural extension of two existing mechanistic population models [24, 27] with emphasis on the climatedriven diapause conditions and stagespecific moderating variables. Although the former models also considered diapause, several issues remained unaddressed, such as identifying the time when diapause eggs are laid and hatched, demarcating the thresholds of diapause and nondiapause periods, and incorporating the mortality rate and the hatching rate of diapause eggs. To remedy these flaws, an improved generic model is proposed, capturing the multifaceted climatedriven mechanism of diapauserelated effects. The formation of the model and the attribution of parameters are finetuned to existing research and field data collected in two Chinese cities over a respective fiveyear period. Overall, the simulation results are relatively compelling and fit the majority of our field observations. The study also confirms the respective as well as the joint effects of model structure and temperaturedriven parameters, corroborating findings from other mechanistic models [23, 24, 27].
Admittedly, the proposed model is methodologically flawed in several aspects. First, the parameter of precipitation was included but was not closely examined in our model. Existing studies have identified either a positive [41], a negative [12], or no influence of precipitation on population abundance [43]. The actual effect of precipitation should be carefully weighted through rigorous field observations before any attempt to quantify the variable in the model. Secondly, in the study only the Container Index was used for model validation and could not fully represent the multifaceted population growth of Ae. albopictus. Future research should explore other monitoring indices, such as the Breteuil Index (BI) and the Housing Index (HI) to generate a solid and robust conclusion. Thirdly, the sensitivity analysis on the nonclimatic variables (Table 3) is exclusive of other variables, while the joint effect of multiple variables is not evaluated. In the future the extension of the assessment should include multivariate analysis methods, such as the Fourier amplitude sensitivity testing [59, 61].
Lastly, a promising direction to extend the study considers the incorporation with other nonmechanistic ecological models, such as the GARP [62] model and the CLIMEX model [63] and the correlation with epidemic models, such as the SIR model [64, 65] to describe the mechanistic transmission and the spatiotemporal distribution of mosquitoborne epidemics. For example, one such an attempt is the work conducted by Erickson et al. [66] that integrated the research on Ae. albopcitus population [23] with an SEIR model [66]. To this end, better understanding the mechanism and variables effecting the Ae. albopictus population growth and improving the model performance will eventually contribute to proactive strategies to predict and prevent contingent mosquitoborne epidemics.
References
 1.
Novak R. The Asian tiger mosquito, Aedes albopictus. Wing Beats. 1992;3(3):5.
 2.
Gratz NG. Critical review of the vector status of Aedes albopictus. Med Vet Entomol. 2004;18(3):215–27.
 3.
Knudsen A. Global distribution and continuing spread of Aedes albopictus. Parassitologia. 1995;37(23):91–7.
 4.
Laird M, Calder L, Thornton RC, Syme R, Holder P, Mogi M. Japanese Aedes albopictus among four mosquito species reaching New Zealand in used tires. J Am Mosq Control Assoc. 1994;10(1):14–23.
 5.
Reiter P, Sprenger D. The used tire trade: a mechanism for the worldwide dispersal of container breeding mosquitoes. J Am Mosq Control Assoc. 1987;3(3):494–501.
 6.
Hales S, De Wet N, Maindonald J, Woodward A. Potential effect of population and climate changes on global distribution of dengue fever: an empirical model. Lancet. 2002;360(9336):830–4.
 7.
Hawley WA. The biology of Aedes albopictus. J Am Mosq Control Assoc Suppl. 1988;1:1–39.
 8.
Juliano SA. Population dynamics. J Am Mosq Control Assoc. 2007;23(2):265–75.
 9.
Kendall BE, Briggs CJ, Murdoch WW, Turchin P, Ellner SP, McCauley E, Nisbet RM, Wood SN. Why do populations cycle? A synthesis of statistical and mechanistic modeling approaches. Ecology. 1999;80(6):1789–805.
 10.
Alto BW, Juliano SA. Temperature effects on the dynamics of Aedes albopictus (Diptera: Culicidae) populations in the laboratory. J Med Entomol. 2001;38(4):548–56.
 11.
Alto BW, Juliano SA. Precipitation and temperature effects on populations of Aedes albopictus (Diptera: Culicidae): implications for range expansion. J Med Entomol. 2001;38(5):646–56.
 12.
Roiz D, Rosa R, Arnoldi D, Rizzoli A. Effects of temperature and rainfall on the activity and dynamics of hostseeking Aedes albopictus females in northern Italy. Vectorborne Zoonot. 2010;10(8):811–6.
 13.
Shaman J, Spiegelman M, Cane M, Stieglitz M. A hydrologically driven model of swamp water mosquito population dynamics. Ecol Model. 2006;194(4):395–404.
 14.
Rueda L, Patel K, Axtell R, Stinner R. Temperaturedependent development and survival rates of Culex quinquefasciatus and Aedes aegypti (Diptera: Culicidae). J Med Entomol. 1990;27(5):892–8.
 15.
Teng HJ, Apperson CS. Development and survival of immature Aedes albopictus and Aedes triseriatus (Diptera: Culicidae) in the laboratory: effects of density, food, and competition on response to temperature. J Med Entomol. 2000;37(1):40–52.
 16.
Ahumada JA, Laoointe 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.
 17.
Depinay JMO, Mbogo CM, Killeen G, Knols B, Beier J, Carlson J, et al. A simulation model of African Anopheles ecology and population dynamics for the analysis of malaria transmission. Malaria J. 2004;3(1):29.
 18.
Focks D, Haile D, Daniels E, Mount G. Dynamic life table model for Aedes aegypti (Diptera: Culicidae): simulation results and validation. J Med Entomol. 1993;30(6):1018–28.
 19.
Focks DA, Haile D, Daniels E, Mount GA. Dynamic life table model for Aedes aegypti (Diptera: Culicidae): analysis of the literature and model development. J Med Entomol. 1993;30(6):1003–17.
 20.
Gong H, DeGaetano AT, Harrington LC. Climatebased models for West Nile Culex mosquito vectors in the Northeastern US. Int J Biometeorol. 2011;55(3):435–46.
 21.
Otero M, Solari HG, Schweigmann N. A stochastic population dynamics model for Aedes aegypti: formulation and application to a city with temperate climate. B Math Biol. 2006;68(8):1945–74.
 22.
Schaeffer B, Mondet B, Touzeau S. Using a climatedependent model to predict mosquito abundance: application to Aedes (Stegomyia) africanus and Aedes (Diceromyia) furcifer (Diptera: Culicidae). Infect Genet Evol. 2008;8(4):422–32.
 23.
Erickson RA, Presley SM, Allen LJ, Long KR, Cox SB. A stagestructured, Aedes albopictus population model. Ecol Model. 2010;221(9):1273–82.
 24.
Cailly P, Tran A, Balenghien T, L’Ambert G, Toty C, Ezanno P. A climatedriven abundance model to assess mosquito control strategies. Ecol Model. 2012;227:7–17.
 25.
Clements A. The biology of mosquitoes: sensory, reception and behaviour. Oxon: CABI Publishing; 1999.
 26.
Clements A. Vol. 1: Development, nutrition and reproduction. London: CABI Publishing; 1992.
 27.
Tran A, L’Ambert G, Lacour G, Benoit R, Demarchi M, Cros M, et al. A rainfall and temperaturedriven abundance model for Aedes albopictus populations. Int J Environ Res Public Health. 2013;10(5):1698–719.
 28.
Chapman RF. The insects: structure and function. Cambridge: Cambridge university press; 1998.
 29.
Tauber MJ, Tauber CA, Masaki S. Seasonal adaptations of insects. Oxford: Oxford University Press; 1986.
 30.
Hanson SM. Field overwinter survivorship of Aedes albopictus eggs in Japan. J Am Mosq Control Assoc. 1995;11(3):354–7.
 31.
Wang R. Observations on the influence of photoperiod on egg diapause in Aedes albopictus Skuse. Acta Entomol Sinica. 1966;15:75–7.
 32.
Imai C, Maeda O. Several factors effecting on hatching of Aedes albopictus eggs. Jpn J Sanit Zool. 1976;27:363–72.
 33.
Hong H, Shim J, Young H. Hibernation studies of forest mosquitoes in Korea. Korean J Entomol. 1971;1:13–6.
 34.
Hanson SM, Craig GB. Cold acclimation, diapause, and geographic origin affect cold hardiness in eggs of Aedes albopictus (Diptera: Culicidae). J Med Entomol. 1994;31(2):192–201.
 35.
Sota T, Mogi M. Survival time and resistance to desiccation of diapause and non‐diapause eggs of temperate Aedes (Stegomyia) mosquitoes. Entomol Exp Appl. 1992;63(2):155–61.
 36.
Mosquito Container Index Database. Guangzhou Center for Disease Control and Prevention, Guangzhou, China. 2016. http://www.gzcdc.org.cn/Health/View.aspx?id=199. Accessed 5 Mar 2016.
 37.
National Oceanic and Atmospheric Administration: Earth System Research Laboratory. http://www.srrb.noaa.gov/highlights/sunrise/calcdetails.html (2016). Accessed 5 Mar 2016.
 38.
EstradaFranco J. Biology, disease relationships, and control of Aedes albopictus. Washington D.C: Pan American Health Organization; 1995.
 39.
Chan K, Ho B, Chan Y. Aedes aegypti (L.) and Aedes albopictus (Skuse) in Singapore City: 2. Larval habitats. B World Health Organ. 1971;44(5):629–33.
 40.
Ho B, Chan K, Chan Y. Aedes aegypti (L.) and Aedes albopictus (Skuse) in Singapore city: 3. Population fluctuations. B World Health Organ. 1971;44(5):635–41.
 41.
LourençodeOliveira R, Castro MG, Braks MA, Lounibos LP. The invasion of urban forest by dengue vectors in Rio de Janeiro. J Vector Ecol. 2004;29:94–100.
 42.
Mori A, Oda T. Studies on the egg diapause and overwintering of Aedes albopictus in Nagasaki. Trop Med. 1981;23(2):79–90.
 43.
Toma L, Severini F, Luca M, Bella A, Romi R. Seasonal patterns of oviposition and egg hatching rate of Aedes albopiczus in Rome. J Am Mosq Control Assoc. 2003;19(1):19–22.
 44.
Soti V, Tran A, Degenne P, Chevalier V, Lo Seen D, Thiongane Y, et al. Combining hydrology and mosquito population models to identify the drivers of Rift Valley fever emergence in semiarid regions of West Africa. PLoS Negl Trop Dis. 2012;6(8):e1795.
 45.
Edgerly JS, Willey MS, Llvdahl TP. The community ecology of Aedes egg hatching: implications for a mosquito invasion. Ecol Entomol. 1993;18(2):123–8.
 46.
Yu X. Comparison studies on the tolerance range of Aedes albopictus and Aedes aegypti to temperature and photoperiod (Master thesis, In Chinese). Beijing: China CDC; 2013.
 47.
Delatte H, Gimonneau G, Triboire A, Fontenille D. Influence of temperature on immature development, survival, longevity, fecundity, and gonotrophic cycles of Aedes albopictus, vector of chikungunya and dengue in the Indian Ocean. J Med Entomol. 2009;46(1):33–41.
 48.
Brady OJ, Johansson MA, Guerra CA, Bhatt S, Golding N, Pigott DM, et al. Modelling adult Aedes aegypti and Aedes albopictus survival at different temperatures in laboratory and field settings. Parasit Vector. 2013;6(1):351–62.
 49.
Martens W. Health impacts of climate change and ozone depletion: an ecoepidemiologic modeling approach. Environ Health Persp. 1998;106 Suppl 1:241–51.
 50.
Craig M, Snow R, Le Sueur D. A climatebased distribution model of malaria transmission in subSaharan Africa. Parasitol Today. 1999;15(3):105–11.
 51.
Ruiz D, Poveda G, Vélez ID, Quiñones ML, Rúa GL, Velásquez LE, et al. Modelling entomologicalclimatic interactions of Plasmodium falciparum malaria transmission in two Colombian endemicregions: contributions to a National Malaria Early Warning System. Malaria J. 2006;5:66.
 52.
Parham PE, Michael E. Modeling the effects of weather and climate change on malaria transmission. Environ Health Persp. 2010;118(5):620–6.
 53.
Yang H, Macoris M, Galvani K, Andrighetti M, Wanderley D. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiol Infect. 2009;137(8):1188–202.
 54.
MathWorks: Matlab. http://cn.mathworks.com/products/matlab/. Accessed 5 Mar 2016.
 55.
Xu R, Cai E, Xu Y, Wang S, Liu W. Surveilance on larva and adult of Aedes albopictus and the study of insecticide in Shanghai. Chin J Vec Biol Contr. 2006;12(5):345–8.
 56.
Wang S, Xu R, Xu Y, Ren W, Jiang L, Hu J, Cai E. Surveillance on dengue vector and breading habit of vector in Shanghai area (in Chinese). Chin J Vector Bio Control. 2003;14(2):124–6.
 57.
Gao J, Liu J, Liu J, Liu S. Research and control strategies for Aedes albopictus in Guangzhou, China, during recent five years (in Chinese). Chin J Vec Biol Contr. 2014;25(5):483–5.
 58.
Liu F, Zhou C, Lin P. Studies on the population ecology of Aedes albopictus: 5. the seasonal abundance of natural population of Aedes albopictus in Guangzhou (in Chinese). Acta Sci Nat Univ Sunyatseni. 1990;29(2):118–22.
 59.
Saltelli A, Ratto M, Andres T, Campolongo F, Cariboni J, Gatelli D, et al. Global sensitivity analysis: the primer. Hoboken: Wiley; 2008.
 60.
Wu H, Liu Q, Liu X, Lu L, Guo Y, Yue Y. Survillence for Aedes albopictus in China, 2006—2013 (in Chinese). Dis Surveill. 2015;30(4):310–5.
 61.
Saltelli A, Tarantola S, Campolongo F. Sensitivity analysis as an ingredient of modeling. Stat Sci. 2000;15(4):377–95.
 62.
Warren DL, Glor RE, Turelli M. Environmental niche equivalency versus conservatism: quantitative approaches to niche evolution. Evolution. 2008;62(11):2868–83.
 63.
Poutsma J, Loomans AJM, Aukema B, Heijerman T. Predicting the potential geographical distribution of the harlequin ladybird, Harmonia axyridis, using the CLIMEX model. BioControl. 2008;53(1):103–25.
 64.
Derouich M, Boutayeb A, Twizell EH. A model of dengue fever. Biomed Eng Online. 2003;2:4.
 65.
Derouich M, Boutayeb A. Dengue fever: mathematical modelling and computer simulation. Appl Math Comput. 2006;177(2):528–44.
 66.
Erickson RA, Presley SM, Allen LJ, Long KR, Cox SB. A dengue model with a dynamic Aedes albopictus vector population. Ecolo Model. 2010;221(24):2899–908.
 67.
Lunde TM, Korecha D, Loha E, Sorteberg A, Lindtjørn B. A dynamic model of some malariatransmitting anopheline mosquitoes of the Afrotropical region. I. Model description and sensitivity analysis. Malar J. 2013;12:28.
Acknowledgements
The research was financially supported by National Key Basic Research and Development Plan (2012CB955501) of China Ministry of Science and Technology.
The research was financially supported by National Natural Science Foundation of China (81273139) of China Ministry of Science and Technology.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
PJ formulated the mechanistic population model, performed the statistical analysis, and drafted the manuscript. LL coordinated the collection of field data of Ae. albopictus in Shanghai and Guangzhou. XC led the research by creating a research agenda and restructuring the manuscript. JC suggested critical directions for incorporating the climatedriven mechanism of diapause in the model. LG suggested directions for formulating the model and helped with testing the variables in the model. XY helped with collecting field data of Ae. albopictus. QL helped with evaluating the model and suggested directions for the involvement of diapauserelated behaviors. All authors read and approved the final version of the manuscript.
Appendix. Formulation of the mortality rate at the larval stage and the pupal stage
Appendix. Formulation of the mortality rate at the larval stage and the pupal stage
To establish Equation 7, we first estimated the daily survival rate at the larval stage and at the pupal stage with the existing dataset (Table 4). Specifically, S _{X}(t) denotes the survival rate at time t at stage X (larval or pupal) and S _{X}(t → t + d) denotes the total survival rate from day t to day (t + d). Supposing S _{X}(t) was stable from day (t) to day (t + d), S _{X}(t) was calculated by Equation 9. Within the time length of S _{X}(t → t + d), the daily survival rate at stage X under different temperatures was calculated (Table 5). According to Lunde et al. [67], we also obtained the daily mortality rate m _{X}(t) at stage X under the corresponding temperature by Equation 10.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Jia, P., Lu, L., Chen, X. et al. A climatedriven mechanistic population model of Aedes albopictus with diapause. Parasites Vectors 9, 175 (2016). https://doi.org/10.1186/s130710161448y
Received:
Accepted:
Published:
Keywords
 Aedes albopictus
 Population dynamics
 Mechanistic model
 Climatedriven
 Diapause