Journal of Function Spaces
Volume 2018 (2018), Article ID 5626039, 10 pages
https://doi.org/10.1155/2018/5626039

Research Article

Hopf Bifurcation of a Delayed Ecoepidemic Model with Ratio-Dependent Transmission Rate

Wanjun Xia,¹ Soumen Kundu,² and Sarit Maitra²

¹School of Management Science and Engineering, Anhui University of Finance and Economics, Bengbu 233030, China
²National Institute of Technology Durgapur, Durgapur, West Bengal 713209, India

Correspondence should be addressed to Wanjun Xia; xiawj0982@sina.com

Received 1 July 2018; Accepted 27 August 2018; Published 12 September 2018

Academic Editor: Lishan Liu

Copyright © 2018 Wanjun Xia 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.

A delayed ecoepidemic model with ratio-dependent transmission rate has been proposed in this paper. Effects of the time delay due to the gestation of the predator are the main focus of our work. Sufficient conditions for local stability and existence of a Hopf bifurcation of the model are derived by regarding the time delay as the bifurcation parameter. Furthermore, properties of the Hopf bifurcation are investigated by using the normal form theory and the center manifold theorem. Finally, numerical simulations are carried out in order to validate our obtained theoretical results.

1. Introduction

In recent years, many dynamical models characterizing the propagation of infectious disease [1–3], spread of computer viruses [4–6], and dynamics of some other systems [7–10] are studied by scholars. Ecoepidemiological research deals with the study of the spread of diseases among interacting populations, where the epidemic and demographic aspects are merged within one model. And they have been investigated by many scholars at home and abroad since the pioneer work of Kermack and McKendrick [11], and the interests in investigating the dynamics of ecoepidemic models will be increasing steadily due to its importance from both the mathematical and the ecological points of view.

Many scholars studied different predator-prey models with disease infection in the prey. Chakraborty et al. [12] studied a ratio-dependent ecoepidemic model with prey harvesting and they assumed that both the susceptible and infected prey are subjected to combined harvesting. Upadhyay and Roy [13] proposed an ecoepidemic model with simple law of mass action and modified Holling type II functional response based on the model in [14]. They analyzed stability (linear and nonlinear) of the model. Zhang et al. [15] proposed a three species ecoepidemic model perturbed by white noise and they studied stochastic stability and longtime behavior of the model. Zhou et al. [16] studied local and global stability of a modified Leslie-Gower predator-prey model with prey infection. Some delayed ecoepidemic models with disease infection in the prey have been proposed, and the effect of the delay on the models has been investigated [17–19]. Similarly, some scholars proposed and investigated the ecoepidemic models with disease in predators. Sarwardi et al. [20] and Shaikh et al. [21] studied a Leslie-Gower Holling type II predator-prey model with disease in predator and Leslie-Gower Holling type III predator-prey model with disease in predator, respectively. Some other ecoepidemic models with disease in predators one can refer to include [22–29].

Clearly, most of the epidemic models above are formulated based on the bilinear transmission rate, which is based on the law of mass action. As stated in [30], transmission rate plays an important role in the modelling of epidemic dynamics and the infection probability per contact is likely influenced by the number of infective individuals. Thus, it can be concluded that nonlinear transmission rate seems more reasonable than the bilinear one. To study the effect of a nonlinear incidence rate on the dynamics of an ecoepidemic model, Maji et al. [31] proposed the following ecoepidemic model based the work of Morozov [32]:

where , , and present the densities of the healthy prey, the infected prey, and the predator population, respectively. More parameters are listed in Table 1. They studied stability and persistence of system (1).

Table 1: Parameters and their meanings in this paper.

Parameter
Description

The carrying capacity of the environment

The maximal per capita growth rate of the healthy prey

The transmission rate in the absence of predator

The predator density mediated additional disease transmission rate

The inhibitory effect

The death rate of the infected prey population

The per capita predator consumption rate

The conversion efficiency of the predator

The encounter rate between the predator and the infected prey

As we know, delay differential equations exhibit much more complicated dynamics than ordinary differential equations, and delays can make a dynamical system lose its stability and can induce various oscillations and periodic solutions [17, 23, 26, 33–38]. It is interesting to study the effect of time delay on system (1). To this end, and considering the time required for the gestation of the predator, we incorporate time delay due to the gestation of the predator into system (1) and get the following delayed ecoepidemic system:

subjected to the initial condition:

where is the time delay due to the gestation of the predator.

This paper is organized as follows. Section 2 deals with local stability and existence of the Hopf bifurcation. In Section 3, direction and stability of the Hopf bifurcation are obtained by using center manifold and normal form theory. In Section 4, some numerical simulations are presented in order to verify the analytical findings. Conclusions and discussions are presented in Section 5.

2. Local Stability of the Positive Equilibrium

By direct computation, we can conclude that if , then system (2) has positive equilibrium , where

where is the positive root of (5)

with

and

The Jacobian matrix of system (2) at is

where

Thus, the characteristic equation of about the positive equilibrium is given by

with

When , (10) becomes

where

Based on the Routh-Hurwitz criterion and the discussion in [31], it follows that the positive equilibrium is locally asymptotically stable if the following condition holds: : , and .

For , let be the root of (10); then

Thus,

where

Suppose that

(15) has at least one positive root .

For , from (14)

Differentiating both sides of (10) with respect to yields

Further, we have

where and , .

Obviously, if the condition holds, then . Therefore, based on the Hopf bifurcation theorem in [39], we can obtain the following results.

Theorem 1. Suppose that the conditions - hold for system (2). The positive equilibrium is locally asymptotically stable when and a Hopf bifurcation occurs at the positive equilibrium when .

3. Property of the Hopf Bifurcation

Let ; then is the Hopf bifurcation value of system (2). Rescaling the time delay , then system (2) can be transformed into a functional differential equation in as

where

and

with

and

with

Thus, there exists a matrix function , such that

In view of (21), we choose

where is the Dirac delta function.

For , define

and

Then system (20) is equivalent to

where for .

For , define

and a bilinear inner product

where . Then and are adjoint operators.

Next, we suppose that is the eigenvector of belonging to and is the eigenvector of belonging to . According to the definition of and , we can obtain

From (32), we can get

such that .

Following the method in [39] and using similar computation process in [40], we can get the following coefficients:

with

where and can be determined by the following two equations:

and

Then, we can get the following coefficients which determine the properties of the Hopf bifurcation:

In conclusion, we have the following results.

Theorem 2. For system (2), If (), then the Hopf bifurcation is supercritical (subcritical). If (), then the bifurcating periodic solutions are stable (unstable). If (), then the bifurcating periodic solutions increase (decrease).

4. Numerical Simulation

We choose the same parameters of system (2) as those in [21]: , , , , , , , , , and , while setting as the bifurcation parameter. Then, we get the specific case of system (2) as follows:

from which we can obtain the unique positive equilibrium . Numerically for we have drawn the figure of Lyapunov exponents (Figure 1). Since all the LEs are negative, the system is stable for . Further, we can obtain and the critical value at which a Hopf bifurcation occurs. As is shown in Figure 2, is locally asymptotically stable when . In this case, the three species in system (40) can coexist in an ideal stable state. However, loses its stability and a family of periodic solutions bifurcate from when , which can be illustrated by Figure 3.

Figure 1: Lyapunov exponents for , depicting a stable system.

Figure 2: is locally asymptotically stable when .

Figure 3: loses its stability when .

On the other hand, by some complex calculations, we can obtain and . And further we have , and . Thus, based on the Theorem 2, we can conclude that the Hopf bifurcation is supercritical and the bifurcating periodic solutions are stable and decrease. Since the bifurcating periodic solutions are stable, the three species in system (40) can coexist in an oscillatory mode under some given conditions. This is valuable from the viewpoint of biology.

5. Conclusions

In the present paper, we propose a delayed ecoepidemic model by incorporating the time delay due to the gestation of the predator in the model studied in [31]. Compared with the work in [31], we mainly consider the effect of the time delay on the stability of system (2). The model investigated in our paper is more general since the time required for the gestation of the predator and the results we obtained are suitable complements to the literature [31]. By regarding the time delay due to the gestation of the predator as the bifurcation parameter, sufficient conditions for the local stability of the model and the critical value at which a Hopf bifurcation occurs are derived. It is found that when the value of the time delay is suitablely small, system (2) is locally asymptotically stable. In this case, the densities of the healthy prey, the infected prey, and the predator population will tend to stabilization. Namely, the densities of the three species will be in ideal stable state and the disease spreading among the prey can be controlled. Once the value of the time delay passes through the critical value , system (2) loses stability and a family of periodic solutions bifurcate from the positive equilibrium , which shows that the delay due to the gestation of the predator plays a very complicated role in destabilizing the stability of system (2). In this case, the densities of the three species may coexist in an oscillatory and the disease spreading among the prey will be out of control. In addition, the explicit formulae determining stability and direction of the Hopf bifurcation are derived by using the normal form theory and then center manifold theorem for the further investigation.

It should be pointed out that predator-prey models involving delays and also spatial diffusion are increasingly applied to the study of a variety of situations. Based on this consideration, we will investigate the dynamics of the ecoepidemic model with diffusion based on the delayed model in our present paper in the near future.

Data Availability

All the data can be accessed in our manuscript in the Numerical Simulation.

Conflicts of Interest

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

Acknowledgments

This work was supported by Project of Support Program for Excellent Youth Talent in Colleges and Universities of Anhui Province (No. gxyqZD2018044) and Natural Science Foundation of Anhui Province (Nos. 1608085QF145, 1608085QF151).

References

Q. Liu, M. Sun, and T. Li, “Analysis of an {SIRS} epidemic model with time delay on heterogeneous network,” Advances in Difference Equations, Paper No. 309, 13 pages, 2017.

J. Liu, “Bifurcation of a delayed {SEIS} epidemic model with a changing delitescence and nonlinear incidence rate,” Discrete Dynamics in Nature and Society, Art. ID 2340549, 9 pages, 2017.

X. Chen, J. Cao, J. H. Park, and J. Qiu, “Stability analysis and estimation of domain of attraction for the endemic equilibrium of an {SEIQ} epidemic model,” Nonlinear Dynamics, vol. 87, no. 2, pp. 975–985, 2017.

L.-X. Yang, X. Yang, and Y. Yan Tang, “A bi-virus competing spreading model with generic infection rates,” IEEE Transactions on Network Science and Engineering, no. 99, 2017.

F. Wang, W. Huang, Y. Shen, and C. Wang, “Analysis of SVEIR worm attack model with saturated incidence and partial immunization,” Journal of Communications and Information Networks, vol. 1, no. 4, pp. 105–115, 2016.

L.-X. Yang, X. Yang, and Y. Wu, “The impact of patch forwarding on the prevalence of computer virus: a theoretical assessment approach,” Applied Mathematical Modelling, vol. 43, pp. 110–125, 2017.

Y. Sun and F. W. Meng, “Reachable Set Estimation for a Class of Nonlinear Time-Varying Systems,” Complexity, vol. 2017, Article ID 5876371, 6 pages, 2017.

M. Han, L. Sheng, and X. Zhang, “Bifurcation theory for finitely smooth planar autonomous differential systems,” Journal of Differential Equations, vol. 264, no. 5, pp. 3596–3618, 2018.

X. Hao, H. Wang, L. Liu, and Y. Cui, “Positive solutions for a system of nonlinear fractional nonlocal boundary value problems with parameters and p-Laplacian operator,” Boundary Value Problems, vol. 2017, no. 1, Article ID 182, 2017.

J. Jiang and L. Liu, “Existence of solutions for a sequential fractional differential system with coupled boundary conditions,” Boundary Value Problems, vol. 2016, no. 159, 15 pages, 2016.

W. O. Kermack and A. G. McKendrick, “A contribution to the mathematical theory of epidemics,” Proceedings of the Royal Society A Mathematical, Physical and Engineering Sciences, vol. 115, no. 772, pp. 700–721, 1927.

S. Chakraborty, S. Pal, and N. Bairagi, “Dynamics of a ratio-dependent eco-epidemiological system with prey harvesting,” Nonlinear Analysis: Real World Applications, vol. 11, no. 3, pp. 1862–1877, 2010.

R. K. Upadhyay and P. Roy, “Spread of a disease and its effect on population dynamics in an eco-epidemiological system,” Communications in Nonlinear Science and Numerical Simulation, vol. 19, no. 12, pp. 4170–4184, 2014.

R. K. Upadhyay, N. Bairagi, K. Kundu, and J. Chattopadhyay, “Chaos in eco-epidemiological problem of the Salton SEA and its possible control,” Applied Mathematics and Computation, vol. 196, no. 1, pp. 392–401, 2008.

Q. Zhang, D. Jiang, Z. Liu, and D. O'Regan, “Asymptotic behavior of a three species eco-epidemiological model perturbed by white noise,” Journal of Mathematical Analysis and Applications, vol. 433, no. 1, pp. 121–148, 2016.

X. Zhou, J. Cui, X. Shi, and X. Song, “A modified Leslie-Gower predator-prey model with prey infection,” Applied Mathematics and Computation, vol. 33, no. 1-2, pp. 471–487, 2010.

G.-P. Hu and X.-L. Li, “Stability and Hopf bifurcation for a delayed predator-prey model with disease in the prey,” Chaos, Solitons & Fractals, vol. 45, no. 3, pp. 229–237, 2012.

X. Zhou, X. Shi, and X. Song, “The dynamics of an eco-epidemiological model with distributed delay,” Nonlinear Analysis: Hybrid Systems, vol. 3, no. 4, pp. 685–699, 2009.

Z. Zhang and H. Yang, “Hopf bifurcation control in a delayed predator-prey system with prey infection and modified Leslie-Gower scheme,” Abstract and Applied Analysis, Art. ID 704320, 11 pages, 2013.

S. Sarwardi, M. Haque, and E. Venturino, “Global stability and persistence in LG-Holling type II diseased predator ecosystems,” Journal of Biological Physics, vol. 37, no. 1, pp. 91–106, 2011.

A. A. Shaikh, H. Das, and N. Ali, “Study of {LG}-Holling type {III} predator--prey model with disease in predator,” Applied Mathematics and Computation, vol. 58, no. 1-2, pp. 235–255, 2018.

E. Francomano, F. M. Hilker, M. Paliaga, and E. Venturino, “Separatrix reconstruction to identify tipping points in an eco-epidemiological model,” Applied Mathematics and Computation, vol. 318, pp. 80–91, 2018.

L. Wang, R. Xu, and G. Feng, “Modelling and analysis of an eco-epidemiological model with time delay and stage structure,” Applied Mathematics and Computation, vol. 50, no. 1-2, pp. 175–197, 2016.

F. M. Hilker, M. Paliga, and E. Venturino, “Diseased social predators,” Bulletin of Mathematical Biology, vol. 79, no. 10, pp. 2175–2196, 2017.

M. Haque, “A predator-prey model with disease in the predator species only,” Nonlinear Analysis: Real World Applications, vol. 11, no. 4, pp. 2224–2236, 2010.

J. Zhang, W. Li, and X. Yan, “Hopf bifurcation and stability of periodic solutions in a delayed eco-epidemiological system,” Applied Mathematics and Computation, vol. 198, no. 2, pp. 865–876, 2008.

J. S. Zhang and S. L. Sun, “Analysis of an eco-epidemiological model with epidemic in the predator,” Journal of Biomathematics. Shengwu Shuxue Xuebao, vol. 20, no. 2, pp. 157–164, 2005.

Z. K. Guo, W. L. Li, L. H. Cheng, and Z. Z. Li, “Eco-epidemiological model with epidemic and response function in the predator,” Journal of Lanzhou University. Natural Sciences. Lanzhou Daxue Xuebao. Ziran Kexue Ban, vol. 45, no. 3, pp. 117–121, 2009.

S. Kant and V. Kumar, “Stability analysis of predator-prey system with migrating prey and disease infection in both species,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 42, pp. 509–539, 2017.

R. Xu and Z. Ma, “Stability of a delayed SIRS epidemic model with a nonlinear incidence rate,” Chaos, Solitons & Fractals, vol. 41, no. 5, pp. 2319–2325, 2009.

C. Maji, D. Kesh, and D. Mukherjee, “The Effect of Predator Density Dependent Transmission Rate in an Eco-Epidemic Model,” Differential Equations and Dynamical Systems, vol. 3, pp. 1–15, 2017.

A. Y. Morozov, “Revealing the role of predator-dependent disease transmission in the epidemiology of a wildlife infection: A model study,” Theoretical Ecology, vol. 5, no. 4, pp. 517–532, 2012.

J. Liu, “Dynamical analysis of a delayed predator-prey system with modified Leslie-Gower and Beddington-DeAngelis functional response,” Advances in Difference Equations, 2014:314, 19 pages, 2014.

D. Jana, R. Agrawal, and R. K. Upadhyay, “Top-predator interference and gestation delay as determinants of the dynamics of a realistic model food chain,” Chaos, Solitons & Fractals, vol. 69, pp. 50–63, 2014.

J. G. Ren, X. F. Yang, L. X. Yang, Y. Xu, and F. Yang, “A delayed computer virus propagation model and its dynamics,” Chaos, Solitons & Fractals, vol. 45, no. 1, pp. 74–79, 2012.

X. Y. Meng, H. F. Huo, X. B. Zhang, and H. Xiang, “Stability and Hopf bifurcation in a three-species system with feedback delays,” Nonlinear Dynamics, vol. 64, no. 4, pp. 349–364, 2011.

F. Wei and Q. Fu, “Hopf bifurcation and stability for predator-prey systems with Beddington-DeAngelis type functional response and stage structure for prey incorporating refuge,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 40, no. 1, pp. 126–134, 2016.

X. Meng, H. Huo, and H. Xiang, “Hopf bifurcation in a three-species system with delays,” Applied Mathematics and Computation, vol. 35, no. 1-2, pp. 635–661, 2011.

B. D. Hassard, N. D. Kazarinoff, and Y.-H. Wan, Theory and Applications of Hopf Bifurcation, Cambridge University Press, 1981.

C. Xu, “Delay-Induced Oscillations in a Competitor-Competitor-Mutualist Lotka-Volterra Model,” Complexity, vol. 2017, Article ID 2578043, 12 pages, 2017.