# Codynamics of Four Variables Involved in Dengue Transmission and Its Control by Community Intervention: A System of Four Difference Equations

Authors: T. Awerbuch-Friedlander; Richard Levins; M. Predescu
Journal: Discrete Dynamics in Nature and Society
 (2014)

---

## Abstract

In the case of Dengue transmission and control, the interaction of nature and society is captured by a system of difference equations. For the purpose of studying the dynamics of these interactions, four variables involved in a Dengue epidemic, proportion of infected people (P), number of mosquitoes involved in transmission (M), mosquito habitats (H), and population awareness (A), are linked in a system of difference equations: Pn+1=aPn+1-e-iMn1-Pn, Mn+1=lMne-An+bHn1-e-Mn, Hn+1=cHn/(1+pAn)+1/(1+qAn), and An+1=rAn+fPn, n=0,1,…. The constraints have socioecological meaning. The initial conditions are such that 0≤P0≤1,(M0,H0,A0)≥(0,0,0), the parameters l,a,c,r∈(0,1), and the parameters f,  i,  b, and p are positive. The paper is concerned with the analysis of solutions of the above system for p=q. We studied the global asymptotic stability of the degenerate equilibrium. We also propose extensions of the above model and some open problems. We explored the role of memory in community awareness by numerical simulations. When the memory parameter is large, the proportion of infected people decreases and stabilizes at zero. Below a critical point we observe periodic oscillations.


---

## Body

## 1. Introduction

The response to an epidemic is triggered by awareness of a coming epidemic or by an existing one. The response is aimed at reducing the incidence of the actual disease. In the case of Dengue fever, the disease is caused by a virus that is transmitted by the bite of the mosquito, usuallyAedes aegypti. The mosquitoes deposit eggs in small containers of water. These hatch to produce larvae. Some transform into pupae and then adult mosquitoes. The breeding sites may be ephemeral, such as water in an empty beer can or used tire, an animal drinking trough near a human habitation, or in-doors stored water in large containers [1, 2].The information about a Dengue epidemic can come from the number of reported cases of Dengue, the abundance of mosquitoes, or the numbers of breeding sites for mosquitoes, or some other indicator such as rainfall that predicts breeding sites. The information triggers consciousness, and the response can be either individual and/or one at the community level. In previous work, we studied the dynamics of a discrete time system in which we modeled the awareness as a factor that is triggered by the formation of potential breeding sites and the response was aimed at eliminating them. The system was studied by a pair of two difference equations [3].By expanding the model to introduce an ongoing educational program, our new model predicted that high consciousness over time kept the number of breeding sites low [4]. In a study with three difference equations, we study a system in which the information is related to the number of adult mosquitoes. The more mosquitoes, the greater the awareness of the population, and this awareness leads to action to reduce the mosquito population by controlling breeding sites [5]. This population awareness is prompted and dissipates at a rate determined by the abundance of mosquitoes, similar to a birth and death process. The dynamics then is that mosquitoes are produced when adult females locate breeding sites and deposit eggs which develop into adult mosquitoes, and mosquitoes die at a rate depending on their own biology and environmental conditions as a result of control measures implemented as awareness rises. Thus the pair of variables, mosquitoes and awareness, are linked in a negative feedback loop in a system of equations where decay due to control was modeled with a rational fractional term at the environmental level. With another system of three difference equations we have explored an intervention by spraying mosquitoes [6]. The change in the spraying parameter resulted in almost periodic behavior and fluctuations in the populations of mosquitoes. Simulations show that alertness in consciousness, by keeping the memory parameter of previous week high, has an impact on the behavior of solutions and implicitly on the number of mosquitoes. When the memory parameter is high, there will be a steady decrease in the number of mosquitoes. The present study builds upon the previous models. We present a system of four difference equations, with the proportion of infected people as an additional variable that prompts consciousness:
(1)Pn+1=aPn+1-e-iMn1-Pn,Mn+1=lMne-gAn+bHn1-e-sMn,Hn+1=cHn1+pAn+d1+pAn,An+1=rAn+fPn,mmmmmmmmmnmmmn=0,1,….
This discrete system links the proportion of the infected people (Pn), mosquitoes (Mn), habitats (Hn), and awareness An. The initial conditions are such that 0≤P0≤1,  (M0,H0,A0)≥(0,0,0), the parameters l,a,c,r∈(0,1), and the parameters f,  i,  b, and p  are positive. The current system represents a modification of the system in [5].The first equation describes the proportion of infected people (between0 and 1). They prompt consciousness, while the intervention is against mosquitoes and perhaps habitats. In the relationships among variables, the awareness is prompted by the proportion of sick people. The control of both adult mosquitoes by spraying and habitats is carried out by community intervention.The parameteri is related to the behavior of infected mosquitoes, and it can be viewed as a transmission rate. An explanation of the term (1-e-iMn) goes as follows. If Q represents the probability that a mosquito transmits the infection, then 1-Q is the probability that it does not transmit the infection. Therefore, (1-Q)Mn will be the probability that Mn mosquitoes do not transmit the infection. One can rewrite
(2)1-QMn=eln⁡(1-Q)Mn=eMnln⁡(1-Q).
We denote i=-ln⁡(1-Q)>0.One can observe that if0≤P0≤1 then P1≤1. This is true because
(3)P1=aP0+1-e-iM01-P0≤aP0+1-P0≤1.
It follows by induction that 0≤Pn≤1. Also, if (M0,H0,A0)≥(0,0,0) then (Mn,Hn,An)≥(0,0,0). Thus, we have that (P0,M0,H0,A0)≥(0,0,0,0) then (Pn,Mn,Hn,An)≥(0,0,0,0).By using a series of transformations, one can rescale the parametersg, s. and d in (1). We use the following changes of variables, Mn=(1/s)mn (in the second and first equation), An=(1/g)an (in the third and fourth equation), and Hn=dhn. These transformations will not change the nature of parameters a, c, l, and r, as these remain between 0 and 1. Thus, after relabeling the variables and parameters, one can work with a simplified system of equations as below (it is this system that will get analyzed in the next sections):
(4)Pn+1=aPn+1-e-iMn1-Pn,Mn+1=lMne-An+bHn1-e-Mn,Hn+1=cHn1+pAn+11+pAn,An+1=rAn+fPn,mmmmmmmmmmmin=0,1,….
In the sequel, we look at boundedness properties, local and global asymptotic stability of equilibria. Numerical simulations, open problems, and further directions of improvement will be mentioned.


## 2. Boundedness of Solutions

Lemma 1.
Let{Pn,Mn,Hn,An}n≥0 be a positive solution of system (4). Parameters are such that 0<l<1, 0<a<1,  0<c<1, and 0<r<1. Then limsup⁡n→∞Pn≤1/(1-a), limsup⁡n→∞Mn≤b/(1-l)(1-c), limsup⁡n→∞Hn≤1/(1-c), and limsup⁡n→∞An≤f/(1-a)(1-r).Proof.
First equation of system (4) gives Pn+1≤aPn+(1-Pn)≤aPn+1. Thus limsup⁡n→∞Pn≤1/(1-a) and then for any positive number ϵp, there exists N sufficiently large, such that, for all n≥N,
(5)Pn+1<11-a+ϵp.
Making use of (5) in the fourth equation, we get
(6)An+1≤rAn+f1-a+ϵp.
Since 0<r<1 we obtain limsup⁡n→∞An≤f/(1-a)(1-r) and then for any positive number ϵa, there exists N sufficiently large, such that, for all n≥N,
(7)An+1<f1-a1-r+ϵa.
The third equation of (4) yields Hn+1≤cHn+1 which combined with 0<c<1 gives limsup⁡n→∞Hn≤1/(1-c). Thus, for any positive number ϵh, there exists N sufficiently large, such that, for all n≥N,
(8)Hn+1<11-c+ϵh.
Finally, (8) and Mn+1≤lMn+bHn≤lMn+b/(1-c)+ϵh produce limsup⁡n→∞Mn≤b/(1-l)(1-c). Thus for any positive number ϵm, there exists N sufficiently large, such that, for all n≥N,
(9)Mn+1<b1-l1-c+ϵm.Some notations that will be used throughout the paper are, in order,(10)limsup⁡n→∞Pn=SP,liminf⁡n→∞Pn=IP,limsup⁡n→∞Hn=SH,liminf⁡n→∞Hn=IH,limsup⁡n→∞Mn=SM,liminf⁡n→∞Mn=IM,limsup⁡n→∞An=SA,liminf⁡n→∞An=IA.


## 3. Equilibria

Clearly,(11)0,0,11-c,0
is an equilibrium point of system (4) for all the values of the parameters.Lemma 2.
(1) Assume thatb≤(1-c)(1-l). Then the degenerate equilibrium (0,0,1/(1-c),0) is the only equilibrium point.
(2) Assume thatb>(1-c)(1-l); then there are two equilibrium points, namely, the degenerate one and a positive one denoted by (P¯,M¯,H¯,A¯). The positive equilibrium can take the form
(12)∑f1-e-iM¯1-r2-a-e-iM¯1-e-iM¯2-a-e-iM¯,M¯,11-c+pf1-e-iM¯/1-r2-a-e-iM¯,f1-e-iM¯1-r2-a-e-iM¯∑f1-e-iM¯1-r2-a-e-iM¯.Proof.
The equilibrium solutions verify the system(13)P¯=aP¯+1-e-iM¯1-P¯,M¯=lM¯e-A¯+bH¯1-e-M¯,H¯=cH¯1+pA¯+11+pA¯,A¯=rA¯+fP¯.
The fourth equation in the above system gives
(14)P¯=(1-r)A¯f.
Solving for H¯ in the third equation yields
(15)H¯=11-c+pA¯.
Combining (15) with the second equation of system (13) produces
(16)1-le-A¯M¯=b1-e-M¯1-c+pA¯.
Replacing (14) in first system equation and multiplying by f to both sides,
(17)1-a1-rA¯=1-e-iM¯f-1-rA¯.
Since (1-c+pA¯)≠0 and (1-le-A¯)≠0, (16) can be written in the form
(18)M¯1-e-M¯=b1-c+pA¯1-le-A¯.
Equation (17) gives
(19)A¯=f1-e-iM¯1-r2-a-e-iM¯.
Notice that (2-a-e-iM¯)>0. Set
(20)w(M¯)=f(1-e-iM¯)(1-r)(2-a-e-iM¯).
Notice that
(21)A¯=wM¯,
where function w(M) has the property that it is an increasing function, first order derivative
(22)w′(M)=f(1-a)ie-iM(1-r)2-a-e-iM2>0
for M∈(0,∞). Set the real valued functions
(23)Φ1M=M1-e-M,gA=b1-c+pA1-le-A.
We have thatΦ1(M) is an increasing function in M, Φ1(0+)=1 and Φ1(M¯)=g(A¯).
From the above,(24)gA¯=gwM¯=g∘wM¯,
where we denote Φ2 as
(25)Φ2(M)=(g∘w)(M)=g(w(M)).
Function Φ2 is decreasing. Let M1<M2. Since function w is increasing, one has w(M1)<w(M2). But g is a decreasing function and
(26)Φ2M1=g∘wM1=gwM1>gwM2=Φ2M2.
Using that w(0+)=0, we have that
(27)Φ20+=b1-c1-l.
For (18) to have a unique solution (and thus system to have a unique solution), one must have Φ1(0+)<Φ2(0+) or equivalently 1<b/(1-c)(1-l) and the proof ends.


## 4. Stability of Equilibrium Points

Next we are concerned with the local and global asymptotic stability of equilibrium points. Notations for our map are as follows:(28)Pn+1=ΘPn,Mn,Hn,AnwithΘP,M,H,A=aP+1-e-iM1-P,Mn+1=gPn,Mn,Hn,Anwithg(P,M,H,A)=lMe-A+bH1-e-M,Hn+1=hPn,Mn,Hn,AnwithhP,M,H,A=cH1+pA+11+pA,An+1=ΦPn,Mn,Hn,AnwithΦP,M,H,A=rA+fP.
The Jacobian evaluated at the equilibrium point (P¯,M¯,H¯,A¯) has the form
(29)JP¯,M¯,H¯,A¯=a-1-e-iM¯1-P¯ie-iM¯000le-A¯+bH¯e-M¯b-be-M¯-lM¯e-A¯00c1+pA¯-p1+cH¯1+pA¯2f00r.
Using the third equilibrium equation, 1+cH¯=H¯(1+pA¯). Thus, -p(1+cH¯)/(1+pA¯)2=-pH¯/(1+pA¯). The characteristic equation associated with (P¯,M¯,H¯,A¯) is given by the fourth order polynomial:
(30)a-1-e-iM¯-λle-A¯+bH¯e-M¯-λ×c1+pA¯-λ][r-λ-1-P¯ie-iM¯f-pH¯b1-e-M¯1+pA¯mmmmmmmmmm+lM¯e-A¯c1+pA¯-λ-pH¯b1-e-M¯1+pA¯=0.
One can look at the characteristic equation in the form
(31)λ4-(A1+A2+A3+A4)λ3+(A1A2+A1A3+A1A4+A2A3+A2A4+A3A4)λ2-A1A2A3+A1A2A4mmm+A1A3A4+A2A3A4+A5A7λ+A1A2A3A4+A5A6+A3A5A7=0,
where
(32)A1=a-1-e-iM¯,A2=le-A¯+bH¯e-M¯,A3=c1+pA¯,A4=r,A5=-1-P¯ife-iM¯,A6=-pH¯b1-e-M¯1+pA¯,A7=lM¯e-A¯.
In the region of existence of positive equilibrium point, b>(1-c)(1-l), the values of parameters for which the roots of the fourth order polynomial are inside unit disc generate a locally asymptotically stable equilibrium point. The positive equilibrium point is not always locally asymptotically stable in the region b>(1-c)(1-l) (see Figure 3).The following theorem about the degenerate equilibrium point (Figure1) (0,0,1/(1-c),0) holds.Figure 1
The above graph is generated with parameter valuesa=0.5,  p=5,  b=0.48,  c=0.04,  l=0.5,  i=0.5,  r=0.5, and f=20 (the parameters are placed in the region where b=(1-c)(1-l)). One can see that the solutions converge to the degenerate equilibrium point.Theorem 3.
Assume thatb<(1-c)(1-l). Then (0,0,1/(1-c),0) is globally asymptotically stable.Proof.
WhenP¯=0, M¯=0, H¯=1/(1-c), and A¯=0, the Jacobian becomes
(33)J0,0,11-c,0=ai000l+b1-c0000c-cp1-c-pf00r
with the characteristic equation a polynomial that factors into
(34)a-λl+b1-c-λc-λr-λ=0.
Three of the roots, namely, λ1=a, λ2=c, and λ4=r are less than 1 and if l+b/(1-c)<1 (or b<(1-c)(1-l)) then the degenerate equilibrium is a sink and thus locally asymptotically stable. It remains to be shown that this equilibrium is a global attractor. We offer a proof by contradiction as in [5]. Let us suppose SP>0 and SM>0. Then using the last equation in the system, we conclude
(35)SA≤rSA+f1-a.
Using that 1-e-Mn<Mn in the second equation of the reduced system yields
(36)Mn≤lMn+bHnMn.
Thus
(37)SM≤lSM+bSHSM.
Dividing by SM>0 to both sides one obtains
(38)1-lb≤SH≤11-c
which implies that (1-l)(1-c)≤b (hence the contradiction). Thus SM=0.
First equation in the reduced system yields the inequality(39)Pn+1≤aPn+iMn1-Pn
or further Pn+1≤aPn+iMn. Passing to the limit one has
(40)SP≤aSP+iSM=aSP.
Dividing by (1-a)>0 the above yields SP≤0 which in combination with SP≥0 gives SP=0.
Using the inequality in the third equation(41)IH≥cIH1+pSA+11+pSA≥cIH+1.
It follows IH≥1/(1-c). But SH≤1/(1-c)≤IH and thus SH=IH=1/(1-c).
S
A
=
0 follows easily.


## 5. Conclusions and Open Problems

The global asymptotic stability of the degenerate equilibrium was investigated (but the global asymptotic stability of the positive equilibrium remains an open problem that is worth investigating mathematically). An interesting result pertains to the role that the memory plays in controlling the epidemic. We observed oscillatory behavior for marginally low memory parameter values (r=0.5, Figure 3), meaning that the population might recover only for a short period of time and then get periodically infected. Computer simulations indicate that high awareness (r=0.97, Figure 2) leads to a complete decrease in the proportion of infected people and the solutions stabilize.Figure 2
The above graph is generated with parameter valuesa=0.5,  p=0.5,  b=5,  c=0.6,  l=0.5,  i=0.5,  r=0.97, and f=20 (the parameters are placed in the region where b>(1-c)(1-l)). One can see that the solutions display convergence to the equilibrium for big values of r.Figure 3
The above graph is generated with parameter valuesa=0.5,  p=0.5,  b=5,  c=0.6,  l=0.5,  i=0.5,  r=0.5, and f=20 (the parameters are placed in the region where b>(1-c)(1-l)). One can see that the solutions display oscillatory behavior for smaller values of r.Simulations done with various parameter values seem to suggest that the memory parameter has a threshold below which there are oscillations and above which it exhibits the equilibria, leading to the extinction of the infection. This is consistent with other findings from studies specifically designed to discover thresholds (see [7]). In [7] the authors considered the rate of contact between susceptible people and infectious vectors, a component captured in our system in the first equation by the term (1-e-iMn)(1-Pn). Their study reports that they were surprised to discover that the size of the viral introduction“was not seen to significantly influence the magnitude of the threshold.” In our future study, we shall focus on finding the memory parameter threshold value that leads to the extinction of the infection and to testing whether changing the initial conditions of the proportion of infected people P0 has an impact on the threshold value or not.The average number of mosquitoes per breeding site (parameterb) was estimated to be 9.5, ranging from 3 to 30, in field studies (see [2]). We used b=5 (see Figures 2 and 3), a value within the range suggested by field studies in the aforementioned reference. Computer simulations on system (1) indicate that it is possible that for large values of parameter d (high pollution level such as new empty cans and tires that collect water), the memory parameter r alone may not be sufficiently strong enough to eliminate the infection from the population, and the infection might equilibrate at levels higher than zero. In future work we shall explore the relationship between environmental pollution and the memory that creates awareness in the community.In this section we also want to bring attention to some extensions and open problems related to system (1). An interesting question to be analytically investigated in a further study is the global asymptotic stability of nondegenerate equilibrium of system (1) especially in the case when the system incorporates different parameters that measure the sensitivity of surviving habitats to communal awareness and individual awareness (hence p≠q). Thus, in this case the third equation reads
(42)Hn+1=cHnh1pAn+dh2qAn.
Based on biological considerations, one can take h1(·) and h2(·) as decreasing functions, h1,h2∈C1((0,∞)→(0,1]) with properties (i) h1(0)=1 and h2(0)=1 and (ii) lim⁡y→∞h1(y)=0 and lim⁡y→∞h2(y)=0. Two most used examples of such functions (used in the previous work, [3]) are for instance h1(y)=1/(1+py) and h2(y)=1/(1+qy). Thus, an open problem that we want to pose here refers to the study of the existence and global asymptotic stability of the positive equilibrium of the general system:
(43)Pn+1=aPn+1-e-iMnPn1-Pn,Mn+1=lMne-gAn+bHn1-e-sMn,Hn+1=cHnh1(pAn)+dh2qAn,An+1=rAn+fPn,mmmmmmmiimmmmin=0,1,….
Mathematical models may serve at designing policy interventions and provide a better understanding of phenomena at study [3]. Because at times interventions are implemented when consciousness is prompted by an increase in the incidence of sick people, one can work with the original system in a form as such:
(44)Pn+1=aPn+1-e-iMnPn1-Pn,Mn+1=lMne-gAn+bHn1-e-sMn,Hn+1=cHn1+pAn+d1+qAn,An+1=rAn+fPn,mmmmmmmmiimmmn=0,1,….
The first equation describes the proportion of infected people in the population (between 0 and 1). Proportion of sick people is assumed to prompt consciousness, while the intervention is against mosquitoes and perhaps habitats. The control of both, adult mosquitoes (M) and habitats (H) where mosquitoes lay their eggs, is carried out by spraying and community intervention by reducing breeding sites. One may use this system (system (44)) to compare a few control strategies, where increase in the proportion of infected people is linked to consciousness. Insecticide spraying is a common method in mosquito control despite its many disadvantages; and new ones are continuously being developed and tested [8–10]. In the long run the mosquitoes become resistant and the insecticide ineffective [11]; it poses serious risks to humans and the environment [10, 12–14]. In order to assess the effect of insecticide spraying without habitat management, the equations are modified so that we eliminate the rational control on Hn, and keep the population control on Mn. To assess the effect of habitat control only, through citizens intervention, the equation will keep its intervention parameters as such, Hn+1=cHn/(1+pAn)+d/(1+qAn), for example. We believe that system (44) is not only useful biologically but also interesting mathematically. Both systems (43) and (44) possess bounded solutions.



---

*Source: 101965-2014-11-17.xml*