- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Discrete Dynamics in Nature and Society
Volume 2014 (2014), Article ID 402194, 10 pages
Dynamical Analysis of a Plateau Pika with Cross-Diffusion under Contraception Control
1School of Information and Communication Engineering, North University of China, Shanxi, Taiyuan 030051, China
2Department of Applied Mathematics, Yuncheng University, Shanxi, Yuncheng 044000, China
Received 6 July 2013; Revised 25 December 2013; Accepted 26 December 2013; Published 18 February 2014
Academic Editor: Seenith Sivasundaram
Copyright © 2014 Xiaoyan Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
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 is small enough. However, when is large enough, the model shows Turing bifurcation if . Furthermore, it is proved that if, 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.
Plateau pika (Ochotona curzoniae) is a species of small mammals in the Ochotonidae. According to the studies [1–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–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.  show that family structure has a strong influence on pika gene dynamics even though mating patterns of plateau pikas are variable. Smith and Cheeseman  show that a permanent contraceptive can give the chances of disease eradication which were very similar to that of lethal control. Burkey and Stenseth  discuss the effects of resource patchiness and the effects of a seasonal environment on the population dynamics of the herbivore species. Liu et al.  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  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.  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 , 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  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  consider the positive solution for the population cross-diffusion. Cross-diffusion inducing the Turing bifurcation can also be studied by [16–20]. One can refer to [13–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.
2. The Model Formulation
For the single species model under contraceptive control and lethal control, it is mentioned in the introduction . 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, holds in the procession. They obtain the following lemma. Define ,, and then it is easy to compute .
Lemma 1. When and , the positive equilibrium is globally asymptotically stable.
In this paper, we are concerned with the effect of the cross-diffusion under the contraceptive control and lethal control. The transmission dynamics of infection are governed by the following differential equations: where fertile female plateau pika () densities and sterile female plateau pika () densities are denoted at location , and time , respectively. is a fix bounded domain with smooth boundary , is the outward unit normal vector of the boundary , is the Laplacian operator, and positive constants , are the diffusive coefficients which denote the nature dispersive force of movement of the fertile female, and sterile female, respectively. The nonnegative constants and which mean the tendency that the sterile female (fertile female) keeps away from the fertile female (sterile female) are usually referred as cross-diffusion pressure. is the remove rate from free virus. For , with homogeneous Neumann boundary conditions and initial conditions 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 . Again by the maximum principle, and are bounded on . Therefore, the time existence interval can extend to by the standard theory for semilinear parabolic systems. The other parameters are similarly well defined as . 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.
Lemma 2 (Harnack inequality ). Assume that and let be a positive solution to Then there exists a positive constant such that .
Lemma 3 (maximum principle ). Suppose that . Assume that satisfies If , then .
3. Stability of Steady States
In this section, first we cite Lemma 1 and find that (3) still has a positive steady state if and . Second we discuss the linearized problem of (3) at . For simplicity, we denote Then problem (3) can be written as Let be the eigenvalues of the operator on with the homogeneous Neumann boundary conditions, and let be the eigenspace corresponding to in . Let , be an orthonormal basis of , and . Then
Let where The linearization of system (3) at is of the form , where . For each , is invariant under the operator , and is an eigenvalue of if and only if it is an eigenvalue of the matrix for some , in which case, there is an eigenvector in .
By a direct calculation, we have with thus the characteristic equation of the matrix is Note, when the system (3) has no diffusion effects, the stability of coexistence equilibrium was discussed by . There exists a sufficient condition for proving the existence and the global stability of the coexistence steady state.
Case 1. Consider .
In this case, it is easy to get the characteristic equation of the coexistence steady state as follows: where is the linear operator when and It is obviously that and . Thus, the real part of the characteristic roots of is less than 0. Therefore, we have the following theorem.
Case 2. Consider.
From the computation and notations, the trace of is less than 0, , and . Since , a necessary condition for of (3) becoming unstable is that for some is less than 0. Since is always nonnegative, the symmetric axis of the curve must be at the right half plane, so must be positive; that is, . At the same time, the equation of must have two positive roots; then must be satisfied for the instability of the positive coexistence steady state. Hence, has two positive roots
It is easy to know that, if for some , then . Therefore, we have the following theorem.
Theorem 5. Assume , and hold. If , and , then the positive coexistent steady state is unstable and (3) exhibits Turing bifurcation provided that, for some .
Remark 6. From the definition of , it is easy to get the following.(i)When , , . Then Turing bifurcation does not appear in (3).(ii)When , or ,, and is large enough, which make and , there exists some , such that . 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 is large enough.
4. 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 .
Theorem 7. Assume that ,,,are bounded; then there exist positive constants , such that the positive solution of (19) satisfies
Proof. For using the Lemma 3, we set , then elliptic equation (19) becomes
Let such that . By the first equation of (19) and Lemma 3, we have , and . Then, by the denotation of , we have Similarly, let such that . By the second equation of (19) and Lemma 3, we have ,. Then, by the denotation of , we have The proof is completely finished.
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.
Theorem 8. Assume that ,, are bounded; then there exist positive constants , such that the positive solution of (19) satisfies
By Lemma 2, there exists a constantsuch that . Hence, we have
By integrating the second equation of (19) in , we have
which implies that there exists , such that . Then . By combining with (26) yields that
Now we need to estimate the positive lower bound of . Suppose, on the contrary, that does not hold. Define . Then there exists a sequence with such that the corresponding positive solutions of (19) satisfy Using and the second equation (21) yields By Lemma 2, there exists a constantsuch that Similarly, we have the following relation between and It follows from (30) and (33) that There exists a subsequence of , which is denoted by , and two nonnegative functions , such that as . Combining (30) and (34), one gets ,. By the regularity theory of elliptic equations and the Sobolev embedding theorem, for some , the norm of is uniformly bounded with respect to . Then, taking the limits for the first equation of (29), we have Therefore, ; otherwise on , which is obviously impossible. So, .
Taking the limit and integrating the second equation of (29), we have Therefore, ; that is, , which is a contradiction with the existence of the positive coexistence steady state.
By Theorems 7 and 8, we can define a set where is defined by Theorems 7 and 8. It is easy to get the set which is invariant set for (3). We can discuss the dynamics of (3) in this set .
5. 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 with respect to elliptic equation (19). We obtain the following theorem.
Proof. Let for any . We denote the Lyapunov function as follows: We derivativewith along the system (3) and get We choose and substitute it into (39); thus Let be the largest invariant set of , that, is invariant with respect to (3). It is easy to get that if and only if . We show that consists of only the coexistence equilibrium . Since is invariant, we have and 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 on , then the Leray-Schauder degree 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 , we have the following lemma.
Lemma 10. If , then the number of negative eigenvalues of on is odd if and only if , for each . In addition,
Recall that , and so do . Then the sign of is determined by which is equivalent to since the matrix is . From the remark of Theorem 5, if is large enough, the determination may be less than 0.
In the following, we consider how the cross-diffusion coefficient impacts the determination. By the definitions of ,, and in Section 3, and , one obtains The two roots of are defined by and , respectively. Combining the definition of the upper limits, we have the following lemma.
Lemma 11. Assume ,. Then the two roots and of are real and satisfy the following limits for some such that, for all the ,
Next we will discuss the existence of nonconstant positive solution of (19) with respect to the cross-diffusion coefficient .
Theorem 12. Assume that the other parameters in (19) are fixed except , and , hold. is given by the limit in (51). If for some , and the sum is odd, then there exists a positive constant , such that, if , there exists at least one nonconstant positive solution which bifurcates from the coexistence steady state of (19).
Proof. Due to Lemma 11, there exists a positive , if ; then there exists two positive roots for and
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.
For , define and consider the problem Then is a nonconstant positive solution of (54) for any . By the definition of , for any , is a positive solution of (30) if and only if
It is easy to check that , Theorem 9 indicates that has only the positive solution in . Dividing with respect to , we have In particular, where ,. Then .
By the definition of , and Lemma 10, for , we have And note that is not an eigenvalue of the matrix for all , and If is odd, it follows from Lemma 10 that Note that when for Theorem 9. Thus we have
In addition, on for all by Theorems 7 and 8, while we use the homotopy invariance of the topological degree to obtain It is easy to know that is the root of both and in . Combining (60) and (61), we have which is contradiction to (62). This completes the proof.
6. Simulation and Discussion
We have formulated a plateau pika model with cross-diffusion. By Lemma 1, there is a unique positive coexistent steady state of (3) when , and . As the sterile female cross-diffusion is small enough, (3) does not exhibit Turing bifurcation. However, if the cross-diffusion rate is large enough, Turing bifurcation emerges in coexistence steady state. The system (3) is persistence and bound under some boundary conditions. If the cross-diffusion rates , 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 ; thus . Hence the protection of environment for the plateau pika is also important.
We can illustrate the theoretical results through two simulatinon examples. Set ,,,,,,,,. It is easy to get that and satisfy the conditions of Lemma 1. System (3) with above coefficients has a unique positive steady state . By Theorem 9, we see that the coexistence steady state of (3) is globally stable.
Numerical simulation illustrates our result (see Figure 1). Then we let ,, 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 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 cross-diffusion 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 ,. In addition, we do not study the dynamics of (3) when . In this case , the system may exhibit steady state bifurcation.
Conflict of Interests
The authors declare that they have no conflict of interests.
This study is supported by the NSF of China (61203228, 61275120, and 11371313), Mathematical Tianyuan Foundation (11226258), the Young Sciences Foundation of Shanxi (2011021001-1), the Natural Science Foundation of Shanxi Province (2013011002-5), program for the Outstanding Innovative Teams of Higher Learning Institutions of Shanxi, and the “131” talents.
- A. T. Smith and J. M. Foggin, “The plateau pika (Ochotona curzoniae) is a keystone species for biodiversity on the Tibetan plateau,” Animal Conservation, vol. 2, no. 4, pp. 235–240, 1999.
- C. H. Lai and A. T. Smith, “Keystone status of plateau pikas (Ochotona curzoniae): effect of control on biodiversity of native birds,” Biodiversity and Conservation, vol. 12, no. 9, pp. 1901–1912, 2003.
- C. Wang, H. He, M. Li et al., “Parasite species associated with wild plateau pika (Ochotona curzoniae) in southeastern Qinghai Province, China,” Journal of Wildlife Diseases, vol. 45, no. 2, pp. 288–294, 2009.
- Z. L. Wang, Y. Chen, J. Yang, W. J. Chen, Y. M. Zhang, and X. Q. Zhao, “Behavioral mechanisms of male sterilization on plateau pika in the Qinghai-Tibet plateau,” Behavioural Processes, vol. 89, no. 3, pp. 278–285, 2012.
- F. S. Dobson, A. T. Smith, and W. X. Gao, “The mating system and gene dynamics of plateau pikas,” Behavioural Processes, vol. 51, no. 1–3, pp. 101–110, 2000.
- G. C. Smith and C. L. Cheeseman, “A mathematical model for the control of diseases in wildlife populations: culling, vaccination and fertility control,” Ecological Modelling, vol. 150, no. 1-2, pp. 45–53, 2002.
- T. V. Burkey and N. C. Stenseth, “Population dynamics of territorial species in seasonal and patchy environments,” Oikos, vol. 69, no. 1, pp. 47–53, 1994.
- H. W. Liu, Z. Jin, Y. M. Chen, and F. Q. Zhang, “Population dynamics of plateau pika under lethalcontrol and contraception control,” Advances in Difference Equations, vol. 2012, article 29, 13 pages, 2012.
- H. W. Liu and Q. Y. Li, “Model of single-species population under contraceptive control and lethal control,” Mathematics in Practice and Theory, vol. 39, no. 15, pp. 104–107, 2009.
- H. W. Liu, L. Zhou, W. Liu, and H. K. Zhou, “Using a cellular-automata model to investigate the effects of grazing on plateau pika population dynamics,” International Journal of Biomathematics, vol. 4, no. 3, pp. 275–287, 2011.
- M. Z. Jessie and S. H. David, “Polygynandry and even-sexed dispersal in a population of collared pikas, Ochotona collaris,” Animal Behaviour, vol. 83, no. 4, pp. 1075–1082, 2012.
- F. S. Dobson, A. T. Smith, and W. X. Gao, “Social and ecological influences on dispersal and philopatry in the plateau pika,” Behavioral Ecology, vol. 9, no. 6, pp. 622–635, 1998.
- Y.-X. Wang and W.-T. Li, “Effects of cross-diffusion and heterogeneous environment on positive steady states of a prey-predator system,” Nonlinear Analysis: Real World Applications, vol. 14, no. 2, pp. 1235–1246, 2013.
- Y.-X. Wang and W.-T. Li, “Effect of cross-diffusion on the stationary problem of a diffusive competition model with a protection zone,” Nonlinear Analysis: Real World Applications, vol. 14, no. 1, pp. 224–245, 2013.
- P. Y. H. Pang and M. X. Wang, “Existence of global solutions for a three-species predator-prey model with cross-diffusion,” Mathematische Nachrichten, vol. 281, no. 4, pp. 555–560, 2008.
- J. Liu, H. Zhou, and L. Zhang, “Cross-diffusion induced Turing patterns in a sex-structured predator-prey model,” International Journal of Biomathematics, vol. 5, no. 4, Article ID 1250016, 23 pages, 2012.
- G.-P. Hu and X.-L. Li, “Turing patterns of a predator-prey model with a modified Leslie-Gower term and cross-diffusions,” International Journal of Biomathematics, vol. 5, no. 6, Article ID 1250060, 17 pages, 2012.
- C. R. Tian, Z. G. Lin, and M. Pedersen, “Instability induced by cross-diffusion in reaction-diffusion systems,” Nonlinear Analysis: Real World Applications, vol. 11, no. 2, pp. 1036–1045, 2010.
- G. Gambino, M. C. Lombardo, and M. Sammartino, “Turing instability and traveling fronts for a nonlinear reaction-diffusion system with cross-diffusion,” Mathematics and Computers in Simulation, vol. 82, no. 6, pp. 1112–1132, 2012.
- Z. Xie, “Cross-diffusion induced Turing instability for a three species food chain model,” Journal of Mathematical Analysis and Applications, vol. 388, no. 1, pp. 539–547, 2012.
- R. Ruiz-Baier and C. Tian, “Mathematical analysis and numerical simulation of pattern formation under cross-diffusion,” Nonlinear Analysis: Real World Applications, vol. 14, no. 1, pp. 601–612, 2013.
- W. Ko and K. Ryu, “On a predator-prey system with cross diffusion representing the tendency of predators in the presence of prey species,” Journal of Mathematical Analysis and Applications, vol. 341, no. 2, pp. 1133–1142, 2008.
- W.-M. Ni, “Diffusion, cross-diffusion, and their spike-layer steady states,” Notices of the American Mathematical Society, vol. 45, no. 1, pp. 9–18, 1998.
- T. Kolokolnikov and J. Wei, “Stability of spiky solutions in a competition model with cross-diffusion,” SIAM Journal on Applied Mathematics, vol. 71, no. 4, pp. 1428–1457, 2011.
- Y. Lou and W. M. Ni, “Diffusion, self-diffusion and cross-diffusion,” Journal of Differential Equations, vol. 131, no. 1, pp. 79–131, 1996.
- Y. Lou and W.-M. Ni, “Diffusion vs cross-diffusion: an elliptic approach,” Journal of Differential Equations, vol. 154, no. 1, pp. 157–190, 1999.
- Y. Lou, W.-M. Ni, and Y. Wu, “On the global existence of a cross-diffusion system,” Discrete and Continuous Dynamical Systems, vol. 4, no. 2, pp. 193–203, 1998.
- P. V. Tuôc, “Global existence of solutions to Shigesada-Kawasaki-Teramoto cross-diffusion systems on domains of arbitrary dimensions,” Proceedings of the American Mathematical Society, vol. 135, no. 12, pp. 3933–3941, 2007.
- X. Zeng, “Non-constant positive steady states of a prey-predator system with cross-diffusions,” Journal of Mathematical Analysis and Applications, vol. 332, no. 2, pp. 989–1009, 2007.
- L. S. Lin, W. M. Ni, and I. Takagi, “Large amplitude stationary solutions to a chemotaxis system,” Journal of Differential Equations, vol. 72, no. 1, pp. 1–27, 1988.