Scientific Research

An Academic Publisher

Epidemiological Study and Optimal Control for Lumpy Skin Disease (LSD) in Ethiopia

**Author(s)**Leave a comment

*v*≤ 0.85. It is shown that the vaccination as a control strategy significantly reduced the effects of LSD on the cattle population if properly managed and that an optimal performance of the control strategy adopted hererin is achieved at an approximate value of

*v*= 0.6.

KEYWORDS

1. Introduction

LSD has an economical importance because of permanent damage to hides. It has been reported ( [1] [2] ) that 10% of animals exposed to LSD experience abortion. According to [3] , lumpy skin disease (LSD) is an exhaustive viral disease and a pox disease of cattle. It is characterized by fever, nodules on the skin, mucous membranes in internal organs, enlarged lymph nodes and odema of the skin ( [3] [4] ). Due to chronic debility in affected animals high economic losses occur. The clinical syndrome of LSD was first described in Zambia in 1929. More cases also occurred between 1943 and 1945 in Botswana, Zimbabwe and the republic of South Africa ( [1] [3] [5] [6] ). Transmission of LSD is by insect vectors (through biting flies) found among cattle sharing similar grazing and watering areas. And those that congregate in the same barn may be contaminated by the viruses in saliva of the infected animals or ingestion of already contaminated food or by teratogenic agents and suckling calves ( [7] [8] ). The incidence of LSD occurrence is high during wet seasons [3] .

LSD was first observed in Ethiopia in the northwestern part of the country (southwest of Lake Tana) in 1983, [1] . It occurs in all agro-climatic conditions and has the potential to extend its boundaries. Good understanding of epidemiology, economic significance and control mechanisms of the disease will enhance suitable control measures. The disease could be diagnosed using appropriate serological and molecular techniques. Effective control measures can be achieved through mass vaccination, separation and culling of infected animals, ( [9] [10] [11] ). It is therefore important to formulate an LSD model with optimal control policy to eradicate LSD epidemic and minimize the cost associated for vaccination strategy [3] .

Most studies of epidemic control of LSD focus on increasing immunization coverage (vaccination coverage) in a population to control the disease, but they do not consider how this parameter affects the strategy over a period of time. Though some of these studies have considered vaccination strategy at different levels of immunization, none of them has studied optimal vaccination strategy for LSD disease in Ethiopia. Hence the general objective of this study aims at addressing these shortfalls in order to make the number of infected animals as small as possible.

2. Epidemiological Model Formulation

SEIR Model and Equation Formulation

All Epidemiological ODE models of disease transmission are continuous time compartment level models [12] . The model considered herein is compartmentalized into:

S: susceptible portion of population,

E: exposed population.

I: the infected population.

R: the recovered population.

Newly born population merge into the S compartment at the rate of alpha ( $\alpha $ ), which represents the birth rate. The infected population are those individuals who can contract the disease. However for the SEIR model, we also consider the exposed population that is those individuals who have the disease but are not yet infectious. After treatment the infected population will move into R-Compartment, at the rate of omega $\omega $ i.e. the recovery rate. Cattles can lose their immunity after some time. Let us assume a rate of $\theta $ . As a result, that population moves from R into the S compartment. Hence the following S-E-I-R-S Epidemiological model without control is shown in Figure 1 as:

The system is represented by the temporal vectors S(t), E(t), I(t), R(t); and is characterized by the following assumptions

・ The population is homogeneous.

・ No demographics,(immigration and emigration) consideration

・ Birth rate and natural death rate are constant. Mathematically this means

$S\left(t\right)+E\left(t\right)+I\left(t\right)+R\left(t\right)=N$ (1)

where N is the total population. This facilitates the expression of one of the variables in terms of the others and the conserved constant N. Table 1 shows the description of the model parameters and variables applied in this study.

$\frac{\text{d}S}{\text{d}t}=\alpha N-\beta SI+\theta R-\varphi S$ (2)

$\frac{\text{d}E}{\text{d}t}=\beta SI-\gamma E-\varphi E$ (3)

$\frac{\text{d}I}{\text{d}t}=\gamma E-\omega I-\varphi I$ (4)

Figure 1. SEIRS Model without Control.

Table 1. Descriptions for variables for the model.

$\frac{\text{d}R}{\text{d}t}=\omega I-\theta R-\varphi R$ (5)

Note that the natural death rate on each of the compartment for this paper is the same, which is

$\varphi S+\varphi E+\varphi I+\varphi R=\varphi \left(S+E+I+R\right)=\varphi N$ (6)

3. Optimal Control

Optimal control deals with the problem of finding a control law for a given system such that a certain optimal criterion is achieved. A control problem includes a function of state and control variable. An optimal control is a set of differential equations describing the paths of the control variables that minimize the cost function. The control can be derived using Pontryagins maximum principles (a necessary condition also known as Pontryagins minimum principle or simply Pontryagins Principle), or by solving the Hamilton-Jacobi-Bellamn Equations (a sufficient condition), [13] . This boundary-value problem actually has a special structure because it arises from taking the derivative of a Hamiltonian. Thus, the resulting dynamical system is a Hamiltonian system of the form

$\stackrel{\dot{}}{x}=\frac{\partial H}{\partial \lambda}$ (7)

$\stackrel{\dot{}}{\lambda}=-\frac{\partial H}{\partial x}$ (8)

where

$H=L+{\lambda}^{\text{T}}a-{\mu}^{\text{T}}b$ (9)

is an augmented Hamiltonian [14] . For this study, vaccination is applied as an optimal control and is denoted as u.

Numerical Solutions of Optimal Control Problems

Consider the optimal control problem

$\mathrm{max}U{\displaystyle \underset{{t}_{0}}{\overset{{t}_{1}}{\int}}f\left(t,x\left(t\right),u\left(t\right)\right)\text{d}t}$ (10)

$\text{subject}\text{\hspace{0.17em}}\text{to}:\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}^{\prime}=g\left(t,x\left(t\right),u\left(t\right)\right)$ (11)

$x\left({t}_{0}\right)={x}_{0},\text{\hspace{0.17em}}x\left({t}_{1}\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{free}$ (12)

Any solution to the above optimal control problem must also satisfy

${x}^{\prime}\left(t\right)=g\left(t,x\left(t\right),u\left(t\right)\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\left({t}_{0}\right)={x}_{0}$ (13)

$\lambda \left(t\right)=-\frac{\delta H}{\delta x}=-\left({f}_{x}\left(t,x,u\right)+\lambda \left(t\right){g}_{x}\left(t,x,u\right)\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\lambda \left({t}_{1}\right)=0$ (14)

$0=\frac{\delta H}{\delta x}=-\left(f\left(t,x,u\right)+\lambda \left(t\right){g}_{u}\left(t,x,u\right)\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{at}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{u}^{*}$ (15)

Equation (6), the optimal condition, can usually be manipulated to find a representation of u* in terms of t, x, and $\lambda $ .

An optimal control with bounded control can be written in the form

${\mathrm{max}}_{u}J\left(x\left(t\right),\left(u\right)\right)={\displaystyle {\int}_{{t}_{0}}^{{t}_{f}}f\left(t,x\left(t\right),u\left(t\right)\right)\text{d}t}$ (16)

Subject to $\{\begin{array}{l}{x}^{\prime}=g\left(t,x\left(t\right),u\left(t\right)\right)\\ x\left({t}_{0}\right)<{x}_{0}\\ a\le u\left(t\right)\le b\end{array}$ (17)

where a, b are fixed real constants and a < b.

Proposition (Necessary conditions):- If ${u}^{*}\left(t\right)$ and ${x}^{*}\left(t\right)$ are optimal for Equations (16, 17), then there exists a piecewise differentiable adjoint variable $\lambda \left(t\right)$ such that

$H\left(t,{x}^{*}\left(t\right),u\left(t\right),\lambda \left(t\right)\right)\le H\left(t,x\left(t\right),{u}^{*}\left(t\right),\lambda \left(t\right)\right)$ (18)

For all controls u at each time t, where H is the Hamiltonian, [13] .

${\lambda}^{\prime}\left(t\right)=-\frac{\delta H\left({x}^{*}\left(t\right),{u}^{*}\left(t\right),\lambda \left(t\right)\right)}{\delta x}$ (adjoint condition) (19)

$\lambda \left({f}_{t}\right)=0$ (transversality condition) (20)

By an adaptation of the Pontryagin’s Maximum Principle (PMP), the OC must satisfy (optimal condition):

${u}^{*}=\{\begin{array}{l}a\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\frac{\delta H}{\delta u}<0,\\ a<\frac{\delta H}{\delta u}<b,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\frac{\delta H}{\delta u}=0,\\ b<\frac{\delta H}{\delta u}<a,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\frac{\delta H}{\delta u}>0,\end{array}$ (21)

If we have a minimization problem, then ${u}^{*}$ is chosen to minimize H point wise.

If $\frac{{\partial}^{2}H}{\partial {u}^{2}}<0$ at ${u}^{*}$ then u is the maximization of the problem, while $\frac{{\partial}^{2}H}{\partial {u}^{2}}<0$ goes with minimization.

The optimal control $\frac{\delta H}{\delta u}=0$ obtained without truncation, and is bounded by a and b such that:

${u}^{*}\left(t\right)=\mathrm{min}\left\{a,\mathrm{max}\left\{b,\frac{\delta H}{\delta u}\right\}\right\}$ (22)

Runge-Kutta (RK4) routine is applied for the numerical solution. It is required that:

$\underset{i=1}{\overset{n}{\sum}}\left|u\left(i\right)-{u}_{old}\left(i\right)\right|$ be small, in addition to guaranteeing that the relative error

between the current and the previous computed scalars be negligible, i.e.

$\frac{\left|u-{u}_{old}\right|}{\left|u\right|}\le \delta $ (23)

where $\delta $ is the accepted tolerance? The optimal control problem can now be written in a standard form:

$\mathrm{max}{\displaystyle \underset{0}{\overset{1}{\int}}Ax\left(t\right)-B{u}^{2}\left(t\right)\text{d}t}$ (24)

Subject to:- ${x}^{\prime}\left(t\right)=-\frac{1}{2}{x}^{2}\left(t\right)+Cu\left(t\right)$ (25)

$x\left(0\right)={x}_{0}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{fixed}\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\left(1\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{free}$

In order for this to be a maximization problem, we require B > 0. The Hamiltonian given by

$H=Ax-B{u}^{2}-\frac{1}{2}+C\lambda u$ (26)

By applying the optimal condition

$0=\frac{\delta H}{\delta u}=-2Bu+C\lambda \Rightarrow {u}^{*}=\frac{C\lambda}{2B}$ (27)

${x}^{\prime}\left(t\right)=-\frac{1}{2}{x}^{2}+Cu,\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\left(0\right)={x}_{0}$ (28)

${\lambda}^{\prime}\left(t\right)=-A+x\lambda ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\lambda \left(1\right)=0$ (29)

The Runge-Kutta forward sweep yields the transient profile of the dependent variable. The new ${u}^{*}$ is calculated by setting

${u}_{1}=C\ast \frac{\lambda}{2\ast B}$ (30)

$u=0.5\ast \left({u}_{1}+{u}_{old}\right)$ (31)

4. ODE Formulation for Optimal Control

4.1. Modification of SEIRS Model

A typical optimal control problem requires a performance index or cost function $\left(J\left[x\left(t\right),u\left(t\right)\right]\right)$ , a set of state variables $\left(x\left(t\right)\in X\right)$ and a set of control variables $\left(u\left(t\right)\in U\right)$ in a time t, ${t}_{0}\le t\le {t}_{f}$ . A major goal of the work reported herein, is to find a piecewise continuous control u(t) (vaccination) and the associated state variable x(t) to minimize a given objective functional. This strategy will directly impact on the susceptible “S” and the recovered “R” compartments in the overall mathematical formulation as shown below.

$\frac{\text{d}S}{\text{d}t}=\alpha \left(1-{v}_{1}\right)N-\beta SI-{v}_{2}S+\theta R-\varphi S$ (32)

$\frac{\text{d}E}{\text{d}t}=\beta SI-\gamma E-\varphi E$ (33)

$\frac{\text{d}I}{\text{d}t}=\gamma E-\omega I-\varphi I$ (34)

$\frac{\text{d}R}{\text{d}t}=\alpha {v}_{1}N+{v}_{2}S+\omega I-\theta R-\varphi R$ (35)

In the above equations, vaccination is applied during g birth ${v}_{1}$ , and for populations in the S-compartment ${v}_{2}$ , thus vaccinated populations moved to recovered compartment as $\alpha {v}_{1}N$ and ${v}_{2}S$ . However, for this paper, we applied vaccination on S-compartment only as ( $vS$ ), and the formulation is given as below.

$\frac{\text{d}S}{\text{d}t}=\alpha N-\beta SI-vS+\theta R-\varphi S$ (36)

$\frac{\text{d}E}{\text{d}t}=\beta SI-\gamma E-\varphi E$ (37)

$\frac{\text{d}I}{\text{d}t}=\gamma E-\omega I-\varphi I$ (38)

$\frac{\text{d}R}{\text{d}t}=vS+\omega I-\theta R-\varphi R$ (39)

With initial conditions,

$\begin{array}{l}S\left(0\right)={S}_{0}\ge 0,\text{\hspace{0.17em}}E\left(0\right)={E}_{0}\ge 0\\ I\left(0\right)={I}_{0}\ge 0,\text{\hspace{0.17em}}R\left(0\right)={R}_{0}\ge 0\end{array}$

The above equations are illustrated graphically in Figure 2 as:

The optimal problem (which is the main objectives of this paper) is to minimize the objective (cost) functional (J) by considering the costs of vaccination for susceptible population given by:

$J\left(u\right)={\displaystyle \underset{0}{\overset{t}{\int}}\left(AI+\frac{1}{2}B{u}^{2}\right)\text{d}t}$ (40)

where:-

t: represents the vaccination period

A: a parameter that balancing cost factors due to the size of ineffective

B: represents the weight attached on the cost of vaccination

A and B can also be weight parameters describing the comparative importance of the two terms (disease burden and cost) in the functional respectively.

A high value of A means that it is more important to reduce the disease burden than to reduce the vaccination costs. A less value of B means that it is more important to reduce the vaccination costs than to reduce the disease burden. The control u is the percentage of the susceptible that is vaccinated per unit time. Thus, u lies between 0 and 1 while u_{max} will depend on the amount of resources available to implement each of the control measures.

If u = 0, then no vaccination is done which is uncontrolled. If u = 1, then all susceptible population is vaccinated.

The rate of vaccination is assumed to take values in [0, 0.9] instead of [0, 1] to eliminate the case where the entire susceptible population is vaccinated.

The vaccination cost could include the cost of the vaccine, cost of syringes, cost of safety boxes, the vaccine storage cost, and other related overheads.

Figure 2. Optimal control with vaccination.

Hence, the Optimal Control problem of the nonlinear dynamics of SEIRS LSD epidemic model studied herein is given by;

${\mathrm{min}}_{u}J\left(u\right)=\mathrm{min}{\displaystyle \underset{0}{\overset{t}{\int}}\left(AI+\frac{1}{2}B{u}^{2}\right)\text{d}t}$ (41)

where $U=u\left(t\right):0\le u\le {u}_{\mathrm{max}}<1,t\in \left[0,t\right]$ .

For the Hamiltonian function:

$\begin{array}{l}H\left(t,S,E,I,R,{\lambda}_{1},{\lambda}_{2},{\lambda}_{3},{\lambda}_{4}\right)\\ =f\left(t,S,E,I,R,u,v\right)+{\displaystyle \underset{i=1}{\overset{4}{\sum}}{\lambda}_{i}{g}_{i}\left(t,S,E,I,R,u,v\right)}\\ =f\left(t,S,E,I,R,u,v\right)+{\lambda}_{1}\frac{\text{d}S}{\text{d}t}+{\lambda}_{2}\frac{\text{d}E}{\text{d}t}+{\lambda}_{3}\frac{\text{d}I}{\text{d}t}+{\lambda}_{4}\frac{\text{d}R}{\text{d}t}\end{array}$ (42)

Then by substituting each of the derivatives in Equation (41) into Equation (42) and by applying PMP to the Hamiltonian, we have:

$\begin{array}{c}H=\left(AI+\frac{1}{2}B{u}^{2}\right)+{\lambda}_{1}\left(\alpha N-\beta SI-uS+\theta R-\varphi S\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{\lambda}_{2}\left(\beta SI-\gamma E-\varphi E\right)+{\lambda}_{3}\left(\gamma E-\omega I-\varphi I\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{\lambda}_{4}\left(uS+\omega I-\theta R-\varphi R\right)\end{array}$ (43)

The Adjoint functions with respect to the dependent variables are:

For S:

$\begin{array}{l}{\lambda}_{1}^{*}=-\frac{\delta H}{\delta S}\\ \Rightarrow {\lambda}_{1}^{*}=-\left[{\lambda}_{1}\left(\beta I-u-\varphi \right)+{\lambda}_{2}\beta I+{\lambda}_{4}u\right]\\ \Rightarrow {\lambda}_{1}^{*}={\lambda}_{1}\left(-\beta I+u+\varphi \right)-{\lambda}_{2}\beta I-{\lambda}_{4}u\end{array}$ (44)

Adjoint function with respect to E

$\begin{array}{l}{\lambda}_{2}^{*}=-\frac{\delta H}{\delta E}\\ \Rightarrow {\lambda}_{2}^{*}=-\left[{\lambda}_{2}\left(-\gamma -\varphi \right)+{\lambda}_{3}\gamma \right]\\ \Rightarrow {\lambda}_{2}^{*}={\lambda}_{2}\left(\gamma +\varphi \right)-{\lambda}_{3}\gamma \end{array}$ (45)

Adjoint function with respect to I

$\begin{array}{l}{\lambda}_{3}^{*}=-\frac{\delta H}{\delta I}\\ \Rightarrow {\lambda}_{3}^{*}=-\left[A+{\lambda}_{1}\left(-\beta S\right)+{\lambda}_{2}\beta S+{\lambda}_{3}\left(-\varphi -\omega \right)+{\lambda}_{4}\omega \right]\\ \Rightarrow {\lambda}_{3}^{*}=-A+{\lambda}_{1}\beta S-{\lambda}_{2}\beta S+{\lambda}_{3}\left(\varphi +\omega \right)-{\lambda}_{4}\omega \end{array}$ (46)

Adjoint function with respect to R

$\begin{array}{l}{\lambda}_{4}^{*}=-\frac{\delta H}{\delta R}\\ \Rightarrow {\lambda}_{4}^{*}=-\left[{\lambda}_{1}\left(\theta \right)+{\lambda}_{4}\left(-\theta -\varphi \right)\right]\\ \Rightarrow {\lambda}_{4}^{*}=-{\lambda}_{1}\theta +{\lambda}_{4}\left(\theta +\varphi \right)\end{array}$ (47)

From the Equation (43 to 47); we obtain:

$\begin{array}{l}{\lambda}_{1}^{*}={\lambda}_{1}\left(-\beta I+u+\varphi \right)-{\lambda}_{2}\beta I-{\lambda}_{4}u\\ {\lambda}_{2}^{*}={\lambda}_{2}\left(\gamma +\varphi \right)-{\lambda}_{3}\gamma \\ {\lambda}_{3}^{*}=-A+{\lambda}_{1}\beta S-{\lambda}_{2}\beta S+{\lambda}_{3}\left(\varphi +\omega \right)-{\lambda}_{4}\omega \\ {\lambda}_{4}^{*}=-{\lambda}_{1}\theta +{\lambda}_{4}\left(\theta +\varphi \right)\end{array}$ (48)

With transversality conditions we have ${\lambda}_{i}\left(T\right)=0,i=1,2,3,4$

Next, we consider an optimal condition for determining the characterization of the optimal control ${u}^{*}$ . The optimal condition by differentiating the Hamiltonian H with respect to u based on the Equation (27) for each control ${u}^{*}$ . We have the following optimal system:

$\frac{\delta H}{\delta u}=0$ at $u={u}^{*}$ $\Rightarrow Bu+{\lambda}_{1}\left(-S\right)$ , at $u={u}^{*}$

then we have ${u}^{*}\left(t\right)={\lambda}_{1}\frac{S}{B}$ taking into account the bound ${u}^{*}$ and its characterization obtain;

${u}^{*}=\{\begin{array}{l}0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\frac{\delta H}{\delta u}<0,\\ {\lambda}_{1}\frac{S}{B}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\frac{\delta H}{\delta u}=0,\\ 0\text{\hspace{0.17em}}.9\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}\frac{\delta H}{\delta u}>0\end{array}$ (49)

the optimal control ${u}^{*}$ can now be put in a compact form:

${u}^{*}=\mathrm{max}\left(\mathrm{min}\left({\lambda}_{1}\frac{S}{B},0.9\right),0\right)$ (50)

Hence using the characterization of the optimal control, we now have following optimal system

$\begin{array}{l}{S}^{\prime}=\alpha N-\beta SI-uS+\theta R-\varphi S\\ {E}^{\prime}=\beta SI-\gamma E-\varphi E\\ {I}^{\prime}=\gamma E-\omega I-\varphi I\\ {R}^{\prime}=uS+\omega I-\theta R-\varphi R\end{array}$ (51)

$S\left(0\right)={S}_{0}\ge 0,\text{\hspace{0.17em}}E\left(0\right)={E}_{0}\ge 0,\text{\hspace{0.17em}}I\left(0\right)={I}_{0}\ge 0,\text{\hspace{0.17em}}R\left(0\right)={R}_{0}\ge 0$

To determine the optimal vaccination, we use a forward-backward sweep numerical method. The model solves the SEIRS Equations (36 to 39) using fourth order Runge-Kutta (RK4) method. The values of all the accompanying parameters necessary for the successful numerical solution of the governing equations are given in Table 2.

4.2. Numerical Results and Discussion

In Figure 3, no optimal control was applied. The number of populations in S is dramatically decreased. Whereas the number of Exposed populations exhibits an increase with time (t). After some time, the population shows an increase in E-compartment. As time progresses, the number of population in I compartment increases, while the number of population in R-decrease. This suggests an endemic situation.

Unlike the previous case, Figure 4 displays the effect of optimal control (vaccination) on the susceptible population. As a consequence of vaccination, there is a movement of the vaccinated population from the S compartment to the recovery compartment. This is confirmed by an increase in the recovery population in the graph. We hasten to comment that irrespective of any control strategy, the population in the S compartment displays a decline as time progresses, but the rate of decrease when vaccination is applied is slower. With the population in the S-compartment vaccinated, recovery is enhanced. This shows that LSD disease can be controlled with vaccination. In addition, Figure 4 displays the response of the population in the S-compartment, for optimal control and non optimal control applications. For both cases, the size of the population decreases at different rates. The profiles for both cases clearly demonstrate the effect of the application of this treatment strategy for the susceptible population.

Table 2. Assigned values for model parameters.

Figure 3. SEIRS without Optimal control.

Figure 4. S-compartment with optimal control and without optimal control.

Figure 5 shows the effect of vaccination as a control measure on the infected population. The sharp decline in the infected cattle population demonstrates that if vaccination is properly applied it sure will turn out as a utilitarian strategy to control LSD spread in cattle, The same can also be said of the “E” or the exposed population (see Figure 6). There is less population in the exposed compartment as a result of the application of vaccination. Both cases demonstrate the beneficial effects of vaccination with time.

Figure 5. I-compartment with optimal control and without optimal control.

Figure 6. E-compartment with optimal control and without optimal control.

Figure 7 shows a dramatic rise in the “R” population resulting from the application of vaccination. This reflects an inflow from virtually all the compartments of the model into the “R” population when control is applied. This again shows that vaccination as intervention strategy has a great impact on controlling LSD among cattle.

For demonstration purposes, we present the effects of different vaccination rates $\left(0\le \nu \le 0.85\right)$ on the SEIR profiles in order to assess the impact on the overall control strategy adopted for this study.

Figures 8-13 show a continuous decrease in the “S” population for increasing vaccination rates. This is solely because vaccination plays a big role in decreasing the susceptible population for any successful control strategy. As the susceptible population decreases, they move to both the recovered (“R” population) and the exposed (the “E” i.e. the exposed or the latent population). There is a gradual rise in the “E” population as indicated in Figure 8. Initially it exhibits a higher profile than the “R” population (see Figure 9 and Figure 10). Figures 11-13 show that it finally catches up with and surpasses the exposed population. The actual point of transition happens around $\nu =0.6$ . Given the parameters adopted in this study this can be considered as an optimal point of application and the decision as to whether to continue increasing the vaccination rate will depend on other consideration among which will be that availability of funds. From Figures 8-13, it can also be observed that the profile of the infected population (the “I” compartment) remains the lowest when compared to the dynamic profiles exhibited by other compartments.

Figure 7. R-compartment with OC and without OC.

Figure 8. SEIR-compartment with optimal control at u [0, 0.2].

Figure 9. SEIR-compartment with optimal control at u [0, 0.3].

Figure 10. SEIR-compartment with optimal control at u [0, 0.5].

Figure 11. SEIR-compartment with optimal control at u [0, 0.6].

Figure 12. SEIR-Compartment with Optimal Control at u [0, 0.75].

Figure 13. SEIR-Compartment with Optimal Control at u [0, 0.85].

5. Conclusions

The work reported herein demonstrates the application of vaccination as an optimal control strategy for LSD cattle disease in Ethiopia in addition to the optimization of the cost of the protective strategy adopted. To implement this, we used the Pontryagins maximum principle to characterize the controls and derive the optimality system. For the numerical simulation, we adopted an algorithm based on the Runge-Kutta forward-backward sweep method. The state equations were solved forward in time, while the adjoint equations were simultaneously solved backward in time.

Both the controlled and uncontrolled scenarios were demonstrated graphically in order to identify the optimal vaccination schedule for the treatment and assessment of the impact of LSD on the cattle population. The results of our numerical experiments (Figures 8-13) clearly demonstrate that the applied optimality measures and intervention protocol have a very desirable effect in controlling and reducing the number of infected cattle population during the control period as long as they are properly managed. In addition, as the control period increases the number of cattle population saved from LSD gradually increased.

Based on our LSD model and parameters adopted for this study, the effects of the optimal vaccination strategy can be assessed for the vaccination rates in the range of $\left(0\le \nu \le 0.85\right)$ . It is observed that the control measures adopted herein achieves optimal operation at approximately $\nu =0.6$ .

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

*Advances in Infectious Diseases*,

**9**, 8-24. doi: 10.4236/aid.2019.91002.

[1] | Abera, Z., Degefu, H., Gari, G. and Ayana, Z. (2015) Review on Epidemiology and Economic Importance of Lumpy Skin Disease. International Journal of Basic and Applied Virology, 4, 8-21. |

[2] | Birhanu, H. and Gezahign, A. (2015) Epidemiology, Economic Importance and Control Techniques of Lumpy Skin Disease: A Review. International Journal of Agricultural Research and Review, 3, 197-205. |

[3] | Regbe, H. (2012) Lumpy Skin Disease (LSD): Outbreak Investigation, Isolation and Molecular Detection of LSD in Selected Areas of Eastern Shewa, Ethiopia. Master Thesis, The School of Graduate Study of Addis Ababa University. |

[4] | Alaa, A., Hatem, M. and Khaled, A. (2008) Polymerase Chain Reaction for Rapid Diagnosis of a Recent Lumpy Skin Disease Virus Incursion to Egypt. Arab Journal of Biotechnology, 11, 293-302. |

[5] | Hawsar, Y.A., Heshu, S.R., Hiewa, O.D. and Hemn, H.O. (2016) Lumpy Skin Disease. Reproductive Immunology, 1, 25. |

[6] |
Birhanu, H., Gezahign, A. and Nuru, S. (2015) Techniques of Lumpy Skin Disease. Animal. https://doi.org/10.11648/j.avs.20150302.15 |

[7] |
Davies, G. (1991) Lumpy Skin Disease, an African Capripox Virus Disease of Cattle. The Veterinary Journal, 147, 489-502.
https://doi.org/10.1016/0007-1935(91)90019-J |

[8] | OIE (2010) OIE Manual of Diagnostic Tests and Vaccines for Terrestrial Animals. Lumpy Skin Disease, 14, 768-778. |

[9] |
Tuppurainen, E.S.M. and Oura, C.A.L. (2012) Review: Lumpy Skin Disease: An Emerging Threat to Europe, the Middle East and Asia. Transboundary and Emerging Diseases, 59, 4048. https://doi.org/10.1111/j.1865-1682.2011.01242.x |

[10] |
Chihota, C., Rennie, L.S., Kitching, R.P. and Mellor, P.S. (2001) Mechanical Transmission of Lumpy Skin Disease Virus by Aedes aegypti (Diptera: Culicidae). Epidemiology & Infection, 126, 317-321. https://doi.org/10.1017/S0950268801005179 |

[11] |
Gari, G., Waret-Szkuta, A., Grosbois, V., Jacquite, P. and Roger, F. (2010) Risk Factors Associated with Observed Clinical Lumpy Skin Disease in Ethiopia. Epidemiology & Infection, 138, 1657-1666. https://doi.org/10.1017/S0950268810000506 |

[12] |
Lenhart, S., Bodine, E., Zhong, P. and Joshi, H. (2013) Illustrating Optimal Control Applications with Discrete and Continuous Features. Fields Institute Communications, 66, 209-238. https://doi.org/10.1007/978-1-4614-5389-5_9 |

[13] | Lenhart, S. and Workman, J.T. (2007) Optimal Control Applied to Biological Models. Mathematical and Computational Biology Series, Chapman and Hall/CRC. |

[14] | Bubniakov, L. (2000) The Mathematics of Infectious Diseases. SIAM, 42, 599-653. |

Copyright © 2020 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.