Discrete Dynamics in Nature and Society
Volume 2015 (2015), Article ID 573146, 6 pages
http://dx.doi.org/10.1155/2015/573146

Research Article

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

Steady Mushayabasa

Department of Mathematics, University of Zimbabwe, P.O. Box MP 167, Harare, Zimbabwe

Received 17 March 2015; Accepted 21 June 2015

Academic Editor: Mauro Sodini

Copyright © 2015 Steady Mushayabasa. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

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 this model fully hinges on whether 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.

1. Introduction

Anthrax is an acute zoonotic disease caused by the spore-forming 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.

2. 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 anthrax-induced 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; 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 = 0.000014 day⁻¹.

For biological reasons we will study the dynamics of system (1) in the closed set

where denotes the nonnegative cone of including its lower dimensional faces and . It can easily be verified that is positively invariant with respect to (1).

2.1. Equilibria and the Reproductive Number

In the absence of anthrax infection in the community, system (1) admits an infection-free equilibrium given by . The Jacobian matrix of system (1) evaluated about is

It is clear that is an eigenvalue of matrix . The other two are determined from matrix . Consider

The trace and determinant of are, respectively, given by

Let the reproductive number of system (1) be

Since and (when ) it follows that the equilibrium point is locally asymptotically stable if and only if .

Biologically, the term 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 and the endemic equilibrium point , given by

From (7) it is evident that makes biological sense whenever .

2.2. 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 () the average number of new anthrax infections generated through spores ingestion is modeled by the threshold quantity 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 free-living 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.

Figure 1: Graphical illustration showing the strength of environmental decontamination on controlling the spread of anthrax. We consider day⁻¹ [5] and day⁻¹ [4].

Numerical results in Figure 1 suggest that environmental decontamination rate 0.05 day⁻¹ 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⁻¹ can be effective in attaining 100% control on the spread of anthrax in the community.

2.3. 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:

System (10) satisfies the following initial conditions:

where for denote the nonnegative continuous functions on . 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 , 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 as in system (1). The Jacobian matrix of system (10) about is

From (13) it follows that the characteristic equation is

Since is a root of (14), we only need to consider

where

Equation (14) with is

If , we have already deduced that all roots of system (15) have negative real parts when , and only one root of (15) has positive real part when .

For , we assume that () is a root of (15). This is the case if and only if

Separating the real and imaginary parts yields

Adding up the squares of both equations, we obtain

Let , so that (21) reduces to

Substituting for , , and gives

Substituting into (22) gives

Whenever , (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 of Ruan and Wei [6] that the real parts of all eigenvalues of characteristic equation (14) are negative for all values of the delay . Thus, we have the following theorem.

Theorem 1. The equilibrium point of system (10) is locally asymptotically stable for all time delay if .

In the next theorem we establish the global stability of the infection-free equilibrium for system (10).

Theorem 2. The equilibrium point of system (10) is globally asymptotically stable when for all .

Proof. We denote by the translation of the solution of system (10); that is, , where . Let us define a Lyapunov functional

The derivative of along the solutions of (10) is given by

Therefore, holds for all and . Furthermore, if and only if , or . Therefore, the largest compact invariant set in , when , is the singleton . Thus, LaSalle’s invariance principle [7] implies that is globally asymptotically stable in . This proves the theorem.

2.3.2. Endemic Equilibrium Point and Its Stability

Theorem 3. Suppose ; then the equilibrium point is unstable.

Proof. From the discussion in Section 2.3.1, we have deduced that the characteristic equation associated with is given by (14), and we only need to consider (15). From the above computations it can easily be verified that if ,  (15) has a positive root when . Now, let () be a root of (15). Then solving (23) gives

For 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 .
Further, let be the root of (14) satisfying and .
Differentiating both sides of (14) with respect to we get

This gives

Thus

By applying Lemma in Ruan and Wei [6] and observing that (15) has a positive real root when , we obtain that characteristic equation (14) has a positive root at least for all . Therefore, the equilibrium point of system (10) is unstable when . This completes the proof.

Now, we investigate the effect of the time delay on the local stability of . The characteristic equation of system (10) at the endemic equilibrium takes the following form:

where

When , (32) becomes

By the Hurwitz criteria, all the roots of (34) have only negative real parts. Thus is locally asymptotically stable when . Now, we consider the case . If () is a solution of (32), we have

Separating the real and imaginary parts we have

Squaring and adding the two equations of (36) yields

where

Hence, if , (32) has no positive roots. Accordingly, if , the endemic equilibrium exists and is locally asymptotically stable for all . Further, using a Lyapunov functional, we can obtain the following theorem.

Theorem 4. Whenever , then the unique equilibrium of system (10) is globally asymptotically stable in for all .

Proof. We consider the following Lyapunov functional:

Differentiating along the solution of system (10) and using the identities

one gets

Here,

for all , because the arithmetic mean is greater than or equal to the geometric mean. Further, note that for any , , and ; and the equality is satisfied if and only if , and, consequently, . Moreover, the largest invariant set of is a singleton where , , and . By the Lyapunov-LaSalle invariance principle [7], we obtain global asymptotic stability of the endemic equilibrium under the condition .

3. 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.

Conflict of Interests

The author declares no conflict of interests.

Acknowledgments

The author is grateful to the anonymous referee and handling editor for their valuable comments and suggestions.

References

CDC, “Use of anthrax vaccine in the United States: recommendations of the Advisory Committee on Immunization Practices (ACIP), 2009,” Morbidity and Mortality Weekly Report, vol. 59, no. 6, pp. 1–30, 2010, http://www.cdc.gov/mmwr/

A. N. Survely, B. Kvasnicka, and R. Torell, “Anthrax: a guide for livestock producers,” Cattle Producer's Library CL613, Western Beef Resource Committe, 2006.

A. Friedman and A. A. Yakubu, “Anthrax epizootic and migration: persistence or extinction,” Mathematical Biosciences, vol. 241, no. 1, pp. 137–144, 2013.

B. D. Hahn and P. R. Furniss, “A deterministic model of an anthrax epizootic: threshold results,” Ecological Modelling, vol. 20, no. 2-3, pp. 233–241, 1983.

S. S. Lewerin, M. Elvander, T. Westermark et al., “Anthrax outbreak in a Swedish beef cattle herd—1st case in 27 years: case report,” Acta Veterinaria Scandinavica, vol. 52, no. 1, article 7, 2010.

S. Ruan and J. Wei, “On the zeros of transcendental functions with applications to stability of delay differential equations with twon delays,” Dynamics of Continuous, Discrete and Impulsive Systems, vol. 10, pp. 863–874, 2003.

J. S. LaSalle, The Stability of Dynamical Systems, vol. 25 of CBMS-NSF Regional Conference Series in Applied Mathematics, SIAM, Philadelphia, Pa, USA, 1976.