Study on a Stochastic Avian Influenza Epidemic Model with Generalized Incidence Rate ()
1. Introduction
Avian Influenza, an acute infectious disease caused by influenza A virus, is a complicated disease that can not only infect poultry but also infect humans who have direct exposure to infected birds or contaminated environments. Because humans generally lack immunity to avian influenza virus, the disease has a high mortality rate. For instance, an outbreak of the avian influenza AH7N9 in China caused 134 cases with 45 deaths, from March 31 to August 31, 2013 [1]. According to the data reported by the World Health Organization (WHO), 860 human infections have been reported worldwide since 2003, with about half of those people dying [2]. Hence, the WHO considers the disease to be one of the most significant potential threats to human beings. To better analyze the spread of avian influenza, an increasing amount of research has been studied from the viewpoint of epidemiology and biomathematics. Alexander et al. [3] introduced the pathogenicity of four avian influenza viruses for chickens, turkeys and ducks. Grais et al. [4] and Ferguson et al. [5] took humans as the core and proposed some strategies to retard the spread of avian influenza in the population. Menach [6] regarded poultry farm as a unit to establish a model of infectious disease and focused on the transmission to humans. For more mathematical epidemic models on avian influenza, one can refer to [7] [8] [9] [10] and the references therein.
Recently, under the assumption that avian influenza virus does not spread from person to person and mutate, Liu et al. [1] proposed the following avian influenza bird-to-human transmission model
(1)
where
,
denote the number of the susceptible and infective avian population (i.e. birds) at time t, respectively,
is the transmission rate from infective birds to susceptible birds,
is the natural death rate of the birds,
is the disease-related death rate of the infected birds;
,
and
denote the number of the susceptible human, infective human and recovered human at time t, respectively,
is the transmission rate from the infective birds to the susceptible human,
is the natural death rate of the human population,
is the disease-related death rate of the infected human;
is the recovery rate of the infective human. The function
describes the net growth rate of the avian population, and two cases are studied in the reference [1]. In the present article, we pick the case that the avian population is subject to the logistic growth, then
, where
is the intrinsic growth rate,
is the maximal carrying capacity of the avian population.
The incidence rate of epidemic model plays a quite important role in describing the evolution of infectious disease, and system (1) chose bilinear ones, that is,
and
. However, bilinear incidence rate becomes unrealistic when the number of the infective individual is getting larger [11]. Hence a large amount of research selects nonlinear incidence rates, such as saturated incidence rate of the form
[11], and incidence rate with media coverage effect of the form
[12]. In practice, the incidence function is frequently difficult to obtain because the details of disease transmission vary in different conditions. Therefore choosing generalized incidence rates may allow epidemic models to be more flexible in handling realistic data. For model (1), we replace
with
to denote the generalized incidence rate between susceptible and infective birds, and replace
with
to denote the generalized incidence rate between susceptible human and infective birds, then the model becomes
(2)
where the logistic function
is plugged in, and the recovered population is removed since it has no effect on the dynamics of susceptible population and infective population. Furthermore, the transmission of influenza is disturbed by various noises in the environment, such as the unpredictable contact with infected ones, population mobility and meteorological factors. It is shown that environment fluctuations have a important effect on the development of infectious disease [13] [14] [15]. For instance, Meng et al. [14] showed that a large stochastic disturbance can cause infectious diseases to go to extinction, and Li et al. [15] also found that the average number of infected individuals always with the increase of noise intensity.
These observations imply that stochastic disturbance is conductive to epidemic diseases control. Many researches choose white noise as an appropriate representation of environmental random fluctuations and study the effect of stochastic disturbance on the dynamics of epidemic models. Motivated by the approach in [16], we introduce Gaussian white noise which is directly proportional to
and
, respectively, and finally arrive at the following stochastic model
(3)
where
, are independent standard Brownian motions with
, and
denote the intensities of white noise. Moreover, let (
be a complete probability space with a filtration
satisfying the usual conditions (i.e., it is increasing and right continuous while
contains all
-null sets), then
are defined on this complete probability space. We also introduce the notations:
,
,
.
Throughout the paper, we further assume that
(H1)
and
,
hold for all
;
(H2)
is Lipschitz on
, namely, there exists a constant
, such that
for any
, where
;
(H3)
, where
denote the derivative of the function
at
, respectively.
Our generalized incidence rates can be applied to some specific forms that have been frequently used, such as
(i) Linear type:
;
(ii) Saturated incidence rate:
.
The article is dedicated to investigating the dynamics of the stochastic avian influenza epidemic model (3). The rest of the paper is organized as follows. In the next section, the existence and uniqueness of positive solution is proved for system (3). In Section 3, we discuss the dynamics of the avian-only subsystem, and obtain the sufficient conditions for the extinction of the disease as well as the existence of an ergodic stationary distribution. In Section 4, the dynamics of the avian-human system are discussed, and the asymptotic behavior of system (3) around the unique endemic equilibrium of system (2) is also investigated.
2. Dynamics of the Stochastic Avian-Only Subsystem
Since
and
in system (3) denote the number of individuals, they should be nonnegative from the viewpoint of biology. We first introduce some basic definitions that will be used in the reminder of the article. In general, let X(t) be a homogeneous Markov process in the d-dimension Euclidean space Rddescribed by the stochastic differential equation
(4)
then the diffusion matrix is defined as follows
Furthermore, the differential operator L is defined by
where V(x) is an arbitrary twice continuously differential real-value function.
To begin with, we present the following fundamental theorem which guarantees the existence and uniqueness of positive solution for system (3).
Theorem 2.1. For any initial value
, there is a unique positive solution
of system (3) on
and the solution will remain in
with probability one, namely,
for all
almost surely.
The proof is similar to those of [17] and hence is omitted.
Consider the stochastic avian-only subsystem as follows
(5)
It is obviously observed that the avian system is independent of the human system. Therefore in this section we will focus on the dynamics of the stochastic system (5), and our main goal is to discuss the extinction of disease and the existence of stationary distribution.
The threshold parameter
of the deterministic system with respect to (5) can be computed by application of the next-generation matrix approach as
Motivated by the approach in [7], we introduce Gaussian white noise to
. Let
Theorem 2.2. Let
be a solution of system (5) with any initial value
. If
, and
then
,
namely,
will tend to zero exponentially a.s. and the disease will tend to extinction with probability one.
The proof is similar to those of [17] and hence is omitted.
On the other hand, if we define
,
where
,
, then we will arrive at the following theorem which describes the extinction of the disease.
Theorem 2.3. If
, then the solution
of system (5) satisfies
and the distribution of
converges weakly to the measure which has the density
,
,
where Q is a constant such that
.
Proof. Consider the following auxiliary logistic equation with random perturbation
(6)
set
,
,
, then we can obtain
Therefore,
It is clear that
(7)
where
.
Consequently, the condition of Theorem 1.16 in [18] follows from (7). Thus system (6) has the ergodic property, and the invariant density is given by
,
,
where Q is a constant such that
. From the ergodic theorem it follows that
(8)
Let x(t) be the solution of (6) with the initial value
, then the comparison theorem of stochastic differential equation [19] yields
(9)
On the other hand, according to (9), we get
(10)
Take the superior limit on both sides of (10), and note that
, then
Which shows that
a.s.. Hence for any small
, there exists a
and a set
Such that
and
, for
and
. Hereon,
it follows that the distribution of the process
converges weakly to the measure with the density
. This completes the proof.
We now concentrate on verifying the existence of an ergodic stationary distribution for system (5). The following lemma is fundamental in the paper.
Lemma 2.4. ( [20]) The Markov process
, the solution of system (4), has a unique ergodic stationary distribution
, if there exists a bounded domain
with regular boundary Γ and (B.1) there is a positive number M such that
,
,
, (B.2) there exists a nonnegative
-function V such that LV is negative for any
. Then
for all
, where
is a function integrable with respect to the measure
.
Theorem 2.5. If
and
hold, then for any initial value
, there exists a unique stationary distribution for system (5) and it has the ergodic property.
Proof. We have obtained that for any initial value
, there exists a unique global solution
. In order to prove Theorem 2.5, it suffices to verify the conditions (B.1) and (B.2) in Lemma 2.4. Firstly, we are trying to verify that there exists a neighborhood
and a nonnegative
-function V such that for any
, LV is negative.
Define
(11)
where
(12)
And
It is standard to verify that
has a unique minimum value point
in
.
Define
(13)
where
Let
, due to the hypotheses (H2) and (H3), there exists a
constant
, such that
,
namely,
(14)
According to inequality (14), we can obtain
Moreover,
Thus
Choose a sufficiently small
such that
(15)
(16)
To confirm condition (B.2) of Lemma 2.4, we consider the bounded open set
Define
Obviously,
. Next we will show that
on
.
Case 1. If
, then
, and we have
where
. According to the definition of
and
Inequality (15), we can obtain that
on
.
Case 2. If
, then
, and we have
where
. By virtute of the definition of
and Inequality (15), we can obtain that
on
. Case 3. If
, then we have
by Young inequality.
Hence,
where
,
According to Inequality (16) we can obtain that
on
.
Case 4. If
, similar to the case 3, it is easy to get
According to inequality (16), we obtain that
on
. That is, the condition (B.2) holds.
The diffusion matrix of system (5) is given by
Choose
, we have
Then the condition (B.1) in Lemma 2.4 holds. Consequently, system (5) has a unique stationary distribution and it is ergodic. The proof is completed.
From the biological perspective, stationary distribution can be considered as a weak stability of the system, and the ergodicity tells us that it is persistent in the time average. Moreover, under the assumption
, the condition in the above theorem implies
. It should be mentioned that the long time behavior is remained to be unknown if
. The observation also indicates that it is difficult to obtain a threshold number for system (5), because of the logistic growth rate in the avian population.
3. The Asymptotic Behavior of Stochastic Full System
For any initial value
, system (3) has a unique positive solution
on
, and the solution will remain in
with probability one. In this section, we will investigate the dynamical behavior of the stochastic full system. From analysis of Section 3, we can easily deduce the following result.
Theorem 3.1. If
, then for any initial value
, the disease of system (3) will tend to extinction, almost surely.
The proof is similar to those of theorem 2.3 and hence is omitted.
If the deterministic system (2) has an endemic equilibrium, then it means the disease will persist in the long term. Since stochastic system does not exist endemic equilibrium, it is interesting to investigate the asymptotic behavior of global positive solution of system (3) around endemic equilibrium.
Lemma 3.2. Suppose that the function
is monotone increasing on
, then there exists a unique endemic equilibrium for system (2) if
.
Proof. An endemic equilibrium
of system (2) satisfie
Furthermore,
and
Set
Recall that
, and
, then
if
. Since the function
is monotone increasing on
,
has a unique positive zero point
if
. Hence, system (2) has a unique endemic equilibrium
. The proof is completed.
Theorem 3.3. Suppose the conditions in Lemma 3.2 hold, then system (2) possesses a unique endemic equilibrium
. Let
be the solution of system (3) with any initial value
. If
,
,
and
hold, then
where
Proof. Define
functions as follows
Using the Itô’s formula, we have
(17)
(18)
where the basic inequality
is used. Similarly,
(19)
And
due to the hypotheses (H1)
(20)
Now, define a
function
as follows
Combining (17), (18), (19) and (20), we derive that
where
It then follows from Theorem 6 in [21] that
4. Conclusion
Most systems in the real world are disturbed by various stochastic factors, such as population mobility and meteorological factors including humidity, temperature and precipitation. Hence the effects of environmental fluctuation on the transmission of infectious diseases cannot be neglected. In this paper, we studied a stochastic avian-human influenza epidemic model with logistic growth for birds. To begin with, we proved the existence and uniqueness of global positive solution. For the stochastic avian-only system, set
, then the disease will tend to extinction if
, while system (5) will exist a unique ergodic stationary distribution if
and
. Hence it is interesting to note that a threshold number is difficult to obtain because of the logistic growth rate. For the full stochastic system, the disease will tend to extinction if
. From the viewpoint of biology,
is a proper threshold parameter and the noise intensity in the infected avian population plays a key role. Moreover, we also discussed the asymptotic behavior and proved that the solution of the system (3) oscillates around corresponding endemic equilibrium under some conditions.
Some interesting topics deserve further consideration. For instance, it has not been confirmed that avian influenza virus does not spread from person to person and mutate. Furthermore, the seasonal effect for the transmission of avian influenza is neglected in the present model. We hope to study the comprehensive impacts of seasonal variation and environmental noises in the future.
Acknowledgements
We are very grateful to the editor and anonymous referees for their careful and valuable comments, which greatly improved the presentation of the paper. This work was supported by the National Natural Science Foundations of China (No. 11501364) and the Development Project of Science and Technology of University of Shanghai for Science and Technology (No. 2018KJFZ148 and No. 2019KJFZ171).