Dynamical Analysis of a Plateau Pika with Cross-Diffusion under Contraception

A plateau pika model with spatial cross-diffusion is investigated. By analyzing the corresponding characteristic equations, the local stability of an coexistence steady state is discussed when d 21 is small enough. However, when d 21 is large enough, the model shows Turing bifurcation if B2 −4AC > 0. Furthermore, it is proved that if, R > R 0 , βK > d and cross-diffusion rates are zero, the positive coexistence steady state is globally asymptotically stable. A nonconstant positive solution bifurcates from the coexistent steady state by the Leray-Schauder degree theory. Numerical simulations are carried out to illustrate the main results.


Introduction
Plateau pika (Ochotona curzoniae) is a species of small mammals in the Ochotonidae.According to the studies [1][2][3][4], it is found that the plateau pika's burrows are the primary homes to a wide variety of small birds and lizards since Qinghai-Tibet area owns the treeless environment and also minimize soil erosion for the suitable number of plateau pikas, enhance the ability of soil to absorb precipitation, contribute to nutrient cycling, and create microhabitats resulting in increased plant species richness; it serves as the main prey for most of the predatory animals on the plateau.Hence, the plateau pika is a keystone species in Qinghai-Tibet plateau.Since people overgraze and climate becomes warm, which makes the alpine meadow degrade seriously and provides the suitable environment for the plateau pika, multiplication of plateau pika is overabundance.According to the statistics data [1][2][3][4], there exists the huge number of the plateau pikas in this areas, and it threatens the forage that could otherwise be utilized by the livestock and the huge number of burrows results in soil erosion, vegetation disturbances, and the individuals who ride horses across the grasslands hazard.Hence plateau pika has been listed as a pest for the Chinese Government.
During the recent years, the various methods are applied to manage the number of plateau pikas.In the first stage, since people are lacking enough understanding about the plateau pika, they directly poisoned the plateau pike with botulin of models C and D. The people considered the plateau pika as rodent; in fact, pikas belong to lagomorphs and not to rodents.From the first paragraph, plateau pike is keystone species, whose elimination or major decimation from an ecosystem would have a greater average effect on other species' populations or ecosystem processes.Hence, plateau pika contributes to maintaining the biodiversity in Qinghai-Tibet plateau [1,2].Obviously it does not work for the people to kill the plateau pika with any chemical sterilant.For the balance of the ecology and long-term development, immunocontraception has been used to manage plateau pika.Immunocontraception is different from the traditional method to manage the plateau pika.It aims to reduce the birth rate of the plateau pika and increase the death rate of the plateau pika.This control method can make the plateau pika be maintained at a low or suitable level for a long period.If we carry out the control method at the breeding season, the number of plateau pikas may be controlled at a desirable level.Several mathematical models describing the plateau pika under sterile control have been employed and analyzed.Dobson et al. [5] show that family structure has a strong influence on pika gene dynamics even though mating patterns of plateau pikas are variable.Smith and Cheeseman 2 Discrete Dynamics in Nature and Society [6] show that a permanent contraceptive can give the chances of disease eradication which were very similar to that of lethal control.Burkey and Stenseth [7] discuss the effects of resource patchiness and the effects of a seasonal environment on the population dynamics of the herbivore species.Liu et al. [8] proposed a discrete model incorporating the sterile and the lethal control.In this model, the reproductive number plays a key role in impacting the dynamics of the model.Liu and Li [9] build a mathematical model to discuss the single species population under contraceptive control and lethal control.They studied the stability of the equilibrium of the model.Liu et al. [10] use a cellular-automata model to investigate the influence of grazing on dynamics of plateau pika population.The dispersal dimension depends on the meadow degradation.
Animal populations may not be distributed randomly in space but exhibit spatial patterns, and plateau pika is not exception.Plateau pika spatial distribution in the plateau is directly affected by individual interactions and interspecific interactions.Plateau pika dispersing is a natural phenomenon.Pikas from the different plateaus often contact and move towards favorable habitats due to climate, food, predators, and so on.This phenomena can employ random diffusion combined with directed movement upward along environmental gradients.Many interesting papers investigating species diffusion have been discussed by many authors.In [11], there is no difference in the average distance dispersed by the male and female sexes, since the costs and benefits of dispersal do not differ between males and females in the species.Dispersal movements in [12] result in equalization of density among pika families, if competition for environmental resources influenced dispersal.Wang and Li [13,14] discuss the positive stationary problems in population cross-diffusion using the bifurcation theory.Pang and Wang [15] consider the positive solution for the population crossdiffusion.Cross-diffusion inducing the Turing bifurcation can also be studied by [16][17][18][19][20].One can refer to [13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29] and references therein for more details.Although there are many references which discussed the population diffusion and cross-diffusion, plateau pika population diffusion or cross-diffusion under contraception still remains a problem.At the best of our knowledge as the authors known, we do not find any relative report or study until now.Hence pikas diffusion is very interesting and is worth studying.
This paper is organized as follows.Section 2 introduces a plateau pika model under contraceptive control and lethal control.In Section 3, local stability of the positive steady state is studied.In Section 4, we obtain persistence of the system.In Section 5, global stability and existence of the nonconstant steady state are investigated.The conclusion is given and some simulations to illustrate the theoretical results.

The Model Formulation
For the single species model under contraceptive control and lethal control, it is mentioned in the introduction [9].The total population is divided into two classes: female plateau pika sterile female plateau pika.The numbers of individuals in these classes at time  are denoted, respectively, by () and ().The population density decreases due to the control; a compensatory reaction from all or part of the population occurs due to increased access to resources.In this model, a fraction of the population is implemented by immunocontraception; they separate the female population into two contributions explicitly; one (the birth rate ) is norm and the other (the death rate ) is immunocontraception.They assume that the immunocontraception just impacts the birth rate .Hence birth rate  and death rate  can be expressed as ( This model is described as follows: Here,  is the birth rate,  is mortality rate, the intrinsic growth rate  =  − ,  is the death rate due to chemical sterilant,  denotes the carrying capacity, and  is efficient transmitted rate from female plateau pika to sterile female plateau pika though the contact.For the biological meaning of the model,  −  > 0 holds in the procession.They obtain the following lemma.Define  = −−,  0 = (−)(+)/, and then it is easy to compute  * * = ( * ,  * ) = (( + )/, ( −  0 )/( −  + )).Lemma 1.When  >  0 and  > , the positive equilibrium  * * is globally asymptotically stable.
In this paper, we are concerned with the effect of the crossdiffusion under the contraceptive control and lethal control.The transmission dynamics of infection are governed by the following differential equations: The boundary conditions in (4) imply that the populations do not move across the boundary Ω.By the maximum principle, both (, ) and (, ) are positive for  ∈ Ω and  ∈ [0, max ).Again by the maximum principle,  and  are bounded on Ω × [0,  max ).Therefore, the time existence interval can extend to ∞ by the standard theory for semilinear parabolic systems.The other parameters are similarly well defined as [9].In addition, the model ( 3) is well posed that is, there exists a unique classical solution of (3).
For the convenience to study the dynamics of (3), we first give the classical results for the elliptic equation.

Stability of Steady States
In this section, first we cite Lemma 1 and find that (3) still has a positive steady state  * = ( * ,  * ) if  >  0 and  > .Second we discuss the linearized problem of (3) at  * = ( * ,  * ).For simplicity, we denote Then problem (3) can be written as Let 0 =  1 <  2 < ⋅ ⋅ ⋅ be the eigenvalues of the operator −Δ on Ω with the homogeneous Neumann boundary conditions, and let (  ) be the eigenspace corresponding to   in  1 (Ω).
Let X = ( 1 (Ω)) 2 , {  ,  = 1, . . ., dim (  )} be an orthonormal basis of (  ), and Let where The linearization of system (3) at  * is of the form   = L, where  = (   ).For each  ≥ 1, X  is invariant under the operator L, and  is an eigenvalue of L if and only if it is an eigenvalue of the matrix −  Φ  ( * ) + Ψ( * ) for some  ≥ 1, in which case, there is an eigenvector in X  .
By a direct calculation, we have with thus the characteristic equation of the matrix L is Note, when the system (3) has no diffusion effects, the stability of coexistence equilibrium was discussed by [9].There exists a sufficient condition  >  0 ,  >  for proving the existence and the global stability of the coexistence steady state.
In this case, it is easy to get the characteristic equation of the coexistence steady state  * as follows: where L 1 is the linear operator L when  21 =  12 = 0 and It is obviously that trace(L 1 ) < 0 and det(L 1 ) > 0. Thus, the real part of the characteristic roots of L 1 is less than 0. Therefore, we have the following theorem.
Case 2. Consider  12  21 ̸ = 0. From the computation and notations, the trace of L is less than 0,  > 0, and  > 0. Since   ≥ 0, a necessary condition for  * of (3) becoming unstable is that det L for some   is less than 0. Since   is always nonnegative, the symmetric axis of the curve  2 −  +  must be at the right half plane, so  must be positive; that is,  21 >  12  * +  2 ( − )/.At the same time, the equation of  2 −  +  = 0 must have two positive roots; then  2 − 4 > 0 must be satisfied for the instability of the positive coexistence steady state.Hence,  2   −   +  = 0 has two positive roots It is easy to know that, if  1 <   <  2 for some   ≥ 0, then det{L} < 0. Therefore, we have the following theorem.(ii) When  12 = 0,  21 ̸ = 0 or  12 ̸ = 0,  21 ̸ = 0, and  21 is large enough, which make  > 0 and  2 − 4 > 0, there exists some   > 0, such that  1 <   <  2 .Under this situation,  * is unstable while  * is stable for the ODE situation.Hence, Turing bifurcation may emerge in (3) if the cross-diffusion coefficients of the sterile female plateau pika  21 is large enough.

A Priori Estimation
In this section, we will be interested primarily in steady state of (3).Now we first give a priori estimate for positive solutions to the elliptic system which is the corresponding steady state of (3).
In the following, we will discuss the upper bounds and the lower bounds by using Lemma 2 and Lemma 3. By the regularity theory for elliptic equations, it is easy to know that the positive solution of ( 19) is in  2 (Ω) ×  2 (Ω).
From Theorem 7, the population of plateau pika for (3) is boundary.Since the maximum capability carrying has limitation, the population of the plateau pika for (3) can not tend infinity.For proving the persistence of the plateau pika, we can obtain the following theorem.
Taking the limit and integrating the second equation of ( 29), we have Therefore,  = (+)/ = /(−); that is,  =  0 , which is a contradiction with the existence of the positive coexistence steady state.By Theorems 7 and 8, we can define a set where   ,   = 1, 2 is defined by Theorems 7 and 8.It is easy to get the set B() which is invariant set for (3).We can discuss the dynamics of (3) in this set B.

Global Stability of the Steady State and the Existence of Nonconstant Positive Solution
In order to discuss the existence of nonconstant positive solution of (3), we first prove the global stability of the coexistence equilibrium  * , when the cross-diffusion coefficients are zero.Then there is no nonconstant solution respect to elliptic equation (19).We obtain the following theorem.
Proof.Let  * = 1/|Ω| ∫ Ω   for any  ∈  1 (Ω).We denote the Lyapunov function as follows: We derivative  with  along the system (3) and get We choose  = (( − ) + )/ and substitute it into (39); thus Let  be the largest invariant set of {()/ = 0}, that, is invariant with respect to (3).It is easy to get that ()/ = 0 if and only if  =  * .We show that  consists of only the coexistence equilibrium  * .Since  is invariant, we have ()/ = 0 and Δ = 0 in .Then from the first equation of ( 3), it holds that We must have (, ) =  * for all , .Thus, the set  consists of only the endemic equilibrium  * .This completes the proof.
In the following, we consider that the cross-diffusion induces the nonconstant solution for (19).The Leray-Schauder degree theory is a good method to establish the existence of stationary patterns of (19).For the convenience, we define a positive cone as In order to use Leray-Schauder degree, we must revise system (19) as a fixed problem.From the definition of elliptic operator Δ, we have Hence, if (⋅) ̸ = 0 on B, then the Leray-Schauder degree deg((⋅), 0, B) is well defined.By a direct computation, it is easy to obtain Thus, the index of the operator  at coexistence steady state is defined as the following formula: if   ( * ) is invertible  which is considered as the total number of eigenvalues of   ( * ) with negative real parts (counting multiplicities).Since   is the characteristic eigenvalue of the operator −Δ, the eigenvalue of   ( * ) is determined by the following matrix: For the denotational convenience, we denote According to the references [15], we have the following lemma.
In the following, we consider how the cross-diffusion coefficient  21 impacts the determination.By the definitions of , , and  in Section 3, and A(,  21 ), one obtains lim The two roots of A(,  21 ) = 0 are defined by μ1 ( 21 ) and μ2 21 )(Re( μ1 ( 21 )) ≤ Re ( μ2 ( 21 ))), respectively.Combining the definition of the upper limits, we have the following lemma.
Lemma 11.Assume  >  0 ,  > .Then the two roots μ1 ( 21 ) and μ2 ( 21 ) of A(,  21 ) are real and satisfy the following limits for some  * 12 such that, for all the  12 >  * 12 , Next we will discuss the existence of nonconstant positive solution of (19) In order to prove the result, we use the contradict method.We assert that the result is not true.By using the homotopy argument, we can derive a contradiction.(63) which is contradiction to (62).This completes the proof.

Simulation and Discussion
We have formulated a plateau pika model with crossdiffusion.By Lemma 1, there is a unique positive coexistent steady state of (3) when  >  0 , and  > .
As the sterile female cross-diffusion is small enough, (3) does not exhibit Turing bifurcation.However, if the crossdiffusion rate is large enough, Turing bifurcation emerges in coexistence steady state.The system (3) is persistence and bound under some boundary conditions.If the crossdiffusion rates  12 ,  21 are both zeros, the coexistence steady state is globally stable by using the Lyapunov functions.
As the pressure from the sterile female plateau pika to the female plateau pika is large enough, a nonconstant steady state solution will bifurcate from the positive steady state by using Leray-Schauder degree theory.From the biological view, the sterile female plateau pika giving the pressure to the the female plateau pika plays an important role in controlling the population of plateau pika in Qinghai-Tibet area, since the plateau pika is one keystone species of plateau.Blending the population of plateau pika to eradicate is not realistic and may disturb the biological balance in the plateau.We should carry out the measures to control the population of the plateau pika tending some level, which is not harmful to the plateau, and is also benefitcial to persisting the biological balance.From the mathematical view, it is better for  * ≤  where  is the level for the persisting ecological cycle.Note that  * =  * +  * ; it is easy to get the following formula though a direct computation: To understand how the total population of the plateau pika can be lowered, we consider its elasticity with respect to the key parameters that can be affected by control measures.In particular, we define elasticity of quantity  with respect to parameter ] as We can compute the following elasticity quantities: Assume that  > ,  >  satisfies the existence condition of the positive coexistence steady state.We note that    * ,  = , ,  is positive, while    * is always negative.Hence, if we want to control the population of plateau pika, we should decrease the birth rate, death rate, and capability, while increasing the transmitted rate for themselves from the theoretical view.However, the death rate  and capability carrying  are harder controlled than the birth rate  and transmitted rate .For the specific control measure, increasing the transmitted rate and decreasing the birth rate are possible.Hence the contraception is a better method for this.For the formula of the elasticity of  * with respect to  and , if the impact of birth rate is larger than the transmitted rate,  must satisfy  2 −  +  < 0; thus  2 > 4.Hence the protection of environment for the plateau pika is also important.
Numerical simulation illustrates our result (see Figure 1).Then we let  12 = 0.01,  21 = 0, and other parameters be the same as Figure 1.System (3) with above coefficients still has a unique positive steady state  * which is the same as Figure 1.By Theorem 5, we see that the coexistence steady state  * of problem (3) is locally stable.Numerical simulation illustrates our results (see Figure 2).We can also know that the pressure from the female plateau pika to sterile female plateau pika diffusion  12 does not change the stability of system (3).However, we also find that the diffusion can advance the time for both the female plateau pika and the sterile female plateau pika to arrive at the coexistence steady state.
For the management of the application, larger  resulted in small number of the plateau pikas.Hence we should try our best to combine a powerful contraception method with some additional measures to control the plateau pikas.In addition, the diffusion can advance the time for both the female plateau pika and the sterile female plateau pika to arrive at the coexistence steady state.Increasing the crossdiffusion coefficients can rich the dynamics of the model and enhance the difficulty of the control.The system may show the Turing bifurcation.Unfortunately, we are unable yet to study global stability of coexistence steady state  * when  12 ̸ = 0,  21 = 0.In addition, we do not study the dynamics of (3) when  2 − 4 = 0.In this case det L = 0, the system may exhibit steady state bifurcation.