![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
ABSTRACT
It is widely recognized that the criminal act of poaching has brought tremendous damage to biodiversity. This paper employs a stochastic single-species model with regime switching to investigate the impact of poaching. We first carry out the survival analysis and obtain sufficient conditions for the extinction and persistence in mean of the single-species population. Then, we show that the model is positive recurrent by constructing suitable Lyapunov function. Finally, numerical simulations are carried out to support our theoretical results. It is found that: (i) As the intensity of poaching increases, the odds of being at risk of extinction increases for the single-species population. (ii) The regime switching can suppress the extinction of the single-species population. (iii) The white noise is detrimental to the survival of the single-species population. (iv) Increasing the criminal cost of poaching and establishing animal sanctuaries are important ways to protect biodiversity.
1. Introduction
Poaching has led to the reduction and disappearance of many species around the world and has a major impact on the structure and function of ecosystems. For example, the tibetan antelopes are hunted for their extremely soft, light and warm underfur which is usually obtained after death to make luxury items at high prices. However, because rarity makes the living species attractive, their overexploitation can remain profitable, rendering such species even rarer, or driving them to extinction. This phenomenon is called the anthropogenic Allee effect where species rarity increases price, and therefore, the incentive to harvest, driving small populations extinct [Citation9, Citation13]. In addition, Kenney et al. [Citation19] demonstrated that a critical zone exists in which a small, incremental increase in poaching greatly increases the probability of extinction, by using an individual-based, stochastic spatial model based on extensive field data. Hence, to prevent the biological resources from destruction and protect the ecological environment, various measures have been proposed. For instance, many countries have clearly regarded poaching as a crime (for example, in China, poaching may face more than 10 years in prison); at the beginning of the twentieth century, forest and ungulate species were protected and either recovered naturally or by means of planting or reintroduction. Nowadays, large carnivores are strictly protected or managed by quotas in many countries of the world, and great efforts are undertaken to protect and conserve their recovering populations [Citation26]. These efforts led to successful conservation measures for some species, e.g. wolf [Citation6] and Eurasian Lynx [Citation7]. However, other species still face strong population declines, because the illegal harvest (i.e. poaching) of wild resources is an open access system and is driven by economic interests, it is difficult for the law to exert sufficient protection [Citation4, Citation5, Citation8], e.g. lion[Citation3] and tiger [Citation27].
Among the effective measures, the approach of providing protected areas has become most popular over the past decades. The impact of protected areas in renewing biological resources and protecting the population was quantitatively investigated by Kar and Matsuda [Citation17]. They found that protected patches are an effective means of conserving resource populations. Dubey et al. [Citation10] investigated the dynamics of a fishery resource model with migrations between two patches and showed that even if fishery is exploited continuously in the unreserved zone, fish populations can be maintained at an appropriate equilibrium level in the habitat. In addition, as shown in May [Citation24], Gard [Citation11, Citation12], the growth rates in population systems should be stochastic owing to random noises, and as a result, the solution of the system will not tend to a steady positive point but fluctuate around some average values. Hu et al. [Citation14–16] pointed out that fluctuating environmental conditions, temperature and rainfall, in particular, influence the growth of populations strongly: Temperature affects the mortality, while rainfall supplies water in breeding sites in the wild area for the survival. Therefore, it is important to investigate the effect of random noises in the environment on population dynamics. For example, Zou and Wang [Citation35] studied the dynamics of stochastic model for a single species with protection zone and showed that the establishment of nature reserves can protect endangered species effectively in good conditions. In particular, Wei and Wang [Citation31] considered the following stochastic model for a single species with migrations between two different patches:
(1)
(1)
where
,
are the densities of population in the unprotected environment and the protected one at time t, respectively. a is the intrinsic growth rate,
is the carrying capacity of the environment,
is the migration rate to the nature reserve from the unprotected area,
is the migration rate to the unprotected area from the nature reserve,
,
are the intensity of poaching in the unprotected area and the nature reserve, respectively, and assume that illegal hunting takes place within two patches for pursuing high prices with restriction
.
is a standard Brownian motion.
is the intensity of white noise. In [Citation31], it is proved that for any initial value
, there exists a unique global solution to (Equation1
(1)
(1) ) that remains in
almost surely, and then they provided sufficient conditions for the persistence in mean and extinction of the species.
As shown in [Citation29], the usual state in nature is that the population fluctuates around a trend or a stable average. However, many studies have shown that smooth change can be interrupted by a sudden shift to a completely different regime. For example, after a long period in which vegetation cover fluctuated around a gradually declining trend of vegetation cover, the Sahara region collapsed suddenly into a desert in ancient times [Citation28]; an abrupt climate change, whether warming or cooling, wetting or drying, could have lasting and profound impacts on natural ecosystems [Citation1]. Just as Assaf et al. [Citation2] pointed out that the white noise cannot describe the phenomena that the population may suffer sudden catastrophic shocks or alien species invasion in nature. Therefore, it is necessary to introduce telegraph noise, which can be illustrated as a switching between two or more regimes of environment [Citation30]. Frequently, the switching among different environments is memoryless and the waiting time for the next switch is exponentially distributed [Citation20–22, Citation33, Citation34, Citation36]. Thus we can model the regime switching by a continuous-time Markov chain with values in a finite state space, which drives the changes of the parameters of population models with state switchings of the Markov chain. Inspired by the above facts, in the present paper, we model the random by a continuous-time Markov chain in a finite state space
with generator
given by
where
for
with
and
for each
Assume Markov chain
is irreducible and independent of Brownian motion. Therefore, there exists a unique stationary distribution
of
which satisfies
and
. Then we can further obtain the following model with regime switching:
(2)
(2)
where
,
,
,
,
,
, and
are all positive constants, and
for any
.
This article formulates a stochastic single species model with regime switching and the impact of poaching. Our main idea is to investigate the effect of poaching and the telegraph noise on the dynamics of single species. The rest of this article is organized as follows: We first present some preliminaries and basic results of the model in Section 2. Then in Section 3, we present the main results of system (1.2). In Section 4, we present the detailed proof of the theoretical results. In Section 5, some numerical simulations are carried out to illustrate the obtained results. Finally, in Section 6, we give a brief discussion and the summary of the main results.
2. Preliminaries
Throughout this paper, unless otherwise specified, we let be a complete probability space with a filtration
satisfying the usual conditions (i.e. it is right continuous and
contains all P-null sets).
For the convenience of discussion, we let ,
,
,
.
Next, we give some results on the stationary distribution for stochastic differential equations under regime switching. For more details, we can refer the readers to [Citation32]. Consider the following equation:
(3)
(3)
where
is the d-dimensional Brownian motion, and
,
satisfying
. Let
be a solution of the model (Equation3
(3)
(3) ) for any given initial value
. Let
be any twice continuously differentiable function for any given
, the operator L of system (Equation3
(3)
(3) ) can be defined by
In view of the generalized It
's formula [Citation25], one can get that
In the sequel, we shall introduce the definition of the positive recurrent and give some fundamental results.
Definition 2.1
[Citation32]
If for any
, where
and
is the hitting time of E for
i.e. the first time that the process
enters the set E, or
, then process
is said to be positive recurrent with respect to some nonempty bounded open subset E of
, that is,
for any
.
Lemma 2.1
[Citation32]
Model (Equation3(3)
(3) ) is positive recurrent if the following condition is satisfied: there exists a bounded open subset D of
with a regular (i.e. smooth) boundary satisfying that, for each
there exists a nonnegative function
such that
is twice continuously differentiable and that for some
,
Remark 2.1
According to Theorem 4.3 in [Citation32], one can see that the positive recurrent process has a unique stationary distribution.
Lemma 2.2
[Citation18]
Consider the following linear system of equations
(4)
(4)
where
. Suppose that
is irreducible. Then Equation (Equation4
(4)
(4) ) has a solution if and only if
.
3. Main results
To investigate the dynamical behaviour of stochastic model (Equation2(2)
(2) ), the first concerning thing is whether the solution is global and positive.
Theorem 3.1
For any initial value , there is a unique solution
of model (Equation2
(2)
(2) ) on
, and the solution will remain in
with probability 1.
Corollary 3.1
The solution established in Theorem 3.1 satisfies
Moreover, the solution of system (Equation2
(2)
(2) ) has the properties that
(5)
(5)
Next, we will study some asymptotic properties of model (Equation2(2)
(2) ), such as the extinction, persistence of the single-species population in the natural environment and the nature reserve. We establish the following theorem.
Theorem 3.2
Let be a solution of model (Equation2
(2)
(2) ) with any initial value
. Then we have the following results:
,
, i.e. species in the nature reserve and the natural environment are always extinct, if
.
,
, i.e. species in the nature reserve and the natural environment are always persistent in mean, if
and
.
Now, we study the existence of a stationary distribution of the positive solutions to model (Equation2(2)
(2) ). The existence of a stationary distribution can be regarded as the stochastic weak stability of the model. For the sake of convenience and simplicity, we introduce the following notations:
Theorem 3.3
If , then for any initial value
, the solution
of model (Equation2
(2)
(2) ) is positive recurrent with respect to the domain
, where
is a sufficiently small constant. Moreover, model (Equation2
(2)
(2) ) has a unique stationary distribution
.
It is noteworthy to mention that if for
, that is, model (Equation2
(2)
(2) ) without protection zone, which simply becomes the following stochastic model
(6)
(6)
Similarly, we can obtain the following result.
Corollary 3.2
For model (Equation6(6)
(6) ), the following assertions hold.
(i) | For any initial value | ||||
(ii) |
| ||||
(iii) |
| ||||
(iv) | If |
Remark 3.1
The proof of Theorem 3.1 is very similar to Theorem 2.1 in Wei and Wang [Citation31] and is omitted. The proof of Corollary 3.1 is provided in Appendix A.
4. Proofs of main results
In this section, we provide the detailed proofs of the main results illustrated in Section 3.
4.1. Proof of theorem 3.2
Proof.
By the generalized It's formula, one can see that
(7)
(7)
here we have applied the hypotheses
. Integrating both sides from 0 to t and dividing by t, one can see that
(8)
(8)
where
. Note that
is a local martingales, whose quadratic variation is
. It follows from the strong law of large numbers for martingales (see [Citation23]) that
(9)
(9)
Equations (Equation8
(8)
(8) ) and (Equation9
(9)
(9) ) and the ergodic theorem of Markov chain
give
That is, when
, then
, which implies that
This completes the proof of (i).
Now we prove (ii). By using generalized It's formula to system (Equation2
(2)
(2) ) results in
(10)
(10)
Integrating from 0 to t and dividing by t on both sides of (Equation10
(10)
(10) ) yields
(11)
(11)
Equations (Equation5
(5)
(5) ), (Equation9
(9)
(9) ) and (Equation11
(11)
(11) ) and the ergodic theorem of Markov chain
give
Similarly, we can have that
This completes the proof of (ii).
4.2. Proof of theorem 3.3
Proof.
According to Theorem 3.1, we can derive that for any initial value , the solution of system (Equation2
(2)
(2) ) is regular.
Note that , and
then
where
,
. It follows from Lemma 2.2 that the linear equation
(12)
(12)
has a solution
. From (Equation12
(12)
(12) ), it is easy to find that
(13)
(13)
Then, define a
-function
by
(14)
(14)
where
. Using generalized It
's formula to (Equation13
(13)
(13) ) yields
Now define a bounded closed set
by
where
is a sufficiently small constant such that
For convenience, we divide
into four domains
Clearly,
. Next, we will prove
for any
.
Case 1. If , we have
Case 2. In domain
,
Case 3. In domain
,
Case 4. In domain
,
Consequently
It then follows from Lemma 2.1 that the solution of model (Equation2
(2)
(2) ) is positive recurrent with respect to the domain
. Moreover, according to Remark 2.1, we know that model (Equation2
(2)
(2) ) admits a unique stationary distribution.
5. Numerical simulations
In this section, we will numerically illustrate the effects of the Markov chain, poaching and the migration rates ,
on survival of species by several examples. Here, we only focus on the intrinsic growth rate of model (Equation2
(2)
(2) ) is disturbed by a random switching because in reality it is more sensitive to environmental fluctuations than other parameters for the single species population, i.e. we assume that
(15)
(15)
In addition, we always assume that
is a right-continuous Markov chain taking values in a finite state space
with the generator
and the initial value of system (Equation2
(2)
(2) ) is
, except for the other specification. Then the unique stationary distribution of
is given by
. We then divide our simulations into six examples:
Example 5.1
Motivated by Wei et al. [15], the parameter values in model (Equation2(2)
(2) ) are chosen as follows:
(16)
(16)
By computing, we have
,
, and
. It then follows from Theorems 3.2 (ii) and 3.3 that the single species in the nature reserve and the natural environment are always persistent in mean, and the hybrid system (Equation2
(2)
(2) ) has a unique stationary distribution. As seen in Figure , our simulation supports these results. Moreover, the sample mean and variance of
are 0.3360 and 0.0237, and the sample mean and variance of
are 0.8856 and 0.1022, respectively.
Figure 1. (a) The paths of and
for stochastic model (Equation2
(2)
(2) ) with initial values
. (b) The probability density function (PDF) of
and
, respectively. The parameter values are given in (Equation15
(15)
(15) ) and (Equation16
(16)
(16) ).
![Figure 1. (a) The paths of xt and yt for stochastic model (Equation2(2) dxt=xt(ar(t)−br(t)xt)−m1r(t)xt+m2r(t)yt−E1r(t)xtdt+σr(t)xtdBt,dyt=yt(ar(t)−br(t)yt)−m2r(t)yt+m1r(t)xt−E2r(t)ytdt+σr(t)ytdBt,(2) ) with initial values (x0,y0,r(0))=(1,2,1). (b) The probability density function (PDF) of xt and yt, respectively. The parameter values are given in (Equation15(15) b1=b2,m11=m12,m21=m22,E11=E12,E21=E22,σ1=σ2.(15) ) and (Equation16(16) a1=0.9,a2=0.3,b1=0.8,m11=0.4,m21=0.1,E11=0.3,E21=0.01,σ1=0.3.(16) ).](/cms/asset/c3ed3fbf-961c-4265-9bcc-1532c2afff9a/tjbd_a_1914757_f0001_oc.jpg)
As shown in [Citation31], poaching is a major threat to the long-term persistence of the endangered species. In fact, if below some critical population density, someone is always willing to pay enough to offset the costs of finding and hunting the species, hunting further reduces the population, which increases the price of a unit harvest, which in turn increases the hunting effort, and the species is driven into an extinction vortex. Next we focus on the effects of the intensity of poaching on the single species population dynamics for model (Equation2(2)
(2) ). we carry out the following example.
Example 5.2
Based on Example 5.1, we use different vector fields values of and
to investigate the effect of the poaching on the survival of population. We then divide our simulations into two cases:
Case 1. |
| ||||
Case 2. |
|
Figure 2. (a) The paths of and
for stochastic model (Equation2
(2)
(2) ) with initial values
. (b) The probability density function (PDF) of
and
, respectively. Here
,
, other parameter values are given in (Equation15
(15)
(15) ) and (Equation16
(16)
(16) ).
![Figure 2. (a) The paths of xt and yt for stochastic model (Equation2(2) dxt=xt(ar(t)−br(t)xt)−m1r(t)xt+m2r(t)yt−E1r(t)xtdt+σr(t)xtdBt,dyt=yt(ar(t)−br(t)yt)−m2r(t)yt+m1r(t)xt−E2r(t)ytdt+σr(t)ytdBt,(2) ) with initial values (x0,y0,r(0))=(1,2,1). (b) The probability density function (PDF) of xt and yt, respectively. Here E11=E12=0.32, E21=E22=0.3, other parameter values are given in (Equation15(15) b1=b2,m11=m12,m21=m22,E11=E12,E21=E22,σ1=σ2.(15) ) and (Equation16(16) a1=0.9,a2=0.3,b1=0.8,m11=0.4,m21=0.1,E11=0.3,E21=0.01,σ1=0.3.(16) ).](/cms/asset/5e898132-4cb1-4f77-8949-7578d734d4f0/tjbd_a_1914757_f0002_oc.jpg)
Figure 3. (a) The red line, green line and blue line represent the paths of for hybrid system, state 1 and state 2, respectively. (b) The red line, green line and blue line represent the paths of
for hybrid system, state 1 and state 2, respectively. Here the parameters take the same values as in Figure .
![Figure 3. (a) The red line, green line and blue line represent the paths of xt for hybrid system, state 1 and state 2, respectively. (b) The red line, green line and blue line represent the paths of yt for hybrid system, state 1 and state 2, respectively. Here the parameters take the same values as in Figure 2.](/cms/asset/9bc66c8f-5b48-40a3-a352-cdb267aa02c9/tjbd_a_1914757_f0003_oc.jpg)
Figure 4. The evolution of and
of the stochastic model (Equation2
(2)
(2) ) is graphed for vield fields
,
,
. Other parameter values are given in (Equation15
(15)
(15) ) and (Equation16
(16)
(16) ).
![Figure 4. The evolution of xt and yt of the stochastic model (Equation2(2) dxt=xt(ar(t)−br(t)xt)−m1r(t)xt+m2r(t)yt−E1r(t)xtdt+σr(t)xtdBt,dyt=yt(ar(t)−br(t)yt)−m2r(t)yt+m1r(t)xt−E2r(t)ytdt+σr(t)ytdBt,(2) ) is graphed for vield fields (E1,E2)=(0.3,0.01), (E1,E2)=(0.32,0.3), (E1,E2)=(0.82,0.8). Other parameter values are given in (Equation15(15) b1=b2,m11=m12,m21=m22,E11=E12,E21=E22,σ1=σ2.(15) ) and (Equation16(16) a1=0.9,a2=0.3,b1=0.8,m11=0.4,m21=0.1,E11=0.3,E21=0.01,σ1=0.3.(16) ).](/cms/asset/85a13270-4881-4d9c-88b3-928e7b8d0cc8/tjbd_a_1914757_f0004_oc.jpg)
Figure 5. The paths of and
for stochastic model (Equation2
(2)
(2) ) with initial values
. Here
,
, other parameter values are given in (Equation15
(15)
(15) ) and (Equation16
(16)
(16) ).
![Figure 5. The paths of xt and yt for stochastic model (Equation2(2) dxt=xt(ar(t)−br(t)xt)−m1r(t)xt+m2r(t)yt−E1r(t)xtdt+σr(t)xtdBt,dyt=yt(ar(t)−br(t)yt)−m2r(t)yt+m1r(t)xt−E2r(t)ytdt+σr(t)ytdBt,(2) ) with initial values (x0,y0,r(0))=(1,2,1). Here E11=E12=0.82, E21=E22=0.8, other parameter values are given in (Equation15(15) b1=b2,m11=m12,m21=m22,E11=E12,E21=E22,σ1=σ2.(15) ) and (Equation16(16) a1=0.9,a2=0.3,b1=0.8,m11=0.4,m21=0.1,E11=0.3,E21=0.01,σ1=0.3.(16) ).](/cms/asset/cef009e8-b034-4a8f-8e0e-b50365b4527d/tjbd_a_1914757_f0005_oc.jpg)
As the intensity of poaching and
increases, the hybrid system changes from stochastic permanence to extinction (see Figure ). In other words, poaching has great influence on the survival of species.
Next, we will investigate the effective ways mentioned in [Citation31], which is increasing the migration rate and decreasing migration rate
simultaneously, to protect biodiversity.
Example 5.3
Based on Example 5.1, we use different vector fields values of and
to investigate the effect of the poaching on the survival of population. We then divide our simulations into two cases:
Case 1. | Without protection zone, i.e. | ||||
Case 2. |
|
Figure 6. The evolution of the total population of the stochastic model (Equation2
(2)
(2) ) is graphed for vield fields
,
,
. Other parameter values are given in (Equation15
(15)
(15) ) and (Equation16
(16)
(16) ).
![Figure 6. The evolution of the total population xt+yt of the stochastic model (Equation2(2) dxt=xt(ar(t)−br(t)xt)−m1r(t)xt+m2r(t)yt−E1r(t)xtdt+σr(t)xtdBt,dyt=yt(ar(t)−br(t)yt)−m2r(t)yt+m1r(t)xt−E2r(t)ytdt+σr(t)ytdBt,(2) ) is graphed for vield fields (m1,m2)=(0.4,0.1), (m1,m2)=(0,0), (m1,m2)=(0.5,0.05). Other parameter values are given in (Equation15(15) b1=b2,m11=m12,m21=m22,E11=E12,E21=E22,σ1=σ2.(15) ) and (Equation16(16) a1=0.9,a2=0.3,b1=0.8,m11=0.4,m21=0.1,E11=0.3,E21=0.01,σ1=0.3.(16) ).](/cms/asset/02553a50-ea9c-4e37-9a59-3b60d8003541/tjbd_a_1914757_f0006_oc.jpg)
Example 5.4
Based on Case 1 in Example 5.2, let the generator of be
Similarly, we have
and
. By Theorem 3.2 (i), one can see that the single species in the nature reserve and the natural environment are always extinct (see Figure ). For comparison, we plot the component-wise sample paths of
and
in Figure . From Figures and show that an interesting fact: if one of the subsystems is extinctive, another is stochastic persistent, then the hybrid system may be stochastic persistent if the Markov chain spends enough time in the stochastically permanent state, or vice versa. This means that Markov chain will balance the density of the population under different regimes. From biological point of view, the Markov chain
is conductive to the survival of species. That is to say, the Markov chain can suppress the extinction of the single species.
Figure 7. (a) The red line, green line and blue line represent the paths of for hybrid system, state 1 and state 2, respectively. (b) The red line, green line and blue line represent the paths of
for hybrid system, state 1 and state 2, respectively. Here
other parameters take the same values as in Figure .
![Figure 7. (a) The red line, green line and blue line represent the paths of xt for hybrid system, state 1 and state 2, respectively. (b) The red line, green line and blue line represent the paths of yt for hybrid system, state 1 and state 2, respectively. Here Γ=(−33833827−27), other parameters take the same values as in Figure 2.](/cms/asset/634f121d-60c6-4483-9e7c-23e4781ac3a5/tjbd_a_1914757_f0007_oc.jpg)
We next consider the effect of the regime switching on the survival of population by use two different switching rates.
Example 5.5
To see the effect of the regime switching, we consider two different switching rates, from slow to fast. The switching rates are ,
for the slow switching rate, and
,
for the fast switching rate. In addition, we assume that
and other parameter sets are the same as in the Example 5.1. We have seen that for each set of parameters we get an ordinary differential equations that has a stable positive equilibrium. The upper panel shows an example of relatively slow switching. In each environmental state, the population tends to a fixed point, specific to the environment. It then fluctuates about this fixed point. The PDF is bimodal. As the switching rate is increased (lower row), the stochastic dynamics spends more time in between the two fixed points and the bimodality of the stationary distribution is lost. At fast switching, the stationary distribution becomes unimodal, peaked at a value between the two fixed points. The system spends most of its time fluctuating about a point in the interior of phase space, away from either of the two fixed points. The result is shown in Figure .
Figure 8. Sample paths of the dynamics and the PDF of model (Equation2(2)
(2) ) with initial values
.
![Figure 8. Sample paths of the dynamics and the PDF of model (Equation2(2) dxt=xt(ar(t)−br(t)xt)−m1r(t)xt+m2r(t)yt−E1r(t)xtdt+σr(t)xtdBt,dyt=yt(ar(t)−br(t)yt)−m2r(t)yt+m1r(t)xt−E2r(t)ytdt+σr(t)ytdBt,(2) ) with initial values (x0,y0,r(0))=(1,2,1).](/cms/asset/81ee3f83-4f09-4f45-a9ae-338b4bc6bfb7/tjbd_a_1914757_f0008_oc.jpg)
It is well recognized that the environment is always subject to random disturbances, so it is important to take the impact of stochastic perturbations on the evolution of the species into consideration. Next, we simulate the solution model (Equation2(2)
(2) ) with different values of σ to investigate the effect of the white noise on the survival of the species.
Example 5.6
To see the effect of the white noise intensity, we consider three different fields values of σ: ,
and
. Other parameters take the same values as in Figure . By computing, we have
,
in the first two cases; and
in the last case. It then follows from Theorem 3.2 (ii) that the single species in the nature reserve and the natural environment are always persistent in mean, in the first two cases (
and
), while the single species in the nature reserve and the natural environment are always extinct in the last case (see Theorem 3.2 (i)). The computer simulations shown in Figure clearly support these results. Here we can conclude that the small noise does not change the stability of the equilibrium state of model (Equation2
(2)
(2) ), but with the intensity of white noise σ increasing, the volatility of
and
are getting larger. However when the environmental noises are relatively big, the persistence and extinction of the stochastic model (Equation2
(2)
(2) ) may be affected significantly.
Figure 9. The evolution of infected individuals of stochastic model (Equation2(2)
(2) ) is graphed for
,
and
. Other parameters take the same values as in Figure .
![Figure 9. The evolution of infected individuals of stochastic model (Equation2(2) dxt=xt(ar(t)−br(t)xt)−m1r(t)xt+m2r(t)yt−E1r(t)xtdt+σr(t)xtdBt,dyt=yt(ar(t)−br(t)yt)−m2r(t)yt+m1r(t)xt−E2r(t)ytdt+σr(t)ytdBt,(2) ) is graphed for σ= (0.1,0.1), (0.3,0.3) and (1.25,1.25). Other parameters take the same values as in Figure 1.](/cms/asset/48d40d1b-0c78-48b6-8431-acb80a507b33/tjbd_a_1914757_f0009_oc.jpg)
6. Discussion and conclusion
This article formulates a stochastic single-species model with regime switching and the impact of poaching. We first proved the existence of global positive solution. Then we investigated the extinction, persistence in mean. Moreover, we discussed the positive recurrence. Finally, numerical simulations were given in Section 5. Our investigation shows that the poaching, Markov chain, and white noise may cause great influence on the survival of species as follows:
Poaching had a large impact on the single-species viability. As the intensity and extent of poaching increases, the odds of being at risk of extinction increases for the single species. Increasing the criminal cost of poaching is an important way to effectively protect biodiversity.
The Markov chain
can suppress the extinction of the single species. Due to the hybrid system switching among different states, the species can get more chance to survival. That is to say, the larger probability of hybrid system switching to one good state, the greater the probability of survival. On the other hand, if only part of state is extinctive, the hybrid system still may be stochastically permanent, provided the Markov chain
spend enough time in good states (see Figures and ). From the point of view of ecology, the Markov chain
is beneficial for the survival of species. In other words, the Markov chain can suppress the extinction of the single species.
The white noise, if which is relatively small, cannot change the persistence of the biological populations. However, when the white noise is relatively big, the persistence and extinction of stochastic model (Equation2
(2)
(2) ) may be affected significantly.
From a management perspective, we can protect the biodiversity of a certain area through the following strategies:
Reduce the intensity of poaching
and
by increasing the criminal cost of poaching.
Establishing key protected areas and enlarging the size of protection zone, to increase the migration rate
and decrease the migration rate
.
At the same time, some interesting questions deserve further investigation. One may propose more realistic but complex models, such as incorporating the time delay into the system due to the time-delay environmentally fluctuating stochastic single species system exhibits an altogether different dynamical picture from model (Equation2(2)
(2) ). We leave this for future consideration.
Acknowledgements
This project was supported by the National Natural Science Foundation of China (No. 12001186), the Science Research Fund of Hunan Province Education Department (No. 20B158) and the Postgraduate Scientific Research Innovation Project of Hunan Province (No. CX20200466).
Disclosure statement
No potential conflict of interest was reported by the author(s).
Additional information
Funding
References
- R. Alley, J. Marotzke and W. Nordhaus, et al. Abrupt climate change, Science 299 (2003), pp. 2005–2010.
- M. Assaf, A. Kamenev and B. Meerson, Population extinction risk in the aftermath of a catastrophic event, Phys. Rev. E79 (2009), pp. 011127.
- H. Bauer, G. Chapron, K. Nowell and et al, Lion (Pantheraleo) populations are declining rapidly across Africa except in intensively managed areas, Proc. Natl. Acad. Sci. 112 (2015), pp. 14894–14899.
- P. Berck, Open access and extinction, Econometrica 47 (1979), pp. 877–882.
- F. Berkes, D. Feeny and B. Mccay, et al. The benefits of the commons, Nature 340 (1989), pp. 91–93.
- J. Blanco, S. Reig and L. De-La-Cuesta, Distribution status and conservation problems of the wolf Canis lupus in Spain, Biol. Conserv. 60 (1992), pp. 73–80.
- U. Breitenmoser, C. Breitenmoser-Würsten and S. Capt, et al. Conservation of the lynx Lynx lynx in the Swiss Jura Mountains, Wildlife Biol. 13 (2007), pp. 340–355.
- E. Bulte, H. Folmer and W. Heijman, Open access common property and scarcity rent in fisheries, Environ. Resour. Econ. 6 (1995), pp. 309–320.
- F. Courchamp, E. Angulo and P. Rivalan, et al. Rarity value and species extinction: the anthropogenic Allee effect, PLoS Biol. 4 (2006), pp. e415.
- B. Dubey, P. Chandra and P. Sinha, A model for fishery resource with reserve area, Nonlin. Anal. Real 4 (2003), pp. 625–637.
- T. Gard, Persistence in stochastic food web models, Bull.Math. Biol. 46 (1984), pp. 357–370.
- T. Gard, Stability for multispecies population models in random environments, Nonlin. Anal. 10 (1986), pp. 1411–1419.
- R. Hall, E. Milner-Gulland and F. Courchamp, Endangering the endangered: the effects of perceived rarity on species exploitation, Conserv. Lett. 1 (2008), pp. 75–81.
- L. Hu, M. Huang and M. Tang, et al. Wolbachia spread dynamics in stochastic environments, Theor. Popul. Biol. 106 (2015), pp. 32–44.
- L. Hu, M. Huang and M. Tang, et al. Wolbachia spread dynamics in multi-regimes of environmental conditions, J. Theor. Biol. 462 (2019), pp. 247–258.
- L. Hu, M. Tang and Z. Wu, et al. The threshold infection level for wolbachia invasion in random environments, J. Differ. Equ. 266 (2019), pp. 4377–4393.
- T. Kar and H. Matsuda, A bioeconomic model of a single-species fishery with a marine reserve, J. Environ. Manag. 86 (2008), pp. 171–180.
- R. Khasminskii, C. Zhu and G. Yin, Stability of regime-switching diffusions, Stoch. Proc. Appl. 117 (2007), pp. 1037–1051.
- J. Kenney, J. Smith and A. Starfield, et al. The long-term effects of tiger poaching on population viability, Conserv. Biol. 9 (1995), pp. 1127–1133.
- G. Lan, Z. Lin, C. Wei and S. Zhang, A stochastic SIRS epidemic model with non-monotone incidence rate under regime-switching, J. Frankl. Inst. 356 (2019), pp. 9844–9866.
- M. Liu, X. He and J. Yu, Dynamics of a stochastic regime-switching predator-prey model with harvesting and distributed delays, Nonlin. Anal. Hybrid. 28 (2018), pp. 87–104.
- M. Liu, W. Li and K. Wang, Persistence and extinction of a stochastic delay logistic equation under regime switching, Appl. Math. Lett. 26 (2013), pp. 140–144.
- X. Mao, Stochastic Differential Equations and Applications, Horwood Publishing, Chichester, 1997.
- R. May, Stability and Complexity in Model Ecosystems, Princeton University Press, Princeton, 2001.
- X. Mao and C. Yuan, Stochastic Differential Equations with Markovian Switching, Imperial College Press, London, 2006.
- W. Ripple, J. Estes and R. Beschta, et al. Status and ecological effects of the world's largest carnivores, Science343 (2014), pp. 1241484.
- J. Seidensticker, Saving wild tigers: a case study in biodiversity loss and challenges to be met for recovery beyond 2010, Integr. Zool. 5 (2010), pp. 285–299.
- M. Scheffer and S. Carpenter, Catastrophic regime shifts in ecosystems: linking theory to observation, Trends Ecol. Evol. 12 (2003), pp. 648–656.
- M. Scheffer, S. Carpenter and J. Foley, et al. Catastrophic shifts in ecosystems, Nature 413 (2001), pp. 591–596.
- M. Slatkin, The dynamics of a population in a Markovian environment, Ecol. 59 (1978), pp. 249–256.
- F. Wei and C. Wang, Survival analysis of a single-species population model with fluctuations and migrations between patches, Appl. Math. Model. 81 (2020), pp. 113–127.
- G. Yin and C. Zhu, Asymptotic properties of hybrid diffusion systems, SIAM J. Control Optim. 46 (2007), pp. 1155–1179.
- X. Yu, S. Yuan and T. Zhang, Persistence and ergodicity of a stochastic single species model with allee effect under regime switching, Commun. Nonlin. Sci. Numer. Simul. 59 (2018), pp. 359–374.
- Y. Zhao, S. Yuan and T. Zhang, The stationary distribution and ergodicity of a stochastic phytoplankton allelopathy model under regime switching, Commun. Nonlin. Sci. Numer. Simul. 37 (2016), pp. 131–142.
- X. Zou and K. Wang, A robustness analysis of biological population models with protection zone, Appl. Math. Model. 35 (2011), pp. 5553–5563.
- X. Zou and K. Wang, Dynamical properties of biological population with a protected area under ecological uncertainty, Appl. Math. Model. 39 (2015), pp. 6273–6284.
The proof of Corollary 3.1
Proof.
Denote
then from system (Equation2
(2)
(2) ), we can obtain
where
for
. In addition, we have applied the hypotheses
for
in here. Now consider the following equation
By Theorem 3.9 in [Citation23], we have
almost surely. Then according to the comparison theorem for stochastic differential equations, this yields
According to Theorem 3.1, for any initial value
, we have
which implies that
This completes the proof.