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

Lumpy skin disease (LSD) is an infectious, fatal skin disease of cattle caused by a virus of the family Poxviridae (genus Capripox). In addition, severely affected animals suffer from reduced weight, cessation of milk production and infertility. The aim of this paper is to computationally apply epidemiological (SEIR) and optimal control (OC) techniques to study the transmission and the impact of vaccination on LSD. Based on our numerical experiments, we were able to deduce the overall impact of the optimal strategy adopted for this study on the cattle population for vaccination rates within the range of 0 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 0.6 ν = .


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.

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 ( α ), which represents the birth rate.The infected population are those individ- uals 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 ω i.e. the recovery rate.Cattles can lose their immunity after some time.Let us assume a rate of θ .As a result, that popula- tion moves from R into the S compartment.Hence the following S-E-I-R-S Epidemiological model without control is shown in Figure 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.
d d Figure 1.SEIRS Model without Control.
Table 1.Descriptions for variables for the model.The rate of population exposed to the disease.

γ
The rate which will be infected from the population of exposed 9 ω The rate of getting immunity after sort of treatment in move to the recovered compartment from infected compartment.

θ
The rate of populations loss their immunity and moved to the susceptible compartment.

φ
The rate of natural death Advances in Infectious Diseases Note that the natural death rate on each of the compartment for this paper is the same, which is ( )

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 where 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 ( ) ( ) ( ) Any solution to the above optimal control problem must also satisfy Equation ( 6), the optimal condition, can usually be manipulated to find a re-presentation of u* in terms of t, x, and λ .
An optimal control with bounded control can be written in the form where a, b are fixed real constants and a < b.
For all controls u at each time t, where H is the Hamiltonian, [13].
By an adaptation of the Pontryagin's Maximum Principle (PMP), the OC must satisfy (optimal condition): u then u is the maximization of the problem, while goes with minimization.
The optimal control 0 H u δ δ = obtained without truncation, and is bounded by a and b such that: Runge-Kutta (RK4) routine is applied for the numerical solution.It is required that: ( ) ( ) ∑ be small, in addition to guaranteeing that the relative error between the current and the previous computed scalars be negligible, i.e.where δ is the accepted tolerance?The optimal control problem can now be written in a standard form: ( ) ( ) Subject to:- In order for this to be a maximization problem, we require B > 0. The Hamiltonian given by By applying the optimal condition * 0 2 2 The Runge-Kutta forward sweep yields the transient profile of the dependent variable.The new * u is calculated by setting

Modification of SEIRS Model
A typical optimal control problem requires a performance index or cost function x t X ∈ and a set of control variables ( ) ( ) 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.
( ) In the above equations, vaccination is applied during g birth 1 v , and for pop- ulations in the S-compartment 2 v , thus vaccinated populations moved to re- covered compartment as α and 2 v S .However, for this paper, we applied vaccination on S-compartment only as ( vS ), and the formulation is given as below.
d d With initial conditions, 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: ( ) 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.( ) where For the Hamiltonian function: ( Then by substituting each of the derivatives in Equation (41) into Equation (42) and by applying PMP to the Hamiltonian, we have: ) The Adjoint functions with respect to the dependent variables are: For S: Adjoint function with respect to I ( ) ( ) ( ) Adjoint function with respect to R ( ) ( ) ( ) From the Equation (43 to 47); we obtain: With transversality conditions we have Next, we consider an optimal condition for determining the characterization of the optimal control * u .The optimal condition by differentiating the Hamil- tonian H with respect to u based on the Equation ( 27) for each control * u .We have the following optimal system: taking into account the bound * u and its characte- rization obtain; the optimal control * u can now be put in a compact form: Hence using the characterization of the optimal control, we now have follow- 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.

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.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 ( ) 0 0.85 ν ≤ ≤ 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 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.

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 ( ) It is observed that the control measures adopted herein achieves optimal operation at approximately 0.6 ν = .

2 S 3 E 4 I 5 R 6 α
Susceptible to be expose for infection those of who interact with the disease Exposed population that are infected but have not yet become infectious Infected by the disease who are capable of transmitting the disease to any susceptible individuals.Individuals recovered from infection who are temporarily immune from the infection Birth rate from the total population N 7 β

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

Figure 7 .
Figure 7. R-compartment with OC and without OC.
are optimal for Equations (16, 17), then there exists a piecewise differentiable adjoint variable * u t and ( ) * x t

Table 2 .
Assigned values for model parameters.