Next Article in Journal
An Improved Tikhonov-Regularized Variable Projection Algorithm for Separable Nonlinear Least Squares
Next Article in Special Issue
Degenerated and Competing Dirichlet Problems with Weights and Convection
Previous Article in Journal
Differential Evolution with Shadowed and General Type-2 Fuzzy Systems for Dynamic Parameter Adaptation in Optimal Design of Fuzzy Controllers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Directed Dynamic Behaviors of a Lotka-Volterra Competition-Diffusion-Advection System

College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, China
*
Author to whom correspondence should be addressed.
Axioms 2021, 10(3), 195; https://doi.org/10.3390/axioms10030195
Submission received: 7 July 2021 / Revised: 17 August 2021 / Accepted: 18 August 2021 / Published: 20 August 2021
(This article belongs to the Special Issue Nonlinear Dynamical Systems with Applications)

Abstract

:
This paper investigates the problem of the global directed dynamic behaviors of a Lotka-Volterra competition-diffusion-advection system between two organisms in heterogeneous environments. The two organisms not only compete for different basic resources, but also the advection and diffusion strategies follow the dispersal towards a positive distribution. By virtue of the principal eigenvalue theory, the linear stability of the co-existing steady state is established. Furthermore, the classification of dynamical behaviors is shown by utilizing the monotone dynamical system theory. This work can be seen as a further development of a competition-diffusion system.

1. Introduction

In the past few decades, the dynamic behaviors of competition-diffusion systems (see [1]) in homogeneous or heterogeneous environments have been extensively studied. Until 2017, He and Ni [2,3] studied the dynamics of two organisms by changing their random diffusion coefficients, resource carrying capacity and competitiveness, and they also described the global dynamics of two organisms. Their research has made outstanding contributions to the competition-diffusion systems. For the competition model of two organisms, either both organisms survive or win with the extinction of the other organisms (see [4,5,6]). In 2019, Du et al. [7,8] studied a Lotka-Volterra competition system with periodic habitat advection. From a biological point of view, this pulsating travel front provided a way for two competing species to interact in heterogeneous habitats. Based on the assumption that the resource function in spatial variables is decreasing, Lou et al. [9] described the competition between two aquatic organisms with different diffusion strategies for the same resource in the Lotka-Volterra reaction-diffusion-advection system in 2019. Md. Kamrujjaman [10] studied the impact of diffusion strategies on the outcome of competition between two populations while the species are distributed according to their respective carrying capacities in competition-diffusion systems. However, in the competition-diffusion-advection systems, the study of different species with different distribution functions will be more complex. Tang and Chen [11] and Xu et al. [12] studied the population dynamics of competition between two organisms from the perspective of river ecology in 2020. One interesting feature of their system was that the boundary conditions at the upstream end and downstream end can represent the net loss of individuals. In some cases, both organisms leave the site of competition, neither coexisting nor becoming extinct. Such an environment is important enough to demonstrate how organisms change their density and survival time in competition (see [13]). In 2021, Ma and Guo [14] described the feature of the coincidence of bifurcating coexistence steady-state solution branches and the effect of advection on the stability of the bifurcating solution. However, it is worthwhile to point out that all the aforementioned works focus on the global dynamic behaviors of competition-diffusion systems (see [10,15,16]) or advection systems (see [17,18]), in which the diffusion rates and spatial carrying capacity are changed, or the periodic habitat of advection systems is studied, or the upstream and downstream boundary conditions are changed.
Motivated by the effort of the aforementioned studies, we investigate the problem of the global directed dynamic behaviors of a Lotka-Volterra advection system between two organisms in heterogeneous environments, where two organisms are competing for different fundamental resources, their advection and diffusion strategies follow the dispersal towards a positive distribution, and the functions of inter-specific competition ability are variable.
Hence, we discuss the following global dynamics of the advection system:
U t = · [ κ 1 ( x ) ( U Q ( x ) ) μ 1 ( x ) U Q ( x ) ω 1 ( x ) ] + U [ r 1 ( x ) U ρ 2 ( x ) V ] , i n   Ω × R + , V t = · [ κ 2 ( x ) ( V Q ( x ) ) μ 2 ( x ) V Q ( x ) ω 2 ( x ) ] + V [ r 2 ( x ) ρ 1 ( x ) U V ] , i n   Ω × R + , κ 1 ( x ) n ( U Q ) μ 1 ( x ) U Q ω 1 ( x ) n = 0 , o n   Ω × R + , κ 2 ( x ) n ( V Q ) μ 2 ( x ) V Q ω 2 ( x ) n = 0 , o n   Ω × R + , U ( x , 0 ) = U 0 ( x ) , 0 , i n   Ω , V ( x , 0 ) = V 0 ( x ) , 0 , i n   Ω ,
where U ( x , t ) and V ( x , t ) are the population densities of biological organisms, location x   i n   Ω and time t > 0 , which are supposed to be nonnegative; κ 1 ( x ) , κ 2 ( x ) > 0 correspond to the dispersal rates of two competing organisms U and V, respectively. ∇ is the gradient operator. μ 1 ( x ) , μ 2 ( x ) > 0 correspond to the advection rates of two competing organisms U and V, and ω 1 ( x ) , ω 2 ( x ) C 2 ( Ω ¯ ) are the nonconstant functions and represent the advective direction. The intrinsic growth rates of the two competing organisms are bounded functions r 1 ( x ) and r 2 ( x ) , respectively, two positive distributions are Q ( x ) . ρ 1 ( x ) , ρ 2 ( x ) > 0 show the inter-specific competition ability. The habitat Ω is a bounded smooth domain in R N , 1 N Z ; n denotes the outward unit normal vector on the boundary Ω . For the sake of simplicity, we can suppose the initial data U 0 and V 0 not identically zero. The system (1) satisfies no-flux boundary conditions.
When κ 1 ( x ) = κ 1 , κ 2 ( x ) = κ 2 , μ 1 ( x ) = μ 1 , μ 2 ( x ) = μ 2 , Q ( x ) 1 , ρ 1 ( x ) = ρ 1 , ρ 2 ( x ) = ρ 2 , ω 1 ( x ) = ω 2 ( x ) , the system (1) becomes the advection system studied by Zhou and Xiao [19]:
U t = κ 1 Δ U μ 1 · [ U ω ( x ) ] + U [ r 1 ( x ) U ρ 2 V ] , i n   Ω × R + , V t = κ 2 Δ V μ 2 · [ V ω ( x ) ] + V [ r 2 ( x ) ρ 1 U V ] , i n   Ω × R + , κ 1 U n μ 1 U ω ( x ) n = 0 , o n   Ω × R + , κ 2 V n μ 2 V ω ( x ) n = 0 , o n   Ω × R + , U ( x , 0 ) = U 0 ( x ) , 0 , i n   Ω , V ( x , 0 ) = V 0 ( x ) , 0 , i n   Ω ,
where κ 1 , κ 2 , μ 1 , μ 2 , ρ 1 and ρ 2 are positive constants. Δ = i = 1 N 2 x i 2 is the usual Laplace operator. If μ 1 , μ 2 > 0 , readers can take a look at the relevant literature [20] and for the case μ 1 > 0 = μ 2 , please see the references [21,22,23,24,25].
If μ 1 = μ 2 = 0 , the system (2) becomes a diffusion model (see [2,3,5,26]):
U t = κ 1 Δ U + U [ r 1 ( x ) U V ] , i n   Ω × R + , V t = κ 2 Δ V + V [ r 2 ( x ) U V ] , i n   Ω × R + , U n = V n = 0 , o n   Ω × R + , U ( x , 0 ) = U 0 ( x ) , 0 , V ( x , 0 ) = V 0 ( x ) , 0 , i n   Ω .
According to the research of the above models, the purpose of our paper is to deal with a more broader model (1) in a high spatial dimensions. In this system, we consider that the diffusion and advection strategies follow the dispersal towards a positive distribution, growth rates and competitiveness of the two organisms are different. Thus, we have the following basic assumptions in this paper.
( A 1 )
μ 1 ( x ) κ 1 ( x ) ω 1 ( x ) μ 2 ( x ) κ 2 ( x ) ω 2 ( x ) : = ζ 1 ω 1 ( x ) ζ 2 ω 2 ( x ) 0 , Λ : = min x Ω ¯ e ζ 2 ω 2 ( x ) ζ 1 ω 1 ( x ) , ζ 1 and ζ 2 are positive constants;
( A 2 )
( ρ 1 ( x ) , ρ 2 ( x ) ) Π Λ : = { ( ρ 1 ( x ) , ρ 2 ( x ) ) : ρ 1 ( x ) , ρ 2 ( x ) > 0 , ρ 1 ( x ) ρ 2 ( x ) Λ } ;
( A 3 )
r 1 ( x ) > 0 , r 2 ( x ) > 0 in L ( Ω ) ;
( A 4 )
Q ( x ) > 0 is nonconstant, and Q ( x ) r 1 ( x ) , Q ( x ) r 2 ( x ) are also nonconstant.
Conditions ( A 3 ) ( A 4 ) ensure that the distribution of resources is heterogeneous for two species and the positivity is imposed here to guarantee the existence of two semi-trivial steady states for later discussion convenience. Under the conditions of ( A 1 ) ( A 4 ) , we show a complete classification of the global dynamics of the system (1). The rest of this paper is arranged as follows. In Section 2, we mainly do some preparatory work. Some related properties of the system (1) are deduced from the properties of a single organisms model (4). Besides, some lemmas are proved. In Section 3, we investigate our main results. By using principal eigenvalue theory, we obtain the linear stability of coexisting steady states (see Theorem 2). Then, the most important thing is that in virtue of the monotone dynamical system theory (see [4]), we show the classification of global dynamic behaviors (see Theorem 3). A discussion on the main results and problems that deserve future investigation is presented in Section 4.

2. Preliminaries

In order to describe our main results, we show a competition-diffusion-advection system for a single organisms as follows:
U t = · [ κ ( x ) ( U Q ( x ) ) μ ( x ) U Q ( x ) ω ( x ) ] + U [ r ( x ) U ] , i n   Ω × R + , κ ( x ) n ( U Q ) μ ( x ) U Q ω ( x ) n = 0 , o n   Ω × R + , U ( x , 0 ) = U 0 ( x ) , 0 , i n   Ω ,
where κ ( x ) > 0 , μ ( x ) > 0 , Q ( x ) > 0 and r ( x ) > 0 , r ( x ) is bounded. According to the relevant description in [27] and the case that r ( x ) > 0 , there is a unique positive steady state θ d , Q , μ , r in the system (4). If we apply this result to the system (1) and the conditions ( A 3 ) ( A 4 ) , there are two semi-trivial steady states ( θ κ 1 , Q , μ 1 , r 1 , 0 ) and ( 0 , θ κ 2 , Q , μ 2 , r 2 ) , respectively.
Lemma 1.
Assume that κ ( x ) > 0 , μ ( x ) > 0 , Q ( x ) > 0 and r ( x ) > 0 , r ( x ) is bounded. The elliptic boundary value Problem:
· κ ( x ) ( θ Q ( x ) ) μ ( x ) θ Q ( x ) ω ( x ) + θ [ r ( x ) θ ] = 0 , i n   Ω , κ ( x ) n ( θ Q ) μ ( x ) θ Q ω ( x ) n = 0 , o n   Ω ,
has a unique positive solution denoted by θ.
Proof. 
It is known in [27] that the problem (5) admits a solution and the solution is positive, denoted by θ , owning to the positivity of κ ( x ) , μ ( x ) , Q ( x ) , r ( x ) . Next, assume that θ 1 , θ 2 are any two positive solutions of (5) and 0 < θ 1 θ 2 . It is not difficult to see that
κ ( θ Q ) μ θ Q ω = κ e μ κ ω ( e μ κ ω θ Q ) .
Then
  · κ e μ κ ω ( e μ κ ω θ 1 Q ) ( e μ κ ω θ 2 Q )   d x = κ e μ κ ω ( e μ κ ω θ 1 Q ) ( e μ κ ω θ 2 Q )   d x = [ r θ 1 ] e μ κ ω θ 1 θ 2 Q   d x = [ r θ 2 ] e μ κ ω θ 1 θ 2 Q   d x .
We deduce
[ θ 1 θ 2 ] e μ κ ω θ 1 θ 2 Q   d x = 0 .
Therefore, θ 1 = θ 2 . □
To give a complete classification of the global dynamic system (1), we define
( κ 1 , κ 2 , μ 1 , μ 2 ) Γ : = R + × R + × R + × R + .
Based on the approach in [2], we define
Σ U : = { ( κ 1 , κ 2 , μ 1 , μ 2 ) Γ : ( θ κ 1 , Q , μ 1 , r 1 , 0 )   i s   l i n e a r l y   s t a b l e } ; Σ V : = { ( κ 1 , κ 2 , μ 1 , μ 2 ) Γ : ( 0 , θ κ 2 , Q , μ 2 , r 2 )   i s   l i n e a r l y   s t a b l e } ; Σ : = { ( κ 1 , κ 2 , μ 1 , μ 2 ) Γ : ( θ κ 1 , Q , μ 1 , r 1 , 0 )   a n d ( 0 , θ κ 2 , Q , μ 2 , r 2 )   a r e   l i n e a r l y   u n s t a b l e } .
We first recall the well-known Krein-Rutman Theorem:
Theorem 1
(Krein-Rutman Theorem [28]). Let X be a Banach space, K X a total cone and T : X X a compact linear operator that is positive (i.e., T ( K ) K ) with positive spectral radius r ( T ) . Then r ( T ) is an eigenvalue with an eigenvector u K \ 0 : T u = r ( T ) u . Moreover, r ( T * ) = r ( T ) is an eigenvalue of T * with an eigenvector u * K * .
In order to better describe the linear stability of semi-trivial steady states, we give the definition of elliptic eigenvalue problem:
· [ κ ( x ) ( ϕ Q ) μ ( x ) ϕ Q ω ( x ) ] + h ( x ) ϕ + σ ϕ = 0 , i n   Ω , κ ( x ) n ( ϕ Q ) μ ( x ) ϕ Q ω ( x ) n = 0 , o n   Ω ,
where κ ( x ) > 0 , μ ( x ) > 0 , Q ( x ) > 0 and h ( x ) L ( Ω ) . Let
A ϕ = · [ κ ( x ) ( ϕ Q ) μ ( x ) ϕ Q ω ( x ) ] + h ( x ) ϕ .
Since A is uniformly strongly elliptic operator, we declare that the operator A satisfies the conditions in Theorem 1. An eigenvalue σ 1 of the problem (8) is called a principal eigenvalue if σ 1 R and for any eigenvalue σ with σ σ 1 , we have R e   σ > σ 1 . Hence, the problem (8) has a principal eigenvalue, denoted by σ 1 ( κ , Q , μ , h ) , and its corresponding eigenfuntion ϕ ( κ , Q , μ , h ) > 0 in Ω . The principal eigenvalue is expressed as
σ 1 ( κ , Q , μ , h ) = inf 0 ϕ H 1 ( Ω ) κ e μ κ ω [ ( e μ κ ω ϕ Q ) ] 2   d x h · e μ κ ω ϕ 2 Q   d x e μ d ω ϕ 2 Q   d x .
Next, we give a useful lemma related to eigenvalue comparison results, which is used for Lemma 3 and Theorem 3.
Lemma 2
([5]). If h 1 ( x ) h 2 ( x ) within Ω, then σ 1 ( κ , Q , μ , h 1 ) σ 1 ( κ , Q , μ , h 2 ) and the equality holds if and only if h 1 ( x ) h 2 ( x ) in Ω.
According to the description of theory of monotone semi-flow in the literature [6], let X denote the standard Banach space consisting of all continuous functions from Ω ¯ to R , i.e., X : = C ( Ω ¯ ) , and X + be the set of all non-negative continuous functions from Ω ¯ to R + 0 . Define K : = X + × ( X + ) as the usual cone for the study of competitive systems with nonempty interior. Then we define the notion of linear stability of a given steady state ( U , V ) . Linearizing the steady state problem of (1) at ( U , V ) , we obtain
· [ κ 1 ( x ) ( φ Q ( x ) ) μ 1 ( x ) φ Q ( x ) ω 1 ( x ) ] + [ r 1 ( x ) U ρ 2 ( x ) V ] φ U [ φ + ρ 2 ( x ) ψ ] + λ φ = 0 , i n   Ω , · [ κ 2 ( x ) ( ψ Q ( x ) ) μ 2 ( x ) ψ Q ( x ) ω 2 ( x ) ] + [ r 2 ( x ) ρ 1 ( x ) U V ] ψ V [ ρ 1 ( x ) φ + ψ ] + λ ψ = 0 , i n   Ω , κ 1 ( x ) n ( φ Q ) μ 1 ( x ) φ Q ω 1 ( x ) n = 0 , o n   Ω , κ 2 ( x ) n ( ψ Q ) μ 2 ( x ) ψ Q ω 2 ( x ) n = 0 , o n   Ω .
Similar to the scalar problem (8), we can define the principal eigenvalue for the system (10), that is, an eigenvalue λ 1 of the problem (10) is called a principal eigenvalue if λ 1 R and for any eigenvalue λ with λ λ 1 , we have R e   λ > λ 1 . Based on the approach in [6], by using Theorem 1, the problem (10) has a principal eigenvalue λ 1 R . In fact, we can select the corresponding eigenfunction ( φ 1 , ψ 1 ) , which satisfies φ 1 > 0 > ψ 1 in Ω ¯ . Here, for the convenience of readers to better understand the problem (10), we provide a simple illustration. Let us do this simple transformation
Φ = e μ 1 κ 1 ω 1 φ   a n d   Ψ = e μ 2 κ 2 ω 2 ψ ,
then the problem (10) can be changed to
[ κ 1 ( x ) ( Φ Q ( x ) ) ] + μ 1 ( x ) ω 1 ( x ) · [ Φ Q ( x ) ] + [ r 1 ( x ) 2 U ρ 2 ( x ) V ] Φ + ρ 2 ( x ) U e ( ζ 2 ω 2 ( x ) ζ 1 ω 1 ( x ) ) Ψ + λ Φ = 0 , i n   Ω , [ κ 2 ( x ) ( Ψ Q ( x ) ) ] + μ 2 ( x ) ω 2 ( x ) · [ Ψ Q ( x ) ] + ρ 1 ( x ) V e ( ζ 1 ω 1 ( x ) ζ 2 ω 2 ( x ) ) Φ + [ r 2 ( x ) ρ 1 ( x ) U 2 V ] Ψ + λ Ψ = 0 , i n   Ω , n ( Φ Q ) = n ( Ψ Q ) = 0 , o n   Ω ,
which is a linear cooperative elliptic system. Suppose now L is the elliptic operator, let
L Φ = [ κ 1 ( x ) ( Φ Q ( x ) ) ] + μ 1 ( x ) ω 1 ( x ) · [ Φ Q ( x ) ] + [ r 1 ( x ) 2 U ρ 2 ( x ) V ] Φ , L Ψ = [ κ 2 ( x ) ( Ψ Q ( x ) ) ] + μ 2 ( x ) ω 2 ( x ) · [ Ψ Q ( x ) ] + [ r 2 ( x ) ρ 1 ( x ) U 2 V ] Ψ .
According to [28,29], the problem (11) has C α ( Ω ¯ ) coefficients and is strictly uniformly elliptic in the bounded domain Ω which has C 2 , α boundary. Let K be the positive cone in X : = C 0 1 , α ( Ω ¯ ) consisting of nonnegative functions. For any Φ 1 , Ψ 1 X , then we can deduce that T : X X defined by T ( Φ 1 , Ψ 1 ) = ( Φ , Ψ ) is a positive compact linear operator. By applying Theorem 1 for positive compact linear operators and the Neumann type boundary condition, the problem (11) admits a principal eigenvalue λ 1 R , and the corresponding eigenfunction ( Φ , Ψ ) can be chosen to satisfy Φ 0 and Ψ 0 in Ω ¯ . Notice that ( Φ , Ψ ) is the solution of the problem (11). Moreover, since the off-diagonal elements ρ 2 ( x ) U e ( ζ 2 ω 2 ( x ) ζ 1 ω 1 ( x ) ) and ρ 1 ( x ) V e ( ζ 1 ω 1 ( x ) ζ 2 ω 2 ( x ) ) are strictly positive in Ω , it can be further concluded that λ 1 is simple and it is the unique eigenvalue corresponding to a pair of strictly positive eigenfunctions, i.e., Φ > 0 and Ψ > 0 in Ω . In fact, we have Φ > 0 and Ψ > 0 in Ω ¯ due to Hopf boundary lemma, which in turn allows us to choose φ > 0 > ψ in Ω ¯ . See [30] using semi-group theory and [31] using maximum principle, [1,6] for detailed explanation. For the principal eigenvalue theory of general linear cooperative elliptic systems, we refer the interested readers to [29]. If λ λ 1 is an eigenvalue of (10) and the boundary condition is Neumann type, then R e   λ > λ 1 in the coexistence case.
Based on [26], (Corollary 2.10), the following lemma is about the linear stability of ( θ κ 1 , Q , μ 1 , r 1 , 0 ) and ( 0 , θ κ 2 , Q , μ 2 , r 2 ) .
Lemma 3.
The linear stability of ( θ κ 1 , Q , μ 1 , r 1 , 0 ) , ( 0 , θ κ 2 , Q , μ 2 , r 2 ) and ( 0 , 0 ) in the system (1) are determined by the sign of m i n σ 1 ( κ 1 , Q , μ 1 , r 1 ) , σ 1 ( κ 2 , Q , μ 2 , r 2 ) , σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) and σ 1 ( κ 1 , Q , μ 1 , r 1 ρ 2 θ κ 2 , Q , μ 2 , r 2 ) .
Proof. 
For the linear stability of ( θ κ 1 , Q , μ 1 , r 1 , 0 ) , when ( U , V ) = ( θ κ 1 , Q , μ 1 , r 1 , 0 ) in (10), we have
· [ κ 1 ( x ) ( φ Q ( x ) ) μ 1 ( x ) φ Q ( x ) ω 1 ( x ) ] + [ r 1 ( x ) 2 θ κ 1 , Q , μ 1 , r 1 ] φ + λ φ = θ κ 1 , Q , μ 1 , r 1 ρ 2 ( x ) ψ , i n   Ω , · [ κ 2 ( x ) ( ψ Q ( x ) ) μ 2 ( x ) ψ Q ( x ) ω 2 ( x ) ] + [ r 2 ( x ) ρ 1 ( x ) θ κ 1 , Q , μ 1 , r 1 ] ψ + λ ψ = 0 , i n   Ω , κ 1 ( x ) n ( φ Q ) μ 1 ( x ) φ Q ω 1 ( x ) n = 0 , o n   Ω , κ 2 ( x ) n ( ψ Q ) μ 2 ( x ) ψ Q ω 2 ( x ) n = 0 , o n   Ω .
Let λ be an principal eigenvalue of (12) with the eigenfunction ( φ , ψ ) . We get
λ = m i n σ 1 ( κ 1 , Q , μ 1 , r 1 2 θ κ 1 , Q , μ 1 , r 1 ) , σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) .
If ψ 0 , then λ belonging to an eigenvalue of the second equation in (12), is real and the inequality λ σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) holds. Perhaps, if ψ = 0 , then φ 0 and λ is an eigenvalue of the first equation, we get -4.6cm0cm
· [ κ 1 ( x ) ( φ Q ( x ) ) μ 1 ( x ) φ Q ( x ) ω 1 ( x ) ] + [ r 1 ( x ) 2 θ κ 1 , Q , μ 1 , r 1 ] φ + λ φ = 0 , i n   Ω , κ 1 ( x ) n ( φ Q ) μ 1 ( x ) φ Q ω 1 ( x ) n = 0 , o n   Ω .
Due to the fact that λ is real and satisfies λ σ 1 ( κ 1 , Q , μ 1 , r 1 2 θ κ 1 , Q , μ 1 , r 1 ) . It follows
λ m i n σ 1 ( κ 1 , Q , μ 1 , r 1 2 θ κ 1 , Q , μ 1 , r 1 ) , σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) .
If now σ 1 ( κ 1 , Q , μ 1 , r 1 2 θ κ 1 , Q , μ 1 , r 1 ) σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) , letting φ be the first eigenfunction corresponding to σ 1 ( κ 1 , Q , μ 1 , r 1 2 θ κ 1 , Q , μ 1 , r 1 ) , then σ 1 ( κ 1 , Q , μ 1 , r 1 2 θ κ 1 , Q , μ 1 , r 1 ) is an eigenvalue of (12) with the eigenfunction ( φ 1 , 0 ) , which deduces λ = σ 1 ( κ 1 , Q , μ 1 , r 1 2 θ κ 1 , Q , μ 1 , r 1 ) .
Suppose that σ 1 ( κ 1 , Q , μ 1 , r 1 2 θ κ 1 , Q , μ 1 , r 1 ) > σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) . Let ψ be the first eigenfunction corresponding to σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) , then σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) is an eigenvalue of (12) with the eigenfunction ( φ , ψ ) = ( φ * , ψ ) , that means λ = σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) . Here φ * satisfies
· [ κ 1 ( x ) ( φ * Q ( x ) ) μ 1 ( x ) φ * Q ( x ) ω 1 ( x ) ] + [ r 1 ( x ) 2 θ κ 1 , Q , μ 1 , r 1 ] φ * + σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) φ * = θ κ 1 , Q , μ 1 , r 1 ρ 2 ( x ) ψ , i n   Ω , κ 1 ( x ) n ( φ * Q ) μ 1 ( x ) φ * Q ω 1 ( x ) n = 0 , o n   Ω .
The existence of φ * is inferred from
σ 1 ( κ 1 , Q , μ 1 , r 1 2 θ κ 1 , Q , μ 1 , r 1 + σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) ) = σ 1 ( κ 1 , Q , μ 1 , r 1 2 θ κ 1 , Q , μ 1 , r 1 ) σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) > 0 .
So our claim is right. Owing to (6) and (9), it is inferred that σ 1 ( κ 1 , Q , μ 1 , r 1 θ κ 1 , Q , μ 1 , r 1 ) = 0 . Hence, according to Lemma 2, we gain
σ 1 ( κ 1 , Q , μ 1 , r 1 2 θ κ 1 , Q , μ 1 , r 1 ) > σ 1 ( κ 1 , Q , μ 1 , r 1 θ κ 1 , Q , μ 1 , r 1 ) = 0 ,
then λ has the same sign as the first eigenvalue σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) . Applying the definition of λ and linear stability, we deduce that the linear stability of ( θ κ 1 , Q , μ 1 , r 1 , 0 ) is determined by the sign of σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) .
Through completely similar arguments, we demonstrate that the stability of ( 0 , 0 ) and ( 0 , θ κ 2 , Q , μ 2 , r 2 ) , is determined by m i n σ 1 ( κ 1 , Q , μ 1 , r 1 ) , σ 1 ( κ 2 , Q , μ 2 , r 2 ) , σ 1 ( κ 1 , Q , μ 1 , r 1 ρ 2 θ κ 2 , Q , μ 2 , r 2 ) respectively. □
Remark 1.
From the variational characteristics of the first eigenvalue, we can see that ( 0 , 0 ) is linearly unstable for any κ 1 ( x ) , κ 2 ( x ) , μ 1 ( x ) , μ 2 ( x ) , ρ 2 ( x ) , ρ 1 ( x ) > 0 .
Therefore, we give equivalent descriptions of (7) below:
Σ U : = { ( κ 1 , κ 2 , μ 1 , μ 2 ) Γ : σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) > 0 } ; Σ V : = { ( κ 1 , κ 2 , μ 1 , μ 2 ) Γ : σ 1 ( κ 1 , Q , μ 1 , r 1 ρ 2 θ κ 2 , Q , μ 2 , r 2 ) > 0 } ; Σ : = { ( κ 1 , κ 2 , μ 1 , μ 2 ) Γ : σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) < 0   a n d σ 1 ( κ 1 , Q , μ 1 , r 1 ρ 2 θ κ 2 , Q , μ 2 , r 2 ) < 0 } .
The neutrally stable case is defined as follows
Σ U , 0 : = { ( κ 1 , κ 2 , μ 1 , μ 2 ) Γ : σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) = 0 } ; Σ V , 0 : = { ( κ 1 , κ 2 , μ 1 , μ 2 ) Γ : σ 1 ( κ 1 , Q , μ 1 , r 1 ρ 2 θ κ 2 , Q , μ 2 , r 2 ) = 0 } ; Σ 0 , 0 : = { ( κ 1 , κ 2 , μ 1 , μ 2 ) Γ : σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) = σ 1 ( κ 1 , Q , μ 1 , r 1 ρ 2 θ κ 2 , Q , μ 2 , r 2 ) = 0 } .
By the definition, it is easy to see Σ 0 , 0 = Σ U , 0 Σ V , 0 .
In the following, “ g . a . s ” is used to mean that the steady state is globally asymptotically stable among all non-negative and not identically zero initial conditions.
Lemma 4
([5]). For any κ 1 , κ 2 , μ 1 , μ 2 > 0 , assume that ( A 1 ) ( A 4 ) hold and every coexistence steady state of the system (1), if it exists, is asymptotically stable. Then one of the following alternatives holds:
(i)
There exists a unique coexistence steady state of (1) that is g . a . s .
(ii)
The system (1) has no coexistence steady state and either one of ( θ κ 1 , Q , μ 1 , r 1 , 0 ) or ( 0 , θ κ 2 , Q , μ 2 , r 2 ) is g . a . s , while the other is unstable.

3. Main Results

In this section, we present the results which are related to the co-existence steady state and the classification of global dynamic behaviors of the system (1).
Theorem 2.
Suppose that ( A 1 ) ( A 4 ) hold. For any ( κ 1 , κ 2 , μ 1 , μ 2 ) Γ \ Σ 0 , 0 , then every co-existence steady state of the system ( 1 ) , if exists, is linearly stable, i.e., λ 1 > 0 .
Theorem 3.
Suppose that ( A 1 ) ( A 4 ) hold. Then we have the mutually disjoint decomposition of Γ:
Γ = ( Σ U Σ U , 0 \ Σ 0 , 0 ) ( Σ V Σ V , 0 \ Σ 0 , 0 ) Σ Σ 0 , 0 .
Moreover, the following statements hold for the system ( 1 ) :
( i )
For all ( κ 1 , κ 2 , μ 1 , μ 2 ) ( Σ U Σ U , 0 \ Σ 0 , 0 ) , ( θ κ 1 , Q , μ 1 , r 1 , 0 ) is g . a . s ;
( i i )
For all ( κ 1 , κ 2 , μ 1 , μ 2 ) ( Σ V Σ V , 0 \ Σ 0 , 0 ) , ( 0 , θ κ 2 , Q , μ 2 , r 2 ) is g . a . s ;
( i i i )
For all ( κ 1 , κ 2 , μ 1 , μ 2 ) Σ , the system ( 1 ) has a unique coexistence steady state that is g . a . s ;
( i v )
For all ( κ 1 , κ 2 , μ 1 , μ 2 ) Σ 0 , 0 , θ κ 1 , Q , μ 1 , r 1 ρ 2 ( x ) θ κ 2 , Q , μ 2 , r 2 in Ω ¯ and the system ( 1 ) has a compact global attractor consisting of a continuum of steady states
{ η ( x ) θ κ 1 , Q , μ 1 , r 1 , ( 1 η ( x ) ) θ κ 1 , Q , μ 1 , r 1 ρ 2 ( x ) : η ( x ) [ 0 , 1 ] }
connecting with two semi-trivial steady states.

3.1. Co-Existence Steady State

In order to prove Theorem 2, we assume that ( U , V ) is the co-existence steady state of the following system (1):
· [ κ 1 ( x ) ( U Q ( x ) ) μ 1 ( x ) U Q ( x ) ω 1 ( x ) ] + U [ r 1 ( x ) U ρ 2 ( x ) V ] = 0 , i n   Ω , · [ κ 2 ( x ) ( V Q ( x ) ) μ 2 ( x ) V Q ( x ) ω 2 ( x ) ] + V [ r 2 ( x ) ρ 1 ( x ) U V ] = 0 , i n   Ω , κ 1 ( x ) n ( U Q ) μ 1 ( x ) U Q ω 1 ( x ) n = 0 , o n   Ω , κ 2 ( x ) n ( V Q ) μ 2 ( x ) V Q ω 2 ( x ) n = 0 , o n   Ω .
Similar to the problem (10), then we get the linear eigenvalue model by linearize system (1) at ( U , V ) ,
· [ κ 1 ( x ) ( φ Q ( x ) ) μ 1 ( x ) φ Q ( x ) ω 1 ( x ) ] + [ r 1 ( x ) U ρ 2 ( x ) V ] φ U [ φ + ρ 2 ( x ) ψ ] + λ φ = 0 , i n   Ω , · [ κ 2 ( x ) ( ψ Q ( x ) ) μ 2 ( x ) ψ Q ( x ) ω 2 ( x ) ] + [ r 2 ( x ) ρ 1 ( x ) U V ] ψ V [ ρ 1 ( x ) φ + ψ ] + λ ψ = 0 , i n   Ω , κ 1 ( x ) n ( φ Q ) μ 1 ( x ) φ Q ω 1 ( x ) n = 0 , o n   Ω , κ 2 ( x ) n ( ψ Q ) μ 2 ( x ) ψ Q ω 2 ( x ) n = 0 , o n   Ω .
According to the problem (8) and using Theorem 1, we can deduce that the problem (18) has a principal eigenvalue λ 1 . Moreover, we can choose the corresponding eigenfunction ( φ , ψ ) , it satisfies φ > 0 > ψ in Ω ¯ .
Now, we are ready to discuss Theorem 2.
Proof of Theorem 2.
Obviously, as long as we can obtain λ 1 > 0 when ( κ 1 , κ 2 , μ 1 , μ 2 ) Γ \ Σ 0 , 0 . Multiplying the first equation in ( 18 ) and ( 17 ) by U Q ( x ) and φ Q ( x ) , respectively, and subtracting the obtained equations, we obtain
κ 1 ( φ Q ) μ 1 φ Q ω 1 U Q κ 1 ( U Q ) μ 1 U Q ω 1 φ Q U 2 Q [ φ + ρ 2 ( x ) ψ ] = λ 1 U φ Q .
In the similar way, it can be derived from the second equation in ( 17 ) and ( 18 ) that
κ 2 ( ψ Q ) μ 2 ψ Q ω 2 V Q κ 2 ( V Q ) μ 2 V Q ω 2 ψ Q V 2 Q [ ρ 1 ( x ) φ + ψ ] = λ 1 V ψ Q .
Furthermore, multiplying ( 19 ) by e μ 1 κ 1 ω 1 · φ 2 U 2 , then we integrate over Ω and deduce (for simplicity, we replace Ω with ∫)
λ 1 e μ 1 κ 1 ω 1 φ 3 U Q   d x = κ 1 ( φ Q ) μ 1 φ Q ω 1 e μ 1 κ 1 ω 1 φ 2 U Q   d x + κ 1 ( U Q ) μ 1 U Q ω 1 e μ 1 κ 1 ω 1 φ 3 U 2 Q   d x + e μ 1 κ 1 ω 1 φ 2 Q [ φ + ρ 2 ( x ) ψ ]   d x = κ 1 ( φ Q ) μ 1 φ Q ω 1 · ( e μ 1 κ 1 ω 1 φ 2 U Q )   d x κ 1 ( U Q ) μ 1 U Q ω 1 · ( e μ 1 κ 1 ω 1 φ 3 U 2 Q )   d x + e μ 1 κ 1 ω 1 φ 2 Q [ φ + ρ 2 ( x ) ψ ]   d x : = I 1 I 2 + I 3 .
By using the similarly method for the Equation (20), we get
  λ 1 e μ 2 κ 2 ω 2 ψ 3 V Q   d x = κ 2 ( ψ Q ) μ 2 ψ Q ω 2 e μ 2 κ 2 ω 2 ψ 2 V Q   d x   + κ 2 ( V Q ) μ 2 V Q ω 2 e μ 2 κ 2 ω 2 ψ 3 V 2 Q   d x   + e μ 2 κ 2 ω 2 ψ 2 Q [ ρ 1 ( x ) φ + ψ ]   d x = κ 2 ( ψ Q ) μ 2 ψ Q ω 2 · ( e μ 2 κ 2 ω 2 ψ 2 V Q )   d x   κ 2 ( V Q ) μ 2 V Q ω 2 · ( e μ 2 κ 2 ω 2 ψ 3 V 2 Q )   d x   + e μ 2 κ 2 ω 2 ψ 2 Q [ ρ 1 ( x ) φ + ψ ]   d x : = J 1 J 2 + J 3 .
We now simplify the formulas I 1 , I 2 , J 1 and J 2 . Then we find
I 1 : = κ 1 ( φ Q ) μ 1 φ Q ω 1 · ( e μ 1 κ 1 ω 1 φ 2 U Q )   d x = e μ 1 κ 1 ω 1 [ κ 1 ( φ Q ) μ 1 φ Q ω 1 ] [ 2 φ U φ φ 2 U U 2 Q ]   d x   μ 1 κ 1 φ 2 U Q e μ 1 κ 1 ω 1 [ κ 1 ( φ Q ) μ 1 φ Q ω 1 ] · [ ω 1 ]   d x   φ 2 U Q 2 e μ 1 κ 1 ω 1 [ κ 1 ( φ Q ) μ 1 φ Q ω 1 ] [ Q ]   d x = e μ 1 κ 1 ω 1 φ 3 U Q 2 [ κ 1 ( φ Q ) Q φ μ 1 ω 1 ] [ 2 φ φ U U ]   d x   μ 1 κ 1 φ 3 U Q 2 e μ 1 κ 1 ω 1 [ κ 1 ( φ Q ) Q φ μ 1 ω 1 ] · [ ω 1 ]   d x   φ 3 U Q 3 e μ 1 κ 1 ω 1 [ κ 1 ( φ Q ) Q φ μ 1 ω 1 ] [ Q ]   d x : = Δ 1 Δ 2 Δ 3 ,
and
I 2 : = κ 1 ( U Q ) μ 1 U Q ω 1 · ( e μ 1 κ 1 ω 1 φ 3 U 2 Q )   d x = e μ 1 κ 1 ω 1 [ κ 1 ( U Q ) μ 1 U Q ω 1 ] [ 3 φ 2 U 2 φ φ 3 2 U U U 4 Q ]   d x   μ 1 κ 1 φ 3 U 2 Q e μ 1 κ 1 ω 1 [ κ 1 ( U Q ) μ 1 U Q ω 1 ] · [ ω 1 ]   d x   φ 3 U 2 Q 2 e μ 1 κ 1 ω 1 [ κ 1 ( U Q ) μ 1 U Q ω 1 ] [ Q ]   d x
= e μ 1 κ 1 ω 1 φ 3 U Q 2 [ κ 1 ( U Q ) Q U μ 1 ω 1 ] [ 3 φ φ 2 U U ]   d x   μ 1 κ 1 φ 3 U Q 2 e μ 1 κ 1 ω 1 [ κ 1 ( U Q ) Q U μ 1 ω 1 ] · [ ω 1 ]   d x   φ 3 U Q 3 e μ 1 κ 1 ω 1 [ κ 1 ( U Q ) Q U μ 1 ω 1 ] [ Q ]   d x : = Δ 4 Δ 5 Δ 6 .
Next, we have
  Δ 1 Δ 4 = e μ 1 κ 1 ω 1 φ 3 U Q 2 [ κ 1 ( φ Q ) Q φ μ 1 ω 1 ] [ 2 φ φ U U ]   d x   e μ 1 κ 1 ω 1 φ 3 U Q 2 [ κ 1 ( U Q ) Q U μ 1 ω 1 ] [ 3 φ φ 2 U U ]   d x = κ 1 e μ 1 κ 1 ω 1 φ 3 U Q 2 { [ ( φ Q ) Q φ ] [ 2 φ φ U U ] [ ( U Q ) Q U ] [ 3 φ φ 2 U U ] }   d x   + e μ 1 κ 1 ω 1 φ 3 U Q 2 [ μ 1 ω 1 ] [ φ φ U U ]   d x = κ 1 e μ 1 κ 1 ω 1 φ 3 U Q 2 { [ φ Q φ φ Q 2 ] Q φ [ 2 φ φ U U ] [ U Q U Q Q 2 ] Q U [ 3 φ φ 2 U U ] }   d x + e μ 1 κ 1 ω 1 φ 3 U Q 2 [ μ 1 ω 1 ] [ φ φ U U ]   d x = 2 κ 1 e μ 1 κ 1 ω 1 φ 3 U Q 2 ( φ φ U U ) 2   d x   + κ 1 e μ 1 κ 1 ω 1 φ 3 U Q 3 [ φ φ U U ] [ Q ]   d x   + e μ 1 κ 1 ω 1 φ 3 U Q 2 [ μ 1 ω 1 ] [ φ φ U U ]   d x .
By the similar method, we deduce that
Δ 2 Δ 5 = e μ 1 κ 1 ω 1 φ 3 U Q 2 [ μ 1 ω 1 ] [ φ φ U U ]   d x ,
and
Δ 3 Δ 6 = κ 1 e μ 1 κ 1 ω 1 φ 3 U Q 3 [ φ φ U U ] [ Q ]   d x .
Thus
I 1 I 2 = [ Δ 1 Δ 4 ] [ Δ 2 Δ 5 ] [ Δ 3 Δ 6 ] = 2 κ 1 e μ 1 κ 1 ω 1 φ 3 U Q 2 ( φ φ U U ) 2   d x .
By a similar method, one obtains
J 1 J 2 = 2 κ 2 e μ 2 κ 2 ω 2 ψ 3 V Q 2 ( ψ ψ V V ) 2   d x .
Replace (21) and (22) with (26) and (27), respectively. Multiplying (22) by ρ 2 ( x ) 3 and subtracting it from (21), we can obtain
  λ 1 [ e μ 1 κ 1 ω 1 φ 3 U Q e μ 2 κ 2 ω 2 ρ 2 ( x ) 3 ψ 3 V Q ]   d x = 2 κ 1 e μ 1 κ 1 ω 1 φ 3 U Q 2 ( φ φ U U ) 2   d x + e μ 1 κ 1 ω 1 φ 2 Q [ φ + ρ 2 ( x ) ψ ]   d x   2 κ 2 e μ 2 κ 2 ω 2 ρ 2 ( x ) 3 ψ 3 V Q 2 ( ψ ψ V V ) 2 d x   e μ 2 κ 2 ω 2 ρ 2 ( x ) 2 ψ 2 Q [ ρ 1 ( x ) ρ 2 ( x ) φ + ρ 2 ( x ) ψ ]   d x e μ 1 κ 1 ω 1 φ 2 Q [ φ + ρ 2 ( x ) ψ ]   d x e μ 2 κ 2 ω 2 ρ 2 ( x ) 3 ψ 3 Q   d x   e μ 1 κ 1 ω 1 ρ 2 ( x ) 2 ψ 2 φ Q [ ρ 1 ( x ) ρ 2 ( x ) e ( μ 1 κ 1 ω 1 μ 2 κ 2 ω 2 ) ]   d x e μ 1 κ 1 ω 1 φ 2 Q [ φ + ρ 2 ( x ) ψ ]   d x e μ 1 κ 1 ω 1 ρ 2 ( x ) 3 ψ 3 Q   d x   e μ 1 κ 1 ω 1 ρ 2 ( x ) 2 ψ 2 φ Q   d x = e μ 1 κ 1 ω 1 1 Q [ φ + ρ 2 ( x ) ψ ] 2 [ φ ρ 2 ( x ) ψ ]   d x 0 ,
where we use the facts ( A 1 ) , ( A 2 ) and φ > 0 > ψ in Ω ¯ . So, λ 1 0 .
Next, we will prove that λ 1 = 0 can not happen. According to (28), we infer that λ 1 = 0 if and only if
ρ 1 ( x ) ρ 2 ( x ) = 1 ,   μ 1 κ 1 ω 1 = μ 2 κ 2 ω 2 ,   φ φ = U U ,   ψ ψ = V V ,   φ = ρ 2 ( x ) ψ ,
which means that
U U = V V ,
i.e.,
[ l n U ] = [ l n V ] .
Then, one obtains
U = d V for   some   constant d > 0 .
In addition, by applying (30) to (17), and the uniqueness of the positive steady state of the system (4), it can be concluded that
( 1 + ρ 2 ( x ) d ) U = θ κ 1 , Q , μ 1 , r 1 and ( ρ 1 ( x ) d + 1 ) V = θ κ 2 , Q , μ 2 , r 2 .
Noting that ρ 1 ( x ) ρ 2 ( x ) = 1 , we deduce
θ κ 1 , Q , μ 1 , r 1 θ κ 2 , Q , μ 2 , r 2 = ( 1 + ρ 2 ( x ) d ) U ( ρ 1 ( x ) d + 1 ) V = U + ρ 2 ( x ) V ρ 1 ( x ) U + V = ρ 2 ( x ) = 1 ρ 1 ( x ) .
Based on (31), one can easily check
σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) = σ 1 ( κ 2 , Q , μ 2 , r 2 θ κ 2 , Q , μ 2 , r 2 ) = 0 ,
and
σ 1 ( κ 1 , Q , μ 1 , r 1 θ κ 1 , Q , μ 1 , r 1 ) = σ 1 ( κ 1 , Q , μ 1 , r 1 ρ 2 θ κ 2 , Q , μ 2 , r 2 ) = 0 .
According to the assumption ( κ 1 , κ 2 , μ 1 , μ 2 ) Σ 0 , 0 , we get λ 1 0 . Therefore, λ 1 > 0 . □

3.2. Classification of the Global Dynamics

In this subsection, we are ready to discuss the complete global dynamic behaviors of the system (1).
Proof of Theorem 3.
According to the proof procedure (see [2,5]), this proof can be divided into two steps.
Step 1. On the proof of the disjoint decomposition in (16).
Obviously, we can get the decomposition in (16). According to the relevant conclusion and definitions, we only need to prove
( Σ U Σ U , 0 \ Σ 0 , 0 ) ( Σ V Σ V , 0 \ Σ 0 , 0 ) = .
By Lemma 3, the linear stability of ( θ κ 1 , Q , μ 1 , r 1 , 0 ) , ( 0 , θ κ 2 , Q , μ 2 , r 2 ) can be determined by the sign of σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) , σ 1 ( κ 1 , Q , μ 1 , r 1 ρ 2 θ κ 2 , Q , μ 2 , r 2 ) respectively. For the sake of convenience of in writing, let
θ κ 1 , Q , μ 1 , r 1 : = θ 1 ,   σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) : = σ 1 ( θ 1 * ) , θ κ 2 , Q , μ 2 , r 2 : = θ 2 ,   σ 1 ( κ 1 , Q , μ 1 , r 1 ρ 2 θ κ 2 , Q , μ 2 , r 2 ) : = σ 1 ( θ 2 * ) .
According to the properties of the variational characterization and (6), we obtain
σ 1 ( θ 2 * ) = inf 0 ϕ H 1 ( Ω ) κ 1 e μ 1 κ 1 ω 1 [ ( e μ 1 κ 1 ω 1 ϕ Q ) ] 2   d x [ r 1 ρ 2 ( x ) θ 2 ] e μ 1 κ 1 ω 1 ϕ 2 Q   d x e μ 1 κ 1 ω 1 ϕ 2 Q   d x κ 1 e μ 1 κ 1 ω 1 [ ( e μ 1 κ 1 ω 1 θ 1 Q ) ] 2   d x [ r 1 ρ 2 ( x ) θ 2 ] e μ 1 κ 1 ω 1 θ 1 2 Q   d x e μ 1 κ 1 ω 1 θ 1 2 Q   d x = e μ 1 κ 1 ω 1 θ 1 2 Q [ ρ 2 ( x ) θ 2 θ 1 ]   d x e μ 1 κ 1 ω 1 θ 1 2 Q   d x ,
and
σ 1 ( θ 1 * ) = inf 0 ϕ H 1 ( Ω ) κ 2 e μ 2 κ 2 ω 2 [ ( e μ 2 κ 2 ω 2 ϕ Q ) ] 2   d x [ r 2 ρ 1 ( x ) θ 1 ] e μ 2 κ 2 ω 2 ϕ 2 Q   d x e μ 2 κ 2 ω 2 ϕ 2 Q   d x κ 2 e μ 2 κ 2 ω 2 [ ( e μ 2 κ 2 ω 2 θ 2 Q ) ] 2   d x [ r 2 ρ 1 ( x ) θ 1 ] e μ 2 κ 2 ω 2 θ 2 2 Q d x e μ 2 κ 2 ω 2 θ 2 2 Q   d x = e μ 2 κ 2 ω 2 θ 2 2 Q [ ρ 1 ( x ) θ 1 θ 2 ]   d x e μ 2 κ 2 ω 2 θ 2 2 Q   d x .
Since 0 < ρ 1 ( x ) ρ 2 ( x ) 1 , combining with (33) and (34) together, we have
  σ 1 ( θ 2 * ) · e μ 1 κ 1 ω 1 θ 1 2 Q   d x + σ 1 ( θ 1 * ) · e μ 2 κ 2 ω 2 ρ 2 ( x ) 3 θ 2 2 Q   d x e μ 1 κ 1 ω 1 θ 1 2 Q [ ρ 2 ( x ) θ 2 θ 1 ]   d x e μ 2 κ 2 ω 2 ρ 2 ( x ) 3 θ 2 3 Q   d x   + e μ 1 κ 1 ω 1 ρ 2 ( x ) 2 θ 2 2 θ 1 Q [ ρ 1 ( x ) ρ 2 ( x ) e ( μ 1 κ 1 ω 1 μ 2 κ 2 ω 2 ) ]   d x e μ 1 κ 1 ω 1 θ 1 2 Q [ ρ 2 ( x ) θ 2 θ 1 ]   d x e μ 1 κ 1 ω 1 ρ 2 ( x ) 3 θ 2 3 Q   d x   + e μ 1 κ 1 ω 1 ρ 2 ( x ) 2 θ 2 2 θ 1 Q   d x = e μ 1 κ 1 ω 1 1 Q [ ρ 2 ( x ) θ 2 θ 1 ] 2 [ θ 1 + ρ 2 ( x ) θ 2 ]   d x 0 ,
where all the inequalities become equalities if and only if
ρ 1 ( x ) ρ 2 ( x ) = 1 ,   μ 1 κ 1 ω 1 = μ 2 κ 2 ω 2   a n d   θ κ 1 , Q , μ 1 , r 1 ρ 2 ( x ) θ κ 2 , Q , μ 2 , r 2   i n   Ω ¯ .
It follows from (35) that the conclusion (32) holds.
Step 2. On the proof of the statements ( i ) ( i v ) .
Firstly, we will prove the statements ( i ) ( i i i ) hold. In consideration of (16) in Theorem 3 and (35) in step 1, we see that for any ( κ 1 , κ 2 , μ 1 , μ 2 ) ( Γ \ Σ 0 , 0 ) , there are five possibilities as follows:
( b 1 )
( κ 1 , κ 2 , μ 1 , μ 2 ) Σ U , i.e., ( θ κ 1 , Q , μ 1 , r 1 , 0 ) is linearly stable, ( 0 , θ κ 2 , Q , μ 2 , r 2 ) is linearly unstable;
( b 2 )
( κ 1 , κ 2 , μ 1 , μ 2 ) Σ V , i.e., ( θ κ 1 , Q , μ 1 , r 1 , 0 ) is linearly unstable, ( 0 , θ κ 2 , Q , μ 2 , r 2 ) is linearly stable;
( b 3 )
( κ 1 , κ 2 , μ 1 , μ 2 ) Σ , i.e., both ( θ κ 1 , Q , μ 1 , r 1 , 0 ) and ( 0 , θ κ 2 , Q , μ 2 , r 2 ) are linearly unstable;
( b 4 )
( κ 1 , κ 2 , μ 1 , μ 2 ) Σ U , 0 \ Σ 0 , 0 , i.e., ( θ κ 1 , Q , μ 1 , r 1 , 0 ) is neutrally stable, ( 0 , θ κ 2 , Q , μ 2 , r 2 ) is linearly unstable;
( b 5 )
( κ 1 , κ 2 , μ 1 , μ 2 ) Σ V , 0 \ Σ 0 , 0 , i.e., ( 0 , θ κ 2 , Q , μ 2 , r 2 ) is neutrally stable, ( θ κ 1 , Q , μ 1 , r 1 , 0 ) is linearly unstable.
By Lemma 4, we immediately deduce the following conclusion:
( θ κ 1 , Q , μ 1 , r 1 , 0 ) and ( 0 , θ κ 2 , Q , μ 2 , r 2 ) are g . a . s based on the assumptions ( b 1 ) and ( b 2 ) , respectively, and there is a unique co-existence steady state under the condition ( b 3 ) .
We now claim that there is no coexistence steady state under the condition ( b 4 ) or the condition ( b 5 ) . Then we can infer that ( θ κ 1 , Q , μ 1 , r 1 , 0 ) and ( 0 , θ κ 2 , Q , μ 2 , r 2 ) are also g . a . s based on the assumptions ( b 4 ) and ( b 5 ) , respectively, from Lemma 4.
We only need to verify the above statement for the case ( b 4 ). Indeed, if the system (1) has a co-existence steady state ( U ˜ , V ˜ ) for some ( κ 1 , κ 2 , μ 1 , μ 2 ) = ( κ 1 ˜ , κ 2 ˜ , μ 1 ˜ , μ 2 ˜ ) Σ U , 0 \ Σ 0 , 0 and ( U ˜ , V ˜ ) satisfies
· [ κ 1 ˜ ( x ) ( U ˜ Q ˜ ( x ) ) μ 1 ˜ ( x ) U ˜ Q ˜ ( x ) ω 1 ( x ) ] + U ˜ [ r 1 ( x ) U ˜ ρ 2 ( x ) V ˜ ] = 0 , i n   Ω , · [ κ 2 ˜ ( x ) ( V ˜ Q ˜ ( x ) ) μ 2 ˜ ( x ) V ˜ Q ˜ ( x ) ω 2 ( x ) ] + V ˜ [ r 2 ( x ) ρ 1 ( x ) U ˜ V ˜ ] = 0 , i n   Ω , [ κ 1 ˜ ( x ) n ( U ˜ Q ˜ ) μ 1 ˜ ( x ) U ˜ Q ˜ ω 1 ( x ) n ] Ω = [ κ 2 ˜ ( x ) n ( V ˜ Q ˜ ) μ 2 ˜ ( x ) V ˜ Q ˜ ω 2 ( x ) n ] Ω = 0 .
We have
σ 1 ( κ 2 ˜ , Q ˜ , μ 2 ˜ , r 2 ρ 1 θ κ 1 ˜ , Q ˜ , μ 1 ˜ , r 1 ) = 0   a n d   σ 1 ( κ 1 ˜ , Q ˜ , μ 1 ˜ , r 1 ρ 2 θ κ 2 ˜ , Q ˜ , μ 2 ˜ , r 2 ) < 0 .
Define the operator G : R + × R + × M 1 × M 1 M 2 × M 2
G : ( ϱ 1 , ϱ 2 , w 1 , w 2 ) κ 1 ˜ Δ ( w 1 Q ) μ 1 ˜ · [ w 1 Q ω 1 ] + w 1 [ r 1 w 1 ϱ 2 w 2 ] κ 2 ˜ Δ ( w 2 Q ) μ 2 ˜ · [ w 2 Q ω 2 ] + w 2 [ r 2 ϱ 1 w 1 w 2 ] ,
with
M 1 : = { φ W 2 , p ( Ω ) : [ κ 1 ˜ n ( φ Q ) μ 1 ˜ ( φ Q ) ω 1 Q ] Ω = 0 } ,   M 2 : = L p ( Ω ) , p > N .
From (36), G ( ρ 1 ( x ) , ρ 2 ( x ) , U ˜ , V ˜ ) = 0 and Theorem 2, it yields that
G ( ϱ 1 , ϱ 2 , w 1 , w 2 ) ( w 1 , w 2 ) ( ϱ 1 , ϱ 2 , w 1 , w 2 ) = ( ρ 1 ( x ) , ρ 2 ( x ) , U ˜ , V ˜ )   i s   i n v e r t i b l e .
Applying implicit function theorem, one gets ( ϱ 1 , ϱ 2 ) is closed to ( ρ 1 ( x ) , ρ 2 ( x ) ) . We have a positive solution ( U ˜ , V ˜ ) to the equation G ( ϱ 1 , ϱ 2 , w 1 , w 2 ) = 0 . Let us choose ( ϱ 1 ¯ , ϱ 2 ¯ ) , which implies (the solution corresponding to ( ϱ 1 ¯ , ϱ 2 ¯ ) is denoted by ( U ¯ , V ¯ ) )
ϱ 1 ¯ > ρ 1 ( x ) ,   0 < ϱ 2 ¯ < ρ 2 ( x )   a n d   ϱ 1 ¯ · ϱ 2 ¯ Λ .
Let us see the following auxiliary problem
U t = · [ κ 1 ˜ ( x ) ( U Q ) μ 1 ˜ ( x ) U Q ω 1 ( x ) ] + U [ r 1 ( x ) U ϱ 2 ¯ V ] = 0 , i n   Ω × R + , V t = · [ κ 2 ˜ ( x ) ( V Q ) μ 2 ˜ ( x ) V Q ω 2 ( x ) ] + V [ r 2 ( x ) ϱ 1 ¯ U V ] = 0 , i n   Ω × R + , [ κ 1 ˜ ( x ) n ( U Q ) μ 1 ˜ ( x ) U Q ω 1 ( x ) n ] Ω = 0 , [ κ 2 ˜ ( x ) n ( V Q ) μ 2 ˜ ( x ) V Q ω 2 ( x ) n ] Ω = 0 ,
which has the same semi-trivial steady states ( θ κ 1 ˜ , Q ˜ , μ 1 ˜ , r 1 , 0 ) and ( 0 , θ κ 2 ˜ , Q ˜ , μ 2 ˜ , r 2 ) . From (37), (38) and Lemma 2, it then follows that
σ 1 ( κ 2 ˜ , Q ˜ , μ 2 ˜ , r 2 ϱ 1 ¯ θ κ 1 ˜ , Q ˜ , μ 1 ˜ , r 1 ) > 0   a n d   σ 1 ( κ 1 ˜ , Q ˜ , μ 1 ˜ , r 1 ϱ 2 ¯ θ κ 2 ˜ , Q ˜ , μ 2 ˜ , r 2 ) < 0 .
According to the case ( b 1 ) , ( θ κ 1 , Q , μ 1 , r 1 , 0 ) is also g . a . s in the system (39) which contradicts with the existence of ( U ¯ , V ¯ ) . Therefore, there is no coexistence steady state under the condition ( b 4 ) . Similarly, we can get the conclusion that there is also no coexistence steady state under the condition ( b 5 ) . The above descriptions of the cases ( b 1 ) ( b 5 ) represent the expected results described in the statements ( i ) ( i i i ) .
Secondly, we prove the statement ( i v ) . We will show
Σ 0 , 0 = Σ : = { ( κ 1 , κ 2 , μ 1 , μ 2 ) Γ : ρ 1 ( x ) ρ 2 ( x ) = 1 , θ κ 1 , Q , μ 1 , r 1 ρ 2 ( x ) θ κ 2 , Q , μ 2 , r 2   i n   Ω ¯ } .
It makes the same description of Σ 0 , 0 , which means the expected result in the statement ( i v ) .
Let ( κ 1 , κ 2 , μ 1 , μ 2 ) Σ , then
ρ 1 ( x ) ρ 2 ( x ) = 1   a n d   θ κ 1 , Q , μ 1 , r 1 ρ 2 ( x ) θ κ 2 , Q , μ 2 , r 2   i n   Ω ¯ .
Based on the proof of Theorem 2, we get
σ 1 ( κ 2 , Q , μ 2 , r 2 ρ 1 θ κ 1 , Q , μ 1 , r 1 ) = σ 1 ( κ 1 , Q , μ 1 , r 1 ρ 2 θ κ 2 , Q , μ 2 , r 2 ) = 0 ,
which implies Σ Σ 0 , 0 . When (42) holds, the last three inequalities in (35) become equalities, we have
ρ 1 ( x ) ρ 2 ( x ) = 1 , θ κ 1 , Q , μ 1 , r 1 ρ 2 ( x ) θ κ 2 , Q , μ 2 , r 2   i n   Ω ¯ ,
which shows Σ 0 , 0 Σ . Hence, the equality (41) is confirmed.
Let ( κ 1 , κ 2 , μ 1 , μ 2 ) Σ 0 , 0 and ( U , V ) be the corresponding coexistence steady state of (1). We claim that U V constant. Let λ 1 be a principal eigenvalue for ( U , V ) . Moreover, we choose the corresponding eigenfunction ( φ , ψ ) , which satisfies φ > 0 > ψ in Ω ¯ and φ 2 2 + ψ 2 2 = 1 . In order to prove it, it is enough to show that (29) holds. Suppose that (29) is not true. Then (28) means λ 1 > 0 . Similar to the proof of the case ( b 4 ) , we get (29) holds, i.e., U V constant. This yields that
( U , V ) = η ( x ) θ κ 1 , Q , μ 1 , r 1 , ( 1 η ( x ) ) θ κ 1 , Q , μ 1 , r 1 ρ 2 ( x ) : η ( x ) [ 0 , 1 ] .
Therefore, we conclude that for any ( κ 1 , κ 2 , μ 1 , μ 2 ) Σ 0 , 0 , the set of equilibria of (1) is
{ ( 0 , 0 ) } { η ( x ) θ κ 1 , Q , μ 1 , r 1 , ( 1 η ( x ) ) θ κ 1 , Q , μ 1 , r 1 ρ 2 ( x ) : η ( x ) [ 0 , 1 ] } ,
where ( 0 , 0 ) is a repeller by Remark 1. Then each solution of (1) converges to a single equilibrium { η ( x ) θ κ 1 , Q , μ 1 , r 1 , ( 1 η ( x ) ) θ κ 1 , Q , μ 1 , r 1 ρ 2 ( x ) : η ( x ) [ 0 , 1 ] } . □

4. Discussion

In this paper, by using principal eigenvalue theory and monotone dynamical system theory, we mainly analyzed the global directed dynamic behaviors of a Lotka-Volterra competition-diffusion-advection system between two organisms in heterogeneous environments. The two organisms compete for different fundamental resources, their advection and diffusion strategies follow a positive diffusion distribution, the functions of inter-specific competition ability are variable. Our work can be seen as a further development of Wang [5] for the competition-diffusion system, where we bring new ingredients in the arguments to overcome the difficulty caused by the involvement of advection.
In the future, exploring the global directed dynamic behaviors under the condition of cross-diffusion may be an interesting research point. We leave this challenge to future investigations.

Author Contributions

All authors contributed equally and significantly in writing this article. Conceptualization, L.C.; Formal analysis, Y.Z.; Funding acquisition, L.C.; Writing—original draft, S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Shandong Provincial Natural Science Foundation under grant ZR2020MA006 and the Introduction and Cultivation Project of Young and Innovative Talents in Universities of Shandong Province.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to express our thanks to the anonymous referees and the editor for their constructive comments and suggestions, which greatly improved this article

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Smith, H. Monotone Dynamical System. An Introduction to the Theory of Competitive and Cooperative Systems; American Mathematical Society: Providence, RI, USA, 1995; Volume 41. [Google Scholar]
  2. He, X.; Ni, W. Global dynamics of the Lotka-Volterra competition-diffusion system: Diffusion and spatial heterogeneity I. Commun. Pure Appl. Math. 2016, 69, 981–1014. [Google Scholar] [CrossRef]
  3. He, X.; Ni, W. Global dynamics of the Lotka-Volterra competition-diffusion system with equal amount of total resources, II. Calc. Var. Partial Differ Equ. 2016, 55, 25. [Google Scholar] [CrossRef]
  4. Hess, P. Periodic-Parabolic Boundary Value Problems and Positivity; Wiley: New York, NY, USA, 1991. [Google Scholar]
  5. Wang, Q. Global directed dynamics of a Lotka-Volterra competition-diffusion system. Nonlinear Anal. Real World Appl. 2020, 55, 103144. [Google Scholar] [CrossRef]
  6. Zhou, P.; Tang, D.; Xiao, D. On Lotka-Volterra competitive parabolic systems: Exclusion, coexistence and bistability. J. Differ. Equ. 2021, 282, 596–625. [Google Scholar] [CrossRef]
  7. Du, L.; Li, W.; Wu, S. Propagation phenomena for a bistable Lotka-Volterra competition system with advection in a periodic habitat. Z. Angew. Math. Phys. 2019, 71, 2402–2435. [Google Scholar] [CrossRef]
  8. Du, L.; Li, W.; Wu, S. Pulsating fronts and front-like entire solutions for a reaction-advection-diffusion competition model in a periodic habitat. J. Differ. Equ. 2019, 266, 8419–8458. [Google Scholar] [CrossRef]
  9. Lou, Y.; Zhao, X.; Zhou, P. Global dynamics of a Lotka-Volterra competition-diffusion-advection system in heterogeneous environments. J. Math. Pures Appl. 2019, 121, 47–82. [Google Scholar] [CrossRef]
  10. Kamrujjaman, M. Interplay of resource distributions and diffusion strategies for spatially heterogeneous populations. J. Math. Model. 2019, 7, 175–198. [Google Scholar]
  11. Tang, D.; Chen, Y. Global dynamics of a Lotka-Volterra competition-diffusion system in advective homogeneous environments. J. Differ. Equ. 2020, 269, 1465–1483. [Google Scholar] [CrossRef]
  12. Xu, F.; Gan, W.; Tang, D. Global dynamics of a Lotka-Volterra competitive system from river ecology: General boundary conditions. Nonlinearity 2020, 33, 1528–1541. [Google Scholar] [CrossRef]
  13. Tang, D.; Zhou, P. On a Lotka-Volterra competition-diffusion-advection system: Homogeneity vs heterogeneity. J. Differ. Equ. 2020, 268, 1570–1599. [Google Scholar] [CrossRef]
  14. Ma, L.; Guo, S. Bifurcation and stability of a two-species reaction–diffusion–advection competition model. Nonlinear Anal. Real World Appl. 2021, 59, 103241. [Google Scholar] [CrossRef]
  15. Guo, Q.; He, X.; Ni, W. Global dynamics of a general Lotka-Volterra competition-diffusion system in heterogeneous environments. Discret. Contin. Dyn. Syst. 2020, 40, 6547–6573. [Google Scholar] [CrossRef]
  16. Korobenko, L.; Braverman, E. On evolutionary stability of carrying capacity driven dispersal in competition with regularly diffusing populations. J. Math. Biol. 2014, 69, 1181–1206. [Google Scholar] [CrossRef]
  17. Lou, Y.; Zhou, P. Evolution of dispersal in advective homogeneous environment: The effect of boundary conditions. J. Differ. Equ. 2015, 259, 141–171. [Google Scholar] [CrossRef]
  18. Ma, L.; Tang, D. Existence and Stability of Stationary States of a Reaction-Diffusion-Advection Model for Two Competing Species. Int. J. Bifurc. Chaos 2020, 5, 2050065. [Google Scholar] [CrossRef]
  19. Zhou, P.; Xiao, D. Global dynamics of a classical Lotka-Volterra competition-diffusion-advection system. J. Funct. Anal. 2018, 275, 356–380. [Google Scholar] [CrossRef]
  20. Averill, I.; Lam, K.; Lou, Y. The Role of Advection in a Two-Species Competition Model: A Bifurcation Approach; American Mathematical Society: Providence, RI, USA, 2017; Volume 245. [Google Scholar]
  21. Cantrell, R.; Cosner, C.; Lou, Y. Movement toward better environments and the evolution of rapid diffusion. Math. Biosci. 2006, 204, 199–214. [Google Scholar] [CrossRef]
  22. Cantrell, R.; Cosner, C.; Lou, Y. Advection-mediated coexistence of competing species. Proc. Roy. Soc. Edinb. Sect. A 2007, 137, 497–518. [Google Scholar] [CrossRef]
  23. Chen, X.; Lou, Y. Effects of diffusion and advection on the smallest eigenvalue of an elliptic operator and their applications. Indiana Univ. Math. J. 2012, 61, 45–80. [Google Scholar] [CrossRef]
  24. Lam, K. Concentration phenomena of a semilinear elliptic equation with large advection in an ecological model. J. Differ. Equ. 2011, 250, 161–181. [Google Scholar] [CrossRef] [Green Version]
  25. Lam, K. Limiting profiles of semilinear elliptic equations with large advection in population dynamics II. SIAM J. Math. Anal. 2012, 44, 1808–1830. [Google Scholar] [CrossRef] [Green Version]
  26. Lam, K.; Ni, W. Uniqueness and complete dynamics in heterogeneous competition-diffusion systems. SIAM J. Appl. Math. 2012, 72, 1695–1712. [Google Scholar] [CrossRef] [Green Version]
  27. Cantrell, R.; Cosner, C. Spatial Ecology via Reaction-Diffusion Equations; Wiley: New York, NY, USA, 2003. [Google Scholar]
  28. Krein, M.; Rutman, M. Linear operators leaving invariant a cone in a Banach space. Uspekhi Mat. Nauk. 1948, 3, 3–95. [Google Scholar]
  29. Sweers, G. Strong positivity in C(Ω¯) for elliptic systems. Math. Z. 1992, 209, 251–271. [Google Scholar] [CrossRef]
  30. Nagel, R. Operator matrices and reaction-diffusion system. Rend. Semin. Mat. Fis. Milano 1989, 59, 185–196. [Google Scholar] [CrossRef]
  31. De Figueiredo, D.; Mitidieri, E. Maximum principles for linear elliptic systems. Rend. Istit. Mat. Univ. Trieste 1990, 22, 36–66. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, L.; Lin, S.; Zhao, Y. Global Directed Dynamic Behaviors of a Lotka-Volterra Competition-Diffusion-Advection System. Axioms 2021, 10, 195. https://doi.org/10.3390/axioms10030195

AMA Style

Chen L, Lin S, Zhao Y. Global Directed Dynamic Behaviors of a Lotka-Volterra Competition-Diffusion-Advection System. Axioms. 2021; 10(3):195. https://doi.org/10.3390/axioms10030195

Chicago/Turabian Style

Chen, Lili, Shilei Lin, and Yanfeng Zhao. 2021. "Global Directed Dynamic Behaviors of a Lotka-Volterra Competition-Diffusion-Advection System" Axioms 10, no. 3: 195. https://doi.org/10.3390/axioms10030195

APA Style

Chen, L., Lin, S., & Zhao, Y. (2021). Global Directed Dynamic Behaviors of a Lotka-Volterra Competition-Diffusion-Advection System. Axioms, 10(3), 195. https://doi.org/10.3390/axioms10030195

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop