Global Stability of an Anthrax Model with Environmental Decontamination and Time Delay Steady

Anthrax occurs worldwide and is associated with sudden death of cattle and sheep. This paper considers an epidemic model of anthrax in animal population, only. The susceptible animal is assumed to be infected, only, through ingestion of the disease causing pathogens. The proposed model incorporates time delay and environmental decontamination by humans. The time delay represents the period an infected animal needs to succumb to anthrax-related death. By constructing suitable Lyapunov functionals, we demonstrate that the global dynamics of thismodel fully hinges onwhether the associated reproductive number is greater or less than unity. The effectiveness of environmental decontamination on eradication of anthrax in the community is explored through the reproductive number.


Introduction
Anthrax is an acute zoonotic disease caused by the sporeforming bacterium Bacillus anthracis [1].The disease can infect all warm-blooded animals, including humans.However, ruminants, particularly cattle and sheep, are more susceptible [2].The disease is associated with death of cattle and sheep, such that very few livestock producers or veterinaries have witnessed the disease or its signs [2].Outbreaks often occur when livestock are grazing on neutral or slightly alkaline soils and have been exposed to the spores through ingestion of contaminated soil in endemic areas or forage is sparse because of overgrazing or drought or when soil has been disturbed by digging or other human activities [2].Anthrax is one of the most dramatic diseases affecting wild animal in Africa.For instance, within Etosha National Park, anthrax has been attributed to the death of a variety of herbivores, from elephants to ostriches [3].
The use of mathematical models to explore the spread and control of infectious diseases has proved to be an important tool for the scientists and epidemiologists.In 1983, Hahn and Furniss [4] constructed a deterministic model to explore the impact of environmental contamination and animal carcasses on driving anthrax epidemic.Within their framework they assumed that upon infection infected animals die immediately and there are no births or deaths from causes other than anthrax related [4].More, recently, Friedman and Yakubu [3] extended the model of Hahn and Furniss [4] to study the effects of anthrax transmission, carcass ingestion, carcass induced environmental contamination, and migration rates on the persistence or extinction of animal population.Their work revealed among others that decreasing the levels of carcass ingestion by removal of carcasses in game reserves may not always lead to a reduction in anthrax cases.
Although anthrax infection in cattle is regarded often as a fatal disease with no signs [2], biologically the infection should take a period before an infected animal succumbs to anthrax-induced mortality, and this period may play an important role in controlling anthrax outbreaks.So incorporating this reason we extend the work of Hahn and Furniss [4] to include a fixed delay.Apart from the delay, our new model incorporates the role of human effort on decontamination of the environment through the destruction of animal carcass and disinfection of the ground or area that contained the carcass by adding lime.The main goal of this paper is to explore the effectiveness of environmental decontamination on controlling anthrax outbreaks and to study the global stability of the anthrax model with time delay.

Model Formulation and Analysis
A framework to assess the dynamics of anthrax in livestock is proposed.Hahn and Furniss's [4] deterministic model of anthrax dynamics provides the starting point for our discussion on the effectiveness of environmental decontamination on controlling the spread of anthrax.The framework is governed by the following system of nonlinear ordinary differential equations: ( The first equation describes the dynamics of the susceptible animals.The parameter  denotes the entrance of new animals through birth and they are assumed to be susceptible to the disease;  is the disease transmission rate which occurs when a susceptible animal ingests the free-living spores while grazing in an endemic area;  denote the permanent exit of the animals due to natural causes or other reasons not related to anthrax infection.Upon infection, we assume that an infected animal dies without displaying clinical signs of the disease.Thus, the second equation captures the dynamics of the carcasses of animals that may have succumbed to anthraxinduced death;  denote the decay rate of carcasses which is assumed to be constant;  denote the rate of disinfection or decontamination of infected areas through the removal or destruction of animal carcasses and adding of lime to the ground where a decomposing animal carcass would have been identified.The third equation describes dynamics of free-living spores or pathogens;  denote the rate at which carcasses of animals that would have died of anthrax and not properly destroyed shed the bacteria into the environment;  −1 is the average life span of free-living pathogens.The length of survival of anthrax free-living spores in the environment is estimated to be around 200 years [5]; thus  = 1/(365×200) = 0.000014 day −1 .
For biological reasons we will study the dynamics of system (1) in the closed set where R 3 + denotes the nonnegative cone of R 3 including its lower dimensional faces and  ≤ min(, , ).It can easily be verified that Γ is positively invariant with respect to (1).

Equilibria and the Reproductive Number.
In the absence of anthrax infection in the community, system (1) admits an infection-free equilibrium given by E 0 = (/, 0, 0).The Jacobian matrix of system (1) evaluated about It is clear that − is an eigenvalue of matrix .The other two are determined from matrix .
The trace and determinant of  are, respectively, given by tr () = − (2 +  + ) , ] . ( Let the reproductive number of system (1) be Since tr() < 0 and det() > 0 (when R  < 1) it follows that the equilibrium point E 0 is locally asymptotically stable if and only if R  < 1.
Biologically, the term R  represents the average number of new anthrax cases generated when susceptible animals ingest disease causing pathogens.It gives a measure of the power of anthrax to invade the cattle population in the presence of environmental decontamination.
Direct calculation can easily show that system (1) has two possible equilibrium points, namely, the infection-free E 0 and the endemic equilibrium point E * , given by From ( 7) it is evident that E * makes biological sense whenever R  > 1.

Effectiveness of Environmental Decontamination.
In this section we explore the strength of environmental decontamination on reducing or eliminating new anthrax infections.
In the absence of environmental decontamination in the community ( = 0) the average number of new anthrax infections generated through spores ingestion is modeled by the threshold quantity R 0 which is given by To explore the effectiveness of environmental decontamination on controlling new anthrax infections we define the efficacy function From ( 9) it is clear that the effectiveness of environmental decontamination depends on the length of survival of freeliving pathogens in the environment and the decay rate of undestroyed carcasses of animals that have succumbed to anthrax infection.In Figure 1 we numerically explore the effectiveness of different environmental decontamination levels.Numerical results in Figure 1 suggest that environmental decontamination rate 0.05 day −1 can be effective in reducing the generation of new anthrax infections by 50%.Further, we note that environmental decontamination rate greater than or equal to 0.01 day −1 can be effective in attaining 100% control on the spread of anthrax in the community.

Anthrax Model with Time Delay.
Although anthrax infection in cattle is often a fatal disease with no clinical signs displayed by an infected animal, it is worth noting that the infection should take a period before an infected animal succumbs to anthrax-induced death and the size of this period may play an important role in controlling the outbreak of this disease.So incorporating this reason we introduce a time delay into system (1) to represent the aforementioned period: (10) System (10) satisfies the following initial conditions: where   () for  = 1, 2, 3 denote the nonnegative continuous functions on  ∈ [−, 0].All the parameters in system (10) are the same as in system (1) except for the positive constant  which represents the length of the delay.It can be verified that Γ is positively invariant with respect to system (10).We denote by Γ and Γ the boundary and the interior of Γ in R 3 + , respectively.With the same motivation as before, we introduce the reproductive number of differential-delay model (10) which is given by a similar expression 2.3.1.Anthrax-Free Equilibrium Point and Its Stability.Direct calculation shows that system (10) has the same disease-free equilibrium E 0 = (/, 0, 0) as in system (1).The Jacobian matrix of system (10) about E 0 is From (13) it follows that the characteristic equation is Since  = − < 0 is a root of ( 14), we only need to consider where Equation ( 14) with  = 0 is If  = 0, we have already deduced that all roots of system (15) have negative real parts when R  < 1, and only one root of (15) has positive real part when R  > 1.
For  > 0, we assume that  =  ( > 0) is a root of (15).This is the case if and only if −  2 +  11  +  12 +  11 (cos  −  sin ) = 0. (18) Separating the real and imaginary parts yields Adding up the squares of both equations, we obtain Let  =  2 , so that (21) reduces to Substituting for  11 ,  12 , and  11 gives Substituting R  into (22) gives Whenever R  < 1, (23) has two roots which have a positive product implying that they are complex or they are real but they have the same sign.In addition, they have negative sum which implies that they are either real and negative or complex conjugates with negative real parts.Consequently, (23) does not have positive real roots which lead to the conclusion that there is no  such that  is a solution of ( 14).Therefore, it follows from Lemma 2.4 of Ruan and Wei [6] that the real parts of all eigenvalues of characteristic equation ( 14) are negative for all values of the delay  ≥ 0. Thus, we have the following theorem.
Theorem 1.The equilibrium point E 0 of system ( 10) is locally asymptotically stable for all time delay  ≥ 0 if R  < 1.
In the next theorem we establish the global stability of the infection-free equilibrium for system (10).

Endemic Equilibrium Point and Its Stability
Theorem 3. Suppose R  > 1; then the equilibrium point E 0 is unstable.
Proof.From the discussion in Section 2.3.1, we have deduced that the characteristic equation associated with E 0 is given by ( 14), and we only need to consider (15).From the above computations it can easily be verified that if  = 0, (15) has a positive root when R  > 1.Now, let  ( > 0) be a root of (15).Then solving (23) gives For R  > 1 it is evident that (26) makes sense while (27) is meaningless.Define Then ( 14) has a pair of purely imaginary roots ± when  =   and has no roots appearing on the imaginary axis when  ̸ =   for  = 0, 1, . . . .Further, let () = () + () be the root of (14) satisfying (  ) = 0 and (  ) =  + .
Differentiating both sides of ( 14) with respect to  we get By applying Lemma 2.4 in Ruan and Wei [6] and observing that (15) has a positive real root when  = 0, we obtain that characteristic equation ( 14) has a positive root at least for all  ≥ 0. Therefore, the equilibrium point E 0 of system (10) is unstable when R  > 1.This completes the proof.Now, we investigate the effect of the time delay on the local stability of E * .The characteristic equation of system (10) at the endemic equilibrium E * takes the following form: where When  = 0, (32) becomes By the Hurwitz criteria, all the roots of (34) have only negative real parts.Thus E * is locally asymptotically stable when  = 0. Now, we consider the case  ̸ = 0.If  ( > 0) is a solution of (32), we have − where Hence, if R  > 1, (32) has no positive roots.Accordingly, if R  > 1, the endemic equilibrium E * exists and is locally asymptotically stable for all  ≥ 0. Further, using a Lyapunov functional, we can obtain the following theorem.
Theorem 4. Whenever R  > 1, then the unique equilibrium E * of system (10) is globally asymptotically stable in Γ for all  ≥ 0.

Concluding Remarks
Anthrax epidemic is now recognized among other factors as the leading cause of species extinctions.In this paper, we propose and analyze an anthrax epidemic model.The model focuses on anthrax transmission in animal population only.Our model is based on the assumption that the disease transmission occurs only when a susceptible animal ingests the disease causing pathogen from contaminated soil in endemic areas when forage is sparse because of overgrazing or drought or when soil has been disturbed by digging or excavations.The proposed model incorporates time delay and environmental decontamination effort.The time delay represents the period that is needed for an infected animal to succumb to anthrax-induced death.We determine the reproductive number and establish that the global dynamics depends on whether the reproductive number is greater than one.