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

In the case of Dengue transmission and control, the interaction of nature and society is captured by a system of diﬀerence 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 diﬀerence equations: . . The constraints have socio-ecological meaning. The initial conditions are such that 0 ≤ P 0 ≤ 1 , ( M 0 , H 0 , A 0 ) ≥ (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


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, usually Aedes 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: This discrete system links the proportion of the infected people (  ), mosquitoes (  ), habitats (  ), and awareness   .The initial conditions are such that 0 ≤  0 ≤ 1, ( 0 ,  0 ,  0 ) ≥ (0, 0, 0), the parameters , , ,  ∈ (0, 1), and the parameters , , , and  are positive.The current system represents a modification of the system in [5].
The first equation describes the proportion of infected people (between 0 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 parameter  is related to the behavior of infected mosquitoes, and it can be viewed as a transmission rate.An explanation of the term (1 −  −  ) goes as follows.If  represents the probability that a mosquito transmits the infection, then 1− is the probability that it does not transmit the infection.Therefore, (1 − )   will be the probability that   mosquitoes do not transmit the infection.One can rewrite We denote  = − ln(1 − ) > 0.
By using a series of transformations, one can rescale the parameters , . and  in (1).We use the following changes of variables,   = (1/)  (in the second and first equation),   = (1/)  (in the third and fourth equation), and   = ℎ  .These transformations will not change the nature of parameters , , , and , 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): 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.
Proof.First equation of system (4) gives and then for any positive number   , there exists  sufficiently large, such that, for all  ≥ , Making use of (5) in the fourth equation, we get Since 0 <  < 1 we obtain lim sup  → ∞   ≤ /(1 − )(1 − ) and then for any positive number   , there exists  sufficiently large, such that, for all  ≥ , The third equation of ( 4) yields  +1 ≤   + 1 which combined with 0 <  < 1 gives lim sup  → ∞   ≤ 1/(1 − ).Thus, for any positive number  ℎ , there exists  sufficiently large, such that, for all  ≥ , Finally, (8) and Thus for any positive number   , there exists  sufficiently large, such that, for all  ≥ , Some notations that will be used throughout the paper are, in order, lim sup

Equilibria
Clearly, is an equilibrium point of system (4) for all the values of the parameters.
(2) Assume that  > (1 − )(1 − ); then there are two equilibrium points, namely, the degenerate one and a positive one denoted by (, , , ).The positive equilibrium can take the form Proof.The equilibrium solutions verify the system The fourth equation in the above system gives Solving for  in the third equation yields Combining (15) with the second equation of system ( 13) produces Replacing (14) in first system equation and multiplying by  to both sides, Since (1 −  + ) ̸ = 0 and (1 −  − ) ̸ = 0, (16) can be written in the form Equation ( 17) gives Notice that where function () has the property that it is an increasing function, first order derivative for  ∈ (0, ∞).Set the real valued functions We have that Φ 1 () is an increasing function in , Φ 1 (0 + ) = 1 and Φ 1 () = ().

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: The Jacobian evaluated at the equilibrium point (, , , ) has the form Using the third equilibrium equation, 1 +  = (1 + ).Thus, −(1 + )/(1 + ) 2 = −/(1 + ).The characteristic equation associated with (, , , ) is given by the fourth order polynomial: One can look at the characteristic equation in the form where In the region of existence of positive equilibrium point,  > (1 − )(1 − ), 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  > (1 − )(1 − ) (see Figure 3).
Then using the last equation in the system, we conclude Using that 1 −  −  <   in the second equation of the reduced system yields Thus Dividing by   > 0 to both sides one obtains which implies that (1− )(1− ) ≤  (hence the contradiction).
Thus   = 0. First equation in the reduced system yields the inequality or further  +1 ≤   +   .Passing to the limit one has Dividing by (1 − ) > 0 the above yields   ≤ 0 which in combination with   ≥ 0 gives   = 0. Using the inequality in the third equation It follows = 0 follows easily.

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 ( = 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 ( = 0.97, Figure 2) leads to a complete decrease in the proportion of infected people and the solutions stabilize.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 −  −  )(1 −   ).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  0 has an impact on the threshold value or not.
The average number of mosquitoes per breeding site (parameter ) was estimated to be 9.5, ranging from 3 to 30, in field studies (see [2]).We used  = 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  (high pollution level such as new empty cans and tires that collect water), the memory parameter  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.
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: 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: 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 () and habitats () 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][9][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][13][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   , and keep the population control on   .To assess the effect of habitat control only, through citizens intervention, the equation will keep its intervention parameters as such,  +1 =   /(1 +   ) + /(1 +   ), for example.We believe that system (44) is not only useful biologically but also interesting mathematically.Both systems (43) and (44) possess bounded solutions.
for their comments that helped improve the quality of their paper.They also thank program Hope at Harvard University for supporting this paper.