A mathematical model of dengue diseases transmission will be discussed in this paper. Various interventions, such as vaccination of adults and newborns, the use of insecticides or fumigation, and also the enforcement of mechanical controls, will be considered when analyzing the best intervention for controlling the spread of dengue. From model analysis, we find three types of equilibrium points which will be built upon the dengue model. In this paper, these points are the mosquito-free equilibrium, disease-free equilibrium (with and without vaccinated compartment), and endemic equilibrium. Basic reproduction number as an endemic indicator has been found analytically. Based on analytical and numerical analysis, insecticide treatment, adult vaccine, and enforcement of mechanical control are the most significant interventions in reducing the spread of dengue disease infection caused by mosquitoes rather than larvicide treatment and vaccination of newborns. From short- and long-term simulation, we find that insecticide treatment is the best strategy to control dengue. We also find that, with periodic intervention, the result is not much significantly different with constant intervention based on reduced number of the infected human population. Therefore, with budget limitations, periodic intervention of insecticide strategy is a good alternative to reduce the spread of dengue.
Ministry of Research and Higher Education1. Introduction
Dengue is the most rapidly growing disease in the world [1]. The disease is spread by Aedes mosquitoes and is therefore often referred to as a mosquito-borne viral disease. The disease has become endemic in more than 100 countries, including the Caribbean, Africa, the Americas, the Pacific, and Asia, including Indonesia [1]. As an endemic disease, dengue occurs regularly in subtropical and tropical regions of the world, and approximately 40% of people live in regions of the world where there is a risk of contracting it [2]. Dengue is a vector-borne disease transmitted from an infected human to a female Aedes aegypti mosquito by a bite. The mosquito, which needs regular meals of blood to mature its eggs, completes the cycle by biting a healthy human, transmitting the disease in one act [3].
Until now, the primary prevention for dengue has been control of mosquitoes, in both larval and adult forms. Larval control is carried out by larvicide treatment using long-lasting chemicals to kill larvae, which sure preferably have WHO clearance for use in drinking water [4]. Mechanical controls are also used to control larvae, with assistance from campaigns and educational programs carried out by governments. In Indonesia, such a program is known as 3M and consists of educating people about the importance of draining and shutting down and burying all tubs, buckets, or containers of water that which be used by female mosquitoes to breed and lay their eggs [5]. The larvae of Aedes aegypti can also grow in used goods that can hold water, and it is therefore recommended to make sure the environment around the house has no space that could allow mosquitoes to breed.
Adult mosquito control is achieved by the use of insecticide. Insecticide fumigation targets the vector Aedes aegypti mosquito as the main control of dengue epidemics. However, the long-term use of insecticides and larvicides poses several risks: one is resistance of the mosquito to the product, reducing its efficacy, while genetic mutation of the mosquito, making it less susceptible to the effects of the product, is another. Such products have also been linked to numerous adverse health effects including the worsening of asthma and respiratory problems [3, 6]. In Surabaya, Indonesia, larval mortality rates of under 80% indicate possible resistance of Aedes aegypti to the insecticide temephos [7], and outside Indonesia, resistance to insecticide has been reported in multiple countries. In recent years, the frequency of kdr mutations associated with pyrethroid resistance has increased rapidly [8]. Pyrethroids have become the most frequently used public health insecticides globally due to their low cost and low toxicity to mammals [9], and they are of considerable concern when kdr is found in wild populations of vector mosquitoes [8]. With the many cases of Aedes aegypti mosquito’s resistance to insecticide, it is therefore necessary to develop alternative strategies to slow its evolution.
Besides controlling dengue via control of mosquitoes population, one of the alternative strategies that is being used is dengue vaccine. In December 2015, the first vaccine against dengue by Sanofi Pasteur, Dengvaxia (CYD-TDV), was approved in three highly endemic countries: Mexico, Philippines, and Brazil [10]. This vaccine is the world’s first dengue vaccine and is already licensed for individuals aged 9–45 years for the prevention of infectious disease caused by four dengue virus serotypes (DEN 1, DEN 2, DEN 3, and DEN 4) [11]. In Indonesia, so far, the government is still conducting clinical trials to determine its effectiveness. However, assessments of the public’s acceptance of the dengue vaccine and its associated factors are widely lacking [12]. A lack of understanding about the importance of vaccination against dengue to the public will be able to reduce the success rate of vaccination interventions in various countries, especially in a country that is less intensive to educate people about the importance of dengue fever vaccination [13]. In 2017, Dengvaxia would have been applied if proven effective and suitable for the dengue serotypes which are pandemic in Indonesia [14].
The earliest mathematical models for dengue disease transmission are developed in [15, 16] which are closely related to the models for the transmission of malaria discussed in [17, 18]. The authors in [16] create the model for two types of viruses by allowing temporary cross immunity and increased susceptibility to the second infection due to the first infection. The intervention has not yet been used into the mathematical model [16]. In [19], the mathematical model with only insecticide campaign intervention is discussed. It has been shown that, with a steady insecticide campaign, it is possible to reduce the number of infected humans and mosquitoes and prevent an outbreak that could transform an epidemiological episode to an endemic disease [19]. A year after the research discussed in [19], it was updated [20], and the mathematical model for dengue was updated continuously with all controls included, that is, (1) proportion of larvicide, (2) proportion of adulticide, and (3) proportion of mechanical control. The results have shown that, even with a low, although continuous, index of control adulticide over time, the results are surprisingly positive [20]. However, it has been stated that to rely only on adulticide is a risky decision [20]. The research in [6, 7, 21] supports this claim, citing the problem of Aedes aegypti mosquito’s resistance to insecticide. Under the new achievement in the field of vaccination technology with the discovery of the first vaccine against dengue by Sanofi Pasteur, the work in [22] devised two models, one assuming that unintentional vaccination increases the infectious period and another assuming that unintentional vaccination leads to the development of symptoms. This argument is also supported by [3], in which the mathematical model is created with the vaccine as the new compartment, arguing that the vaccine must divide the human population into classes, that is, the perfect pediatric vaccine for newborns and perfect adult vaccine (conferring 100% protection throughout life), and also classes of human with imperfect vaccine effect [3]. Other mathematical models with different intervention were also introduced in [23] which discuss the use of mosquito repellent to reduce probability of success of infection in human population and in [4] which discuss the use of sterile mosquito strategy.
According to above explanation, it is important to find the best strategy for controlling dengue spreads for both short-term and long-term interventions. Therefore, a mathematical model of dengue disease transmission by using adult and newborn vaccines with waning immunity, the use of insecticides and larvicides, and mechanical control will be developed in the next section. Equilibrium points will be found, which ensure the existence of local stability. Basic reproduction numbers will be obtained as the main factor in whether the disease will become epidemic in a population or not. Numerical analysis for comparing the dynamic of infected humans and mosquitoes will be used to support the model interpretation.
2. Mathematical Model Construction
To construct our model, firstly we divide the human population into four compartments, that is,
Sh(t): susceptible (individuals who can be infected with dengue);
Vh(t): vaccinated (individuals who have had the vaccine injected into their bodies, making them resistant to infectious disease. However, the use of the vaccine does not provide perfect immunity. There will be a time when the vaccine does not work properly in the body or when the effect of the vaccine has begun to subside [3]);
Inh(t): infected (individuals who are infected with dengue. In this case, the infected human is incapable of transmitting the disease to other humans);
Rh(t): recovered (individuals who have recovered from dengue and have acquired temporal immunity to respective DEN virus).
On the other hand, we divide the mosquito population into three compartments, that is,
Av(t): aquatic phase (the phase that includes the egg, larvae, and pupa stages, which live in water);
Sv(t): susceptible (mosquitoes that are able to infect with dengue);
Inv(t): infected (mosquitoes that have been infected with dengue by an infected human and are capable of transmitting dengue to humans).
Secondly, we have made some assumptions that we will use to describe a dynamic process in our model that we will construct: (1) There is no migration in either human or mosquito population. (2) Humans and mosquitoes are assumed to be born susceptible, there is no natural protection, and dengue is not passed onto the next generation (no vertical transmission) [19]. (3) The transmission process in susceptible and vaccinated humans is simply by the bite of an infected mosquito. Infected humans cannot transmit the virus to other susceptible or vaccinated humans [19]. (4) The death rate is considered to be a natural death rate in both populations. (5) Vaccinated human status is considered temporary because of the ability of the vaccine to subside over time [3]. (6) There is no recovered phase in mosquitoes due to their short lifespan [20]. (7) There is no resistant (immune) effect in mosquitoes to the use of synthetic fumigation, such as insecticides and larvicides [6, 7]; in this article, we assume that there is no resistant (immunity) effect to mosquitoes due to use of synthetic fumigation, such as insecticide and larvicide.
With the assumptions, variables, and transmission diagram given in Figure 1, the model is represented as a seven-dimensional system of differential equations which are given by (1)dShdt=1-u2Ah-bβhShInvNh-u1Sh-μhSh+δhRh+θhVh,dVhdt=u2Ah+u1Sh-bξβhVhInvNh-μhVh-θhVh,dInhdt=bβhInvSh+ξVhNh-γhInh-μhInh,dRhdt=γhInh-μhRh-δhRh,dAvdt=ϕ1-Avu5kNhSv+Inv-u3Av-ηvAv-μAAv,dSvdt=ηvAv-bβvSvInhNh-μvSv-u4Sv,dInvdt=bβvSvInhNh-μvInv-u4Inv,with parameters description being the following:
Nh: total of human population
Ah: human per capita birth rate
b: average number of bites of humans by mosquitoes
βh: average of successful transmission in human
βv: average of successful transmission in mosquitoes
μh: human death rate
μv: mosquito death rate
μA: larval death rate
δh: the rate of change from Rh to Sh because of the disappearance of the temporal natural immunity
θh: the rate of change from Vh to Sh because of the disappearance of the temporal vaccine effect
ξ: reduction of βh because of vaccine
γh: human recovery rate
ϕ: average number of eggs at each deposit
k: ratio for number of larvae per human
ηv: transition rate from Av to Sv
u1: adult vaccination rate
u2: newborn vaccine
u3: larvicide rate
u4: fumigation rate
u5: enforcement of mechanical control proportion to reduce k
Mathematical model using vaccination of adults and newborns, fumigation and larvicide treatment, and enforcement of mechanical control.
In the next section, mathematical model analysis to find equilibrium points and their local stability criteria will be given.
3. Equilibrium Points and Local Stability
From the system of equation (1), we find three types of equilibrium points.
3.1. Mosquito-Free Equilibrium (MFE) Point
It is the equilibrium where the mosquito virus does not exist in the living environment, so any infectious disease never occurs at the MFE point (sterile conditions). This equilibrium is given by(2)Sh1=-Ahμh+θhR2-1μhμh+u1+θh,Vh1=Ahμhu2+u1μhμh+u1+θh,Inh1=0,Rh1=0,Av1=0,Sv1=0,Inv1=0.At the MFE point, it can be seen that Sh1 can have either a positive or a negative value. However, from an epidemiological point of view, a population has biological meaning only if it has a nonnegative value. Therefore, to ensure that the MFE point exists, the positiveness for Sh1 should be made; that is,(3)R2=μhu2μh+θh<1.To guarantee local stability of equilibrium, it is necessary to ensure that all eigenvalues of system (1) in its Jacobian matrix evaluated in the MFE have negative values. The condition is (4)ηvϕμv+u4u3+ηv+μA=R1<1.Please note that R1 is known as the basic offspring number, which will guarantee the existence of mosquitoes in the population. The mosquito population will exist if, and only if, R1≥1 as will be discussed in the next equilibrium point.
3.2. Disease-Free Equilibrium (DFE) Point
DFE is an equilibrium where mosquito and human populations exist in the living environment, but the virus does not occur.
We have two types of DFE equilibrium point, which are DFE-1 and DFE-2. DFE-1 describes a condition where infected, recovered human groups and infected, recovered mosquito groups do not exist, while DFE-2 describes a condition where infected, recovered human groups, infected, recovered mosquito groups, and human vaccinated groups do not exist (special case without intervention of vaccines (u1=u2=0)).
3.2.1. Disease-Free Equilibrium I (DFE-1) Point
This equilibrium is given by(5)Sh2=-Ahμh+θhR2-1μhμh+u1+θh,Vh2=Ahμhu2+u1μhμh+u1+θh,Inh2=0,Rh2=0,Av2=u5kNhR1-1,Sv2=u5kNhηvR1-1μv+u4,Inv2=0,which will exist if and only if (6)R1=ηvϕμv+u4u3+ηv+μA>1,R2=μhu2μh+θh<1.DFE-1 will be locally asymptotically stable if and only if(7)ηvϕμv+u4u3+ηv+μA=R1>1,(8)b2βvu5kβhAhu1ξ+μh+θh+μh+θhξ-1R2R1-1ϕμv+u42γh+μhNhμhμh+u1+θh=R0<1.Based on (8), it is known that
if R1<1, then R0<0 will be obtained;
if R1>1, then R0∈(0,1) or R0∈(1,∞) will be obtained.
3.2.2. Disease-Free Equilibrium II (DFE-2) Point
This equilibrium is given by(9)Sh3=Ahμh,Vh3=0,Inh3=0,Rh3=0,Av3=u5kNhR1-1,Sv3=u5kNhηvR1-1μv+u4,Inv3=0,which only exist if and only if(10)R1=ηvϕμv+u4u3+ηv+μA>1.To guarantee the local stability of equilibrium, it is necessary to ensure that all eigenvalues of system (1), evaluated in the Jacobi matrix on DFE-2 point, are negative. The condition is (11)Ahβhβvb2ku5ϕηvμhNhμv+u4μh+γhϕ+b2u5βhβvkAhu3+μA+ηvμv+u4=R6<1,ηvϕμv+u4u3+ηv+μA=R1>1.
3.3. Endemic Equilibrium (EE)
Endemic equilibrium describes a condition where all compartments, both human and mosquitoes, achieve coexistence. The endemic equilibrium point of system (1) is not in simple way to be written in explicit form. However, the existence of this equilibrium point might be written as equilibrium points that depend on values of Inv and Inh which are given by (12)Sh∗=Inh∗Inv∗bβhξ+μhInh∗Nhθh+μhμh+γh-NhbβhAhu2ξInv∗Nhξu1+bβhInv∗ξ+Nhμh+NhθhbβhInv∗,Vh∗=Nhu2AhbβhInv∗+u1NhγhInh∗+u1NhμhInh∗Nhξu1+bβhInv∗ξ+Nhμh+NhθhbβhInv∗,Rh∗=γhInh∗μh+δh,Av∗=Inv∗μv+u4Nhu4+Nhμv+bβvInh∗ηvbβvInh∗,Sv∗=Inv∗Nhμh+u4bβvInh∗,while Inv∗ and Inh∗ are taken from positive solution of(13)F1=bβvInh∗+μvNh+u4NhbβvInh∗ϕu5kNhηv-muv+u4bkInh∗Nhβvηvu5+bkInh∗NhβvμAu5+bkInh∗Nhβvu3u5+bϕInh∗Inv∗βv+ϕInv∗Nhμv+u4=0,F2=AInv2+BInh+CInh+DInv2+EInv,with (14)A=-b2ξβh2μhδh+γh+μh,B=-bNhβhμhξδhγh+ξδhμh+ξδhu1+ξγhμh+ξγhu1+ξμh2+ξμhu1+δhμh+δhθh+γhμh+γhθh+μh2+μhθh,C=-Nh2μhμh+θh+u1μh+γhδh+μh,D=b2ξAhβh2δh+μh,E=bAhNhβhξμhu2+ξu1-μhu2+μh+θhδh+μh.Substituting all parameters values from Table 1 into above couple of equations will give us existence of Inv and Inh numerically as shown in Figure 2. It can be seen that, as long as the intersection between F1 and F2 is in the first quadrant, we will have a positive endemic equilibrium.
Parameters values.
Parameters
Value
Description
Nh
1000
Total of human population is assumed to be 1000 people.
μh
1/(65×365)
Since human life expectation is approximately 65 years, we have μh=1/(65×365) [23].
Ah
1000/(65×365)
Since in our model the total of the human population is constant, we have Ah=Nhμh=1000/(65×365).
ϕ
300
We assume that each female Aedes aegypti produces 300 eggs at each spawning.
βh,βv
0,1
It is assumed that it needs 10 successful contacts to infect a human/mosquito with dengue [23].
θh
1/60
We assume that the effect of vaccination will have disappeared in 60 days.
γh
1/14
The natural recovery rate for the human population from dengue is 14 days [4].
ξ
0.1
With vaccination, the infection rate from mosquito to human population will be reduced by 90%.
μv
1/30
Life expectation of the mosquito population is 30 days [4].
μA
0.75/21
We assume that there is only a 25% chance that larvae might grow and become adult mosquitoes, with time to transition being 21 days.
ηv
0.25/21
Transition from aquatic phase to adult mosquito [4].
δh
1/30
Short-term immunity of humans to dengue after recovery is 30 days [23].
b
1
Mosquitoes only bite once a day [23].
k
2
We assume that the ratio between human and adult mosquitoes is 2; that is, each human related to 2 adult mosquitoes.
Existence of endemic equilibrium for Inv and Inh depending on couple of polynomial characteristics.
For simple case when no intervention is given into system (1) (u1=u2=u3=u4=0 and u5=1), endemic equilibrium point is given by (15)Sh+,Vh+,Inh+,Rh+,Av+,Sv+,Inv+,where (16)Sh+=-μvμh+γhNhμh2μv+Nhμvδh+γh+bβvAhμh+bβvAhδhR4-1μh2μv+μvδh+γh+bkηvβhμh+μvδhγh+bkηvβhδh+γhμhβvb,Vh+=0,Inh+=bkAhηvβhR3-1μh+δhR4-1μhμh2μv+μvδh+γh+bkηvβhμh+μvδhγh+bkηvβhδh+γh,Rh+=bkAhηvβhγhR3-1R4-1μhμh2μv+μvδh+γh+bkηvβhμh+μvδhγh+bkηvβhδh+γh,Av+=-NhkR5-1,Sv+=-R4-1μh2μv+μvδh+γh+bkηvβhμh+μvδhγh+bkηvβhδh+γhμhNh2βhNhμh2μv+Nhμvδh+γh+bβvAhμh+bβvAhδhb,Inv+=-AhNhbβvηvkR3-1μh+δhμvNhμhδh+γh+μhμv+bβvAhμh+δh.It is seen that Sh+,Vh+,Inh+,Rh+,Av+,Sv+, and Inv+ points can have either a positive or a negative value. But as in the MFE and DFE case, an equilibrium point has biological meaning only if it has positive values. Therefore, to ensure that the EE point exists, the condition for Sh+,Vh+,Inh+,Rh+,Av+,Sv+, and Inv+ which has positive value needs to be made; that is,(17)R3=ϕNhμhμh+γhμv+b2βhβvkAhμA+ηvμvAhb2βhβvηvkϕ<1,R4=bβhkμvμA+ηvδh+γh+μhbβhkδh+γh+μhηv+μvμh+δhμh+γhϕ<1,R5=μvμA+ηvϕηv<1.Numerical simulation, using data parameters in Table 1, is performed to show an example of the stability of endemic equilibrium points. The numerical simulation result of the equilibrium point stability can be seen in Table 2.
Numerical example to show existence and stability of equilibrium points for various values of parameters in Table 1 and u1=0,u3=0,u5=1.
4. Basic Reproduction Number4.1. Construction of Basic Reproduction Number
Basic reproduction number (R0) is defined as the expected number of secondary cases from one primary case in a virgin population during the infection period [24]. R0 can be taken from the spectral radius of the next-generation matrix. Please see [25] for further explanation about the construction of the next-generation matrix of the compartmental model in various disease models.
According to our model in system (1) and evaluating it at disease-free equilibrium (DFE) in (5), our next-generation matrix is given by (18)K=0-bβhNh-μv-u4Ahμhu2+u1ξμhμh+u1+θh-Ahμhu2-μh-θhμhμh+u1+θhbβvu5kμv+u4u3+μA+ηv-ηvϕϕμv+u4-γh-μh0.The element of the next-generation matrix K can be interpreted as follows: the number of new infections in jth column is caused by one infection from ith row of K. Please note that i and j for 1 and 2 represent Inh and Inv group. Therefore, for example, K2,1 represent the case that one infected mosquito will produce bβvu5kμv+u4u3+μA+ηv-ηvϕ/ϕμv+u4-γh-μh number of new infected people. On the other hand, K1,2 represent the case that one infected human will produce -(bβh/Nh-μv-u4)Ahμhu2+u1ξ/μhμh+u1+θh-Ahμhu2-μh-θh/μhμh+u1+θh number of new infected mosquitos. Supported by dengue facts that the new infection in human and mosquito population cannot occur from contact between human and human or mosquito and mosquito, we have that K1,1 and K2,2 are equal to 0.
Finding the spectral radius of (18), our basic reproduction number associated with system (1) is given by(19)R0=b2βvu5kβhAhu1ξ+μh+θh+μh+θhξ-1R2R1-1ϕμv+u42γh+μhNhμhμh+u1+θh,with R1 and R2 already defined in the previous section.
Please note that, according to the previous section, this R0 becomes a threshold number to guarantee the existence and local stability of the disease-free equilibrium point (see (5)) and endemic equilibrium point. We find that the disease-free equilibrium point will be locally asymptotically stable when R0<1. This situation will tend the system to possibility equilibrium, that is, stable in MFE if R1<1 and stable in DFE if R1>1. On the other hand, if R0>1, the DFE or MFE will be unstable and the system will tend to endemic equilibrium point.
4.2. Sensitivity Analysis of Basic Reproduction Number
In this subsection, a sensitivity analysis of the basic reproduction number will be performed to find and compare the most sensitive parameters (ui) to determine the value of R0. In the first subsection, we will compare the sensitivity analysis of intervention in the mosquito population (ui for i=3,4,5) and in the next subsection we will compare the sensitivity of intervention in the human population (ui for i=1,2).
4.2.1. Sensitivity Analysis in Mosquito Population
As already stated in the previous section, we include larvicide, fumigation, and mechanical control in our model as u3,u4, and u5, respectively. To find the sensitivity curve of basic reproduction number as shown in Figure 3 for ui and uj, we input all parameters into R0 except ui and uj and then plot its implicit equation.
Sensitivity analysis of R0 with respect to dengue intervention in mosquito population.
In Figure 3(a), a comparison of the efficacy of larvicide and insecticide is performed and we find that insecticide is much more efficacious in reducing R0 than larvicide. In the next figure, Figure 3(b), we find that intervention using mechanical control is more efficacious in reducing R0. Finally, we compare the efficacy of insecticide and mechanical control, and we find that insecticide is much more efficacious in reducing R0. Therefore, from these three figures, we conclude that insecticide is the best way of controlling dengue spread, followed by mechanical control and larvicide, respectively.
4.2.2. Sensitivity Analysis in Human Population
In this subsection, the same procedure is applied to find the sensitivity of R0 in Figure 4. It can be seen that the larger the intervention of vaccination we give, the smaller R0 will be, and reducing R0 with intervention of adult vaccination (u1) is faster than with vaccination of newborns (u2).
Sensitivity analysis of R0 with respect to dengue intervention in human population.
To back up the result given in the two previous figures in Figures 3 and 4, we will determine the value of each intervention in a single-intervention scenario. This means that if we only use vaccination intervention in the model, we will set the other parameters to 0 for u2,u3,u4 and 1 for u5. With this scenario, we need u1=0.04765,u2=293.04,u3=71.363,u4=0.02436, and u5=0.33364 to reduce R0 to 0.99. It can be seen that u4 is the smallest value of intervention to reduce R0 with respect to controlling the spread of dengue.
In the next section, we present some numerical experiments to show the long- and short-term behavior of our model in (1) with respect to the value of various interventions.
5. Numerical Experiment
To perform a numerical experiment in this section, we use the parameters values given in Table 1 and the initial condition given by (20)Sh0=980;Vh0=0;Inh0=20;Rh0=0;Av0=2000;Sv0=980;Inv0=20.As it can be seen from (20), total numbers of human (Sh+Vv+Inh+Rh) and adult mosquito (Sv+Inv) are 1000 in t=0. The number of infected humans and mosquitoes is small to describe the situation when the infection of dengue has just started. Using a value of ui in the previous section to reduce R0 to 0.99, it can be seen that u2 and u3 will not satisfy the condition that ui should be between 0 and 1. For the next simulation, therefore, we will only perform the dynamic behavior of infected groups (Inh,Inv) with respect to intervention of u1,u4, and u5 as shown in Figures 5 and 6 for Inh and Inv, respectively.
Dynamic of infected humans for short-term (a) and long-term (b) intervention.
Dynamic of infected mosquitoes for short-term (a) and long-term (b) intervention.
Figures 5 and 6 show the dynamic of infected humans and mosquitoes in the short term (t∈[0,100]) and long term (t∈[300,500]). It can be seen that, without intervention, the number of infected humans and infected mosquitoes will tend to endemic equilibrium, since R0=2.99>1. After intervention is given until R0=0.99<1 (we take u1=0.04765 or u4=0.02436 and/or u5=0.33364 to represent each simulation), the number of infected humans and mosquitoes will be decreased and pushed to the disease-free equilibrium point. It can also be seen in Figure 5 that intervention of adult vaccination in short-term simulations is the best way to reduce the number of infected humans to the lowest level, rather than other interventions, following this with fumigation and mechanical control interventions, respectively. Unfortunately, for long-term simulations, intervention by fumigation is the best way to reduce the number of infected humans, rather than an adult vaccination strategy. On the other hand, in both short- and long-term simulations, intervention by fumigation is the best way to reduce the number of infected mosquitoes, as shown in Figure 6.
The next simulation is performed to show the efficacy of u4 as the best strategy for long-term intervention in both human and mosquito populations, as shown in Figure 7. It can be seen that an intervention of u4 gradually from 0 to 0.1 will reduce R0 from 2.99 to 0.18. As a consequence, a smaller R0 will reduce the infected population and delay the outbreak for some time. It can also be seen that it needs a proper value of u4 in order that the dynamic of the infected population will never reach outbreak level.
Sensitivity of u4 with respect to the number of infected mosquitoes (a) and infected humans (b).
The last simulation is performed to show the effect of periodic as opposed to constant fumigation intervention. For this purpose, the fumigation is implemented biweekly as a constant (0.02436), and the effect of fumigation will disappear linearly after two weeks, as illustrated in Figure 8. As a result, although constant intervention is much better at significantly reducing the number of infected and susceptible mosquitoes in the mosquito population, as shown in Figure 9, biweekly intervention is only slightly different from constant intervention in reducing the number of infected humans and increasing the number of susceptible humans, as shown in Figure 10. This result indicates that, rather than implementing fumigation constantly, which is more expensive, it would be better to implement a fumigation strategy periodically, since it involves a lower cost.
Biweekly fumigation strategy.
Susceptible (a) and infected (b) mosquito dynamics with constant (solid curve) and biweekly (dash curve) strategy.
Susceptible (a) and infected (b) human dynamics with constant (solid curve) and biweekly (dash curve) strategy.
6. Conclusions
In this article, we have proposed a mathematical model of dengue spread, with various interventions, such as vaccination of adults and newborns in the human population, and larvicide, insecticide, and mechanical control of the mosquito population. Basic reproduction number and basic offspring as the endemic threshold for disease existence and mosquito existence, respectively, have been shown analytically. We find that disease-free equilibrium will be locally asymptotically stable if, and only if, basic reproduction number is smaller than one and will be unstable otherwise.
From sensitivity analysis and backed up with some numerical simulations, we find that fumigation is the best strategy for long-term intervention to reduce the infected populations of both mosquitoes and humans. But for short-term intervention, vaccination of the adult population is the best way to reduce the number of infected people. From numerical simulation, intervention using fumigation, apart from reducing the outbreak, can also delay an outbreak for some period of time.
We also find that although periodic intervention strategy of fumigation cannot reduce the number of infected humans and mosquitoes as efficiently as constant intervention, it is only slightly different. Therefore, if the government has a limited budget, then periodic intervention could be a good option to implement.
For future research, reconstructing the model in this article as an optimal control problem will be considered to show the effectiveness of interventions, not only based on a reduction of the number of infected humans and mosquitoes but also with the lowest cost for intervention purposes.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research is funded by Ministry of Research and Higher Education (Kemenristek Dikti), with PUPT research grant, 2017.
World Health Organization2009World Health Organization. ISBN 9789241547871Education NatureDengue Transmission 2016. http://www.nature.com/scitable/topicpage/dengue-transmission-22399758RodriguesH. S.MonteiroM. T. T.TorresD. F. M.2013Universidade de Aveiro Departamento de Matematica, PortugalAldilaD.NurainiN.SoewonoE.Mathematical Model in Controlling Dengue Transmission with Sterile Mosquito StrategiesProceedings of the AIP Conference2015IndonesiaKementrian KesehatanR. I.Demam Berdarah DengueAugust 201021114Aguirre-ObandoO. A.PietrobonA. J.BonaA. C. D.Navarro-SilvaM. A.Contrasting patterns of insecticide resistance and knockdown resistance (kdr) in Aedes aegypti populations from Jacarezinho (Brazil) after a Dengue Outbreak2016601941002-s2.0-8495828051610.1016/j.rbe.2015.11.009MulyatnoK. C.YamanakaA.NgadinoKonishiE.Resistance of Aedes aegypti (L.) larvae to temephos in Surabaya, Indonesia201243129332-s2.0-84856182166DemingR.Spatial Variation of Insecticide Resistance In The Dengue Vector Aedes aegypti Presents Unique Vector Control Challenges2016967World Health Organization2013Geneva, SwitzerlandWorld Health OrganizationOliveraG.BaurinN.CoudevilleL.The Role Of Indirect Protection In The Assessment Of Dengue Vaccination Impact2016Evora, PortugalPasteurS.First Dengue Vaccine, Dengvaxia, Awarded the Vaccine Breakthrough Of 2016 at the 6th Annual Biopharma Industry Awards Event in Singapore2016, http://www.dengue.info/#overlay=content/first-dengue-vaccine-awarded-vaccine-breakthrough-2016-singapore-0HarapanH.AnwarS.SetiawanA. M.SasmonoR. T.Dengue vaccine acceptance and associated factors in Indonesia: A community-based cross-sectional survey in Aceh20163432367036752-s2.0-8497165470210.1016/j.vaccine.2016.05.026LeeJ.MogasaleV.LimJ. K.CarabaliM.SirivichayakulC.AnhD. D.LeeK.ThiemV. D.LimkittikulK.ThoL. H.VelezI. D.OsorioJ. E.ChanthavanichP.da SilvaL. J.MaskeryB. A.ShepardD. S.A Multi-country Study of the Household Willingness-to-Pay for Dengue Vaccines: Household Surveys in Vietnam, Thailand, and Colombia201596e000381010.1371/journal.pntd.0003810WahyuniT.Berantas DBD dengan Jumatik dan Vaksin DengvaxiaCNN Indonesia, 2016, http://www.cnnindonesia.com/gaya-hidup/20160112181200-255-103823/berantas-dbd-dengan-jumantik-dan-vaksin-dengvaxia/EstevaL.VargasC.Analysis of a dengue disease transmission model199815021311512-s2.0-003252646910.1016/S0025-5564(98)10003-2ZBL0930.92020Zbl0930.92020FengZ.Velasco-HernandezJ. X.Competitive exclusion in a vector-host model for the dengue fever199735552354410.1007/s002850050064MR1479326BaileyN. T. J.1975London, UKGriffinMR0452809DietzK.LudwigD.Transmission and control of arbovirus diseasesProceedings of the Society for Industrial and Applied Mathematics1974Philadelphia, PA, USA104RodriguesH. S.MonteiroM. T.TorresD. F.ZinoberA.Dengue disease, basic reproduction number and control2012893334346MR287856810.1080/00207160.2011.554540Zbl1237.920422-s2.0-81855166473RodriguesH. S.MonteiroM. T. T.M TorresD. F.Vigo-AguiarJ.Modeling and Optimal Control Applied to a Vector Borne DiseaseJuly 20123Murcia, Spain10631070AldilaD.NurainiN.SoewonoE.SupriatnaA. K.Mathematical Model of Temephos Resistance in Aedes aegypti Mosquito PopulationProceedings of theAIP Conference201310.1063/1.48688432-s2.0-84903649602SupriatnaA. K.SoewonoE.van GilsS. A.A two-age-classes dengue transmission model20082161114121MR247475010.1016/j.mbs.2008.08.0112-s2.0-54449084653AldilaD.GötzT.SoewonoE.An optimal control problem arising from a dengue disease transmission model2013242191610.1016/j.mbs.2012.11.014MR3033464DiekmannO.HeesterbeekJ. A. P.2000Chichester, UKJohn Wiley & SonsMR1882991DiekmannO.HeesterbeekJ. A. P.RobertsM. G.The construction of next-generation matrices for compartmental epidemic models201074787388510.1098/rsif.2009.03862-s2.0-77953160649