# Mathematical Analysis of the Effects of Controls on the Transmission Dynamics of Anthrax in Both Animal and Human Populations

Authors: Elijah B. Baloba; Baba Seidu; Christopher Saaha Bornaa
Journal: Computational and Mathematical Methods in Medicine
 (2020)

---

## Abstract

A nonlinear differential equation model is proposed to study the impact of vaccination on the transmission dynamics of anthrax in both livestock and human populations. The model is shown to exhibit only two equilibria, namely, the disease-free and the endemic equilibrium points, which are proven to be locally stable if the basic reproduction number (ℛ0) is less than unity and greater than unity, respectively. Local sensitivity analysis shows that the infection rate, pathogen-shedding rate, and rate of vaccination of livestock are parameters with the most positive impact on disease spread, whereas the rate of disinfection followed by the rate of vaccination are the parameters with the most negative impact on disease spread. Numerical simulation shows that implementing all control measures (i.e., vaccination, education, disinfection, and treatment) is a most effective strategy to curb disease spread.


---

## Body

## 1. Introduction

Poverty is one of the major challenges facing the developing world, and household income is seen as one of the factors worsening poverty situations for many developing countries, where hunger and extreme poverty are widespread [1]. Livestock production has therefore become a key option for policy makers and stakeholders for increasing income and improving livelihood of rural households [2, 3]. The livestock sector does not only provide manure for crop and vegetable farmers [4] but also provides food and nutritional security, financial security, and employment [3] to improve the standard of living, especially in rural households. Over 100,000 households in Ghana depend on livestock production for their livelihood, especially in the five regions of northern Ghana [5]. Globally, about 70% of the world’s 1.4 billion extremely poor people depend on livestock [6]. One of the setbacks of livestock rearing is the anthrax disease which mainly affects livestock and other herbivores [7]. Anthrax as an infection caused by the bacterium, Bacillus anthracis, is endemic throughout Africa and other parts of the world [8]. Anthrax results in considerable livestock and wildlife losses and also causes severe and fatal infections in the human population [8]. There has not been any reported case of human-to-human anthrax transmission [9], but it is one of the diseases causing poverty and negatively affecting both livelihoods and human health with serious implications [10].Anthrax occurs every year, leading to deaths of many animals ranging from hundreds to thousands and transmission to humans. Humans often contract the disease from an exposure to infected animals or their products [8, 11]. It is estimated that about 2,000 to 20,000 human cases of anthrax occur each year worldwide [11].Ghana has 28 scheduled livestock diseases, of which anthrax is one of the major occurrences in the country and causes severe livestock mortality and deaths in human population. Anthrax occurrence rose from 29 cases in the year 2012 to 1875 cases in 2013 and 943 cases in 2014 [3]. Aside causing human and animal mortality, it has devastating effects on the economy and livelihood systems in Ghana and sub-Sahara Africa since it affects the quality and quantity of meat, milk, and manure output thereby causing food insecurity, especially in poor rural households [3].The exploration and control of infectious diseases through mathematical modeling has proved to be an important tool for researchers as it helps in increasing our understanding of infectious diseases by revealing some of the complex dynamics of the diseases being modeled [12]. For anthrax disease, some mathematical models (e.g., [13–17]) have been proposed in the literature to help improve the understanding of how the disease spreads and can be eradicated.In a study by Hahn and Furniss [13], a deterministic model of an anthrax epizootic was studied with the hypothesis that anthrax is transmitted through environmental contamination. An analytic form of the generation multiplication rate was established as well as an analytical expression for the threshold of the disease.Friedman and Yakubu [14] extended the model of Hahn and Furniss [13] by studying the effects of anthrax transmission, carcass ingestion, carcass-induced environmental contamination and migration rates on the persistence or extinction of animal population. They showed that decreasing the levels of carcass ingestion may not always lead to a reduction in the number of infective animals but increasing the rate of carcass-induced environmental contamination can result in extinction of a persistent animal population.Mushayabasa [15] also extended the model of Hahn and Furniss [13] by including a fixed time delay and environmental decontamination on eradication of anthrax disease. They explored the effectiveness of environmental decontamination on eradication of anthrax disease through the reproduction number.Sinkie and Murthy [16] extended the model of Mushayabasa [15] by assuming that some infective animals can display clinical signs of the disease and can be cured by treatment. Their study showed that the pathogen level decreased and susceptible animals increased when the curing rate is increased through treatment.Osman et al. [18] extended the model of Friedman and Yakubu [14] by constructing a mathematical model of anthrax in human and animal population to investigate the impact of vaccination through sensitivity analysis.Whilst so much research has been done on zoonotic diseases and anthrax in particular, it is worth noting that not much has been done in the use of Mathematics to model the impact of vaccination and education on the transmission dynamics of anthrax disease in both human and animal populations. In this paper, the impact of anthrax on food security, income, and the health of human beings have necessitated investigation of the transmission and control of anthrax disease in livestock and human populations through mathematical modeling by extending the study of Sinkie and Murthy [16]. This was done by incorporating vaccination of susceptible livestock, cross contamination of the human population, and including education of the human population as a control measure.


## 2. Model Formulation

In this section, a deterministic model is proposed for the study of anthrax transmission in both livestock and human populations in the presence of vaccination programme and education. The following assumptions are considered:(i)
Anthrax is spread by environmental contamination(ii)
There is no vertical transmission(iii)
Susceptible livestock get infected only through ingestion of the anthrax pathogens(iv)
Some infected livestock can show symptoms of the disease and can be cured through treatment(v)
Susceptibles, infectives, educated population, and vaccinated livestock can exit naturally(vi)
Newly born and other susceptible livestock are vaccinated(vii)
Susceptible humans get infected by coming into contact with infected livestock, infected carcass, and pathogens(viii)
There is no human-to-human infection.The compartmental diagram in Figure1 shows the transmission of anthrax in both humans and livestock population. The livestock population is subdivided into three mutually exclusive compartments of susceptible livestock SLt, vaccinated VLt, and infected individuals ILt, with adjoining compartments of infected carcass Ct and pathogens PLt shed into the environment. The livestock that are at risk of getting the infection are called the susceptible livestock. Vaccinated livestock are those that received vaccination before the outbreak of the disease. Infected livestock include those that show symptoms of the anthrax disease.Figure 1
Schematic diagram of anthrax epidemic in both livestock and human populations.The spillover compartmental model for the human population consists of the susceptible human compartmentSHt, who are at risk of contracting the anthrax disease upon contact with the infected animal, infected carcass, or pathogen in the environment. A portion Ed of the human population is a compartment of careful susceptible human population who cannot contract the disease due to education or experience. The infected human compartment IHt includes individuals who show symptoms of the disease.Table1 summarizes the definition of parameters and variables used in the model.Table 1
Parameters used in the model.
Par.
Description
Value
References
β
1
Infection rate
0.02
K
The rate of inflow of animals
0.99
[16]
γ
The recovery rate of infected animals
0.0025
[18]
μ
Natural death rate
0.0001
[16]
π
Anthrax-induced death rate
0.5
[16]
ϕ
1
The rate of disinfection of infected carcass
0.6
Assumed
ϕ
2
The rate of disinfection of pathogens in the environment
0.5
[16]
d
Decay rate of infective carcasses
0.8
[16]
α
The rate of shedding pathogens into the ground
0.45
[18]
A
Rate of inflow of livestock
0.75
Assumed
q
Fraction of recruited livestock that are vaccinated
0.7
Assumed
Q
Rate of loss of immunity
0.002
Assumed
a
Recruitment rate of humans
0.92
e
0
Natural death rate in humans
0.0001
β
2
Infection rate of susceptible humans
0.0001
b
The recovering rate of infected humans
0.04
ω
Rate at which humans acquire knowledge to avoid infection
0.50
ρ
Rate of vaccination of susceptibles
0.75
Assumed
S
L
Population size of susceptible livestock
State variable
V
L
Population size of vaccinated livestock
State variable
I
L
Population size of infected livestock
State variable
C
Population size of carcass of livestock
State variable
P
L
Concentration of anthrax pathogens in the environment
State variable
S
H
Population size of susceptible humans
State variable
E
Population size of humans educated enough to avoid infection
State variable
I
H
Population size of infected humans
State variableThe susceptible livestock in compartmentSL gets infected when they ingest or come into contact with pathogens in compartment PL at rate β1. Livestock that show symptoms of the disease are treated to recover without immunity at rate γ to become susceptible.A fraction,q, of newborns (recruits) A are vaccinated and thus are placed in compartment VL. That is, qA livestock enters compartment VL while the rest of newborns 1−qA enter compartment SL. Susceptible livestock are vaccinated at rate ρ to join compartment, VL. It is assumed that the vaccinated individuals have temporary immunity so that Q is the rate at which vaccinated individuals lose immunity and return to compartment SL.Infected carcasses of livestock are placed in compartmentC and are disinfected at a rate, ϕ1. The remaining pathogens from the carcasses get into the environment at a rate α to form compartment PL. The environment is disinfected of pathogens at a rate of ϕ2.Humans can get the disease by consuming infected animal products or coming into contact with anthrax infected carcasses. In addition, it is assumed that human beings can acquire anthrax directly by getting into contact with free living spores or pathogens in the environment. As a means of controlling the disease, education is assumed to lead susceptible humans to gain knowledge at rateω to avoid being infected. The model under discussion is thus represented by the following system of equations:(1)dSLdt=k+1−qA+QVLt+γILt−β1PLt+μ+ρSLt,dVLdt=ρSLt+qA−Q+μVLt,dILdt=β1SLtPLt−γ+π+μILt,dCdt=π+μILt−ϕ1+dCt,dPLdt=αCt−ϕ2PLt,dSHdt=a+bIHt−β2ILt+Ct+PLt+ω+e0SHt,dEddt=ωSHt−e0Edt,dIHdt=β2ILt+Ct+PLtSHt−b+e0IHtwith initial conditionsSL0≥0,VL0≥0,IL0≥0,C0≥0,PL0≥0,SH0≥0,Ed0≥0,IH0≥0.,For convenience, we use the following substitutions:(2)ψ1=k+1−qA,ψ2=μ+ρ,ψ3=γ+π+μ,ψ4=ϕ1+d,ψ5=Q+μ,ψ6=ω+e0,ψ7=π+μ,Γ1=b+e0,δ=β2IL+C+PL.All model parameters are assumed to be nonnegative, and hence, we are only concerned with nonnegative solutions of equation (1).In the next section, the qualitative analysis of the model is carried out.


## 3. Qualitative Analysis

### 3.1. Positivity of Solutions

Theorem 1.
If the initial conditions of the state equations are nonnegative (i.e.,SL0≥0, VL0≥0, IL0≥0, C0≥0, PL0≥0, SH0≥0, Ed0≥0, IH0≥0), then future solutions are also nonnegative (i.e., SLt≥0, VLt≥0, ILt≥0, Ct≥0, PLt≥0, SHt≥0, Edt≥0, and IHt≥0).Proof.
Lett>0 be the first time such that SLt=0. Then, from the first equation, we have dSL/dt=k+1−qA+QVL+γIL≥0. Hence, SLt≥∀t>0.
Similarly, consider a positive timet. If VLt=0, then the second equation becomes dVL/dt=ρSL+qA≥0, which implies that VLt≥0,∀t>0. All the remaining state variables can similarly be shown to be nonnegative for all positive times. This concludes the proof.


### 3.2. Boundedness and Feasible Region of Solutions of Model

The feasible region of the model is the set containing reasonable solutions of the model. In this section, we show that there exists a region within which the model can sufficiently be studied. In this region, the model is epidemiologically and mathematically well posed. This is provided for in the following theorem.Theorem 2.
All solutions of model (1) are uniformly bounded and contained in the feasible region defined by ξ=ξL×ξH,∀t>0, where ξL=SLt,VLt,ILt∈ℝ+3limsupt⟶∞NLt⟶∞≤k+A/μ and ξH=SHt,Edt,IHt∈ℝ+3limsupt⟶∞NHt≤a/e0. That is the region, ξ is positively invariant under the flow defined by model (1).Proof.
Let the total livestock and human populations be denoted byNL and NH.
ThenNL=SL+VL+IL and NH=SH+Ed+IH, respectively. Then we have dNL/dt=k+A−μNL−πIL and dNH/dt=a−e0NH.
FromdNL/dt, we have dNL/dt≤k+A−μNL.
By standard comparison theorems [19], we have NLt≤NL0e−μt+k+A/μ1−e−μt, so that limsupt⟶∞NLt≤k+A/μ.
Also, fromdNH/dt, we have NHt=NH0e−e0t+a/e01−e−e0t.
Hence all solutions are uniformly bounded.
Now, define(3)ξL=SLt,VLt,ILt∈ℝ+3limsupt⟶∞NLt≤k+Aμ,ξH=SHt,Edt,IHt∈ℝ+3limsupt⟶∞NHt≤ae0.
Then, the regionξ=ξL×ξH is positively invariant under model (1). Therefore, the model is mathematically and epidemiologically well posed, and hence, it is sufficient to study the model in ξ. This concludes the proof.


### 3.3. Equilibrium Points

Model (1) has basically two equilibrium points, namely, the disease-free and endemic equilibrium points. The disease-free equilibrium is given by(4)ℰ0=ψ1ψ2+Qρψ1+qAψ2ψ2ψ2ψ5−ρQ,ρψ1+qAψ2ψ2ψ5−ρQ,0,0,0,aψ6,ωaψ6,0.Next, we find the basic reproduction number using the method described in [20]. The basic reproduction number is defined as the average number of secondary anthrax infections a single infected animal causes during its infectious life if it is introduced into an initially anthrax-free human livestock system. The basic reproduction number ℛ0 is the dominant eigenvalue (spectral radius) of the next generation matrix FV−1, where F and V are the transmission and transition matrices [20]. For our model, the matrices are given by(5)V=ψ3000−ψ7ψ4000−αϕ20000Γ1,F=00β1SLt000000000β2SHtβ2SHtβ2SHt0.The basic reproduction number is thus given by(6)ℛ0=β1αψ7SL0ψ3ψ4ϕ2.In the presence of pathogens in the environment, the model exhibits the endemic equilibrium point,ℰ∗=SL∗,VL∗,IL∗,C∗,PL∗,SH∗,Ed∗,IH∗, where(7)SL∗=SL0ℛ0,VL∗=1ψ5ρSL0ℛ0+qA,IL∗=μQ+ψ2SL0ψ52ψ71−1ℛ0,C∗=μQ+ψ2SL0ψ4ψ521−1ℛ0,PL∗=αμQ+ψ2SL0ϕ2ψ4ψ521−1ℛ0SH∗=aΓ1Γ1ψ6+e0β2IL∗+C∗+PL∗,Ed∗=ωe0SH∗,IH∗=β2IL∗+C∗+PL∗SH∗Γ1,IL∗+C∗+PL∗=μQ+ψ2SL0ψ52ϕ2ψ4ψ7ϕ2ψ4+α+ϕ2ψ71−1ℛ0,,Clearly, the endemic equilibrium exists only whenℛ0>0. In the next section, we study the local stability of the equilibrium points of the model.


### 3.4. Stability Analysis of Equilibrium Points

In this section, we study the local stability of the disease-free and endemic equilibrium points using the indirect Lyapunov method. By this method, an equilibrium pointx∗ of a model dx/dt=fx is said to be locally asymptotically stable if all eigenvalues of the Jacobian of the model evaluated at x∗ have negative real parts.The Jacobian of model (1) is given by(8)J=−β1PL−ψ2Qγ0−β1SL000ρ−ψ5000000β1PL0−ψ30β1SL00000ψ7−ψ40000000α−ϕ200000−β2SH−β2SH−β2SH−β2δ−ψ60b00000ω−e0000β2SHβ2SHβ2SHβ2δ0−Γ1.The Jacobian,J, evaluated at ℰ0 has eigenvalues −e0, −ψ6, and −Γ1, and the rest satisfy the following characteristic polynomial equation.(9)λ5+k4λ4+k3λ3+k2λ2+k1λ+k0=0,where(10)k0=μϕ2ψ4ψ3Q+ψ21−ℛ0,k1=ϕ2ψ3ψ4ψ5+ψ21−ℛ0+μψ4ψ3+ϕ2ψ4+ϕ2ψ3Q+ψ2,k2=ϕ2ψ4ψ31−ℛ0+ψ5+ψ2ψ4ψ3+ϕ2ψ4+ψ3+μϕ2+ψ3+ψ4Q+ψ2,k3=ψ2+ψ4+ψ3+ψ5ϕ2+ψ4ψ3+ψ5+ψ2+ψ3ψ5+ψ2+μQ+ψ2,k4=ϕ2+ψ2+ψ3+ψ4+ψ5.Clearly,k3>0,k4>0 and also k0≥0,k1>0,k2>0 whenever ℛ0≤1.By the Descartes rule of signs, it can be shown that equation (9) has all its roots with real parts lying in the left half of the plane whenever ℛ0≤1. The following result is thus established.Theorem 3.
The disease-free equilibrium pointℰ0 is locally asymptotically stable whenever ℛ0≤1 and unstable otherwise.
Also, the Jacobian evaluated at the endemic equilibrium−e0 as an eigenvalue, and the remaining eigenvalues can be shown to satisfy the following characteristic polynomial:(11)λ7+η6λ6+η5λ5+η4λ4+η3λ3+η2λ2+η1λ+η0=0,where(12)η6=Γ1+β1PL∗+β2δ+ϕ2+ψ2+ψ3+ψ4+ψ5+ψ6,η5=β2δ+Γ1+ψ7+ψ5+ψ6+ϕ2+ψ4β1PL∗+ψ5+ψ2+e0+ϕ2+ψ4+ψ3δβ2+ϕ2+ψ6+ψ2+ψ4+ψ3+ψ5Γ1+ψ3+ψ6+ψ2+ψ4+ψ5ϕ2+ψ3ψ5+ψ3ψ2+ψ6ψ4+ψ6ψ3+ψ6ψ5+ψ6ψ2+ψ4ψ2+μQ+ψ2+ψ4ψ5+ψ4ψ3,η4=β1PL∗ψ4+ψ5+ψ7+ϕ2β2δ+Γ1+ψ6+μψ2+Qβ2δ+ψ6+ϕ2+Γ1+ψ3+ψ4+aΓ1SH∗ψ4+ψ5+ϕ2+β1PL∗+ψ2+ψ3+ψ5+ψ2ϕ2+ψ3+ψ4β2δ+Γ1+ψ6+PL∗β1ϕ2ψ5+ψ4ψ5+ϕ2ψ4+ψ4ψ7+ψ5ψ7+ϕ2ψ7+ψ4ψ3+ϕ2ψ4+ϕ2ψ3Γ1+β2δ+ψ5+ψ6+ψ2,η3=aΓ1SH∗ψ4+ψ5+ψ7+ϕ2β1PL∗+ϕ2+ψ3+μ+ψ4ψ2+ψ4+ψ5+ϕ2ψ3+ψ4ϕ2+ψ5+μQ+ϕ2ψ5+ψ4+ψ5+ϕ2ψ7+ϕ2ψ4+ψ5+ψ4ψ5β1PL∗+ψ5+ψ2ψ4ψ3+ϕ2ψ4+ϕ2ψ3+μϕ2+ψ3+ψ4ψ2+Qβ2δ+Γ1+ψ6+ϕ2ψ5+ϕ2ψ4+ψ4ψ5ψ7+ϕ2ψ4ψ5β1PL∗+μψ4ψ3+ϕ2ψ4+ϕ2ψ3ψ2+Q,η2=aΓ1SH∗μϕ2+ψ3+ψ4ψ2+Q+ψ5+ψ2ψ4ψ3+ϕ2ψ4+ϕ2ψ3+ψ4+ψ5+ϕ2ψ7+ψ4+ϕ2ψ5+ϕ2ψ4β1PL∗+β1ϕ2ψ4PL∗ψ5ψ7+β2δ+Γ1+ψ6ϕ2ψ5+ϕ2ψ4+ψ4ψ5ψ7+ϕ2ψ4ψ5β1PL∗+μψ4ψ3+ϕ2ψ4+ϕ2ψ3ψ2+Q,η1=β1PL∗β2δ+ψ6+Γ1ϕ2ψ4ψ5ψ7+aμΓ1SH∗ψ4ψ3+ϕ2ψ4+ϕ2ψ3ψ2+Q+aΓ1SH∗ψ4ψ7+ψ4ψ5+ψ5ψ7ϕ2+ψ5ψ7ψ4β1PL∗,η0=aΓ1ϕ2ψ4SH∗β1ψ5ψ7PL∗+2ψ3μQ+ψ2.
Recall that the endemic equilibrium exists only whenℛ0>1 and when that happens, we have ηi>0,∀i=0,1,2,…,6. By Descartes’s rule of signs, equation (9) has all its roots in the left half of the real plane. This proves that the eigenvalues of the Jacobian of the model evaluated at the endemic equilibrium point (whenever it exists) all have negative real parts. The following results are thus established.Theorem 4.
The endemic equilibrium pointℰ∗ is locally asymptotically stable whenever ℛ0>1 and unstable otherwise.


### 3.5. Local Sensitivity Analysis

In this section, a local sensitivity analysis is conducted on the basic reproduction number using the Normalized Forward Sensitivity Index to determine the impact of parameter values on the basic reproduction numberℛ0 [21, 22]. The sensitivity indexes in Table 2 were evaluated using the parameter values in Table 1.Table 2
Sensitivity indexes.
Parameter
Sensitivity index
ϕ
2
−1.0000
α
+1.0000
β
1
+1.0000
π
+0.9998
ρ
−0.9972
μ
−0.9662
Q
+0.9636
k
+0.5773
d
−0.5333
ϕ
1
−0.4667
A
+0.4227
q
−0.0146
γ
−7.95E − 04


## 3.1. Positivity of Solutions

Theorem 1.
If the initial conditions of the state equations are nonnegative (i.e.,SL0≥0, VL0≥0, IL0≥0, C0≥0, PL0≥0, SH0≥0, Ed0≥0, IH0≥0), then future solutions are also nonnegative (i.e., SLt≥0, VLt≥0, ILt≥0, Ct≥0, PLt≥0, SHt≥0, Edt≥0, and IHt≥0).Proof.
Lett>0 be the first time such that SLt=0. Then, from the first equation, we have dSL/dt=k+1−qA+QVL+γIL≥0. Hence, SLt≥∀t>0.
Similarly, consider a positive timet. If VLt=0, then the second equation becomes dVL/dt=ρSL+qA≥0, which implies that VLt≥0,∀t>0. All the remaining state variables can similarly be shown to be nonnegative for all positive times. This concludes the proof.


## 3.2. Boundedness and Feasible Region of Solutions of Model

The feasible region of the model is the set containing reasonable solutions of the model. In this section, we show that there exists a region within which the model can sufficiently be studied. In this region, the model is epidemiologically and mathematically well posed. This is provided for in the following theorem.Theorem 2.
All solutions of model (1) are uniformly bounded and contained in the feasible region defined by ξ=ξL×ξH,∀t>0, where ξL=SLt,VLt,ILt∈ℝ+3limsupt⟶∞NLt⟶∞≤k+A/μ and ξH=SHt,Edt,IHt∈ℝ+3limsupt⟶∞NHt≤a/e0. That is the region, ξ is positively invariant under the flow defined by model (1).Proof.
Let the total livestock and human populations be denoted byNL and NH.
ThenNL=SL+VL+IL and NH=SH+Ed+IH, respectively. Then we have dNL/dt=k+A−μNL−πIL and dNH/dt=a−e0NH.
FromdNL/dt, we have dNL/dt≤k+A−μNL.
By standard comparison theorems [19], we have NLt≤NL0e−μt+k+A/μ1−e−μt, so that limsupt⟶∞NLt≤k+A/μ.
Also, fromdNH/dt, we have NHt=NH0e−e0t+a/e01−e−e0t.
Hence all solutions are uniformly bounded.
Now, define(3)ξL=SLt,VLt,ILt∈ℝ+3limsupt⟶∞NLt≤k+Aμ,ξH=SHt,Edt,IHt∈ℝ+3limsupt⟶∞NHt≤ae0.
Then, the regionξ=ξL×ξH is positively invariant under model (1). Therefore, the model is mathematically and epidemiologically well posed, and hence, it is sufficient to study the model in ξ. This concludes the proof.


## 3.3. Equilibrium Points

Model (1) has basically two equilibrium points, namely, the disease-free and endemic equilibrium points. The disease-free equilibrium is given by(4)ℰ0=ψ1ψ2+Qρψ1+qAψ2ψ2ψ2ψ5−ρQ,ρψ1+qAψ2ψ2ψ5−ρQ,0,0,0,aψ6,ωaψ6,0.Next, we find the basic reproduction number using the method described in [20]. The basic reproduction number is defined as the average number of secondary anthrax infections a single infected animal causes during its infectious life if it is introduced into an initially anthrax-free human livestock system. The basic reproduction number ℛ0 is the dominant eigenvalue (spectral radius) of the next generation matrix FV−1, where F and V are the transmission and transition matrices [20]. For our model, the matrices are given by(5)V=ψ3000−ψ7ψ4000−αϕ20000Γ1,F=00β1SLt000000000β2SHtβ2SHtβ2SHt0.The basic reproduction number is thus given by(6)ℛ0=β1αψ7SL0ψ3ψ4ϕ2.In the presence of pathogens in the environment, the model exhibits the endemic equilibrium point,ℰ∗=SL∗,VL∗,IL∗,C∗,PL∗,SH∗,Ed∗,IH∗, where(7)SL∗=SL0ℛ0,VL∗=1ψ5ρSL0ℛ0+qA,IL∗=μQ+ψ2SL0ψ52ψ71−1ℛ0,C∗=μQ+ψ2SL0ψ4ψ521−1ℛ0,PL∗=αμQ+ψ2SL0ϕ2ψ4ψ521−1ℛ0SH∗=aΓ1Γ1ψ6+e0β2IL∗+C∗+PL∗,Ed∗=ωe0SH∗,IH∗=β2IL∗+C∗+PL∗SH∗Γ1,IL∗+C∗+PL∗=μQ+ψ2SL0ψ52ϕ2ψ4ψ7ϕ2ψ4+α+ϕ2ψ71−1ℛ0,,Clearly, the endemic equilibrium exists only whenℛ0>0. In the next section, we study the local stability of the equilibrium points of the model.


## 3.4. Stability Analysis of Equilibrium Points

In this section, we study the local stability of the disease-free and endemic equilibrium points using the indirect Lyapunov method. By this method, an equilibrium pointx∗ of a model dx/dt=fx is said to be locally asymptotically stable if all eigenvalues of the Jacobian of the model evaluated at x∗ have negative real parts.The Jacobian of model (1) is given by(8)J=−β1PL−ψ2Qγ0−β1SL000ρ−ψ5000000β1PL0−ψ30β1SL00000ψ7−ψ40000000α−ϕ200000−β2SH−β2SH−β2SH−β2δ−ψ60b00000ω−e0000β2SHβ2SHβ2SHβ2δ0−Γ1.The Jacobian,J, evaluated at ℰ0 has eigenvalues −e0, −ψ6, and −Γ1, and the rest satisfy the following characteristic polynomial equation.(9)λ5+k4λ4+k3λ3+k2λ2+k1λ+k0=0,where(10)k0=μϕ2ψ4ψ3Q+ψ21−ℛ0,k1=ϕ2ψ3ψ4ψ5+ψ21−ℛ0+μψ4ψ3+ϕ2ψ4+ϕ2ψ3Q+ψ2,k2=ϕ2ψ4ψ31−ℛ0+ψ5+ψ2ψ4ψ3+ϕ2ψ4+ψ3+μϕ2+ψ3+ψ4Q+ψ2,k3=ψ2+ψ4+ψ3+ψ5ϕ2+ψ4ψ3+ψ5+ψ2+ψ3ψ5+ψ2+μQ+ψ2,k4=ϕ2+ψ2+ψ3+ψ4+ψ5.Clearly,k3>0,k4>0 and also k0≥0,k1>0,k2>0 whenever ℛ0≤1.By the Descartes rule of signs, it can be shown that equation (9) has all its roots with real parts lying in the left half of the plane whenever ℛ0≤1. The following result is thus established.Theorem 3.
The disease-free equilibrium pointℰ0 is locally asymptotically stable whenever ℛ0≤1 and unstable otherwise.
Also, the Jacobian evaluated at the endemic equilibrium−e0 as an eigenvalue, and the remaining eigenvalues can be shown to satisfy the following characteristic polynomial:(11)λ7+η6λ6+η5λ5+η4λ4+η3λ3+η2λ2+η1λ+η0=0,where(12)η6=Γ1+β1PL∗+β2δ+ϕ2+ψ2+ψ3+ψ4+ψ5+ψ6,η5=β2δ+Γ1+ψ7+ψ5+ψ6+ϕ2+ψ4β1PL∗+ψ5+ψ2+e0+ϕ2+ψ4+ψ3δβ2+ϕ2+ψ6+ψ2+ψ4+ψ3+ψ5Γ1+ψ3+ψ6+ψ2+ψ4+ψ5ϕ2+ψ3ψ5+ψ3ψ2+ψ6ψ4+ψ6ψ3+ψ6ψ5+ψ6ψ2+ψ4ψ2+μQ+ψ2+ψ4ψ5+ψ4ψ3,η4=β1PL∗ψ4+ψ5+ψ7+ϕ2β2δ+Γ1+ψ6+μψ2+Qβ2δ+ψ6+ϕ2+Γ1+ψ3+ψ4+aΓ1SH∗ψ4+ψ5+ϕ2+β1PL∗+ψ2+ψ3+ψ5+ψ2ϕ2+ψ3+ψ4β2δ+Γ1+ψ6+PL∗β1ϕ2ψ5+ψ4ψ5+ϕ2ψ4+ψ4ψ7+ψ5ψ7+ϕ2ψ7+ψ4ψ3+ϕ2ψ4+ϕ2ψ3Γ1+β2δ+ψ5+ψ6+ψ2,η3=aΓ1SH∗ψ4+ψ5+ψ7+ϕ2β1PL∗+ϕ2+ψ3+μ+ψ4ψ2+ψ4+ψ5+ϕ2ψ3+ψ4ϕ2+ψ5+μQ+ϕ2ψ5+ψ4+ψ5+ϕ2ψ7+ϕ2ψ4+ψ5+ψ4ψ5β1PL∗+ψ5+ψ2ψ4ψ3+ϕ2ψ4+ϕ2ψ3+μϕ2+ψ3+ψ4ψ2+Qβ2δ+Γ1+ψ6+ϕ2ψ5+ϕ2ψ4+ψ4ψ5ψ7+ϕ2ψ4ψ5β1PL∗+μψ4ψ3+ϕ2ψ4+ϕ2ψ3ψ2+Q,η2=aΓ1SH∗μϕ2+ψ3+ψ4ψ2+Q+ψ5+ψ2ψ4ψ3+ϕ2ψ4+ϕ2ψ3+ψ4+ψ5+ϕ2ψ7+ψ4+ϕ2ψ5+ϕ2ψ4β1PL∗+β1ϕ2ψ4PL∗ψ5ψ7+β2δ+Γ1+ψ6ϕ2ψ5+ϕ2ψ4+ψ4ψ5ψ7+ϕ2ψ4ψ5β1PL∗+μψ4ψ3+ϕ2ψ4+ϕ2ψ3ψ2+Q,η1=β1PL∗β2δ+ψ6+Γ1ϕ2ψ4ψ5ψ7+aμΓ1SH∗ψ4ψ3+ϕ2ψ4+ϕ2ψ3ψ2+Q+aΓ1SH∗ψ4ψ7+ψ4ψ5+ψ5ψ7ϕ2+ψ5ψ7ψ4β1PL∗,η0=aΓ1ϕ2ψ4SH∗β1ψ5ψ7PL∗+2ψ3μQ+ψ2.
Recall that the endemic equilibrium exists only whenℛ0>1 and when that happens, we have ηi>0,∀i=0,1,2,…,6. By Descartes’s rule of signs, equation (9) has all its roots in the left half of the real plane. This proves that the eigenvalues of the Jacobian of the model evaluated at the endemic equilibrium point (whenever it exists) all have negative real parts. The following results are thus established.Theorem 4.
The endemic equilibrium pointℰ∗ is locally asymptotically stable whenever ℛ0>1 and unstable otherwise.


## 3.5. Local Sensitivity Analysis

In this section, a local sensitivity analysis is conducted on the basic reproduction number using the Normalized Forward Sensitivity Index to determine the impact of parameter values on the basic reproduction numberℛ0 [21, 22]. The sensitivity indexes in Table 2 were evaluated using the parameter values in Table 1.Table 2
Sensitivity indexes.
Parameter
Sensitivity index
ϕ
2
−1.0000
α
+1.0000
β
1
+1.0000
π
+0.9998
ρ
−0.9972
μ
−0.9662
Q
+0.9636
k
+0.5773
d
−0.5333
ϕ
1
−0.4667
A
+0.4227
q
−0.0146
γ
−7.95E − 04


## 4. Numerical Simulation

We use the fourth order Runge-Kutta scheme to perform a numerical simulation of the model using parameter values in Table1. We perform the simulation for different combinations of the controls (vaccination, education, disinfection, and treatment) (see Figures 2–7). We also study the impact of various parameters on model dynamics by running the model for different values of the parameters (see Figures 8–11).Figure 2
Trajectories of model solution without implementing any control.Figure 3
Trajectories of model solution when only vaccination is implemented.Figure 4
Trajectories of model solution when only education is implemented.Figure 5
Trajectories of model solution when vaccination and education are implemented.Figure 6
Trajectories of model solution when all controls are implemented.Figure 7
Trajectories of model solution for varying rate of infection.
(a)
(b)Figure 8
Trajectories of model solution for varying infection rate of livestock.
(a)
(b)Figure 9
Trajectories of vaccinated and educated populations for varying rate of infection.Figure 10
Trajectories of susceptible human population for varying rates vaccination.Figure 11
Impact of education on susceptible and educated humans.


### 4.1. Findings

The following are the findings of this research:(1)
The disease-free and endemic equilibrium points of the model are stable if the reproduction number is less than one (ℛ0<1) and unstable if the reproduction number is greater than one (ℛ0>1).(2)
The sensitivity analysis shows that the most sensitive positive parameters are the pathogen-shedding rate, the anthrax infection rate, anthrax-induced death rate, and rate of livestock that lose immunity(3)
The sensitivity analysis also shows that the most sensitive negative parameters are the rate of pathogen disinfection in the environment, the rate of vaccination, and natural death rate. These parameters are inversely related to the reproduction number and indicate that a smaller amount of increment in these parameters will lead to a greater amount of reduction in the basic reproduction number while a smaller amount of decrement will cause a significant increment in the basic reproduction number.(4)
The numerical simulation of the model shows that in the absence of vaccination, the number of infectives will increase drastically.(5)
Further simulation shows that the pathogen level in the environment increases exponential when disinfection is absent.(6)
The most efficient single implementation of the controls in the model is education of susceptible human population.(7)
The most efficient pair of controls are vaccination of susceptible livestock and education of the human population.(8)
The most efficient triple combination of controls is vaccination of susceptible livestock, education of the human population, and disinfection of infected carcass and the pathogen in the environment.


### 4.2. Results and Discussion

The sensitivity analysis conducted on the basic reproduction number,ℛ0, shows that, the most sensitive positive parameters are the rate of pathogens shed into the environment (α), the rate of anthrax infection (i.e., β1), anthrax-induced death rate (π), and the rate of animals that lose immunity (Q), and these parameters are directly related to the reproduction number.Thus, a 10% increase (or decrease) in the pathogen shedding rate and the rate of infection will result in a 10% increase (or decrease) in the basic reproduction number, respectively, while a 10% increase (or decrease) in the anthrax-induced death rate and in the rate of animals that lose immunity will lead to 9.998% and 9.636% increase (or decrease) in the basic reproduction number, respectively.The rate of pathogen disinfection in the environment (ϕ2), the rate of vaccination (ρ), and the natural death rate (μ) are the most sensitive negative parameters; hence, a 10% increase (or decrease) in the rate of pathogen disinfection in the environment (ϕ2), the rate of vaccination (ρ), and the natural death rate (μ) will cause 10%, 9.972%, and 9.662% decrease (or increase) in the basic reproduction number, respectively.Figure2 shows a simulation of the basic model when all control measures are absent (that is, ρ=0,q=0,ϕ=0, and γ=0), and the number of infectives, the carcass, and the pathogen level increase, leading to more infection thus a reduction in the susceptible population. In Figure 3, the number of infectives decreases resulting in a corresponding decrease in the carcass and pathogen levels as more Susceptibles join the vaccinated class. The vaccination increases the population of careful and susceptible humans and also reduces the infective human population. When only education is implemented, it results in high level of careful susceptible humans and livestock populations as well as the overall educated human population as indicated in Figure 4. However, additional measures are required to further reduce the infective populations and the pathogen level.In Figure5, a combination of vaccination and education is implemented; this results in low level of infectives of both populations and the number of carcasses but the number of susceptible human population increases. However, the level of pathogens in the environment remains high and therefore requires additional controls.The number of human and livestock susceptibles and the careful susceptible humans (educated) as well as the vaccinated livestock increase exponentially when all control parameters are implemented, and the infective population for both livestock and humans, the number of carcasses, and the pathogen level decrease as indicated in Figure6.Increasing the rate of infection as susceptible humans come into contact with infected livestock, infected carcass, and pathogens in the environment lead to an increment in the number of infective humans as indicated in Figure7.It is established from Figure8 that the susceptible population of both human and livestock increase at smaller values of the infection rate and reduces at higher values. The reduction of susceptible might be as a result of movement from susceptible compartment to the infectious compartment as infection rate increases.In Figure9, the vaccinated and educated population increases when the rate of infection is low, and as the rate of infection increases, both vaccinated and educated population reduces. This result shows that it is much safer to vaccinate susceptibles when the infection rate is very lower or if there is not yet an outbreak.Simulation of susceptible humans at various rates of vaccination in Figure10 shows that when the vaccination rate is increased the number of human susceptible increases and when the vaccination rate is decreased the susceptible human population also decreases. The vaccinated population increases as the rate of vaccination increases, and when the rate of vaccination is lowered the population decreases.Finally, the simulation in Figure11 shows that the human susceptible population is directly related to the education rate, and the rate of education is inversely related to the educated population between time 0 and 4 but above time 4 the rate of education and the educated population are directly related. This simulation shows that the anthrax disease is best controlled in the human population when there is a long-term programme of education on the anthrax disease.Table3 shows a summary of the comparison between the controls on each state variable. Here, the control parameters vaccination, education, disinfection, and treatment are ranked as R1 to R4 to determine the most efficient parameter that influence each of the state variables positively. The values from the simulation results are used for this analysis, and the ranks R1, R2, R3, and R4 imply first, second, third, and fourth positions, respectively.Table 3
Single combination of controls.
State variables
Vaccination
Education
Disinfection
Treatment
Susceptible livestock,SL
900
R
3
1100
R
2
1900
R
1
700
R
4
Vaccinated livestock,VL
1001
R
1
10
R
2
10
R
2
10
R
2
Infected livestock,IL
1500
R
1
2700
R
2
3100
R
4
2700
R
2
Infected carcass,C
600
R
1
1100
R
2
1100
R
2
1100
R
2
Pathogen,P
1400
R
2
2300
R
3
600
R
1
2300
R
3
Susceptible humans,SH
700
R
2
1100
R
1
700
R
2
700
R
2
Educated humans,Ed
2100
R
4
3400
R
1
2490
R
3
2600
R
2
Infected humans,IH
300
R
2
400
R
3
600
R
4
200
R
1
Livestock pop. (SL+VL−IL)
401
R
1
−
1590
R
2
−
1190
R
3
−
1990
R
4
Human pop. (SH+Ed−IH)
2500
R
4
4100
R
1
2590
R
3
3100
R
2
Total rank per control
21
19
25
24In Table3, the (1st row, 1st column) is R3 and indicates that the vaccination programme is ranked 3rd in terms of increasing the number of susceptible livestock population when implemented. The (1st row, 3rd column) is ranked R1 over vaccination, education, and treatment, and this implies that disinfection is ranked first in terms of increasing the number of susceptible livestock when implemented.The best single control measure for the fight of anthrax within the livestock population is vaccination programme and this is indicated in the (9th row, 1st column) asR1, and the best single control programme for the human population is education and it is indicated in the (10th row, 2nd column) as R1. The overall best single control for the entire model is education of the human population, and it is indicated in the (last row, 2nd column) as 19 (that is, the value with the least sum of ranks). The vaccination programme is ranked first in the following: increasing the vaccinated livestock, decreasing the number of infective livestock, decreasing the number of carcasses, and reducing the level of pathogen in the environment.Table4 shows the results of pairing the controls in fighting anthrax disease. The double combinations are ranked according to their efficiency in increasing all susceptible and decreasing all infective populations of both livestock and humans as well as decreasing the number of carcasses and pathogen level in the environment. The most efficient pair in controlling the disease within the livestock population is implementing combined vaccination, and the disinfection program is indicated in (9th row, 1st column) by R1, whereas the most efficient combination pair for the human population is combined vaccination and education indicated in (10th row, 3rd column) by R1. The overall best pair in the entire model is implementing combined vaccination and education which is 19 (the least value of the total ranks) in the (last row, 3rd column).Table 4
Double combination of controls. V = vaccination, D = disinfection, T = treatment, and E = education.
State variables
V & D
V & T
V & E
T & D
T & E
E & D
Susceptible livestock,SL
1,100
R
3
900
R
4
900
R
4
1,800
R
1
1,400
R
2
1,800
R
1
Vaccinated livestock,VL
1,500
R
1
1,001
R
2
1,001
R
2
10
R
3
10
R
3
10
R
3
Infected livestock,IL
1,800
R
2
1,500
R
1
1,500
R
1
3,600
R
4
2,700
R
3
3,700
R
5
Infected carcass,C
600
R
2
600
R
2
600
R
2
1,200
R
3
100
R
1
1,200
R
3
Pathogen,P
400
R
1
1,400
R
3
1,400
R
3
600
R
2
2,300
R
4
600
R
2
Susceptible humans,SH
600
R
4
600
R
4
2,600
R
1
600
R
4
1,600
R
3
2,200
R
2
Educated humans,Ed
2,500
R
5
2,500
R
5
4,900
R
1
3,400
R
4
4,400
R
2
4,001
R
3
Infected humans,IH
300
R
2
10
R
1
300
R
2
400
R
3
10
R
1
600
R
4
Livestock pop. (SL+VL−IL)
800
R
1
401
R
2
401
R
2
−
1,790
R
4
−
1,290
R
3
−
1,890
R
5
Human pop. (SH+Ed−IH)
2,800
R
6
3,091
R
5
7,200
R
1
3,600
R
4
5,990
R
2
5,601
R
3
Total rank per control
27
29
19
33
25
32Table5 shows the implementation of three different controls at the same time. The results show that the most efficient triple combination of controls within the livestock population is implementing combined vaccination, education, and disinfection programme, and this is indicated in (9th row, 3rd column) by R1. The most efficient triple combination for the human population is combining vaccination, education and disinfection as indicated in (10th row, 3rd column) by R1. The overall best triple combination in the entire model is implementing vaccination, education, and disinfection indicated in the (last row, 3rd column) by 12 (the least value of the total ranks).Table 5
Triple combination of controls. V = vaccination, D = disinfection, T = treatment, and E = education.
State variables
V, T & E
V, T & D
V, E & D
E, T & D
Susceptible livestock,SL
1,000
R
4
1,010
R
3
1,050
R
2
1,800
R
1
Vaccinated livestock,VL
1,000
R
2
1,500
R
1
1,500
R
1
10
R
3
Infected livestock,IL
1,500
R
1
1,800
R
2
1,800
R
2
3,900
R
3
Infected carcass,C
500
R
1
500
R
1
500
R
1
1,000
R
2
Pathogen,P
1,200
R
3
400
R
1
400
R
1
700
R
2
Susceptible humans,SH
3,100
R
3
600
R
4
4,400
R
1
3,900
R
2
Educated humans,Ed
6,000
R
1
3,010
R
2
6,000
R
1
6,000
R
3
Infected humans,IH
10
R
1
10
R
1
400
R
2
400
R
2
Livestock pop. (SL+VL−IL)
500
R
3
710
R
2
750
R
1
−
2,090
R
4
Human pop. (SH+Ed−IH)
9,090
R
3
3,600
R
4
10,000
R
1
9,500
R
2
Total rank per control
22
21
13
24


## 4.1. Findings

The following are the findings of this research:(1)
The disease-free and endemic equilibrium points of the model are stable if the reproduction number is less than one (ℛ0<1) and unstable if the reproduction number is greater than one (ℛ0>1).(2)
The sensitivity analysis shows that the most sensitive positive parameters are the pathogen-shedding rate, the anthrax infection rate, anthrax-induced death rate, and rate of livestock that lose immunity(3)
The sensitivity analysis also shows that the most sensitive negative parameters are the rate of pathogen disinfection in the environment, the rate of vaccination, and natural death rate. These parameters are inversely related to the reproduction number and indicate that a smaller amount of increment in these parameters will lead to a greater amount of reduction in the basic reproduction number while a smaller amount of decrement will cause a significant increment in the basic reproduction number.(4)
The numerical simulation of the model shows that in the absence of vaccination, the number of infectives will increase drastically.(5)
Further simulation shows that the pathogen level in the environment increases exponential when disinfection is absent.(6)
The most efficient single implementation of the controls in the model is education of susceptible human population.(7)
The most efficient pair of controls are vaccination of susceptible livestock and education of the human population.(8)
The most efficient triple combination of controls is vaccination of susceptible livestock, education of the human population, and disinfection of infected carcass and the pathogen in the environment.


## 4.2. Results and Discussion

The sensitivity analysis conducted on the basic reproduction number,ℛ0, shows that, the most sensitive positive parameters are the rate of pathogens shed into the environment (α), the rate of anthrax infection (i.e., β1), anthrax-induced death rate (π), and the rate of animals that lose immunity (Q), and these parameters are directly related to the reproduction number.Thus, a 10% increase (or decrease) in the pathogen shedding rate and the rate of infection will result in a 10% increase (or decrease) in the basic reproduction number, respectively, while a 10% increase (or decrease) in the anthrax-induced death rate and in the rate of animals that lose immunity will lead to 9.998% and 9.636% increase (or decrease) in the basic reproduction number, respectively.The rate of pathogen disinfection in the environment (ϕ2), the rate of vaccination (ρ), and the natural death rate (μ) are the most sensitive negative parameters; hence, a 10% increase (or decrease) in the rate of pathogen disinfection in the environment (ϕ2), the rate of vaccination (ρ), and the natural death rate (μ) will cause 10%, 9.972%, and 9.662% decrease (or increase) in the basic reproduction number, respectively.Figure2 shows a simulation of the basic model when all control measures are absent (that is, ρ=0,q=0,ϕ=0, and γ=0), and the number of infectives, the carcass, and the pathogen level increase, leading to more infection thus a reduction in the susceptible population. In Figure 3, the number of infectives decreases resulting in a corresponding decrease in the carcass and pathogen levels as more Susceptibles join the vaccinated class. The vaccination increases the population of careful and susceptible humans and also reduces the infective human population. When only education is implemented, it results in high level of careful susceptible humans and livestock populations as well as the overall educated human population as indicated in Figure 4. However, additional measures are required to further reduce the infective populations and the pathogen level.In Figure5, a combination of vaccination and education is implemented; this results in low level of infectives of both populations and the number of carcasses but the number of susceptible human population increases. However, the level of pathogens in the environment remains high and therefore requires additional controls.The number of human and livestock susceptibles and the careful susceptible humans (educated) as well as the vaccinated livestock increase exponentially when all control parameters are implemented, and the infective population for both livestock and humans, the number of carcasses, and the pathogen level decrease as indicated in Figure6.Increasing the rate of infection as susceptible humans come into contact with infected livestock, infected carcass, and pathogens in the environment lead to an increment in the number of infective humans as indicated in Figure7.It is established from Figure8 that the susceptible population of both human and livestock increase at smaller values of the infection rate and reduces at higher values. The reduction of susceptible might be as a result of movement from susceptible compartment to the infectious compartment as infection rate increases.In Figure9, the vaccinated and educated population increases when the rate of infection is low, and as the rate of infection increases, both vaccinated and educated population reduces. This result shows that it is much safer to vaccinate susceptibles when the infection rate is very lower or if there is not yet an outbreak.Simulation of susceptible humans at various rates of vaccination in Figure10 shows that when the vaccination rate is increased the number of human susceptible increases and when the vaccination rate is decreased the susceptible human population also decreases. The vaccinated population increases as the rate of vaccination increases, and when the rate of vaccination is lowered the population decreases.Finally, the simulation in Figure11 shows that the human susceptible population is directly related to the education rate, and the rate of education is inversely related to the educated population between time 0 and 4 but above time 4 the rate of education and the educated population are directly related. This simulation shows that the anthrax disease is best controlled in the human population when there is a long-term programme of education on the anthrax disease.Table3 shows a summary of the comparison between the controls on each state variable. Here, the control parameters vaccination, education, disinfection, and treatment are ranked as R1 to R4 to determine the most efficient parameter that influence each of the state variables positively. The values from the simulation results are used for this analysis, and the ranks R1, R2, R3, and R4 imply first, second, third, and fourth positions, respectively.Table 3
Single combination of controls.
State variables
Vaccination
Education
Disinfection
Treatment
Susceptible livestock,SL
900
R
3
1100
R
2
1900
R
1
700
R
4
Vaccinated livestock,VL
1001
R
1
10
R
2
10
R
2
10
R
2
Infected livestock,IL
1500
R
1
2700
R
2
3100
R
4
2700
R
2
Infected carcass,C
600
R
1
1100
R
2
1100
R
2
1100
R
2
Pathogen,P
1400
R
2
2300
R
3
600
R
1
2300
R
3
Susceptible humans,SH
700
R
2
1100
R
1
700
R
2
700
R
2
Educated humans,Ed
2100
R
4
3400
R
1
2490
R
3
2600
R
2
Infected humans,IH
300
R
2
400
R
3
600
R
4
200
R
1
Livestock pop. (SL+VL−IL)
401
R
1
−
1590
R
2
−
1190
R
3
−
1990
R
4
Human pop. (SH+Ed−IH)
2500
R
4
4100
R
1
2590
R
3
3100
R
2
Total rank per control
21
19
25
24In Table3, the (1st row, 1st column) is R3 and indicates that the vaccination programme is ranked 3rd in terms of increasing the number of susceptible livestock population when implemented. The (1st row, 3rd column) is ranked R1 over vaccination, education, and treatment, and this implies that disinfection is ranked first in terms of increasing the number of susceptible livestock when implemented.The best single control measure for the fight of anthrax within the livestock population is vaccination programme and this is indicated in the (9th row, 1st column) asR1, and the best single control programme for the human population is education and it is indicated in the (10th row, 2nd column) as R1. The overall best single control for the entire model is education of the human population, and it is indicated in the (last row, 2nd column) as 19 (that is, the value with the least sum of ranks). The vaccination programme is ranked first in the following: increasing the vaccinated livestock, decreasing the number of infective livestock, decreasing the number of carcasses, and reducing the level of pathogen in the environment.Table4 shows the results of pairing the controls in fighting anthrax disease. The double combinations are ranked according to their efficiency in increasing all susceptible and decreasing all infective populations of both livestock and humans as well as decreasing the number of carcasses and pathogen level in the environment. The most efficient pair in controlling the disease within the livestock population is implementing combined vaccination, and the disinfection program is indicated in (9th row, 1st column) by R1, whereas the most efficient combination pair for the human population is combined vaccination and education indicated in (10th row, 3rd column) by R1. The overall best pair in the entire model is implementing combined vaccination and education which is 19 (the least value of the total ranks) in the (last row, 3rd column).Table 4
Double combination of controls. V = vaccination, D = disinfection, T = treatment, and E = education.
State variables
V & D
V & T
V & E
T & D
T & E
E & D
Susceptible livestock,SL
1,100
R
3
900
R
4
900
R
4
1,800
R
1
1,400
R
2
1,800
R
1
Vaccinated livestock,VL
1,500
R
1
1,001
R
2
1,001
R
2
10
R
3
10
R
3
10
R
3
Infected livestock,IL
1,800
R
2
1,500
R
1
1,500
R
1
3,600
R
4
2,700
R
3
3,700
R
5
Infected carcass,C
600
R
2
600
R
2
600
R
2
1,200
R
3
100
R
1
1,200
R
3
Pathogen,P
400
R
1
1,400
R
3
1,400
R
3
600
R
2
2,300
R
4
600
R
2
Susceptible humans,SH
600
R
4
600
R
4
2,600
R
1
600
R
4
1,600
R
3
2,200
R
2
Educated humans,Ed
2,500
R
5
2,500
R
5
4,900
R
1
3,400
R
4
4,400
R
2
4,001
R
3
Infected humans,IH
300
R
2
10
R
1
300
R
2
400
R
3
10
R
1
600
R
4
Livestock pop. (SL+VL−IL)
800
R
1
401
R
2
401
R
2
−
1,790
R
4
−
1,290
R
3
−
1,890
R
5
Human pop. (SH+Ed−IH)
2,800
R
6
3,091
R
5
7,200
R
1
3,600
R
4
5,990
R
2
5,601
R
3
Total rank per control
27
29
19
33
25
32Table5 shows the implementation of three different controls at the same time. The results show that the most efficient triple combination of controls within the livestock population is implementing combined vaccination, education, and disinfection programme, and this is indicated in (9th row, 3rd column) by R1. The most efficient triple combination for the human population is combining vaccination, education and disinfection as indicated in (10th row, 3rd column) by R1. The overall best triple combination in the entire model is implementing vaccination, education, and disinfection indicated in the (last row, 3rd column) by 12 (the least value of the total ranks).Table 5
Triple combination of controls. V = vaccination, D = disinfection, T = treatment, and E = education.
State variables
V, T & E
V, T & D
V, E & D
E, T & D
Susceptible livestock,SL
1,000
R
4
1,010
R
3
1,050
R
2
1,800
R
1
Vaccinated livestock,VL
1,000
R
2
1,500
R
1
1,500
R
1
10
R
3
Infected livestock,IL
1,500
R
1
1,800
R
2
1,800
R
2
3,900
R
3
Infected carcass,C
500
R
1
500
R
1
500
R
1
1,000
R
2
Pathogen,P
1,200
R
3
400
R
1
400
R
1
700
R
2
Susceptible humans,SH
3,100
R
3
600
R
4
4,400
R
1
3,900
R
2
Educated humans,Ed
6,000
R
1
3,010
R
2
6,000
R
1
6,000
R
3
Infected humans,IH
10
R
1
10
R
1
400
R
2
400
R
2
Livestock pop. (SL+VL−IL)
500
R
3
710
R
2
750
R
1
−
2,090
R
4
Human pop. (SH+Ed−IH)
9,090
R
3
3,600
R
4
10,000
R
1
9,500
R
2
Total rank per control
22
21
13
24


## 5. Conclusions and Recommendations

### 5.1. Conclusions

In this work, a nonlinear mathematical model is proposed to study the transmission of anthrax in both livestock and human populations in the presence of vaccination of livestock, education of the human population, treatment of infectives, and disinfection of carcasses and pathogens in the environment.In order to control the disease from the populations, the parameters such as the pathogen-shedding rate, the anthrax infection rate, the anthrax-induced death rate, and the rate of animals that lose immunity which have the most sensitive positive indexes should be reduced since a reduction in these parameters will lead to a decrease in the basic reproduction number. In addition, the most sensitive negative parameters such as the rate of disinfection of the pathogen in the environment, the rate of vaccination of susceptible livestock, and the natural death rate of livestock must be increased in other to control the disease, and an increase in these parameters will lead to a decrease in the basic reproduction number. The numerical simulation of the model confirms the results of the sensitivity analysis.The study of different control parameters combined to fight the anthrax disease shows that the most efficient single implementation of the control parameters is education while the combination of vaccination and education is the most efficient double implementation of the controls. In addition, implementing the combination of vaccination, education, and disinfection is the best option for triple combination of controls, and lastly, implementing all the controls is an excellent option of controlling the transmission of anthrax in both livestock and human population, and all these options depend solely on the availability of resources.


### 5.2. Recommendations

The transmission of anthrax can greatly be minimised or possibly be eradicated in the population by ensuring that the reproduction number is less than one. This can be achieved through the following recommendations:(1)
Intensify long-term awareness creation on the disease and the practice of prophylaxes through education of the human population.(2)
Carry out vaccination of newborns and susceptible livestock as well as constant monitoring and continual vaccination of livestock in the vaccinated class, and these will prevent the livestock from losing immunity and becoming susceptible again, thereby reducing the number of susceptible as the vaccinated population increases.(3)
Infected carcass should be thoroughly destroyed and carefully disposed off. In addition, the infected area should be disinfected to kill the pathogens, and these practices will reduce the infected population drastically.(4)
Livestock population should be randomly screened for anthrax infection and the infected animals should be rapidly treated, and this will reduce the pathogen that goes into the environment since anthrax-induced deaths would be reduced.


## 5.1. Conclusions

In this work, a nonlinear mathematical model is proposed to study the transmission of anthrax in both livestock and human populations in the presence of vaccination of livestock, education of the human population, treatment of infectives, and disinfection of carcasses and pathogens in the environment.In order to control the disease from the populations, the parameters such as the pathogen-shedding rate, the anthrax infection rate, the anthrax-induced death rate, and the rate of animals that lose immunity which have the most sensitive positive indexes should be reduced since a reduction in these parameters will lead to a decrease in the basic reproduction number. In addition, the most sensitive negative parameters such as the rate of disinfection of the pathogen in the environment, the rate of vaccination of susceptible livestock, and the natural death rate of livestock must be increased in other to control the disease, and an increase in these parameters will lead to a decrease in the basic reproduction number. The numerical simulation of the model confirms the results of the sensitivity analysis.The study of different control parameters combined to fight the anthrax disease shows that the most efficient single implementation of the control parameters is education while the combination of vaccination and education is the most efficient double implementation of the controls. In addition, implementing the combination of vaccination, education, and disinfection is the best option for triple combination of controls, and lastly, implementing all the controls is an excellent option of controlling the transmission of anthrax in both livestock and human population, and all these options depend solely on the availability of resources.


## 5.2. Recommendations

The transmission of anthrax can greatly be minimised or possibly be eradicated in the population by ensuring that the reproduction number is less than one. This can be achieved through the following recommendations:(1)
Intensify long-term awareness creation on the disease and the practice of prophylaxes through education of the human population.(2)
Carry out vaccination of newborns and susceptible livestock as well as constant monitoring and continual vaccination of livestock in the vaccinated class, and these will prevent the livestock from losing immunity and becoming susceptible again, thereby reducing the number of susceptible as the vaccinated population increases.(3)
Infected carcass should be thoroughly destroyed and carefully disposed off. In addition, the infected area should be disinfected to kill the pathogens, and these practices will reduce the infected population drastically.(4)
Livestock population should be randomly screened for anthrax infection and the infected animals should be rapidly treated, and this will reduce the pathogen that goes into the environment since anthrax-induced deaths would be reduced.



---

*Source: 1581358-2020-02-18.xml*