The Scientific World Journal

Volume 2014 (2014), Article ID 949718, 9 pages

http://dx.doi.org/10.1155/2014/949718

Research Article

Multiple Control Strategies for Prevention of Avian Influenza Pandemic

Roman Ullah,1,2 Gul Zaman,³ and Saeed Islam²

¹Department of Mathematics, Bacha Khan University, Charsadda 24420, Pakistan

²Department of Mathematics, Abdul Wali Khan University Mardan, Mardan 23200, Pakistan

³Department of Mathematics, University of Malakand, Chakdara 23101, Pakistan

Received 28 August 2013; Accepted 31 December 2013; Published 24 February 2014

Academic Editors: P. Crooke and G. Fedele

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

Abstract.

We present the prevention of avian influenza pandemic by adjusting multiple control functions in the human-to-human transmittable avian influenza model. First we show the existence of the optimal control problem; then by using both analytical and numerical techniques, we investigate the cost-effective control effects for the prevention of transmission of disease. To do this, we use three control functions, the effort to reduce the number of contacts with human infected with mutant avian influenza, the antiviral treatment of infected individuals, and the effort to reduce the number of infected birds. We completely characterized the optimal control and compute numerical solution of the optimality system by using an iterative method.

1. Introduction

Highly pathogenic avian influenza is a zoonotic disease caused by H5N1 virus and it has a devastating impact on poultry causing above 90 percent mortality within 48 hours of infection [1]. The disease also spreads to humans through the direct contact with the infected poultry. Despite all advances of medical, avian influenza still poses a significant threat to health and economy of the society. The laboratory confirmed that the fatality rate in human case is about 60 percent [2].

Several reports have so far been made of possible coinfection of humans with an H5N1 strain and a human strain [3, 4]. One of the coinfection reports is of a Turkish family cluster. In late December 2005 and early January 2006, a cluster of 8 confirmed highly pathogenic avian influenza (H5N1) cases was detected in Dogubayazit district in eastern Turkey [5]. All confirmed case patients were hospitalized after onset of symptoms. Four of them were died and the other four recovered. Another cluster of 8 cases of HPAI (H5N1) was detected in the village of Kubu Sambilang, Samatra, Indonesia, in April-May 2006 [6]. Only one of them recovered and the remaining seven members died. Though the human-to-human transmission of the disease is rare, its potential to change into an extremely high virulent human-to-human transmittable pandemic strain is the real danger to human health. Due to this threat, some international organizations have developed and began to implement different strategies (biosecurity, isolation, and antiviral treatment) for forestalling the onset of a pandemic. It has been estimated that a future avian influenza pandemic could cause hundreds of billions of dolors due to lost productivity, associated medical policies, and prevention policies.

Mathematical models have been widely used to evaluate the effects of control measures of avian influenza [7–9]. In 2007, Iwami et al. [10] proposed a mathematical model of the spread of avian influenza from bird population to human population. They discussed that to minimize the spread of the disease in human population, we must take the measures for the infected human with bird flu to quarantine when mutant bird flu has already occurred. Lee et al. [11] show that the optimal control strategies that rely on the use of antiviral treatment or isolation strategies can reduce the number of clinical cases; therein they emphasize isolation strategies in mitigation of pandemic of influenza particularly when the access to antiviral resources is limited. Later on Ullah et al. [12] proposed the control strategies where they focus on minimizing the impact of influenza by minimizing the vaccine wanning.

In this paper, we focus on identifying the optimal control strategies for the model developed in [13]; these controls minimize the impact of avian influenza by isolating the infected individuals, the judicious use of drug supply, and eliminating the infected birds. Recommended medicines like oseltamivir or zanamivir must be taken within two days after the appearance of symptoms. The clinically ill and infectious individuals can be isolated by reasonably effective ways to reduce the transmission of disease, like to educate them to cover their sneeze and cough, not to spit openly, and to avoid the closed contacts with others, sanitizing the rooms or equipments occupied by the patients. The diagnosed infected birds should be culled within a few hours and all the infected dead birds and other contaminated objects (faeces, blood, and feathers) must be destroyed properly (burying, incineration) as soon as possible. We can get the most effective results if we use the isolation, treatment of clinically infectious humans, and elimination of infected birds concurrently.

2. Model Frame Work

In this section, the compartmental model that divides the human and birds populations into two different classes presents the optimal control problem for the transmission dynamics of avian influenza. Our main aim is to show that it is possible to implement the time-dependent antiavian influenza control techniques while minimizing the cost of such measures. In order to do this, we introduce three time-dependent control functions , , and . The control represents the effort to reduce the number of contacts with human infected by mutant avian influenza, represents the fraction of clinically infectious cases treated with antiviral per unit of time, and is the effort to reduce the number of infected birds. Note that the controls are fully effective when for , while there is no control if . We divide the total human population at time denoted by into five distinct epidemiological subclasses which are susceptible , exposed , clinically ill and infectious , treated , and recovered and the birds population into three distinct subclasses which are susceptible , exposed , and infected . Taking into account the assumptions above, the dynamics of the control problem is given by

with the initial conditions

The human population is recruited at a constant birth rate , is the rate of immunity loss, is the progression rate from to , represents the treatment rate of human, and are the natural death rates of humans and birds, respectively, and are disease induced death rates in humans and birds, respectively, represents recruitment rate of birds population, shows recovery due to treatment, is progression rate of birds from to , and , , and are effective contact rates between and , and , and and , respectively.

The objective of our control problem is to minimize the number of clinically infectious humans and infected birds and the cost of implementing the control by using possible minimal control variables for . We use the Lebesgue measurable control and define our objective functional as

where and are positive constants to keep a balance in the size of and . The square of the control variables reflects the severity of the side effects of the controls. , , and are positive weight parameters such that , , and . The objective of the optimal control problem is to seek optimal control functions such that

where the control set is defined as

subject to the system (1) with appropriate initial conditions. Pontryagin’s maximum principle [14] is used to solve this optimal control problem and the derivation of the necessary conditions. First we prove the existence of an optimal control for the system (1) and then derive the optimality system.

3. Existence of Control Problem

In this section, we consider the control system (1) with initial conditions at to show the existence of the control problem. Note that, for the bounded Lebesgue measurable controls and nonnegative initial conditions, nonnegative bounded solutions to the state system exist (see [15]). To prove the existence of the solution of system (1), we can write

where

where denotes derivative of with respect to time . Equation (3) is a nonlinear system with a bounded coefficient. We set

Now,

where

Here the constant is independent of the state variables , , , , and . Also we get , where , which shows that is uniformly Lipschitz continuous; hence by definition, the solution of the system exists (see [15]).

Let us go back to the optimal control problem (1)–(3). In order to find an optimal solution, first we should find the Lagrangian and Hamiltonian for the optimal control problem. The Lagrangian is

We seek for the minimal value of the Lagrangian. To do this, we define the Hamiltonian for the control problem, where , , are the adjoint variables:

For the existence of the control system (1), we state and prove the following theorem.

Theorem 1. There exists an optimal control such that

subject to the control system (1) with the initial conditions (2).

Proof. To prove the existence of an optimal control, we use the result in [16–18]. Note that the control and the state variables are non-negative values. In this minimizing problem, the necessary convexity of the objective functional in , , and is satisfied. The set of all the control variables is also convex and closed by definition. The optimal system is bounded which determines the compactness needed for the existence of an optimal control. In addition the integrand in the functional (3), , is convex on the control set . Also we can see that there exist a constant and positive numbers and such that

because the state variables are bounded, which complete the existence of an optimal control.

In order to derive the necessary conditions, we use Pontryagin’s maximum principle [14] as follows. If is an optimal solution of an optimal control problem, then there exists a nontrivial vector function satisfying the following equations:

We now derive the necessary conditions that optimal control functions and corresponding states must satisfy. In the following theorem, we present the adjoint system and control characterization.

Theorem 2. Let , , , , , , , and be optimal state solutions with associated optimal control variables (, , and ) for the optimal control problem (1)–(3). Then there exist adjoint variables , for , satisfying

with transversality conditions

Furthermore the control functions , , and are given by

Proof. To determine the adjoint equations and the transversality conditions, we use the Hamiltonian in (12). The adjoint system results from Pontryagin’s maximum principle [14]:

with , .

To get the characterization of the optimal control given by (18), solving the equations

on the interior of the control set, and setting the property of the control space , we can derive the desired characterization (18).

4. Numerical Results and Discussion

In this section, we present a semi-implicit finite difference method by discretizing the interval at the points , , where represents the time step such that . We define the state and adjoint variables , , , , , , , , , , , , , , , and and the controls , , and in terms of nodal points , , , , , , , , , , , , , , , , , , and . By combination of forward and backward difference approximation, the method developed by [19] to adopt it in our case is given by

To approximate the time derivative of the adjoint variables by the first-ordered backward difference, we use the appropriate scheme as follows:

The algorithm that describes the approximation method for obtaining the optimal control is as follows.

Algorithm 3. Consider the following.

Step 1. Consider , , , , , , , , , and , .

Step 2. For , do

end for.

Step 3. For , write , , , , , , , , , , and end for.

We have plotted susceptible, exposed, infected, treated, and recovered population of humans and susceptible, exposed, and infected population of birds with and without control by considering parameter values given in Table 1, with initial values , , , , , , , and . In each of the given figures the undashed and the dashed lines represent the individuals without and with control, respectively. The weight constant of the objective functional , , , , and is considered for numerical simulation. Figure 1 shows that the population of the susceptible humans increased after control, in Figures 2 and 3 we see that the number of the exposed and infected humans with control decreased more sharply than that without control, Figure 4 shows that per day clinically reported humans decreased after control, Figure 5 shows that the number of recovered human individuals increased after control, Figure 6 shows that the birds population increased after control, and in Figures 7 and 8 we see that the number of the exposed and infected birds decreased after control. Figure 9 shows the graphs of control variables.

Table 1: Parameter values used for numerical simulations.

Notation
Parameters definition
Value

Recruitment rate of humans
1.7/day

Effective contact rate between and
0.3

Effective contact rate between and
0.15

Effective contact rate between and
0.12

Recruitment rate of birds
5/day

Rate of immunity loss
0.005

Natural death rate of humans
0.1/day

Natural death rate of birds
0.3/day

Recovery due to treatment
0.3

Treatment rate
0.75

Disease induced death rate in humans
0.4/day

Disease induced death rate in birds
3.7/day

Progression rate from to
0.45

Progression rate from to
0.85

Figure 1: The plot represents the population of susceptible humans with and without control.

Figure 2: The plot represents the population of exposed humans with and without control.

Figure 3: The plot represents the population of infected humans with and without control.

Figure 4: The plot represents the population of treated humans with and without control.

Figure 5: The plot represents the population of recovered humans with and without control.

Figure 6: The plot represents the population of susceptible birds with and without control.

Figure 7: The plot represents the population of exposed birds with and without control.

Figure 8: The plot represents the population of infected birds with and without control.

Figure 9: The plot represents the control variables , , and .

5. Conclusion

An optimal control problem of the transmission dynamics of the human-avian influenza disease has been presented. We sought to determine three types of control functions associated with isolation and antiviral treatment of the clinically infectious individuals and elimination of the infected birds. Our control plots indicated that the number of exposed, infected, and hospitalized humans and the number of exposed, infected birds decreased in the optimal system. We developed the necessary conditions for the existence of the optimal control by using Pontryagin’s maximum principle. Using the state and adjoint system together with the characterization of the optimal control, we solved the problem numerically by an efficient numerical method based on optimal control with the estimated parameter values based on avian influenza. The results show that the control practices are very effective in reducing the incidence of infectious population.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

Centers for Disease Control and Prevention, “Key Facts about Avian Influenza (Bird Flu) and Avian Influenza A, (H5N1) Virus,” 2006, http://www.cdc.gov/flu/avian/gen-info/facts.htm

M. Martcheva, B. M. Bolker, and R. D. Holt, “Vaccine-induced pathogen strain replacement: what are the mechanisms?” Journal of the Royal Society Interface, vol. 5, no. 18, pp. 3–13, 2008.

Y. Yang, M. E. Halloran, J. D. Sugimoto, and I. M. Longini Jr., “Detecting human-to-human transmission of avian influenza A (H5N1),” Emerging Infectious Diseases, vol. 13, no. 9, pp. 1348–1353, 2007.

M. Zaman, S. Ashraf, N. A. Dreyer, and S. Toovey, “Human infection with avian influenza virus, Pakistan, 2007,” Emerging Infectious Diseases, vol. 17, no. 6, pp. 1056–1059, 2011.

World Health Organization Europe, “Making preparation count: lessons from the avian influenza outbreak in Turkey,” 2006, http://www.euro.who.int/__data/assets/pdf_file/0018/90513/E89139.pdf

D. Butler, “Family tragedy spotlights flu mutations,” Nature, vol. 442, no. 7099, pp. 114–115, 2006.

N. M. Ferguson, D. A. T. Cummings, S. Cauchemez et al., “Strategies for containing an emerging influenza pandemic in southeast Asia,” Nature, vol. 437, no. 7056, pp. 209–214, 2005.

R. F. Grais, J. H. Ellis, and G. E. Glass, “Assessing the impact of airline travel on the geographic spread of pandemic influenza,” European Journal of Epidemiology, vol. 18, no. 11, pp. 1065–1072, 2003.

I. M. Longini Jr., A. Nizam, S. Xu et al., “Containing pandemic influenza at the source,” Science, vol. 309, no. 5737, pp. 1083–1087, 2005.

S. Iwami, Y. Takeuchi, and X. Liu, “Avian-human influenza epidemic model,” Mathematical Biosciences, vol. 207, no. 1, pp. 1–25, 2007.

S. Lee, G. Chowell, and C. Castillo-Chávez, “Optimal control for pandemic influenza: the role of limited antiviral treatment and isolation,” Journal of Theoretical Biology, vol. 265, no. 2, pp. 136–150, 2010.

R. Ullah, G. Zaman, and S. Islam, “Prevention of influenza pandemic by multiple control strategies,” Journal of Applied Mathematics, vol. 2012, Article ID 294275, 14 pages, 2012.

R. Ullah, G. Zaman, S. Islam, and M. R. Muhyuddin, “Global dynamics of the avian-human influenza with horizontal transmission in human population,” Life Science Journal, vol. 9, no. 4, pp. 5747–5753, 2012.

L. S. Pontryagin, V. G. Boltyanskii, R. V. Gamkrelidze, and E. F. Mishchenko, The Mathematical Theory of Optimal Processes, John Wiley & Sons, Hoboken, NJ, USA, 1962.

G. Birkhoff and G. C. Rotta, Ordinary Differential Equations, John Wiley & Sons, New York, NY, USA, 4th edition, 1989.

A. A. Lashari and G. Zaman, “Optimal control of a vector borne disease with horizontal transmission,” Nonlinear Analysis: Real World Applications, vol. 13, no. 1, pp. 203–212, 2012.

D. L. Lukes, Differntial Equations: Classical to Controlled, vol. 162 of Mathematics in Science and Engineering, Academic Press, New York, NY, USA, 1982.

G. Zaman, Y. Han Kang, and I. H. Jung, “Stability analysis and optimal vaccination of an SIR epidemic model,” BioSystems, vol. 93, no. 3, pp. 240–249, 2008.

A. B. Gumel, P. N. Shivakumar, and B. M. Sahai, “A mathematical model for the dynamics of HIV-1 during the typical course of infection,” in Proceedings of the 3rd World Congress of Non-linear Analysts, vol. 47, pp. 2073–2083, 2001.