![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
We investigate a mosquito population suppression model, which includes the release of Wolbachia-infected males causing incomplete cytoplasmic incompatibility (CI). The model consists of two sub-equations by considering the density-dependent birth rate of wild mosquitoes. By assuming the release waiting period T is larger than the sexual lifespan of Wolbachia-infected males, we derive four thresholds: the CI intensity threshold
, the release amount thresholds
and
, and the waiting period threshold
. From a biological view, we assume
throughout the paper. When
, we prove the origin
is locally asymptotically stable iff
, and the model admits a unique T-periodic solution iff
, which is globally asymptotically stable. When
, we show the origin
is globally asymptotically stable iff
, and the model has a unique T-periodic solution iff
, which is globally asymptotically stable. Our theoretical results are confirmed by numerical simulations.
1. Introduction
Mosquitoes can spread diseases between humans, or from animals to humans. They ingest pathogens from an infected host and then transmit diseases during subsequent bites. Common mosquito-borne diseases include dengue, chikungunya, Zika, etc. [Citation37]. The World Health Organization has reported that an estimated 100–400 million people are infected with dengue each year [Citation38]. In recent years, dengue has spread rapidly in more than 129 countries [Citation3]. Traditional ways to prevent the transmission of dengue are integrated mosquito management which includes spraying insecticide or removing breeding sites [Citation1,Citation36]. However, besides environmental pollution, spraying insecticide has only a short-term effect on eradicating mosquitoes due to the emergence of insecticide resistance [Citation32,Citation34]. At the same time, removing all breeding sites of larva mosquitoes is an impossible mission. Although the dengue vaccine Dengvaxia was approved for use in 2015, the safety of the vaccine cannot be guaranteed due to the phenomenon of antibody-dependent enhancement [Citation9,Citation11].
Nowadays, an innovative biological control strategy involves an intracellular bacterium called Wolbachia, which exists in about 20% arthropod hosts [Citation19,Citation35]. Wolbachia was first identified in 1924 by Hertig and Wolbach, and has gained widespread attention of researchers after its role in cytoplasmic incompatibility (CI) was revealed by Laven [Citation21] in 1956. CI occurs when Wolbachia-infected males mate with Wolbachia-uninfected females, resulting in eggs produced from mated females will not hatch [Citation15,Citation29]. Since then, scientists have made great efforts to inject Wolbachia into Aedes mosquitoes to offer possible biological tools to combat dengue and other mosquito-borne diseases. The groundbreaking work was credited to Xi et al. [Citation39], who successfully established the first stable Wolbachia strain, named as WB1, in Aedes aegypti by embryonic microinjection in 2005. Subsequent experiments showed that WB1 induced complete CI. Owing to this breakthrough, various Wolbachia strains have been stably established in Aedes albopictus [Citation40], which is the sole vector of dengue in Guangzhou. In 2015, Xi's team built the largest mosquito factory of the world in Guangzhou, China, and a three-year field trial in Shazai island and Dadaosha island by releasing factory-reared male mosquitoes eradicated more than 95% of Aedes albopictus mosquito populations [Citation48].
As the successful establishment of various Wolbachia strains in mosquitoes, mathematical modelling and analysing on the interactive dynamics of wild and released mosquitoes has attracted plenty of attention [Citation5–7,Citation22,Citation23,Citation30,Citation33,Citation42,Citation44]. Various mathematical models have been established, including discrete models [Citation4,Citation8,Citation45,Citation47,Citation50], which are suited to cage mosquito populations, ordinary differential equation models [Citation13,Citation20,Citation49] and delay differential equation models [Citation41,Citation43], which consider the overlapping generations in mosquito populations and the maturation delay of mosquitoes from egg to reproductive adults, respectively. Also, stochastic models [Citation16,Citation17] have been developed to include the random switch of environments, and reaction diffusion models [Citation18,Citation28,Citation31] have been investigated to include the effect of mosquito dispersal on the interactive dynamics of wild and released mosquitoes.
Among most of the above models, the authors assumed that CI is complete, i.e. no eggs from the incompatible mating will hatch. Complete CI is the most ideal situation when we aim to sterilize wild females through CI by releasing Wolbachia-infected males. However, experiments showed that CI could be incomplete [Citation10,Citation40] due to different hosts or different Wolbachia strains, which could also be influenced by environmental conditions. In this paper, we aim to develop a mosquito suppression model which includes both the release of male mosquitoes and the incomplete CI mechanism. To this end, we introduce the first parameter to estimate the CI intensity. That is,
is the relative reduction in egg hatch caused by incompatible crosses, and the hatch rate of the compatible case is
. Let
and
be the numbers of wild and Wolbachia-infected mosquitoes at time t, respectively. Then the ratio
represents the probability of wild mosquitoes mating with Wolbachia-infected mosquitoes under the assumption of random mating behaviour [Citation8]. Let a be the total number of offspring per wild mosquito per unit time. With the interference of Wolbachia-infected males, the birth rate of wild mosquitoes is reduced from a to
Meanwhile, mosquitoes go through four stages of complete metamorphosis: egg, pupa, larva, and adult, with the first three stages named as aquatic stages. Based on the experimental observations that the crowding effect basically takes place in aquatic stages [Citation12,Citation14], many authors included the density-dependent effect in the growth of the mosquito population [Citation24,Citation26]. Under this consideration, we also assume that the density dependence (denoted by ξ) of mosquitoes development process takes place in aquatic stages. Thus, the birth rate of wild mosquitoes is
Together with the assumption that the death rate of wild mosquitoes is density-independent, we reach the following model,
(1)
(1) where μ is the death rate of wild mosquitoes, and we assume that
.
We should mention here that the modelling idea in our model which treated the released male mosquitoes as a given function was first proposed by Yu [Citation41]. In [Citation41], Yu developed an effective strategy to suppress wild mosquitoes by releasing Wolbachia-infected male mosquitoes, which is governed by a given function
. This strategy reduces the difficulties in theoretical researches by reducing dimensions of systems. Since then, this novel modelling idea has attracted a great attention of many scholars. Yu introduced a parameter
as the sexual activity lifespan of Wolbachia-infected males. The release function
is then determined by the relation of
with the waiting period T between two consecutive releases. We assume that the Wolbachia-infected mosquitoes are released impulsively and periodically at discrete time points
, and the release amount of each batch is maintained at a constant level c. The simplest case of
happens when
, that is, once the released males lose their mating ability, another batch of infected males is released supplementarily. Thus, we have the release function
and Model (Equation1
(1)
(1) ) takes the form
(2)
(2) We provide a complete dynamics of Model (Equation2
(2)
(2) ) in Section 2, which is the basis for the discussion of
and
. We only focus on the case
in this paper and leave the case
in future work. Thus, the release function
becomes
and Model (Equation1
(1)
(1) ) can be rewritten as two sub-equations
(3)
(3)
(4)
(4) where
.
In this study, we formulate a mosquito population suppression model interfered by Wolbachia-infected males causing incomplete CI. In Section 2, the local and global dynamics of wild mosquitoes are studied completely when . In Section 3, by using Poincaré map, we obtain the sufficient and necessary conditions for the stability of the origin
and the existence and uniqueness of a globally asymptotically stable T-periodic solution for the case
, respectively. Finally, some numerical examples are provided to confirm our main results in Section 4.
2. The dynamical behaviour for the case ![](//:0)
![](//:0)
In this section, we focus on a special release strategy, that is, the waiting period T between two consecutive releases is equal to the sexual activity lifespan . The release amount of each batch is controlled at stable level c, and the loss of released mosquitoes during the short sexual activity lifespan can be ignored. Thus, the total amount of released mosquitoes with sexual activity is maintained at the level c. The dynamics of wild mosquitoes
are mainly governed by (Equation2
(2)
(2) ). Thus, to better understand the dynamical behaviours, we need to count the positive equilibria of (Equation2
(2)
(2) ) and analyse their stability.
We claim that Model (Equation2(2)
(2) ) has no positive equilibrium when the release amount
. In fact, from (Equation2
(2)
(2) ) we easy to find that the part in the first bracket is positive and the part in the second bracket is negative when
. Taking together, we derive that the right part of (Equation2
(2)
(2) ) is negative for all
implying that (Equation2
(2)
(2) ) has no positive equilibrium. In the following, we only need to consider the case that
and analyse the existence of the positive equilibria over the rectangular area in
-c plane, as shown in Figure . We rewrite (Equation2
(2)
(2) ) as
(5)
(5) where
The equilibria of Model (Equation5
(5)
(5) ) satisfy
. It is very obvious that Model (Equation5
(5)
(5) ) has at least one trivial equilibrium
. We then search other equilibria from
. Since
is a quadratic equation, its discriminant is
which is a quadratic polynomial with respect to c. It is easy to find that there are two positive real roots
and
such that
, where
and
are given as
The two roots divide
into three intervals:
and
. We compare
with
and find
which implies that Model (Equation5
(5)
(5) ) has no positive equilibrium when the release amount
. When
, the discriminant
holds, then
is positive for all w>0, thus Model (Equation5
(5)
(5) ) also has no positive equilibrium.
Figure 1. The distribution of positive equilibria N in Model (Equation5(5)
(5) ).
![Figure 1. The distribution of positive equilibria N in Model (Equation5(5) {dwdt=−wf(w)w+c,w(0)=u,(5) ).](/cms/asset/1885bc04-777d-4237-a1e5-e8ee892e9111/tjbd_a_2037760_f0001_oc.jpg)
We only need to find the positive equilibria for . When
, we have
and derive a unique equilibrium of (Equation5
(5)
(5) ), that is,
From the expression of
, it is a function of
when
varies over
. This curve divides the rectangle in Figure into two areas. For the case that the release amount
, the equilibrium
is positive only when
, where
(6)
(6) When
, Model (Equation5
(5)
(5) ) has two equilibria, that is,
Now, we need to determine the signs of the two equilibria. Since
is a convex function, it is easy to get that
has a unique positive root when
, or
and
. We consider
and derive another curve
which is tangent to
at
. The two curves divide the rectangle in Figure into six areas. When
, we have
and Model (Equation5
(5)
(5) ) has a unique equilibrium
. When
, Model (Equation5
(5)
(5) ) has a unique positive equilibrium
if
and has no positive equilibrium if
. By further analysing, we find
when
, and Model (Equation5
(5)
(5) ) has two positive equilibria
and
if
and
, where
, and has no positive equilibrium if
. Now, we can give the main results as follows.
Theorem 2.1
For the case , the following conclusions hold.
If
, then Model (Equation5
(5)
(5) ) has two non-negative equilibria: the origin
and the unique positive equilibrium
, where
is unstable and
is asymptotically stable.
If
, then Model (Equation5
(5)
(5) ) has only one non-negative equilibrium
, and
is globally asymptotically stable.
Theorem 2.2
For the case , the following conclusions hold.
If
, then Model (Equation5
(5)
(5) ) has two non-negative equilibria: the origin
and the unique positive equilibrium
, where
is unstable and
is asymptotically stable.
If
, then Model (Equation5
(5)
(5) ) has three non-negative equilibria: the origin
, two positive equilibria
and
, where both
and
are asymptotically stable, and
is unstable.
If
, then Model (Equation5
(5)
(5) ) has two non-negative equilibria: the origin
and the unique positive equilibrium
, where
is asymptotically stable and
is semi-stable, namely, the right-side of
is stable and the left-side is not.
If
, then Model (Equation5
(5)
(5) ) has only one equilibrium
, and
is globally asymptotically stable.
Proof.
The existence of equilibria in each domain can be derived from the above analysis, and we only need to prove the stability. We consider the most complicated case, that is, and
. In this case, Model (Equation5
(5)
(5) ) has three non-negative equilibria
,
and
, hence (Equation5
(5)
(5) ) can be rewritten as
(7)
(7) The interval
is divided into three intervals:
,
and
. From (Equation7
(7)
(7) ), we have
which implies that
and
are asymptotically stable, and
is unstable. The proof of Theorem 2.2(2) is completed. For other cases, the proof is similar and omitted here.
Furthermore, combining Theorems 2.1 and 2.2, we can obtain the CI intensity threshold and the release amount threshold (see Figure ),
(8)
(8) In fact, what we most like to see is that
is globally asymptotically stable, which means that no matter what the initial number of wild mosquitoes is, the population will be eventually eradicated in real life. The above theorems state that wild mosquitoes will tend to be extinct when the release amount of Wolbachia-infected males is larger than
, otherwise, the mosquito population can only be suppressed locally.
On this basis, we continue to study the case . When
and
, the global dynamics behaviour of Model (Equation1
(1)
(1) ) will be quite different from the case
. For biological relevance, we only consider
in the following.
3. The dynamical behaviour for the case ![](//:0)
![](//:0)
Since the release period T is greater than the sexual activity lifespan of Wolbachia-infected males, we easy to find that
is a periodic piecewise function and jumps between two constants c and 0. The release strategy makes the time evolution of
switch between (Equation3
(3)
(3) ) and (Equation4
(4)
(4) ). Then Model (Equation1
(1)
(1) ) has a complex dynamical behaviour.
3.1. The Poincaré map of model (1)
For the case , the dynamical behaviour of (Equation1
(1)
(1) ) mainly depends on the constantly switching between (Equation3
(3)
(3) ) and (Equation4
(4)
(4) ). According to (Equation3
(3)
(3) ) and (Equation4
(4)
(4) ), Model (Equation1
(1)
(1) ) has only one equilibrium
. Let
be a solution of (Equation1
(1)
(1) ) with the initial condition
. Clearly,
is continuous and piecewise differentiable. Initial values satisfy
on
for (Equation3
(3)
(3) ) and
on
for (Equation4
(4)
(4) ). In reality,
and
are continuously differentiable with respect to u and
, respectively, where
can be regarded as a solution of the initial value problem of (Equation4
(4)
(4) ) with
. The solution can be defined in the same way in other intervals. Furthermore,
is a T-periodic solution if
. For convenience, we define two continuously differentiable functions
throughout this paper.
For Equations (Equation3(3)
(3) ) and (Equation4
(4)
(4) ), we introduce two new thresholds which are the release amount threshold
and the waiting period threshold
, that is,
(9)
(9) and
(10)
(10) Since
is strictly increasing with respect to
and
when
, we have
for
.
Next, we study the initial value problem of Equation (Equation3(3)
(3) ) with
on
and Equation (Equation4
(4)
(4) ) with
on
. We first rewrite (Equation3
(3)
(3) ) as
(11)
(11) where
Equation (Equation11
(11)
(11) ) is equivalent to
that is,
(12)
(12) where
Integrating (Equation12
(12)
(12) ) from 0 to
, we obtain
(13)
(13) We rewrite (Equation4
(4)
(4) ) as
Integrating it from
to T, we have
(14)
(14) where
. Here
is implicitly determined by (Equation13
(13)
(13) ) and (Equation14
(14)
(14) ). Following (Equation13
(13)
(13) ) and (Equation14
(14)
(14) ), we have
(15)
(15) since
and
as
. Thus, we have
(16)
(16) According to the Poincaré map of Equations (Equation3
(3)
(3) ) and (Equation4
(4)
(4) ), we define two function sequences
and
by
By mathematical induction, we have
and
(17)
(17) By using a similar discussion as in [Citation44], we obtain the following lemma.
Lemma 3.1
For any given initial value u>0, the following conclusions hold.
If
, then sequences
and
are both strictly increasing.
If
, then
. Furthermore,
is a T-periodic solution.
If
, then sequences
and
are both strictly decreasing.
3.2. Release amount ![](//:0)
![](//:0)
In this case, the release amount c of mosquitoes is controlled between and
. To study the dynamical behaviour, we establish two lemmas before giving the main result.
Lemma 3.2
If and
hold, then (Equation1
(1)
(1) ) has a unique T-periodic solution.
Proof.
First, we show the existence of T-periodic solution. According to (Equation16(16)
(16) ), we have
Thus, there exists a sufficiently small
such that
(18)
(18) According to (Equation1
(1)
(1) ), we know that the solution
is strictly decreasing on
, that is, we have
, which means that
though
is strictly increasing on
. Hence, there must exist a
such that
(19)
(19) By Lemma 3.1, we have
is a T-periodic solution of (Equation1
(1)
(1) ).
Next, the uniqueness of the T-periodic solution is showed by contradiction. Assume that there exists another T-periodic solution in
such that
(20)
(20) Without loss of the generality, we assume that
and
are the minimal and maximal T-periodic solutions on
, respectively. There must exist a
such that
due to the continuous differentiability of
. From the fact that
, there are two cases:
, or
.
For the case , let
(21)
(21) Differentiating
with respect to u gives
Substituting (Equation21
(21)
(21) ) into (Equation13
(13)
(13) ), we obtain
(22)
(22) Taking the derivative with respect to u in (Equation22
(22)
(22) ), we have
that is,
(23)
(23) At the same time, taking the derivative of Equation (Equation14
(14)
(14) ), we have
By (Equation14
(14)
(14) ) again, we get
From (Equation23
(23)
(23) ), we have
(24)
(24) By (Equation14
(14)
(14) ), we obtain
(25)
(25) Substituting (Equation25
(25)
(25) ) into (Equation24
(24)
(24) ), we get
(26)
(26) According to (Equation19
(19)
(19) ), (Equation20
(20)
(20) ) and (Equation26
(26)
(26) ), we have
which is equivalent to
where
(27)
(27) Set
then
and
. Since
, we have
. Thus,
is a concave quadratic polynomial, that is, E<0 and F>0. However,
which is a contradiction.
For the case , without loss of generality, we assume that
. Then, we have
Let k be a small constant such that k>1, mk<1, and
has three roots
where
(28)
(28) According to (Equation26
(26)
(26) ), we have
Equation (Equation28
(28)
(28) ) is equivalent to
where
Set
then
,
and
. Since
is a quadratic polynomial,
is impossible. Similar to the above proof,
is in contradiction to the concavity of
. The proof is completed.
Lemma 3.3
If and
hold, then (Equation1
(1)
(1) ) has a unique T-periodic solution.
Proof.
To prove the existence of T-periodic solution, similar to the proof of Lemma 3.2, we only need to show (Equation18(18)
(18) ) is true. Applying (Equation13
(13)
(13) ) and (Equation14
(14)
(14) ), we have
Since
when
, we obtain
that is,
(29)
(29) Set
Clearly,
if and only if
where
. The derivative of
is
where
due to
, which means that
(30)
(30) Thus,
is strictly increasing for
, that is,
. Therefore, there exists a
such that
. Then Model (Equation1
(1)
(1) ) has at least one T-periodic solution for
with the fact that
for
.
Now, we show the uniqueness of the T-periodic solution. From the existence of T-periodic solution, we choose a such that
If (Equation1
(1)
(1) ) has another T-periodic solution, then there exists a
such that
According to (Equation30
(30)
(30) ) and the fact
, we obtain
Then, we get
Since
is strictly increasing for u>0, we have
, which is a contradiction.
Theorem 3.1
If and
are satisfied, then we have
The origin
of Model (Equation1
(1)
(1) ) is locally asymptotically stable if and only if
.
Model (Equation1
(1)
(1) ) has a unique globally asymptotically stable T-periodic solution if and only if
.
Proof.
Assume that is locally asymptotically stable. Based on (Equation16
(16)
(16) ), to show
, it suffices to prove that
by contradiction. If
, then according to Lemmas 3.2 and 3.3, there exists a
such that
By Lemma 3.1, we know
is strictly increasing for
, that is, we have
, which is a contradiction. Therefore,
is true.
We assume and show
is locally asymptotically stable. Following (Equation15
(15)
(15) ), we have
Select a
such that
for
, where
. For any
and
, we let
and prove that if
, then we have
(31)
(31) Suppose that (Equation31
(31)
(31) ) is not true. If there is a
such that
then there must exist a nonnegative integer p such that either
, or
, or
.
For the case , we obtain
When
, we have
, which is a contradiction. When
, from (Equation14
(14)
(14) ), we get
Clearly, the minimum of
is
for
. Hence,
is contradictory.
For the case , we know that
is the minimum of
for
, a contradiction.
For the case , we have
due to
. Similarly, by (Equation14
(14)
(14) ),
Since
, we get
, which leads to a contradiction. In summary,
is stable.
Then, we show that is locally attractive, namely
Obviously,
is the maximum of
for
,
, hence we only need to prove that
By (Equation16
(16)
(16) ), we have
, namely, there exists a
such that
(32)
(32) By (Equation32
(32)
(32) ) and Lemma 3.1, we know that the sequence
is strictly decreasing for
. Therefore, the limit
exists. Choosing l = 0, otherwise there will exist a T-periodic solution
. In conclusion,
is locally asymptotically stable.
Now, we give the proof of (2). Assume that Model (Equation1(1)
(1) ) has a unique globally asymptotically stable T-periodic solution. If
, then there is a contradiction with the above conclusion. Hence, we have
.
Inversely, we assume that and show that Model (Equation1
(1)
(1) ) has a unique globally asymptotically stable T-periodic solution. Employing Lemmas 3.2 and 3.3, we know that Model (Equation1
(1)
(1) ) has a unique T-periodic solution. Define the T-periodic solution as
Obviously,
(33)
(33) We first show the stability of the T-periodic solution. For any
,
, select
, it suffices to show that
, that is,
(34)
(34) Assume that (Equation34
(34)
(34) ) is not true. Then there must exist a
such that
Without loss of generality, we assume that
there must be a nonnegative integer p such that
, or
, or
. Assume that
. Otherwise, the solutions
and
can reverse to the point
.
For the case , we get
which means that
From Lemma 3.1, we know that sequence
is strictly increasing and
where
. Hence,
is another T-periodic solution of Model (Equation1
(1)
(1) ), which leads to a contradiction.
For the case , we obtain
By (Equation14
(14)
(14) ), we have
which implies that the sequence
is strictly increasing by Lemma 3.1. Thus, there is another T-periodic solution, a contradiction again.
For the case , we know
. When
, integrating (Equation12
(12)
(12) ) from pT to
, we have
where
is the increasing function defined in (Equation21
(21)
(21) ). Since
, we obtain
that is,
(35)
(35) Similarly, integrating (Equation12
(12)
(12) ) from
to
, we obtain
From (Equation35
(35)
(35) ), we have
which means that
Hence, we get
. Clearly, the sequence
is strictly increasing by Lemma 3.1, which is a contradiction with a unique T-periodic solution.
When , integrating (Equation4
(4)
(4) ) from
to
, we obtain
Similarly, integrating (Equation4
(4)
(4) ) from
to pT, we have
Since
it follows that
that is,
. Obviously, this is a contradiction.
In conclusion, (Equation34(34)
(34) ) is true, which implies that the T-periodic solution
is stable. Next, we show
is globally attractive, which only need to prove that
(36)
(36) According to Lemma 3.1, (Equation17
(17)
(17) ) and (Equation33
(33)
(33) ), the sequence
is strictly monotonous, thus,
which leads to (Equation36
(36)
(36) ). Therefore,
is globally asymptotically stable.
3.3. Release amount ![](//:0)
![](//:0)
Before giving the main result, we first give a lemma for the case .
Lemma 3.4
If and
hold, then
for any u>0.
Proof.
The proof here is divided into two cases: and
. For the case
, by Lemma 3.2 and
, we have
, so
for
. Therefore,
is strictly decreasing. Furthermore, from (Equation29
(29)
(29) ), we obtain
that is,
Thus, we have
for
.
For the case , according to (Equation14
(14)
(14) ) and (Equation16
(16)
(16) ), we see that
and
. Next, we show the result by contradiction. Assume that there exists a
such that
. Then, there exists a
such that
and
. From (Equation26
(26)
(26) ), we get
which leads to
where E, F and G are defined in (Equation27
(27)
(27) ). Rewriting the above equation, we have
Thus, we get
Set
Obviously,
is a convex quadratic polynomial. When
, we obtain
and
. Then, we have
for
, which is in contradiction to
.
In conclusion, with the fact for
,
for u>0 is true under the given conditions. The proof, therefore, is completed.
Theorem 3.2
If and
are satisfied, then we have
(1) | The origin | ||||
(2) | Model (Equation1 |
Proof.
(1) Assume that is globally asymptotically stable. Therefore, to prove
, it suffices to show that
. Similar to the proof of Theorem 3.1(1), the conclusion can be proved by contradiction. Assume that
. The stability can be shown in the same way as Theorem 3.1(1). Using Lemma 3.4, we have
for any u>0, which means that the sequence
is strictly decreasing. Hence, we have
, same as the proof of Theorem 3.1(1), which implies that
is globally attractive. Thus, the origin
is globally asymptotically stable.
(2) Assume that the T-periodic solution of Model (Equation1(1)
(1) ) is unique globally asymptotically stable. If
, according to the above result, we know that
is globally asymptotically stable, which contradicts the assumption. Assume that
. By Lemma 3.2, we have a unique T-periodic solution. Similar to the proof of Theorem 3.1(2), the global asymptotic stability of T-periodic solution can be proved. The proof is completed.
Corollary 3.1
If and
are satisfied, then Model (Equation1
(1)
(1) ) has a unique globally asymptotically stable T-periodic solution if and only if the origin
is unstable.
Proof.
Assume that Model (Equation1(1)
(1) ) has a unique globally asymptotically stable T-periodic solution denoted by
with
. We prove that
is unstable by contradiction. If
is stable, by the definition of stability, we know that the solution
does not approach
, which is a contradiction. In turn, assume that
is unstable. Similar to the proof of Theorems 3.1(2) and 3.2(2), we can show that Model (Equation1
(1)
(1) ) has a unique globally asymptotically stable T-periodic solution. Detailed steps are omitted here.
4. Numerical simulations
According to the data in [Citation2,Citation25,Citation27,Citation46] and the method of choosing data in [Citation49], we get the values of parameters as
The most difficult parameter to determine is the density-dependent death rate ξ. Note that ξ is primarily determined by the environmental conditions of the area studied. Normally, when the environment is suitable for mosquitoes breeding, ξ will be very small. Not to change the dynamics significantly, we can take
as a typical example. We provide two numerical simulations to confirm the analytical results of Theorems 3.1 and 3.2.
Example 4.1
Let the parameters a, μ, ξ, and be given as
From (Equation6
(6)
(6) ), (Equation8
(8)
(8) ), (Equation9
(9)
(9) ) and (Equation10
(10)
(10) ), we have
We choose
and
. The release period is set at
. As shown in Figure (a), when the initial amount is small, the wild mosquitoes can be suppressed finally, implying that the origin
is locally asymptotically stable, as stated by Theorem 3.1(1). If the release period T is set at
or
, the amount of wild mosquitoes oscillates periodically around some level, as shown in Figures (b) and (c), which verifies the result given in Theorem 3.1(2). And we can find the amplitude increases with the release period.
Figure 2. The dynamical behaviours of Model (Equation1(1)
(1) ) with
. (a) The origin
is locally asymptotically stable with
. (b) Model (Equation1
(1)
(1) ) has a unique globally asymptotically stable T-periodic solution with
. (c) Model (Equation1
(1)
(1) ) has a unique globally asymptotically stable T-periodic solution with
.
![Figure 2. The dynamical behaviours of Model (Equation1(1) dw(t)dt=a(1−shg(t)w(t)+g(t))(1−ξ(w(t)+g(t)))w(t)−μw(t),(1) ) with g∗<c<c∗. (a) The origin E0 is locally asymptotically stable with T<T∗. (b) Model (Equation1(1) dw(t)dt=a(1−shg(t)w(t)+g(t))(1−ξ(w(t)+g(t)))w(t)−μw(t),(1) ) has a unique globally asymptotically stable T-periodic solution with T=T∗. (c) Model (Equation1(1) dw(t)dt=a(1−shg(t)w(t)+g(t))(1−ξ(w(t)+g(t)))w(t)−μw(t),(1) ) has a unique globally asymptotically stable T-periodic solution with T>T∗.](/cms/asset/34e330f5-8da2-4c54-9c5a-7dc1e19f52fb/tjbd_a_2037760_f0002_oc.jpg)
Example 4.2
Same as the above example, we select the same parameter values a, μ, ξ and , and get four thresholds
,
,
, and
. Choosing
and
, if the release period
or
, then the origin
is globally asymptotically stable as shown in Figures (a) and (b). If
, then the amount of wild mosquitoes oscillates periodically as shown in Figure (c). Next, we choose
and select
,
or
, respectively. The same conclusions are reached (see Figure (d,e,f)).
Figure 3. The dynamical behaviours of Model (Equation1(1)
(1) ) with
. (a) The origin
is globally asymptotically stable with
and
. (b) Model (Equation1
(1)
(1) ) has a unique globally asymptotically stable T-periodic solution with
and
. (c) Model (Equation1
(1)
(1) ) has a unique globally asymptotically stable T-periodic solution with
and
. (d) The origin
is globally asymptotically stable with
and
. (e) Model (Equation1
(1)
(1) ) has a unique globally asymptotically stable T-periodic solution with
and
. (f) Model (Equation1
(1)
(1) ) has a unique globally asymptotically stable T-periodic solution with
and
.
![Figure 3. The dynamical behaviours of Model (Equation1(1) dw(t)dt=a(1−shg(t)w(t)+g(t))(1−ξ(w(t)+g(t)))w(t)−μw(t),(1) ) with c≥c∗. (a) The origin E0 is globally asymptotically stable with c=c∗ and T<T∗. (b) Model (Equation1(1) dw(t)dt=a(1−shg(t)w(t)+g(t))(1−ξ(w(t)+g(t)))w(t)−μw(t),(1) ) has a unique globally asymptotically stable T-periodic solution with c=c∗ and T=T∗. (c) Model (Equation1(1) dw(t)dt=a(1−shg(t)w(t)+g(t))(1−ξ(w(t)+g(t)))w(t)−μw(t),(1) ) has a unique globally asymptotically stable T-periodic solution with c=c∗ and T>T∗. (d) The origin E0 is globally asymptotically stable with c>c∗ and T<T∗. (e) Model (Equation1(1) dw(t)dt=a(1−shg(t)w(t)+g(t))(1−ξ(w(t)+g(t)))w(t)−μw(t),(1) ) has a unique globally asymptotically stable T-periodic solution with c>c∗ and T=T∗. (f) Model (Equation1(1) dw(t)dt=a(1−shg(t)w(t)+g(t))(1−ξ(w(t)+g(t)))w(t)−μw(t),(1) ) has a unique globally asymptotically stable T-periodic solution with c>c∗ and T>T∗.](/cms/asset/abac75b1-c59e-4367-a729-d11d38685559/tjbd_a_2037760_f0003_oc.jpg)
5. Conclusion
Based on the modelling methods of Yu [Citation41] and Li [Citation24], and the incomplete CI mechanism, we considered a population suppression model consisting of two sub-equations switching each other in this study by releasing Wolbachia-infected male mosquitoes. In modelling, we introduced the CI intensity and the sexual lifespan
of Wolbachia-infected mosquitoes, and ignored the independent dynamical equation for Wolbachia-infected mosquitoes.
We assumed that the release amount of Wolbachia-infected mosquitoes was periodic and impulsive at discrete time points
,
, and focused on the case that the waiting period T was greater than the sexual lifespan
. In this paper, we defined four thresholds: the CI intensity threshold
, the release amount thresholds
and
, and the release waiting period threshold
. We obtained a relatively complete analysis for the dynamics of Model (Equation1
(1)
(1) ) when CI intensity
is greater than its threshold
.
Several important lemmas were established and used to prove the main results. For many dynamical systems, it is challenging to obtain the conditions for the existence and uniqueness of periodic solutions. We established necessary and sufficient conditions for both the stability of the origin and the existence and uniqueness of the globally asymptotically stable T-periodic solution in Theorems 3.1 and 3.2, respectively.
In Theorem 3.1, we showed that when the conditions and
hold, Model (Equation1
(1)
(1) ) has either a locally asymptotically stable origin
or a unique globally asymptotically stable T-periodic solution depending on
or
. That is, if we don't release enough Wolbachia-infected mosquitoes making the release amount be less than the threshold
, then whether the mosquito population is successfully suppressed depends on its initial value even though we release at a sufficient frequency
(see Figure (a)). At the same time, when the mosquitoes are released infrequently, namely,
, then the suppression of mosquito population is unsuccessful (see Figure (b,c)). In Theorem 3.2, we obtained that when the conditions
and
hold, the globally asymptotically stable origin
or the unique globally asymptotically stable T-periodic solution of Model (Equation1
(1)
(1) ) are obtained if
or
. In other words, if we release enough Wolbachia-infected mosquitoes (
) and release them frequently (
), then the mosquito population will be suppressed (see Figure (a,b,e,d)). On the contrary, if the frequency of release is not frequent (
), then the mosquito population is persistent (see Figure (c,f)). Combining Theorems 3.1 and 3.2, under the conditions
and
, the necessary and sufficient condition for the existence of the unique globally asymptotically stable T-periodic solution is that the origin
is unstable.
In summary, we studied the dynamical behaviours of the wild mosquito population when and
, and gave complete conclusions and detailed steps. These studies have reference values for guiding mosquito control. The cases of
or
have not been studied, and these investigations are left to future work.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Additional information
Funding
References
- F. Baldacchino, B. Caputo, F. Chandre, A. Drago, A. della Torre, F. Montarsi, and A. Rizzoli, Control methods against invasive Aedes mosquitoes in Europe: A review, Pest Manag. Sci. 71 (2015), pp. 1471–1485.
- S. Boyer, J. Gilles, D. Merancienne, G. Lempérière, and D. Fontenille, Sexual performance of male mosquito Aedes albopictus, Med. Vet. Entomol. 25 (2011), pp. 454–459.
- O.J. Brady, P.W. Gething, S. Bhatt, J.P. Messina, J.S. Brownstein, A.G. Hoen, C.L. Moyes, A.W.Farlow, T.W. Scott, S.I. Hay, and R. Reithinger, Refining the global spatial limits of dengue virus transmission by evidence-based consensus, Plos Neglect. Trop. D 6(8) (2012), pp. 1–15.
- J.S. Brownstein, E. Hett, and S.L. O'Neill, The potential of virulent Wolbachia to modulate disease transmission by insects, J. Invertebr. Pathol. 84 (2003), pp. 24–29.
- L.M. Cai, S. Ai, and J. Li, Dynamics of mosquitoes populations with different strategies for releasing sterile mosquitoes, SIAM J. Appl. Math. 74(6) (2014), pp. 1786–1809.
- L.M. Cai, S. Ai, and G. Fan, Dynamics of delayed mosquitoes populations models with two different strategies of releasing sterile mosquitoes, Math. Biosci. Eng. 15(5) (2018), pp. 1181–1202.
- L.M. Cai, J. Huang, X.Y. Song, and Y.Y. Zhang, Bifurcation analysis of a mosquito population model for proportional releasing sterile mosquitoes, Discrete Contin. Dyn. Syst. B 24(11) (2019), pp. 6279–6295.
- E. Caspari and G.S. Watson, On the evolutionary importance of cytoplasmic sterility in mosquitoes, Evolution 13(4) (1959), pp. 568–570.
- J. Cohen, Dengue may bring out the worst in Zika, Science 355(6332) (2017), p. 1362.
- J.E. Crawford, D.W. Clarke, V. Criswell, M. Desnoyer, D. Cornel, B. Deegan, K. Gong, K.C. Hopkins, P. Howell, J.S. Hyde, J. Livni, C. Behling, R. Benza, W. Chen, K.L. Dobson, C. Eldershaw, D. Greeley, Y. Han, B. Hughes, E. Kakani, J. Karbowski, A. Kitchell, E. Lee, T. Lin, J. Liu, M. Lozano, W.MacDonald, J.W. Mains, M. Metlitz, S.N. Mitchell, D. Moore, J.R. Ohm, K. Parkes, A. Porshnikoff, C. Robuck, M. Sheridan, R. Sobecki, P. Smith, J. Stevenson, J. Sullivan, B. Wasson, A.M. Weakley, M. Wilhelm, J. Won, A. Yasunaga, W.C. Chan, J. Holeman, N. Snoad, L. Upson, T. Zha, S.L. Dobson, F.S. Mulligan, P. Massaro, and B.J. White, Efficient production of male Wolbachia-infected Aedes aegypti mosquitoes enables large-scale suppression of wild populations, Nat. Biotechnol. 38 (2020), pp. 482–492.
- W. Dejnirattisai, P. Supasa, W. Wongwiwat, A. Rouvinski, G. Barba-Spaeth, T. Duangchinda, A.Sakuntabhai, V.-M. Cao-Lormeau, P. Malasit, F.A. Rey, J. Mongkolsapaya, and G.R. Screaton, Dengue virus sero-cross-reactivity drives antibody-dependent enhancement of infection with Zika virus, Nat. Immunol. 17(9) (2016), pp. 1102–1108.
- C. Dye, Intraspecific competition amongst larval Aedes aegypti: Food exploitation or chemical interference, Ecol. Entomol. 7(1) (1982), pp. 39–46.
- J.Z. Farkas and P. Hinow, Structured and unstructured continuous models for Wolbachia infections, Bull. Math. Biol. 72 (2010), pp. 2067–2088.
- R.M. Gleiser, J. Urrutia, and D.E. Gorla, Effects of crowding on populations of Aedes albifasciatus larvae under laboratory conditions, Entomol. Exp. Appl. 95(2) (2000), pp. 135–140.
- A.A. Hoffmann, B.L. Montgomery, J. Popovici, I. Iturbe-Ormaetxe, P.H. Johnson, F. Muzzi, M.Greenfield, M. Durkan, Y.S. Leong, Y. Dong, H. Cook, J. Axford, A.G. Callahan, N. Kenny, C.Omodei, E.A. McGraw, P.A. Ryan, S.A. Ritchie, M. Turelli, and S.L. O'Neill, Successful establishment of Wolbachia in Aedes populations to suppress dengue transmission, Nature 476 (2011), pp. 454–457.
- L. Hu, M. Huang, M. Tang, J. Yu, and B. Zheng, Wolbachia spread dynamics in stochastic environments, Theor. Popul. Biol. 106 (2015), pp. 32–44.
- L. Hu, M. Tang, Z. Wu, Z. Xi, and J. Yu, The threshold infection level for Wolbachia invasion in random environments, J. Differ. Equ. 266 (2019), pp. 4377–4393.
- M.G. Huang, M.X. Tang, and J.S. Yu, Wolbachia infection dynamics by reaction-diffusion equations, Sci. China Math. 58(1) (2015), pp. 77–96.
- A. Jeyaprakash and M.A. Hoy, Long PCR improves Wolbachia DNA amplification: Wsp sequences found in 76% of sixty-three arthropod species, Insect Mol. Biol. 9 (2000), pp. 393–405.
- M.J. Keeling, F.M. Jiggins, and J.M. Read, The invasion and coexistence of competing Wolbachia strains, Heredity 91 (2003), pp. 382–388.
- H. Laven, Cytoplasmic inheritance in Culex, Nature 177 (1956), pp. 141–142.
- J. Li, Simple mathematical models for interacting wild and transgenic mosquito populations, Math. Biosci. 189 (2004), pp. 39–59.
- J. Li, Differential equations models for interacting wild and transgenic mosquito populations, J. Biol. Dyn. 2(3) (2008), pp. 241–258.
- J. Li, New revised simple models for interactive wild and sterile mosquito populations and their dynamics, J. Biol. Dyn. 11(S2) (2017), pp. 316–333.
- Y. Li, F. Kamara, G. Zhou, S. Puthiyakunnon, C. Li, Y. Liu, Y. Zhou, L. Yao, G. Yan, X.-G. Chen, and P. Kittayapong, Urbanization increases Aedes albopictus larval habitats and accelerates mosquito development and survivorship, PLoS Negl. Trop. Dis. 8 (2014), p. e3301.
- G.H. Lin and Y.X. Hui, Stability analysis in a mosquito population suppression model, J. Biol. Dyn. 14(1) (2020), pp. 578–589.
- F.S. Liu, C.S. Yao, P.Q. Lin, and C.Q. Zhou, Studies on life table of the natural population of Aedes albopictus, Acta Sci. Nat. Univ. Sunyatseni 31(8) (1992), pp. 84–93.
- Y.F. Liu, F. Jiao, and L.C. Hu, Modeling mosquito population control by a coupled system, J. Math. Anal. Appl. 506 (2022), p. 125671.
- C.J. McMeniman, R.V. Lane, B.N. Cass, A.W.C. Fong, M. Sidhu, Y.-F. Wang, and S. L. O'Neill, Stable introduction of a life-shortening Wolbachia infection into the mosquito Aedes aegypti, Science 323 (2009), pp. 141–144.
- J.L. Rasgon and T.W. Scott, Impact of population age structure on Wolbachia transgene driver efficacy: Ecologically complex factors and release of genetically modified mosquitoes, Insect Biochem. Mol. Biol. 34 (2004), pp. 707–713.
- P. Schofield, Spatially explicit models of Turelli-Hoffmann Wolbachia invasive wave fronts, J. Theor. Biol. 215 (2002), pp. 121–131.
- P. Somwang, J. Yanola, W. Suwan, C. Walton, N. Lumjuan, L. Prapanthadara, and P. Somboon, Enzymes-based resistant mechanism in pyrethroid resistant and susceptible Aedes aegypti strains from northern Thailand, Parasitol. Res. 109 (2011), pp. 531–537.
- M. Turelli, Evolution of incompatibility-inducing microbes and their hosts, Evolution 48(5) (1994), pp. 1500–1513.
- Y. Wang, X. Liu, C.L. Li, S. Tianyun, J. Jin, Y. Guo, D. Ren, Z. Yang, Q. Liu, and F. Meng, A survey of insecticide resistance in Aedes albopictus (Diptera: Culicidae) during a 2014 dengue fever outbreak in Guangzhou China, J. Econ. Entomol. 110(1) (2017), pp. 239–244.
- J.H. Werren, D. Windsor, and L. Gao, Distribution of Wolbachia among neotropical arthropods, Proc. R. Soc. Lond. B 262 (1995), pp. 197–204.
- WHO, Dengue: Guidelines for Diagnosis, Treatment, Prevention and Control, World Health Organization, Geneva, 2009.
- WHO, Vector-borne diseases, 2 March 2020. Available at https://www.who.int/news-room/fact-sheets/detail/vector-borne-diseases.
- WHO, Dengue and severe dengue, 19 May 2021. Available at https://www.who.int/news-room/fact-sheets/detail/dengue-and-severe-dengue.
- Z.Y. Xi, C.C.H. Khoo, and S.L. Dobson, Wolbachia establishment and invasion in an Aedes aegypti laboratory population, Science 310 (2005), pp. 326–328.
- Z.Y. Xi, C.C.H. Khoo, and S.L. Dobson, Interspecific transfer of Wolbachia into the mosquito disease vector Aedes albopictus, Proc. R. Soc. B 273 (2006), pp. 1317–1322.
- J.S. Yu, Modelling mosquito population suppression based on delay differential equations, SIAM J. Appl. Math. 78(6) (2018), pp. 3168–3187.
- J.S. Yu, Existence and stability of a unique and exact two periodic orbits for an interactive wild and sterile mosquito model, J. Differ. Equ. 269 (2020), pp. 10395–10415.
- J.S. Yu and J. Li, Dynamics of interactive wild and sterile mosquitoes with time delay, J. Biol. Dyn. 13(1) (2019), pp. 606–620.
- J.S. Yu and J. Li, Global asymptotic stability in an interactive wild and sterile mosquito model, J. Differ. Equ. 269 (2020), pp. 6193–6215.
- J.S. Yu and B. Zheng, Modelling Wolbachia infection in mosquito population via discrete dynamical models, J. Differ. Equ. Appl. 25(11) (2019), pp. 1549–1567.
- D. Zhang, X. Zheng, Z. Xi, K. Bourtzis, J.R.L. Gilles, and R. Cordaux, Combining the sterile insect technique with the incompatible insect technique: I-Impact of Wolbachia infection on the fitness of triple-and double-infected strains of Aedes albopictus, PLoS One 10 (2015), p. e0121126.
- B. Zheng and J.S. Yu, Existence and uniqueness of periodic orbits in a discrete model on Wolbachia infection frequency, Adv. Nonlinear Anal. 11 (2022), pp. 212–224.
- X. Zheng, D. Zhang, Y. Li, C. Yang, Y. Wu, X. Liang, Y. Liang, X. Pan, L. Hu, Q. Sun, X. Wang, Y.Wei, J. Zhu, W. Qian, Z. Yan, A.G. Parker, J.R.L. Gilles, K. Bourtzis, J. Bouyer, M. Tang, B. Zheng, J. Yu, J. Liu, J. Zhuang, Z. Hu, M. Zhang, J.-T. Gong, X.-Y. Hong, Z. Zhang, L. Lin, Q. Liu, Z. Hu, Z. Wu, L.A. Baton, A.A. Hoffmann, and Z. Xi, Incompatible and sterile insect techniques combined eliminate mosquitoes, Nature 572 (2019), pp. 56–61.
- B. Zheng, J.S. Yu, and J. Li, Modeling and analysis of the implementation of the Wolbachia incompatible and sterile insect technique for mosquito population suppression, SIAM J. Appl. Math. 81(2) (2021), pp. 718–740.
- B. Zheng, J. Li, and J.S. Yu, One discrete dynamical model on Wolbachia infection frequency in mosquito populations, Sci. China Math. (2022). https://doi.org/10.1007/s11425-021-1891-7.