Global Dynamics of an Avian Influenza A(H7N9) Epidemic Model with Latent Period and Nonlinear Recovery Rate
An SEIR type of compartmental model with nonlinear incidence and recovery rates was formulated to study the combined impacts of psychological effect and available resources of public health system especially the number of hospital beds on the transmission and control of A(H7N9) virus. Global stability of the disease-free and endemic equilibria is determined by the basic reproduction number as a threshold parameter and is obtained by constructing Lyapunov function and second additive compound matrix. The results obtained reveal that psychological effect and available resources do not change the stability of the steady states but can indeed diminish the peak and the final sizes of the infected. Our studies have practical implications for the transmission and control of A(H7N9) virus.
Avian influenza A(H7N9) is a subtype of influenza viruses that have been detected in birds and confirmed to be low pathogenic among poultry in the past . Human infections by this particular A(H7N9) virus had not previously been reported until it was found in March, 2013 in China (WHO). It appears that A(H7N9) virus has become a highly pathogenic virus for human species who directly or indirectly contacts poultry carrying virus [2, 3]. From September 1, 2016, to April 31, 2017, 643 cases of avian influenza A(H7N9) laboratory-confirmed cases have been reported in Mainland China, including 233 cases that have died (China CDC), which imposes a serious threat to public health.
There are different types of models to analyze the dynamical behavior of avian influenza virus and assess useful control measures. Iwami et al.  showed that when mutant avian influenza had already occurred, reducing the contact rate of susceptible with infectious humans may have a positive effect on preventing the second outbreak. Liu and Fang  formulated a two-host model to investigate the impact of screening and culling of infected poultry. Liu et al.  considered different growth laws of the avian population, to present that the necessary and sufficient condition for periodic solution existing is the Allee effect in avian population. However, most of these models ignore the latent period between infection and symptom onset in human populations, which does exist on the basis of the reported infection cases. Hence, we introduce the incubation period into our model to further study the internal transmission mechanism of A(H7N9) virus.
When a disease breaks out, people’s awareness of its severity can generate a profound psychological impact on the individuals’ behaviors to reduce unnecessary contact with infections . Wang et al.  found that 77% of urban respondents in their investigation reported that they visited live markets less often after influenza A(H7N9) cases were first identified in China in March 2013. Wu et al.  showed that, in the second wave of avian influenza A(H7N9), greater worry among respondents led to changes in protective behaviors such as less visit to live poultry markets and less purchase of live poultry. To model the reduction in contacts due to the psychological effect, various incidence rates were formulated by researchers [10–13]. In this paper, we will modify these functions to investigate the psychological effect on the transmission of A(H7N9) virus.
In previous dynamic models of avian influenza A(H7N9), one usually assumed the recovery rate as a constant, which means that the treatments were always sufficient. But in fact, hospital resources (such as doctors, drugs, hospital bed, and isolation places) are limited to public, especially when a disease breaks out . According to reported cases by CDC, human infections with A(H7N9) virus and common flu virus have similarities in infected time and the early clinical manifestations; therefore, some available hospital resources have already been occupied. Hospital bed-population ratio, the number of available hospital beds per 10,000 population, is widely used by health planners as a method of estimating resource availability to the public . Abdelrazec et al.  established a model for the transmission dynamics of dengue fever, with the recovery rate function reflected by the hospital bed-population ratio and the number of infections. They found that both the oscillations and backward bifurcation occur attributed to limited hospital resource capacity. And that the basic reproduction ratio is not enough to determine whether the disease eliminates or not. In this paper we will introduce the recovery rate function to explore the impact of available hospital resources in spreading A(H7N9) virus.
This paper is organized as follows. In Section 2, we formulate the A(H7N9) transmission model incorporating the combined impact of psychological effect and available hospital resources. The dynamical behavior of the model is investigated in Section 3, including the existence and global stability of equilibria. In Section 4, we carry out numerical simulations to verify the theoretical results, and the conclusions and discussions for further work are presented in Section 5.
2. Model Formulation
Based on information reported, there is no evidence of sustained human-to-human transmission, although there have been two family clusters reported. Thus, we always assume the transmission of A(H7N9) virus is not from person to person. In our model, we divide the poultry into two subclasses: susceptible and infectious , respectively, and the human population is divided into four subclasses: susceptible , latent , infectious , and recovered . Before constructing the model, we make the following assumptions.
(i) Taking into account the factors such as poultry market mobility, environment capacity, and the existing populations, the susceptible poultry is subject to the logistic growth  where and are the intrinsic growth rate and maximal carrying capacity of the poultry.
(ii) Due to psychological effect, the infection force may decrease when the number of infectious individuals increases. Hence, we modify a nonlinear incidence rate proposed by Liu et al.  to describe the transmission of the virus from infected poultry to susceptible individuals, with the following form: where is the transmission coefficient, measures the infection force of the disease, is a nonnegative constant, and measures the inhibition due to the psychological effect.
(iii) We assume that latent humans do not take up the hospital bed resources during the latent period and, meanwhile, consider the impact of hospital resources on the recovery rate, first proposed by Shan and Zhu , which can be expressed in the following formula: with where is the maximum per capita recovery rate due to the sufficient health care resources and few infectious individuals, is the minimum per capita recovery rate due to the basic clinical resources, and is the hospital bed-population ratio which is a nonnegative constant.
Due to the above assumptions, we can formulate the system as follows:
Detailed descriptions of system parameters and their estimated values are listed in Table 1. The variable can be decoupled from the first four equations of system. Hence, we can reduce system (5) to the following system:
For system (6), we first show the following result.
Lemma 1. The set is a positively invariant and attracting region of system (6).
Proof. For system (6) with nonnegative initial conditions, the following holds: hence, the solutions of , , and are nonnegative. Then we get so all solutions of system (6) are nonnegative.
Let , and it follows that which implies that Moreover, if , we have Therefore, each solution of system (6) with nonnegative initial conditions initiating from will remain in for .
3. Analysis of Equilibria
3.1. Existence of Equilibria
Therefore, the coordinates of equilibria are determined by nonnegative solutions of equations (12). Simple calculation yields that system (6) always has two equilibria and , where for all parameter values. We call and disease-free equilibria, which represent the state that there is no infection. Using the method proposed by Diekmann et al.  and van den Driessche and Watmough , the basic reproduction number of system (6), which is the dominant eigenvalue of the next-generation matrix, can be given bywhere is the spectral radius of a matrix.
Next, we discuss the endemic equilibrium denoted by . From a straightforward calculation of the first and second equations of (12), we have Obviously, if , is positive. From the last three equations of (12), the coordinates of point must satisfySubstituting (16) into equation , after some calculations we have the following equation of : where
And , and obviously Intermediate Value Theorem indicates that there exists at least one positive root of (17). In the following, we consider all the situations.
(i) Assuming there are three real roots, Vieta Theorem indicates that From (20) we can see that there are two cases; one is that all three roots are positive and the other is one positive and two negative roots. Equation (21) indicates that at least one root should be negative; hence, in this case, equation (21) has a unique positive root.
(ii) Otherwise, suppose that there are a pair of complex roots and a positive real root, denoted by , , and , where are real numbers. We have Comparing with (17), holds, which contradicts with .
Lemma 2. In system (6) two disease-free equilibria and exist, where if and a unique endemic equilibrium if .
3.2. The Dynamical Behavior of the Poultry-Only Subsystem
In order to better discuss the full system, we first learn the poultry-only subsystem in , Clearly, the poultry-only subsystem (23) is independent of the full system (6). From Lemma 2, we can directly obtain two disease-free equilibria of (23), denoted by and and a unique endemic equilibrium, denoted by if .
Linearizing the subsystem (23) at the equilibria , , and , respectively, we can obtain the Jacobian matrices. For , the characteristic equation always has a positive root . For , the characteristic equation has two negative roots , if . Otherwise it has one positive root. If , exists and the characteristic equation is . All roots of the equation have negative real parts. Hence, we summarize the results as follows.
Lemma 3. The disease-free equilibrium is always unstable. Further, (i) if , the disease-free equilibrium is locally asymptotically stable and (ii) if , the disease-free equilibrium is unstable and the endemic equilibrium exists and is locally asymptotically stable.
The following theorem shows the global stability of the equilibria.
Theorem 4. If , the disease-free equilibrium is globally asymptotically stable in ; if the endemic equilibrium is globally asymptotically stable in .
Proof. If , construct Lyapunov function Calculate the derivative along subsystem (23); it yields The set has a unique point . According to the invariance principle of Lasalle, all solutions of subsystem (23) approach the largest positively invariant subset of the set . Hence, if , is globally asymptotically stable in .
If , consider the Lyapunov function in . Calculate the derivative along subsystem (23); it satisfies The set has a unique point . According to the invariance principle of Lasalle, all solutions of subsystem (23) approach the largest positively invariant subset of the set . Hence, if , is globally asymptotically stable in .
3.3. The Dynamical Behavior of System (6)
In this section, we will discuss the dynamical behavior of system (6) and study the local stability of equilibria , , and . First, we present the following results.
Lemma 5. The disease-free equilibrium is always unstable. Further, (i) if , the disease-free equilibrium is locally asymptotically stable and (ii) if , the disease-free equilibrium is unstable and the endemic equilibrium exists and is locally asymptotically stable.
Proof. (i) The Jacobian matrix at isSince the characteristic equation always has a positive root , is always unstable.
(ii) The Jacobian matrix at is One root of the characteristic equation is , and others are negative roots. Obviously if , disease-free equilibrium is locally asymptotically stable; otherwise is unstable.
(iii) The Jacobian matrix at is The characteristic equation reads where It follows that By the Routh-Hurwitz criterion, the roots of (31) have negative real parts. The next theorem shows the global dynamics of the system.
Theorem 6. If , the disease-free equilibrium is globally asymptotically stable in ; if and , the endemic equilibrium is globally asymptotically stable in .
Proof. If , Theorem 4 indicates that the disease-free equilibrium is globally asymptotically stable in subsystem (23). By calculation, system (6) can be reduced to the following system: From the first two equations of system (34), we can obtain Clearly, we have that further, since that is, is a monotonically decreasing function and . In summary, the disease-free equilibrium is globally asymptotically stable.
If , Theorem 4 indicates that the endemic equilibrium is globally asymptotically stable in subsystem (23). Similarly, we can also simplify system (6) as From the Jacobian matrix of system (38), we can obtain the second additive compound matrix : where We choose the matrix and calculate , which denotes the matrix whose components are , so . Rewrite the matrix in block matrixwhere We consider the norm in as with vector in and denote by the Lozinskil measure with respect to this norm. It follows that where , are matrix norms with respect to the vector norm and denotes the Lozinskil measure with respect to the norm.
We calculate , where Hence, Further, , where and, hence,From the last two equations of system (38), we have Taking into consideration (49), the following holds: and thusWe assume ; henceand then The Bendixson condition is satisfied; then the result follows.
4. Numerical Simulations
In this section, we carry out numerical simulations for system (5) in order to illustrate the influences of the basic reproduction number , psychological effect, and hospital resources on the disease evolution. The lifespan of poultry for chickens is to years under favorable conditions ; thus we assume the poultry can survive years and fix the parameter . People can usually live for years, so the natural death rate of human is . The latent period is about days (China CDC) and . We assume the following parameters: . We choose the initial values as .
Our theoretical results show that the basic reproduction number determines the global dynamics of the system (5). Fix , , and , for ; when equals 1, we obtain . If , the solutions of converge to the disease-free steady state and the disease will finally be extinct (see Figure 1(a)). If , the solutions of converge to the endemic state, which implies that the disease will persist (see Figure 1(b)).
We then use Latin hypercube sampling (LHS)  and partial rank correlation coefficients (PRCCs)  to explore parameter space and find to which parameter the prevalence at endemic equilibrium is sensitive when parameters vary. Due to limited data on the distribution for each parameter, we choose a uniform distribution for all input parameters with the mean value listed in Table 1. PRCC results in Figure 2(a) indicate that the first four parameters with the most significant impact on the equilibrium prevalence are the psychological effect parameter , the hospital bed-population ratio , the minimum recovery rate of human , and maximum recovery rate of human . It is reasonable that the four parameters play important roles in the infections. In fact, a larger psychological effect parameter means that the public improve their awareness of A(H7N9) virus and take more preventive measures, which leads to lower incidence rate and then lower new infections. A larger hospital bed ratio indicates that more sufficient hospital resources and treatments are provided, which then can improve the recovery rate and lead to lower new infections. The results can be seen explicitly from Figure 2(b). When the impact of psychological effect and hospital resources is introduced, the amount of equilibrium prevalence obviously decreases with the parameters and (or) increasing.
To further examine the impact of psychological effect and hospital resources on infections, respectively, we take and with one of parameters and fixed and the other varying. Figure 3(a) shows that slightly increasing parameter can not only diminish the final size of the infected but also result in a much lower peak of the disease. Similar results can be obtained when parameter varies (see Figure 3(b)).
5. Conclusions and Discussions
In this work, in order to evaluate the combined impact of psychological effect and available hospital resources on the transmission of A(H7N9) virus from poultry to humans, we formulated and analyzed a dynamical model with a nonlinear incidence rate and a nonlinear recovery rate. From the mathematical point of view, we obtained the basic reproduction number , which determines the extinction of the avian influenza. Theoretical analysis of system (6) indicates that the disease-free equilibrium is globally asymptotically stable in when the basic reproduction number is less than unity; that is, the avian influenza A(H7N9) will die out (see Figure 1(a)); and the endemic equilibrium is globally asymptotically stable in when the basic reproduction number is larger than unity and . Note that although the global stability of endemic equilibrium is obtained under this specific condition, which may be due to the limitations of the analytical method, numerical simulations show that all solutions can converge to eventually without the specific condition (see Figure 1(b)).
Both the psychological effect and available hospital resources cannot neither change the stability of endemic equilibrium nor alter the basic reproduction number, but they indeed play a significant role in affecting the number of infectious humans, seen from PRCC results (Figure 2(a)) and the impact of parameters and on equilibrium prevalence (Figure 2(b)). Comparing the number of infectious humans with or without psychological effects, that is, parameter or , it can be seen that bigger parameter can significantly decrease the peak of A(H7N9) infections; meanwhile, the final size of the disease can be reduced. However, no matter whether there is psychological effect or not, the disease cannot die out, seen from Figure 3(a). Figure 3(b) indicates that when the available hospital resources are more sufficient, a bigger parameter leads to a smaller size of the outbreak and a lower number of infectious humans. Similarly, the impact of available hospital resources cannot eradicate the disease either.
Different from the previous avian influenza dynamics models, which usually use bilinear and standard incidence rates and constant recovery rate, in this work, incorporating the combined impact of psychological effect and available hospital resources, we formulate A(H7N9) dynamic model with nonlinear incidence rate and nonlinear recovery rate. We introduce the recovery function , where parameter represents hospital bed-population ratio, which reflects the available resources of the health care system to public. The number of hospital beds is a critical index and with the number of infected cases increasing it may become a limiting factor in controlling the spread of A(H7N9) virus. Our results demonstrate that both psychological effect and available hospital resources can dramatically affect the A(H7N9) virus transmission dynamics. This work is an improvement of existing models of the avian influenza A(H7N9) and the results can provide some practical implications for the control of A(H7N9) virus transmission.
Note that, from current data for A(H7N9) infection, there is an incubation period between infection and symptom onset in both avian and human populations . We consider latent class in our model, which is more realistic to exhibit the epidemiology of A(H7N9). Based on this characteristic of A(H7N9) virus, we will incorporate time delay in our model for future study. There have been five seasonal outbreaks of human infection by A(H7N9) virus in China, since the first outbreak was observed in 2013. Except for the first outbreak, others usually started in October, significantly increased in late December, and then peaked in January of the next year . Thus seasonal variation may affect the spread of A(H7N9) virus as one of the important factors. Zhao et al.  presented a model with period parameters to analyze the effect of climate change on the transmission of A(H7N9) and discussed the global stability and threshold conditions. In our future work, we can also consider the incidence rate as a periodic function.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This work is supported by the National Natural Science Foundation of China (NSFC 11401349) and the Foundation for Outstanding Young Scientist in Shandong Province (BS2014SF008).
Study Fingers Chickens, Quail, In Spread Of H7N9 Influenza Virus, 2014, https://www.sciencedaily.com/releases/2014/03/140318093722.htm.
WHO, Dengue and dengue naemorrhagic fever, Fact Sheet 117, 2012.
J. Lucchetti, M. Roy, and M. Martcheva, “An avian influenza model and its fit to human avian influenza cases,” Advances in Disease Epidemiology, pp. 1–30, 2009.View at: Google Scholar