Positive solutions for a three-trophic food chain model with diffusion and Beddington–Deangelis functional response☆
Introduction
Understanding of spatial and temporal behaviors of interaction species in ecological systems is a center issue in population ecology. One aspect of great interest for a model with multispecies interactions is whether the involved species can persist or even stabilize at a coexistence steady state. In the case whether the species are homogeneously distributed, this would be indicated by a constant positive solution of a ordinary differential equation system. In the spatially inhomogeneous case, the existence of a nonconstant time-independent positive solution, also called stationary pattern, which is an indication of the richness of the corresponding partial differential equation dynamics. In recent years, stationary pattern induced by diffusion has been studied extensively, and many important phenomena have been observed.
The role of diffusion in the modeling of many physical, chemical and biological processes has been extensively studied. Starting with Turing’s seminal 1952 paper [1], diffusion and cross-diffusion have been observed as causes of the spontaneous emergence of ordered structures, called patterns, in a variety of non-equilibrium situations. These include the Gierer–Meinhardt model [2], [3], [4], [5], [6], the Sel’kov model [7], [8], the Noyes–Field model for Belousov–Zhabotinskii reaction [9], the chemotactic diffusion model [10], [11], the Brusselator model [12], the competition model [13], [14], [15], [16], [17], the predator–prey model [18], [19], [20], [21], [22], [23], [24], [25], [26], as well as models of semiconductors, plasmas, chemical waves, combustion systems, embryogenesis, etc., see e.g. [27], [28], [29] and references therein. Diffusion-driven instability, also called Turing instability, has also been verified empirically [30], [31].
This paper studies the stationary solution of the parabolic model with initial values in , where is a bounded domain of , , with boundary of class for some , is the outward unit vector on , the parameters are strictly positive, and , and stand for the initial conditions. We denote .
Problem (1.1) models three species food chain with diffusion and Beddington–Deangelis functional response, where are the respective densities of prey, predator, and top predator (cf. [32], [33], [34], and the list of references therein).
In our work here, we are mainly characterizing the existence of steady-states of (1.1), which are the solution of
Owing to the classical theory of parabolic equations, the solutions of (1.1) are globally defined in time and satisfy So, in this paper, we will discuss the nonnegative solutions of (1.2). We have a trivial nonnegative solution . As in [35], [36], the other nonnegative solutions of (1.2) can be classified by three types:
- (i)
nonnegative solutions with exactly two components identically zero
- (ii)
nonnegative solutions with exactly one component identically zero
- (iii)
nonnegative solutions with no component identically zero.
In the present work, we attempt to further understand the influence of diffusion and functional response on pattern formation. As a consequence, the existence and nonexistence results for nonconstant positive steady state solution to (1.2) indicate the stationary pattern arises as the diffusion coefficient enter into certain regions. In other words, diffusion does help to create stationary pattern. On the other hand, our results also demonstrate that diffusion and functional response can become determining factors in the formation pattern.
The distribution of this paper is the following. In Section 2, we give some fundamental theorems. In Section 3, the necessary and sufficient conditions for positive solutions of (1.2) are investigated. To achieve this, some degree theorems are developed, which is playing an important role in the text. Finally, in Section 4, sufficient conditions for the extinction and permanence to the time-dependent system (1.1) are investigated.
Section snippets
Preliminaries
In this section, we give some fundamental theorems, especially some degree theorems, which play an important role in the text.
For each and , denote the principle eigenvalue of by and simply denote by . It is well known that is strictly increasing in the sense that and implies (see Proposition 1.1 of [37]).
Theorem 2.1 For and be a sufficiently large number such [37], [38], [39]
Existence of positive solutions
To give some sufficient conditions for the existence of positive solutions of (1.2) by using fixed point index theory, we need an a priori estimate for nonnegative solutions of (1.2). So, we first give the following theorem: Theorem 3.1 If and , any nonnegative solution of (1.2) has an a priori bounds:whereand Proof From the first equation of (1.2), we have
Asymptotic behavior: extinction and global attractor
In this section, the asymptotic behavior of the time-dependent solutions of (1.1) is considered, that is, sufficient conditions for the extinction and permanence to system (1.1) are investigated.
Theorem 4.1 Let be a positive solution of (1.1), then we have (i) If , then as ; (ii) If and , then as .
Proof (i) First, we observe that any time-dependent solution of (1.1) satisfies
Acknowledgement
The authors would like to thank the referee for the careful reading of this paper and for the valuable suggestions that improved the presentation and style of the paper.
References (54)
- et al.
The stability of spike solutions to the one-dimensional Gierer–Meinhardt model
Physica D
(2001) - et al.
A Spikes for the Gierer–Meinhardt system in two dimensions: the strong coupling case
J. Differential Equations
(2002) Mini-maximizer for reaction–diffusion systems with skew-gradient structure
J. Differential Equations
(2002)Non-constant positive steady states of the Sel’kov model
J. Differential Equations
(2003)- et al.
Positive steady-state solutions of the Noyes–Field model for Belousov–Zhabotinskii reaction
Nonlinear Anal.
(2004) - et al.
Large amplitude stationary solutions to a chemotaxis systems
J. Differential Equations
(1988) - et al.
Pattern formation of a coupled two-cell Brusselator model
J. Math. Anal. Appl.
(2010) - et al.
Stationary patterns created by cross-diffusion for the competitor-mutualist model
J. Math. Anal. Appl.
(2004) - et al.
Dynamics of Lotka–Volterra competition systems with large interaction
J. Differential Equations
(2002) - et al.
Diffusion, self-diffusion and cross-diffusion
J. Differential Equations
(1996)
Diffusion vs cross-diffusion: an elliptic approach
J. Differential Equations
Existence and instability of Neumann layer solutions for a 3-component Lotka–Volterra model with diffusion
J. Math. Anal. Appl.
Stationary patterns of strongly coupled predator–prey models
J. Math. Anal. Appl.
Coexistence states of a Holling type-II predator–prey system
J. Math. Anal. Appl.
A ratio-dependent food chain model and its applications to biological control
Math. Biosci.
Positive solutions for a three-species competition system with diffusion-I. General existence results
Nonlinear Anal.
Positive solutions for a three-species competition system with diffusion-II. The case of equal birth rates
Nonlinear Anal.
Positive solution for Lotka–Volterra systems wiht cross-diffions
Existence and structure of the set of positive solutions of a general class of sublinear elliptic non-classical mixed boundary value problems
Nonlinear Anal.
Properties of the principle eigenvalues of a general class if non-classical mixed boundary value problems
J. Differential Equations
On positive solutions of some pairs of differential equations, II
J. Differential Equations
On the indices of fixed points of mapping in cones and applications
J. Math. Anal. Appl.
Coexistence regions in Lotka–Volterra models with diffusin
Nonlinear Anal.
Quasisolutions and global attractor of reaction–diffusion systems
Nonlinear Anal.
The chemical basis of morphogenesis
Philos. Trans. Roy. Soc. Ser. B
A theory of biological pattern formation
Kybernetik
Asymmetric spike patterns for the one-dimensional Gierer–Meinhardt model: equilibria and stability
European J. Appl. Math.
Cited by (18)
Asymptotic behavior and multiplicity for a diffusive Leslie-Gower predator-prey system with Crowley-Martin functional response
2014, Computers and Mathematics with ApplicationsCitation Excerpt :Therefore, the predator–prey systems with various functional responses and different boundary conditions have been proposed and studied by many ecologists and mathematicians. These studies applied classic types of functional responses including: Lotka–Volterra type [1,2], Holling type [3–7], Beddington–DeAngelis type [8–12], ratio-dependent type [13–16] and so on. The paper is organized as follows.
Global attractivity of equilibrium in Gierer-Meinhardt system with activator production saturation and gene expression time delays
2013, Nonlinear Analysis: Real World ApplicationsCitation Excerpt :Since the pioneering work of Turing [1], Reaction–Diffusion systems have been used to demonstrate morphogenetic pattern formation [2–7].
Positive solutions of a diffusive predator-prey model with modified Leslie-Gower and Holling-type II schemes
2012, Journal of Mathematical Analysis and ApplicationsCoexistence states for a diffusive one-prey and two-predators model with B-D functional response
2012, Journal of Mathematical Analysis and ApplicationsDynamical behaviors of a diffusive predator-prey system with Beddington-DeAngelis functional response
2014, International Journal of BiomathematicsDynamics in a diffusive predator-prey system with double Allee effect and modified Leslie-Gower scheme
2022, International Journal of Biomathematics
- ☆
J. Zhou is supported by the Fundamental Research Funds for the Central Universities (No. XDJK2009C069) and C.L. Mu is supported by NNSF of China (No. 10771226).