# American Institute of Mathematical Sciences

## Stability analysis and optimal control of production-limiting disease in farm with two vaccines

 Department of Mathematics, City University of Hong Kong, Tat Chee Avenue, Kowloon, Hong Kong, China

* Corresponding author: Yue Liu

Received  September 2020 Revised  December 2020 Published  February 2021

Fund Project: This work was supported by CityU Strategic Research Grants (Project Nos. CityU 11303719 and CityU 11301520)

The transmission of production-limiting disease in farm, such as Neosporosis and Johne's disease, has brought a huge loss worldwide due to reproductive failure. This paper aims to provide a modeling framework for controlling the disease and investigating the spread dynamics of Neospora caninum-infected dairy as a case study. In particular, a dynamic model for production-limiting disease transmission in the farm is proposed. It incorporates the vertical and horizontal transmission routes and two vaccines. The threshold parameter, basic reproduction number $\mathcal{R}_0$, is derived and qualitatively used to explore the stability of the equilibria. Global stability of the disease-free and endemic equilibria is investigated using the comparison theorem or geometric approach. On the case study of Neospora caninum-infected dairy in Switzerland, sensitivity analysis of all involved parameters with respect to the basic reproduction number $\mathcal{R}_0$ has been performed. Through Pontryagin's maximum principle, the optimal control problem is discussed to determine the optimal vaccination coverage rate while minimizing the number of infected individuals and control cost at the same time. Moreover, numerical simulations are performed to support the analytical findings. The present study provides useful information on the understanding of production-limiting disease prevention on a farm.

Citation: Yue Liu, Wing-Cheong Lo. Stability analysis and optimal control of production-limiting disease in farm with two vaccines. Discrete & Continuous Dynamical Systems - B, doi: 10.3934/dcdsb.2021058
##### References:

show all references

##### References:
Flow diagram describing the vaccination model in disease-infected dairy cattle. Animals in cattle are categorized into four compartments: susceptible animals ($S$), infected animals ($I$), susceptible vaccinees ($V$), and infected vaccinees ($W$). The dashed line indicates that the buying or selling rate $\Lambda$ is a variable and depending on other compartments. Here the horizontal infection rate $\rho_h = \zeta\frac{I}{N}$
where the criterion for global asymptotic stability of DFE is satisfied. (b) $\alpha = 0.34$, $\beta_{I1} = 0.05$, $\rho_{v2} = 0.5$, $\delta = 0.2$, $r = 0.2$, $p_1 = \phi_1 = \phi_2 = 0.2$, $p_2 = 0.1$ and other parameters not specified are used as given in Table 1. With these parameter values, the criterion for global stability of EE is satisfied">Figure 2.  Solutions with different initial values converge to the (a) DFE and (b) EE. (a) All the parameters, except $\rho_{v1} = 0.5$, are used as given in Table 1 where the criterion for global asymptotic stability of DFE is satisfied. (b) $\alpha = 0.34$, $\beta_{I1} = 0.05$, $\rho_{v2} = 0.5$, $\delta = 0.2$, $r = 0.2$, $p_1 = \phi_1 = \phi_2 = 0.2$, $p_2 = 0.1$ and other parameters not specified are used as given in Table 1. With these parameter values, the criterion for global stability of EE is satisfied
">Figure 3.  Uncertainty analysis of the basic reproduction number (a) $\mathcal{\widehat{R}}_0$ without control and (b) $\mathcal{R}_0$ with vaccination control. Frequency distribution of basic reproduction number $\mathcal{\widehat{R}}_0$ and $\mathcal{R}_0$. All the parameters are used as Table 1
Simulation results of optimal control with low control cost $C_1 = 1$, $C_2 = 5$. (a)–(c) Number of infected individuals ($I$), susceptible vaccinees ($V$), and infected vaccinees ($W$). The dashed and solid lines represent, respectively, the population dynamics with regular control and optimal control. (d) Control profiles of $p_1(t)$ and $p_2(t)$
Simulation results of optimal control with high control cost $C_1 = 20$, $C_2 = 100$. (a)–(c) Number of infected individuals ($I$), susceptible vaccinees ($V$), and infected vaccinees ($W$). The dashed and solid lines represent, respectively, the population dynamics with regular control and optimal control. (d) Control profiles of $p_1(t)$ and $p_2(t)$
Parameters used in the numerical simulations
 Notation Definition Unit Value Reference $\Lambda$ Buying or selling rate $\text{Year}^{-1}$ – – $\alpha$ Pregnant rate $\text{Year}^{-1}$ 0.30 [14,37] $\beta_S$ Abortion rate of susceptible pregnant cow $\text{Unitless}$ 0.02 [14] $\beta_{I1}$ Abortion rate of infected pregnant cow $\text{Unitless}$ 0.08 [29] $\beta_{I2}$ Abortion rate of infected vaccinee $\text{Unitless}$ 0.05 [29] $\rho_{v1}$ Vertical infection rate of infected cow $\text{Unitless}$ 0.95 [13,14] $\rho_{v2}$ Vertical infection rate of infected vaccinee $\text{Unitless}$ 0.40 [13,14] $\zeta$ Prevalence dependent factor $\text{Year}^{-1}$ 0.028 [10,14] $p_1$ Proportion of the vaccine 1 $\text{Year}^{-1}$ 0.50 Assumption $p_2$ Proportion of the vaccine 2 $\text{Year}^{-1}$ 0.50 Assumption $\delta$ Mortality rate $\text{Year}^{-1}$ 0.10 Assumption $\epsilon$ Removal rate $\text{Year}^{-1}$ 0.095 Assumption $\phi_1$ Efficacy of the vaccine 1 Unitless 0.60 Assumption $\phi_2$ Efficacy of the vaccine 2 Unitless 0.60 Assumption $r$ Proportion of cows Unitless 0.50 [23]
 Notation Definition Unit Value Reference $\Lambda$ Buying or selling rate $\text{Year}^{-1}$ – – $\alpha$ Pregnant rate $\text{Year}^{-1}$ 0.30 [14,37] $\beta_S$ Abortion rate of susceptible pregnant cow $\text{Unitless}$ 0.02 [14] $\beta_{I1}$ Abortion rate of infected pregnant cow $\text{Unitless}$ 0.08 [29] $\beta_{I2}$ Abortion rate of infected vaccinee $\text{Unitless}$ 0.05 [29] $\rho_{v1}$ Vertical infection rate of infected cow $\text{Unitless}$ 0.95 [13,14] $\rho_{v2}$ Vertical infection rate of infected vaccinee $\text{Unitless}$ 0.40 [13,14] $\zeta$ Prevalence dependent factor $\text{Year}^{-1}$ 0.028 [10,14] $p_1$ Proportion of the vaccine 1 $\text{Year}^{-1}$ 0.50 Assumption $p_2$ Proportion of the vaccine 2 $\text{Year}^{-1}$ 0.50 Assumption $\delta$ Mortality rate $\text{Year}^{-1}$ 0.10 Assumption $\epsilon$ Removal rate $\text{Year}^{-1}$ 0.095 Assumption $\phi_1$ Efficacy of the vaccine 1 Unitless 0.60 Assumption $\phi_2$ Efficacy of the vaccine 2 Unitless 0.60 Assumption $r$ Proportion of cows Unitless 0.50 [23]
The sensitivity indices of $\mathcal{\widehat{R}}_0$ and $\mathcal{R}_0$ to the parameter $v$
 Parameter ($v$) Value Sensitivity index of $\mathcal{\widehat{R}}_0$ Sensitivity index of $\mathcal{R}_0$ $\alpha$ 0.30 0.8240 0.5225 $\beta_{I1}$ 0.08 -0.0717 -0.0240 $\rho_{v1}$ 0.95 0.8240 0.2751 $\zeta$ 0.028 0.1760 0.0059 $\delta$ 0.10 -0.6893 -0.0311 $\epsilon$ 0.095 -0.3103 -0.2203 $r$ 0.50 0.5137 0.3023 $\beta_{I2}$ 0.05 – -0.0129 $\rho_{v2}$ 0.40 – 0.2475 $p_1$ 0.50 – -0.3219 $p_2$ 0.50 – -0.4267 $\phi_1$ 0.60 – -0.3219 $\phi_2$ 0.60 – -0.4267
 Parameter ($v$) Value Sensitivity index of $\mathcal{\widehat{R}}_0$ Sensitivity index of $\mathcal{R}_0$ $\alpha$ 0.30 0.8240 0.5225 $\beta_{I1}$ 0.08 -0.0717 -0.0240 $\rho_{v1}$ 0.95 0.8240 0.2751 $\zeta$ 0.028 0.1760 0.0059 $\delta$ 0.10 -0.6893 -0.0311 $\epsilon$ 0.095 -0.3103 -0.2203 $r$ 0.50 0.5137 0.3023 $\beta_{I2}$ 0.05 – -0.0129 $\rho_{v2}$ 0.40 – 0.2475 $p_1$ 0.50 – -0.3219 $p_2$ 0.50 – -0.4267 $\phi_1$ 0.60 – -0.3219 $\phi_2$ 0.60 – -0.4267
 [1] Yves Dumont, Frederic Chiroleu. Vector control for the Chikungunya disease. Mathematical Biosciences & Engineering, 2010, 7 (2) : 313-345. doi: 10.3934/mbe.2010.7.313 [2] Bernard Bonnard, Jérémy Rouot. Geometric optimal techniques to control the muscular force response to functional electrical stimulation using a non-isometric force-fatigue model. Journal of Geometric Mechanics, 2021, 13 (1) : 1-23. doi: 10.3934/jgm.2020032 [3] Felix Finster, Jürg Fröhlich, Marco Oppio, Claudio F. Paganini. Causal fermion systems and the ETH approach to quantum theory. Discrete & Continuous Dynamical Systems - S, 2021, 14 (5) : 1717-1746. doi: 10.3934/dcdss.2020451 [4] Miroslav Bulíček, Victoria Patel, Endre Süli, Yasemin Şengül. Existence of large-data global weak solutions to a model of a strain-limiting viscoelastic body. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2021053 [5] Qi Lü, Xu Zhang. A concise introduction to control theory for stochastic partial differential equations. Mathematical Control & Related Fields, 2021  doi: 10.3934/mcrf.2021020 [6] Tobias Geiger, Daniel Wachsmuth, Gerd Wachsmuth. Optimal control of ODEs with state suprema. Mathematical Control & Related Fields, 2021  doi: 10.3934/mcrf.2021012 [7] Rafael Luís, Sandra Mendonça. A note on global stability in the periodic logistic map. Discrete & Continuous Dynamical Systems - B, 2020, 25 (11) : 4211-4220. doi: 10.3934/dcdsb.2020094 [8] Lakmi Niwanthi Wadippuli, Ivan Gudoshnikov, Oleg Makarenkov. Global asymptotic stability of nonconvex sweeping processes. Discrete & Continuous Dynamical Systems - B, 2020, 25 (3) : 1129-1139. doi: 10.3934/dcdsb.2019212 [9] Xu Pan, Liangchen Wang. Boundedness and asymptotic stability in a quasilinear two-species chemotaxis system with nonlinear signal production. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2021064 [10] Diana Keller. Optimal control of a linear stochastic Schrödinger equation. Conference Publications, 2013, 2013 (special) : 437-446. doi: 10.3934/proc.2013.2013.437 [11] Lorenzo Freddi. Optimal control of the transmission rate in compartmental epidemics. Mathematical Control & Related Fields, 2021  doi: 10.3934/mcrf.2021007 [12] Marzia Bisi, Maria Groppi, Giorgio Martalò, Romina Travaglini. Optimal control of leachate recirculation for anaerobic processes in landfills. Discrete & Continuous Dynamical Systems - B, 2021, 26 (6) : 2957-2976. doi: 10.3934/dcdsb.2020215 [13] Simão Correia, Mário Figueira. A generalized complex Ginzburg-Landau equation: Global existence and stability results. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2021056 [14] Sumon Sarkar, Bibhas C. Giri. Optimal lot-sizing policy for a failure prone production system with investment in process quality improvement and lead time variance reduction. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021048 [15] Paula A. González-Parra, Sunmi Lee, Leticia Velázquez, Carlos Castillo-Chavez. A note on the use of optimal control on a discrete time model of influenza dynamics. Mathematical Biosciences & Engineering, 2011, 8 (1) : 183-197. doi: 10.3934/mbe.2011.8.183 [16] Luke Finlay, Vladimir Gaitsgory, Ivan Lebedev. Linear programming solutions of periodic optimization problems: approximation of the optimal control. Journal of Industrial & Management Optimization, 2007, 3 (2) : 399-413. doi: 10.3934/jimo.2007.3.399 [17] Xiaohong Li, Mingxin Sun, Zhaohua Gong, Enmin Feng. Multistage optimal control for microbial fed-batch fermentation process. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021040 [18] John T. Betts, Stephen Campbell, Claire Digirolamo. Examination of solving optimal control problems with delays using GPOPS-Ⅱ. Numerical Algebra, Control & Optimization, 2021, 11 (2) : 283-305. doi: 10.3934/naco.2020026 [19] Livia Betz, Irwin Yousept. Optimal control of elliptic variational inequalities with bounded and unbounded operators. Mathematical Control & Related Fields, 2021  doi: 10.3934/mcrf.2021009 [20] Christian Meyer, Stephan Walther. Optimal control of perfect plasticity part I: Stress tracking. Mathematical Control & Related Fields, 2021  doi: 10.3934/mcrf.2021022

2019 Impact Factor: 1.27

## Tools

Article outline

Figures and Tables