Next Article in Journal
Least-Squares Solutions of Eighth-Order Boundary Value Problems Using the Theory of Functional Connections
Previous Article in Journal
Darboux Associated Curves of a Null Curve on Pseudo-Riemannian Space Forms
Previous Article in Special Issue
Quantitatively Inferring Three Mechanisms from the Spatiotemporal Patterns
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Persistence for a Two-Stage Reaction-Diffusion System

by
Robert Stephen Cantrell
1,*,
Chris Cosner
1 and
Salomé Martínez
2
1
Department of Mathematics, University of Miami, Coral Gables, FL 33124, USA
2
Departamento de Ingeniería Matemática and Centro de Modelamiento Matemático, UMI 2807 CNRS-UChile, Universidad de Chile, 8370456 Santiago, Chile
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(3), 396; https://doi.org/10.3390/math8030396
Submission received: 3 February 2020 / Revised: 2 March 2020 / Accepted: 5 March 2020 / Published: 11 March 2020
(This article belongs to the Special Issue Partial Differential Equations in Ecology: 80 Years and Counting)

Abstract

:
In this article, we study how the rates of diffusion in a reaction-diffusion model for a stage structured population in a heterogeneous environment affect the model’s predictions of persistence or extinction for the population. In the case of a population without stage structure, faster diffusion is typically detrimental. In contrast to that, we find that, in a stage structured population, it can be either detrimental or helpful. If the regions where adults can reproduce are the same as those where juveniles can mature, typically slower diffusion will be favored, but if those regions are separated, then faster diffusion may be favored. Our analysis consists primarily of estimates of principal eigenvalues of the linearized system around ( 0 , 0 ) and results on their asymptotic behavior for large or small diffusion rates. The model we study is not in general a cooperative system, but if adults only compete with other adults and juveniles with other juveniles, then it is. In that case, the general theory of cooperative systems implies that, when the model predicts persistence, it has a unique positive equilibrium. We derive some results on the asymptotic behavior of the positive equilibrium for small diffusion and for large adult reproductive rates in that case.
AMS Classification:
92D40; 92D50; 35P15; 35K57

1. Introduction

The question of how dispersal interacts with spatial heterogeneity to influence population dynamics and species interactions has been studied extensively in recent years, specifically from the viewpoint of reaction-diffusion systems and related models—see, for example, [1,2,3] and the references cited therein. Most work on that topic assumes that each population is structured only by space and has only one mode of dispersal. However, populations are often structured by age, stage, or other attributes, and there may be variation among individuals in their dispersal rates or patterns. Here we will examine how the presence of a stage structure influences how diffusion rates influence population dynamics in a class of reaction-diffusion models for a population with two stages. In the case of a population with logistic growth, without age or stage structure, diffusing in a closed bounded spatially heterogeneous environment that is constant in time, it is well known that reaction-diffusion models predict that slower diffusion rates are advantageous relative to faster diffusion—see [4,5]. The results in [5] also hold for patch models. More broadly, a wide class of models arising in population genetics, population dynamics, and related areas display some version of the reduction principle, which says that dispersal, which causes faster mixing, typically reduces the rate of population growth—see [6]. However, the situation seems to be quite different in the case of stage structured populations. In [7], the authors considered a discrete-time patch model for a structured population and found that, in some cases, there was no selection against faster dispersal. The goal of the present paper is to use a spatially explicit reaction-diffusion model to understand how the spatial distributions of habitats that are favorable for reproduction by adults and those that are favorable for survival and growth by juveniles affect whether faster diffusion is advantageous or harmful for a stage structured population. We will see that the answer depends on the details of the spatial distribution of favorable and unfavorable habitats.
The type of reaction-diffusion model we will study is
u t = d 1 Δ u + r ( x ) v s ( x ) u a ( x ) u b ( x ) u 2 c ( x ) u v in Ω , t > 0 , v t = d 2 Δ v + s ( x ) u e ( x ) v f ( x ) v 2 g ( x ) u v in Ω , t > 0 , u · ν = v · ν = 0 on Ω , t > 0 .
Ω is a bounded domain in R N , and ν is the outward unit normal to Ω , such that the system has Neumann boundary conditions, which are the no-flux boundary conditions for simple diffusion. In this system, u and v represent the population densities of juveniles and individuals that have reached reproductive age, i.e. adults, respectively, of the same species. Thus, the term s ( x ) represents the rate at which juveniles mature into adults, which is determined by the fraction of individuals that reach reproductive age and the rate at which they mature, while r ( x ) accounts for the local fecundity of adults such that r ( x ) v ( x ) describes that rate at which new juveniles are produced by an adult population with density v at location x. The terms a ( x ) , b ( x ) , c ( x ) , e ( x ) , f ( x ) , and g ( x ) account for per-capita death rates and saturation factors due to logistic self-limitation. The diffusion coefficients d 1 and d 2 account for the the dispersal rates of juveniles and adults, respectively. The coefficients are all assumed to be nonnegative and continuous in Ω ¯ . This is the type of model for a stage structured population introduced in [8]. Related models with a different interpretation are discussed in [9,10] and in the references in those papers. The model expressed in Equation (1) is not an explicitly age-structured model. It assumes that individuals in the juvenile stage mature at some spatially dependent rate but does not track the age of individuals within each stage. Explicitly age structured models are considered in [11,12,13]. A different way of modeling an age-structured population, based on delayed reaction diffusion equations, is developed in [14]. Our focus here is on how spatial heterogeneity, dispersal, and stage structure interact, so we chose to use the simplest possible formulation of stage structure. In the case where c = g = 0 , the system is cooperative, and the methods and results of [15,16] would apply to it. The linearization of Equation (1) around ( 0 , 0 ) is cooperative, so the results of [15] apply to it; in particular, with a few technical assumptions, they imply that it has a principal eigenvalue.
The main questions we will address in this work are related to understanding the roles of the different functions and coefficients in Equation (1) in the persistence of the species. For the remainder of the paper we will focus primarily on understanding how the principal eigenvalue of the linearization of Equation (1) around ( 0 , 0 ) depends on the coefficients and what that dependence means biologically. We will see that whether faster diffusion is harmful or helpful for the persistence of the population depends on the details of the distribution of habitats that are favorable for adult reproduction and those that are favorable for juvenile survival and maturation. In some cases slower diffusion is still an advantage, but sometimes faster diffusion turns out to be helpful, and sufficiently fast diffusion may even be necessary for persistence. The spatial distribution of habitats favorable to adult reproduction ( r ( x ) large) relative to those favorable to juvenile development ( s ( x ) large) turns out to be important in some cases. Our analysis here is similar in spirit to the sorts of results obtained for diffusive Lotka-Volterra competition models in [3,17,18,19]. In particular, we will examine the behavior of the system for small, large, and general diffusion rates. Related results for some epidemiological models are derived in [20,21].
The linearization of Equation (1) around ( 0 , 0 ) has a principal eigenvalue whose sign determines whether the model predicts persistence or extinction. Since the sign of the principal eigenvalue of the linearization of Equation (1) around ( 0 , 0 ) determines the fate that Equation (1) predicts for the population it describes, we will study in detail the following problem:
d 1 Δ φ + r ( x ) ψ ( s ( x ) + a ( x ) ) φ = λ φ in Ω , d 2 Δ ψ + s ( x ) φ e ( x ) ψ = λ ψ in Ω , φ · ν = ψ · ν = 0 on Ω .

2. Basic Properties

In this section, we discuss some basic properties of Equation (1). From now on, we assume that r , s , a , b , c , e , f , g C α ( Ω ¯ ) , Ω is of class C 2 , α , and that the following hypotheses hold:
Hypothesis 1 (H1).
r ( x ) , s ( x ) 0 in Ω, with r ( x r ) 0 , s ( x s ) 0 for some x r , x s Ω .
Hypothesis 2 (H2).
a ( x ) , c ( x ) , e ( x ) , g ( x ) 0 in Ω.
Hypothesis 3 (H3).
b ( x ) > 0 , f ( x ) > 0 in Ω ¯ .
The model expressed in Equation (1) has many mathematical features in common with the models discussed in [9,10] for populations where individuals can switch between two different movement modes. A key feature is that the linear part of Equation (1) is cooperative, so it will have a principal eigenvalue which determines the stability of the equilibrium ( 0 , 0 ) and hence the persistence or extinction of the population. Another key feature of Equation (1) is that the nonlinearity is subhomogeneous. The maximum principle and existence of principal eigenvalues for cooperative linear systems such as the linear part of the right side of Equation (1) are derived in [15]. The general theory for systems such as Equation (1) is developed in [16] for the fully cooperative case (where c = g = 0 , such that adults only compete with other adults and juveniles with other juveniles) and in the general case in [8,9,10,22]. As expected, the sign of the principal eigenvalue of the linearization of Equation (1) around ( 0 , 0 ) gives us the relevant information to study the persistence of the species. If it is positive, the population will persist. If it is nonpositive, the population will go extinct. In the case where the coefficients of c and g in Equation (1) are zero such that the system is cooperative, the results and methods of [16] imply that, if the principal eigenvalue of the linear part is positive, then the system has a unique globally attractive equilibrium. If those coefficients are small, the methods of [9,10] can be applied to show that Equation (1) is asymptotically cooperative, and still has a unique globally attractive positive equilibrium. Combining results that are given in [8,9,10,15,16] or that follow directly by the same arguments used in those papers, we have the following:
Lemma 1.
The eigenvalue problem expressed in Equation (2) has a unique principal eigenvalue λ 1 that is characterized by having a positive eigenvector ( φ , ψ ) .
Lemma 2.
If λ 1 > 0 , then the system expressed in Equation (1) is persistent and has at least one positive equilibrium. If λ 1 0 , then ( 0 , 0 ) is globally asymptotically stable in Equation (1).
Lemma 3.
If λ 1 > 0 and c and g are sufficiently small, then the system expressed in Equation (1) has a unique globally attracting positive equilibrium.
Remark 1.
In the case that c = g = 0 , the system expressed in Equation (1) is cooperative and hence generates a monotone semi-flow on appropriate spaces.

3. The Case of d 1 , d 2 Small

Following the approach in [23], we will establish the asymptotic behavior of the principal eigenvalue of Equation (2) when d 1 , d 2 are small and, in the fully cooperative case (where c g 0 ), the profile of the nonnegative solutions of the corresponding steady state system for Equation (1).
d 1 Δ u + r ( x ) v s ( x ) u a ( x ) u b ( x ) u 2 = 0 in Ω , d 2 Δ v + s ( x ) u e ( x ) v f ( x ) v 2 = 0 in Ω , u · ν = v · ν = 0 on Ω ,
as well. Related results are derived in [16]. We observe that the associated kinetic system, which corresponds to Equation (1), is given by
u t = r ( x ) V ( x ) s ( x ) U ( x ) a ( x ) U ( x ) b ( x ) U 2 ( x ) c ( x ) U ( x ) V ( x ) = 0 , v t = s ( x ) U ( x ) e ( x ) V ( x ) f ( x ) V 2 ( x ) g ( x ) U ( x ) V ( x ) = 0 ,
for each x Ω .
For each x, this system shares the same properties as Equation (1) given in Lemmas 1–3, which we state for convenience.
Lemma 4.
Set x Ω . The linearization around (0,0) of Equation (4) has a principal eigenvalue λ 1 ( x ) . Moreover,
(i) 
If λ 1 ( x ) 0 , then all the solutions with nonnegative initial condition of Equation (4) converge to (0,0) when t .
(ii) 
If λ 1 ( x ) > 0 , then Equation (4) is persistent and has at least one positive equilibrium.
(iii) 
If c g 0 and λ 1 ( x ) > 0 , then Equation (4) is cooperative and admits a unique positive equilibrium, which is the global attractor for all nonnegative, non-trivial solutions.
Observe that, when d 1 = d 2 = 0 , the eigenvalues of the linearization around (0,0) of Equation (4) are the roots of det ( A ( x ) λ I ) , with
A ( x ) = ( s ( x ) + a ( x ) ) r ( x ) s ( x ) e ( x )
By a simple computation, we obtain that the maximum eigenvalue is given by
Λ ( x ) = 1 2 ( s ( x ) + a ( x ) + e ( x ) ) + ( s ( x ) + a ( x ) e ( x ) ) 2 + 4 r ( x ) s ( x ) ,
which is positive provided that ( s ( x ) + a ( x ) ) e ( x ) r ( x ) s ( x ) < 0 . Our first result, which is a direct application of Theorem 1.4 of [23] (Theorem A1 in the Appendix A), states that this is indeed the necessary and sufficient condition to have a positive principal eigenvalue when d 1 and d 2 are small.
Proposition 1.
The principal eigenvalue λ 1 of Equation (2) satisfies
λ 1 max x Ω ¯ Λ ( x ) as d 1 , d 2 0 .
Thus, there exists a δ > 0 such that if
min x Ω ¯ ( ( s ( x ) + a ( x ) ) e ( x ) r ( x ) s ( x ) ) < 0 ,
the principal eigenvalue of Equation (2) is positive for all 0 < d 1 , d 2 < δ , while if
min x Ω ¯ ( ( s ( x ) + a ( x ) ) e ( x ) r ( x ) s ( x ) ) > 0 ,
the principal eigenvalue is negative.
As a consequence of this result, if Equation (9) holds, the unique nonnegative equilibrium of the system expressed in Equation (3) is ( 0 , 0 ) , and that equilibrium is globally attracting, whenever d 1 , d 2 are small; if Equation (8) holds, the system expressed in Equation (3) is persistent, and has a positive equilibrium for small d 1 , d 2 . If Equation (8) holds and c g 0 , then by Lemma 3 the system expressed in Equation (3) has a unique globally attracting positive equilibrium, which we denote by ( u d , v d ) with d = ( d 1 , d 2 ) .
Throughout the remainder of this section, we will assume that c g 0 in Ω , in which case, the system expressed in Equation (1) is cooperative.
The next result establishes the convergence of ( u d , v d ) to the unique nonnegative steady state ( U ( x ) , V ( x ) ) of the kinetic system, which satisfies
( U ( x ) , V ( x ) ) is positive where ( s ( x ) + a ( x ) ) e ( x ) r ( x ) s ( x ) < 0 , ( U ( x ) , V ( x ) ) = 0 where ( s ( x ) + a ( x ) ) e ( x ) r ( x ) s ( x ) 0 .
Theorem 1.
Suppose that Equation (8) holds. Then ( u d , v d ) ( U , V ) as d 0 locally uniformly in Ω.
To prove this theorem, we follow the proof of Theorem 1.5 of [23], specifically their Proposition 5.2 in [23] and its hypotheses, which are listed in [23] as (A1)–(A4) and are given in the Appendix A as (L1)-(L4) to avoid confusion with the equation labels there. We should point out that Assumptions (A2) and (A3) of [23] do not hold in our case, so we cannot apply that result directly. The difference is that we allow situations where the kinetic system expressed in Equation (4) has a positive equilibrium for some values of x Ω ¯ but not for others, whereas Condition (A2) requires a positive equilibrium for the kinetic system for all x. For that reason, we need to construct a version of the arguments in [23] that is local in x. Condition (A3) in [23] is used only to prove the existence of a nontrivial subsolution for a system corresponding to Equation (3), which is independent of d 1 , d 2 . We show the existence of the analogous local subsolutions we need in our case in the next lemma.
Lemma 5.
Suppose that x ˜ Ω 0 , where
Ω 0 = { x Ω / ( s ( x ) + a ( x ) ) e ( x ) r ( x ) s ( x ) < 0 } .
Then there exists d 0 > 0 , ρ 0 > 0 , and a function w ̲ 0 > 0 in B ( x ˜ , ρ ) Ω 0 , which is a subsolution of Equation (3) for all 0 < d 1 , d 2 < d 0 .
Proof. 
Let p = ( p 1 , p 2 ) , a positive eigenvector of A ( x ˜ ) with p 1 + p 2 = 1 , associated with its principal eigenvalue σ ˜ > 0 . Set ε > 0 and small. We can choose ρ > 0 such that B ( x ˜ , ρ ) ¯ Ω 0 and
| a ( x ) a ˜ | < ε , | r ( x ) r ˜ | < ε , | s ( x ) s ˜ | < ε and | e ( x ) e ˜ | < ε for all x in B ( x ˜ , ρ ) ¯ ,
where a ˜ = a ( x ˜ ) , r ˜ = r ( x ˜ ) , s ˜ = s ( x ˜ ) , and e ˜ = e ( x ˜ ) . Set η > 0 as the principal eigenfunction associated with λ > 0 , the principal eigenvalue of
Δ η + λ η = 0 in B ( x ˜ , ρ ) , η = 0 on B ( x ˜ , ρ ) ,
with max B ( x ˜ , ρ ) η = 1 . We claim that we can choose δ , ε , ρ , d 0 > 0 such that δ η p is a subsolution of Equation (3) for all d 1 , d 2 < d 0 . For simplicity and to keep the notation consistent with that in [23], we define F ( x , u , v ) = ( F 1 ( x , u , v ) , F 2 ( x , u , v ) ) , with
F 1 ( x , u , v ) = r v s u a u b u 2 and F 2 ( x , u , v ) = s u e v f v 2 = 0 ,
where we have omitted the variable x in a , b , e , r , s , f to shorten the expressions. Observe that
F 1 ( x , δ η p ) = δ η ( σ ˜ p 1 + ( r r ˜ ) p 2 ( s s ˜ ) p 1 ( a a ˜ ) p 1 b δ p 1 2 η ) F 2 ( x , δ η p ) = δ η ( σ ˜ p 2 + ( s s ˜ ) p 1 ( e e ˜ ) p 2 f δ p 2 2 η ) ,
and using Equation (12), we obtain that, if we choose ε > 0 and a small δ , we have that
F 1 ( x , δ η p ) δ η ( σ ˜ p 1 ε p 2 2 ε p 1 b δ p 1 2 η ) > δ η σ ˜ 2 p 1 F 2 ( x , δ η p ) δ η ( σ ˜ p 2 ε p 1 ε p 2 f δ p 2 2 ) > δ η σ ˜ 2 p 2 .
Therefore, replacing these inequalities in Equation (3), we obtain
d 1 δ p 1 Δ η + F 1 ( x , δ η p ) δ η d 1 p 1 λ + σ ˜ 2 p 1 d 2 δ p 2 Δ η + F 2 ( x , δ η p ) δ η d 2 p 2 λ + σ ˜ 2 p 2 ;
hence, if we set d 0 = σ ˜ 2 λ . we obtain the desired result. □
Using this lemma we can follow the proof of Proposition 5.2 in [23]. To facilitate our exposition, we will use the same notation. Set the operators D = diag ( d 1 , d 2 ) , L = diag ( Δ , Δ ) . To prove Theorem 1, we will state the needed lemmas, discussing their relationships with the lemmas in [23] leading to the proof of Proposition 5.2.
Suppose that Equation (8) holds, setting w ̲ 0 = η δ p as in Lemma 5, and w ¯ 0 = M where M > 0 is given in Assumption (A4) such that F 1 ( x , u , v ) c u and F 2 ( x , u , v ) c v for all u , v M and x Ω , with c > 0 fixed. Set K > 0 such that K + u F 1 ( x , u , v ) > 0 and K + v F 2 ( x , u , v ) > 0 for all 0 u , v M , and we define z = w ¯ k as the unique solution of
D L z + K z = K u + F ( x , u ) in Ω , z · ν = 0 on Ω ,
for u = w ¯ k 1 .
Lemma 6.
Suppose that Equation (8) holds. For every k, we have w ̲ 0 < w ¯ k + 1 < w ¯ k < w ¯ 0 , and as k , w ¯ k converges uniformly to the unique positive solution w of Equation (3), which satisfies w ̲ 0 < w < w ¯ k in Ω for all k 0 .
Proof. 
We will prove that w ̲ 0 < w ¯ k by induction. Suppose this is true for k. Observe that w ̲ 0 < w ¯ 0 by construction. In the set, B ( x ˜ , ρ ) Ω 0 as in Lemma 5 w ¯ k + 1 satisfies
D L ( w ¯ k + 1 w ̲ 0 ) + K ( w ¯ k + 1 w ̲ 0 ) = K ( w ¯ k w ̲ 0 ) + F ( x , w ¯ k ) F ( w ̲ 0 ) ,
in B ( x ˜ , ρ ) . By the induction hypothesis w ̲ 0 < w ¯ k , whence K w ¯ k + F ( x , w ¯ k ) > 0 ; hence, by the strong maximum principle applied to each component, we have that w ¯ k + 1 > 0 in Ω ¯ . Thus, we have
D L ( w ¯ k + 1 w ̲ 0 ) + K ( w ¯ k + 1 w ̲ 0 ) > 0 in B ( x ˜ , ρ ) , w ¯ k + 1 w ̲ 0 > 0 in B ( x ˜ , ρ ) ,
such that we have that w ¯ k + 1 w ̲ 0 > 0 in B ( x ˜ , ρ ) ¯ . The remainder of the proof is a standard monotone iteration argument, just as in the proof of Lemma 5.3 of [23]. We observe that
D L ( w ¯ 1 w ¯ 0 ) + K ( w ¯ 1 w ¯ 0 ) = K w ¯ 0 + F ( x , w ¯ 0 ) K w ¯ 0 < 0 in Ω , [ w ¯ 1 w ¯ 0 ] · ν = 0 on Ω ,
Thus, by the strong maximum principle, we have w ¯ 1 < w ¯ 0 .
Similarly, if w ¯ k < w ¯ k 1 , then
D L ( w ¯ k + 1 w ¯ k ) + K ( w ¯ k + 1 w ¯ k ) = K w ¯ k + F ( x , w ¯ k ) K w ¯ k 1 F ( x , w ¯ k 1 ) < 0 in Ω , [ w ¯ k + 1 w ¯ k ] · ν = 0 on Ω .
By induction, the sequence { w ¯ k } is decreasing, and it is bounded below by m a x { 0 , w ̲ 0 ( x ) } , so by standard elliptic theory, it converges to a nonnegative nontrivial solution of Equation (3) as k . Since by Lemma 3 the nontrivial nonnegative solution of Equation (3) is unique, it coincides with the one constructed as the limit of the sequence { w ¯ k } . □
Define w ¯ 0 = w ¯ 0 and w ¯ k + 1 = w ¯ k + F ( x , w ¯ k ) in Ω . Following the proof of Lemmas 5.6 and 5.7 in [23], we can prove the following result.
Lemma 7.
Suppose that Equation (8) holds. For every k, we have
w ̲ 0 < w ¯ k + 1 < w ¯ k .
Then, w ¯ k converges locally uniformly to W as k , with
W = ( U ( x ) , V ( x ) ) in Ω 0 , and W ( x ) = 0 in Ω \ Ω 0 ,
where Ω 0 is given by Equation (11).
Proof. 
Observe that, by Equation (15), the function w ̲ 0 is a subsolution of the kinetic system. Repeating the proof of Lemma 5.6 in [23], or following the arguments leading to the monotonicity of the proof of Lemma 6 above, we have that the sequence { w ¯ k } is monotone decreasing and bounded below by w ̲ 0 . Therefore, w ¯ k W pointwise, which satisfies F ( x , W ) = 0 , i.e. a nonnegative equilibrium of the kinetic system. Therefore, if at some x Ω we have that ( s ( x ) + a ( x ) ) e ( x ) r ( x ) s ( x ) 0 , then W ( x ) = 0 . On the other hand, W ( x ) w ̲ 0 ( x ) > 0 in B ( x ˜ , ρ ) Ω 0 for x Ω 0 . Since x is arbitrary and the sequence does not depend on w ̲ 0 , we obtain that W ( x ) = ( U ( x ) , V ( x ) ) , the unique positive kinetic equilibrium, whenever x Ω 0 . Particularly, W is continuous. Using Theorem 5.8 of [23], we obtain that the convergence is uniform in any compact set of Ω . □
Lemma 8.
For each k, as d 1 , d 2 0 , we have that w ¯ k converges to w ¯ k uniformly in Ω ¯ .
The proof of this result is the same as the one of Lemma 5.5 in [23].
Proof of Theorem 7.
Using a diagonal argument, and Lemma (8), the unique positive solution w of Equation (3) converges to W as d 1 , d 2 . □

4. The Case of d 1 and d 2 Large

We will start by giving a proof of a result that is well known as a “folk theorem.” It is stated in slightly more generality than is needed for the specific application. For i = 1 , , N , let L i denote the operator
L i u = · μ i ( x ) [ u u α i ( x ) ] for x Ω
with no-flux boundary conditions
[ u u α i ] · ν = 0 for x Ω .
Assume that μ i ( x ) μ 0 > 0 on Ω ¯ for all i. Let A = ( a i j ( x ) ) be an N × N irreducible matrix with a i j 0 if i j . Consider the eigenvalue problem
d i L i φ i + j = 1 N a i j φ j = λ φ j , i = 1 N
where d i > 0 for all i and φ i satisfies the boundary condition expressed in Equation (17) for each i. Note that, if we let Φ i = e x p ( α i ( x ) ) φ i , then Φ i satisfies Neumann boundary conditions such that the system expressed in Equation (18) rewritten in terms of the variables Φ i is still cooperative. Because of the classical boundary conditions, the usual results on elliptic regularity and on maximum principles for cooperative systems from [15,23] can be applied to the system for the Φ i values, such that the system and hence Equation (18) will have a principal eigenvalue under suitable conditions on the domain Ω and the coefficients. This idea has been used in models for single populations without an age structure or competing pairs of such populations—see, for example, [2,24].
Furthermore, we have L i ( e x p ( α i ( x ) ) = 0 such that the principal eigenvalue of L i is zero, and the eigenfunction is a multiple of e x p ( α i ) . Let A ¯ be the matrix defined by
A ¯ i j : = Ω a i j e x p ( α i ) d x Ω e x p ( α i ) d x .
Denote the principal eigenvalue of Equation (18) as λ 1 ( d ) where d = ( d 1 , , d N ) . Denote the principal eigenvalue of A ¯ as Λ ¯ .
Lemma 9.
Suppose that, for some γ ( 0 , 1 ) , the coefficients of Equation (18) satisfy α C 2 , γ ( Ω ¯ ) , μ C 1 , γ ( Ω ¯ ) , and a i j C γ ( Ω ¯ ) for i , j = 1 N , and that Ω is of class C 2 , γ . Suppose further that A ¯ is irreducible. If m i n { d i : i = 1 , N } , then λ 1 ( d ) Λ ¯ .
Proof. 
Choose any sequence d n = ( d 1 n , d N n ) such that m i n { d i n : i = 1 , N } . Choose any subsequence, then renumber it as d n . Let λ n be the principal eigenvalue of Equation (18) corresponding to d n and let φ i n ( x ) > 0 be the ith component of the eigenvector, where the eigenvector is normalized by m a x { φ i n ( x ) : x Ω ¯ , i = 1 , , N } = 1 . Integrating the ith equation of Equation (18) over Ω and summing over i yields
λ n Ω i = 1 N φ i n ( x ) d x = Ω i , j = 1 N a i j ( x ) φ j n ( x ) d x A 1 Ω i = 1 N φ i n ( x ) d x
where A 1 is a constant, depending only on A. It follows that λ n is uniformly bounded from above. Similarly, λ n is uniformly bounded from below. Thus, any subsequence of λ n itself has a convergent subsequence. It then follows from dividing the ith equation of Equation (18) by d i n that L i φ i n is uniformly bounded, and L i φ i n 0 as n . By elliptic regularity, the sequence φ i n is uniformly bounded in W 2 , p ( Ω ) for any p < , then by Sobolev embedding, it has a subsequence that is convergent in C 1 ( Ω ¯ ) and weakly convergent in W 2 , p ( Ω ) . This will be true for any i. Taking a further subsequence if necessary and renumbering again, we obtain a sequence where λ n λ * for some λ * and φ i n φ i * for all i, with L i φ i * = 0 . We then must have φ i * = c i e x p ( α i ) for some nonnegative constant c i , and with m a x { φ i * ( x ) : x Ω ¯ , i = 1 N } = 1 . Integrating Equation (18) over Ω and using the no-flux boundary conditions gives
j = 1 N Ω a i j ( x ) φ j * ( x ) d x = λ * Ω φ i * , i = 1 N ,
such that
j = 1 N Ω a i j ( x ) e x p ( α j ( x ) ) d x Ω e x p ( α i ( x ) ) d x c j = λ * c i , i = 1 N .
It follows that ( c 1 , , c N ) must be a nontrivial nonnegative eigenvector of A ¯ with the normalization prescribed by m a x { c i e x p ( α i ( x ) ) : x Ω ¯ , i = 1 , N } = 1 . These last conditions uniquely determine the limits of the subsequence of the original subsequence { λ ( d n ) , φ n } . Since every subsequence of the original sequence { λ ( d n ) , φ n } has a subsequence converging to the values determined by Equation (21), the same must be true for the original sequence. Since the original sequence of values { d n } could be any increasing sequence that approaches infinity as n , the conclusion of the lemma follows. □
In the specific system expressed in Equation (1) that we consider, L i = Δ , such that α i and μ i are constants. In that case, we have A ¯ i j = A ¯ i j , where a ¯ i j is the average of a i j over Ω . Denote the averages of the coefficients in Equation (1) by r ¯ , s ¯ , etc. Calculations analogous to those in Equation (5), Equation (6), and the related discussion then yield the following:
Corollary 1.
Suppose that the hypotheses of Lemma 9 are satisfied. There exists a D > 0 such that if
e ¯ ( s ¯ + a ¯ ) r ¯ s ¯ < 0 ,
the principal eigenvalue λ 1 of Equation (2) is positive for all d 1 , d 2 > D , while if
e ¯ ( s ¯ + a ¯ ) r ¯ s ¯ > 0 ,
the principal eigenvalue is nonpositive.
Remark 2.
In the ODE system corresponding to Equation (1) with coefficients averaged over Ω, one can compute R 0 as r ¯ s ¯ / [ e ¯ ( s ¯ + a ¯ ) ] via the methods of [25]. The first inequality in Corollary 1 is equivalent to R 0 > 1 , while the second is equivalent to R 0 < 1 . By writing R 0 = [ r ¯ / ( s ¯ + a ¯ ) ] [ s ¯ / e ¯ ] , we can interpret the condition for persistence as saying that the products of the ratios of the growth terms over the loss terms for adults and juveniles should be greater than 1 for persistence.

5. General Diffusion Rates

Case1: Persistence or extinction for all diffusion rates
Proposition 2.
If
Ω r s d x 1 2 Ω ( s + a + e ) d x > 0
then λ 1 > 0 for all positive diffusion rates.
If
m i n x Ω ¯ [ 4 ( s ( x ) + a ( x ) ) e ( x ) ( r ( x ) + s ( x ) ) 2 ] > 0
then λ 1 < 0 for all positive diffusion rates.
Proof. 
If we divide the first equation in Equation (2) by φ and integrate over Ω , using Green’s formula to integrate the term Δ φ / φ , we obtain the inequality
| Ω | λ 1 Ω r ψ φ d x Ω ( s + a ) d x .
Similarly, if we divide the second equation by ψ and integrate we obtain
| Ω | λ 1 Ω s φ ψ d x Ω e d x .
If we add Equations (24) and (25) and divide by 2, we obtain
λ 1 1 2 | Ω | Ω r ψ φ + s φ ψ d x Ω ( s + a + e ) d x .
By Cauchy’s inequality, r z + s z 1 2 r s for all z > 0 , so from Equation (26) we obtain
λ 1 1 | Ω | Ω r s d x 1 2 Ω ( s + a + e ) d x
Thus, λ 1 > 0 if Equation (22) holds, so the first part of Proposition 2 holds. Going in the other direction, if we multiply the first equation of Equation (2) by φ and integrate, using integration by parts on the φ Δ φ term, and similarly multiply the second equation by ψ and integrate, and then add the results, we get
λ 1 Ω ( φ 2 + ψ 2 ) d x Ω [ ( s + a ) φ 2 + ( r + s ) φ ψ e ψ 2 ) d x ] .
The integrand on the right side of Equation (28) is a quadratic form in φ and ψ , which will be negative definite if
4 ( s + a ) e > ( r + s ) 2 ,
so λ 1 < 0 if Equation (23) holds, which proves the second part of Proposition 2. □
Remarks: Note that the first integral in Equation (22) is what appears in the formula for the Bhattacharyya coefficient [26,27], which is used to compare how well probability distributions match each other. Specifically, if two probability distributions P and Q have probability density functions p ( x ) and q ( x ) for x U R n , the Bhattacharyya coefficient is
B C ( P , Q ) = U p ( x ) q ( x ) d x .
For any P and Q, 0 B C ( P , Q ) 1 . If B C ( P , Q ) = 1 , then P and Q are the same, that is, p = q a.e. If B C ( P , Q ) = 0 , then the supports of p and q are disjoint. If we write r ( x ) = r 0 ρ ( x ) and s ( x ) = s 0 σ ( x ) such that Ω ρ ( x ) = Ω σ ( x ) = 1 , then we can compute r 0 = r ¯ | Ω | and s 0 = s ¯ | Ω | . We can treat ρ and σ as if they were probability density functions for distributions R and S. We then have
Ω r s d x = | Ω | r ¯ s ¯ B C ( R , S ) .
The maximum of B C ( R , S ) is 1, corresponding to the case where r and s are multiples of each other, and the minimum is 0, corresponding to the case where the supports of r and s are disjoint. Thus, the degree to which ρ and σ match each other has a strong impact on the estimate for λ in Equation (27).
Using Equation (30) and the fact that B C ( R , S ) 1 in Equation (22) shows that Equation (22) implies 2 r ¯ s ¯ > [ ( s ¯ + a ¯ ) + e ¯ ] . Squaring both sides and using Cauchy’s inequality implies e ¯ ( s ¯ + a ¯ ) r ¯ s ¯ < 0 as in the first case of Corollary 1. Similarly, if Equation (22) holds, then 2 r ( x ) s ( x ) > ( s ( x ) + a ( x ) ) + e ( x ) for some x Ω , and it then follows in the same way that the inequality in the first case of Proposition 1 holds. If Equation (23) holds, then Equation (29) holds, and then by Cauchy’s inequality the second case of Proposition 1 holds. Thus, the conditions expressed in Equation (22) and Equation (23) in Proposition 2, which imply λ 1 > 0 or λ 1 < 0 for all diffusion rates, also imply some of the corresponding conditions we have obtained for either large or small diffusion rates.
In the situation where the spatial distributions of habitat quality r for reproduction by adults and s for survival and maturation of juveniles into adults are perfectly correlated, such that r ( x ) = r 1 s ( x ) for some constant r 1 , the eigenvalue problem expressed in Equation (2) can be rewritten as a weighted symmetric eigenvalue problem by multiplying the second equation in Equation (2) by r 1 , which yields
d 1 Δ φ ( x ) ( s ( x ) + a ( x ) ) φ + r ( x ) ψ = λ φ d 2 r 1 Δ ψ + r ( x ) φ r 1 e ( x ) ψ = λ r 1 ψ .
The principal eigenvalue for Equation (31) has a variational characterization of λ 1 as
λ 1 = max φ , ψ W 1 , 2 ( Ω ) Ω ( d 1 | φ | 2 d 2 r 1 | ψ | 2 ( s + a ) φ 2 + 2 r φ ψ r 1 e ψ 2 ) d x Ω ( φ 2 + r 1 ψ 2 ) d x .
It follows in that case that λ 1 is decreasing in both d 1 and d 2 , such that slower diffusion is advantageous.
Case 2: Asymptotic behavior for large reproductive rates
Suppose that r ( x ) = n r 0 ( x ) and that s ( x ) r 0 ( x ) > 0 for x Ω 0 with Ω 0 ; therefore, there is a region where both the adult reproduction rate and the juvenile maturation rate are positive. The factor n scales the reproductive rate of adults in regions where r 0 ( x ) > 0 . For any fixed diffusion rates, it turns out that, for sufficiently large values of the scaling coefficient n, the principal eigenvalue of Equation (2) is positive, so the system expressed in Equation (1) is persistent. We will characterize the asymptotic behavior of the principal eigenvalue as n . If we make the further assumption that g c 0 , then the system expressed in Equation (1) is cooperative, so for an n large enough that the principal eigenvalue of Equation (2) is positive, Equation (1) has a unique positive equilibrium, and we will characterize the behavior of that equilibrium as n as well in that case.
Let λ 1 n denote the principal eigenvalue for Equation (2) with r ( x ) = n r 0 ( x ) . Observe that, since s ( x ) r 0 ( x ) > 0 for x Ω 0 and Ω 0 , Proposition 2 implies that λ 1 n > 0 for n sufficiently large, and in fact by Equation (27), λ 1 n as n . The following proposition states the asymptotic behavior of λ 1 n as n .
Proposition 3.
If r ( x ) = n r 0 ( x ) and s ( x ) r 0 ( x ) > 0 for x Ω 0 with Ω 0 , then
lim n λ 1 n n max x Ω ¯ ( r 0 ( x ) s ( x ) ) .
Proof. 
We start by noting that, if λ 1 n , then ( φ n , ψ n ) are the principal eigenvalue and corresponding eigenfunction of Equation (2) for r ( x ) = n r 0 ( x ) ; therefore, λ 1 n / n , φ ^ n = φ n , and ψ ^ n = n ψ n are the principal eigenvalue and eigenfunction of the following problem:
d 1 n Δ φ ^ ( s ( x ) + a ( x ) ) n φ ^ + r 0 ( x ) ψ ^ = λ ^ φ ^ in Ω , d 2 n Δ ψ ^ e ( x ) n ψ ^ + s ( x ) φ ^ = λ ^ ψ ^ in Ω , φ ^ · ν = ψ ^ · ν = 0 on Ω .
Considering the elliptic operators L 1 u = d 1 Δ u ( s ( x ) a ( x ) ) u and L 2 v = d 2 Δ v e ( x ) v , D = diag 1 n , 1 n and L = diag ( L 1 , L 2 ) , the system expressed in Equation (33) satisfies the hypothesis of Theorem 1.4 of [23]. Thus, n
λ 1 n n max x Ω ¯ λ ( A ( x ) ) ,
where
A ( x ) = 0 r 0 ( x ) s ( x ) 0 ,
which has eigenvalues ± r 0 ( x ) s ( x ) , from whence the result follows. □
In the case where g c 0 such that Equation (1) is cooperative, Lemma 3 implies that Equation (1) has a unique positive equilibrium if the principal eigenvalue of Equation (2) is positive. The next result states the asymptotic behavior of the unique positive equilibrium of Equation (1) for n large in that case.
Proposition 4.
Suppose the hypotheses of Proposition 3 are satisfied and that g c 0 . Let ( u n , v n ) be the unique positive solution of Equation (3). Then n 2 3 ( u n , v n ) ( U , V ) uniformly in Ω ¯ where
U ( x ) = r 0 ( x ) 2 3 b ( x ) 2 3 s ( x ) 1 3 f ( x ) 1 3 , V ( x ) = s ( x ) 2 3 f ( x ) 2 3 r 0 ( x ) 1 3 b ( x ) 1 3 w h e n s ( x ) r 0 ( x ) > 0 U ( x ) = 0 , V ( x ) = 0 w h e n s ( x ) r 0 ( x ) = 0 .
Proof. 
To prove this result, we use a different scaling. After some simple computations, we find that
( w n , z n ) = u n n 2 3 , v n n 1 3 ,
is the unique positive solution of the scaled system
n 2 3 [ d 1 Δ w ( s ( x ) + a ( x ) ) w ] + r 0 ( x ) z b ( x ) w 2 = 0 in Ω , n 1 3 [ d 2 Δ z e ( x ) z ] + s ( x ) w f ( x ) z 2 = 0 in Ω , w · ν = z · ν = 0 on Ω .
We set the operators
L 1 w = d 1 Δ w ( s ( x ) + a ( x ) ) w and L 2 z = d 2 Δ z e ( x ) z ,
D = diag n 2 3 , n 1 3 and L = ( L 1 , L 2 ) , and
F ( x , w , z ) = ( r 0 ( x ) z b ( x ) w 2 , s ( x ) w f ( x ) z 2 ) .
To prove the proposition, we can follow the same steps as in the proof of Theorem 1. Indeed, the principal eigenvalue of the linearization around ( 0 , 0 ) of the associated kinetic system of Equation (35) is s ( x ) r 0 ( x ) , and when it is positive, the kinetic equilibrium is given by the right hand side of Equation (34). This concludes the proof. □

6. Conclusions

The most fundamental conclusion from our analysis is that reaction-diffusion models for populations with a stage structure in spatially heterogeneous environments do not necessarily predict that slower diffusion is advantageous for persistence. This is in contrast to the case where populations are structured only by spatial location, where a version of the reduction principle [6] applies; in a competition between otherwise identical populations with different diffusion rates, the prediction is that “the slower diffuser wins” [4,5]. The mechanism underlying this observation is that, in our structured model, the regions where it is possible for adults to produce offspring may be separated from those where juveniles can survive and mature into adults. The conditions we find that imply persistence generally require that the product r ( x ) s ( x ) of the reproductive rate of adults and the maturation rate of juveniles be sufficiently large relative to their death rates. For slow diffusion, the condition for persistence is that r ( x ) s ( x ) > e ( x ) ( s ( x ) + ( a ( x ) ) at some point x Ω . For fast diffusion, it is r ¯ s ¯ > e ¯ ( s ¯ + a ¯ ) where r ¯ , s ¯ , e ¯ , and a ¯ are the spatial averages of those quantities. If the spatial distributions of r and s are closely correlated and are large in a few places but small in most, such that the maximum of r s is large but the averages r ¯ and s ¯ are small, the condition for persistence with slow diffusion may be satisfied, while the condition with fast diffusion may fail. In that type of environment, slow diffusion is clearly favored. Furthermore, if r and s are perfectly correlated in the sense that they are multiples of each other, the principal eigenvalue determining the growth rate of the population at low density is decreasing with respect to the diffusion rates, as in the case of unstructured populations in heterogeneous environments. On the other hand, if both r and s are large on some regions but very small outside of them, and the regions where they are large are disjoint (that is, separated from each other), then the product r s could be small everywhere, but the averages r ¯ and s ¯ could be large. In that case, the condition for persistence with small diffusion may fail, but the condition with fast diffusion may be satisfied, such that fast diffusion is favored.
We found that a sufficient condition for persistence for all diffusion rates is
Ω r s d x 1 2 Ω ( s + a + e ) d x > 0 .
The first term can be written as r ¯ s ¯ | Ω | B C ( r ( x ) / r ¯ , s ( x ) / s ¯ ) , where B C denotes the Bhattacharrya coefficient (see [26,27]), which measures how closely probability densities match each other. For distributions that are equal to each other, B C = 1 , but for distributions that are mutually exclusive in the sense that the regions where they are positive do not intersect, B C = 0 . This observation again shows that the degree to which the spatial distributions of r and s match each other is significant in determining the predictions of the model expressed in Equation (1).
Finally, we found that, if we scale the adult reproductive rate as r ( x ) = n r 0 ( x ) and there is some overlap between the distributions of r and s such that s ( x ) r 0 ( x ) > 0 on some subset of Ω with positive measure, then for any fixed diffusion rates the system expressed in Equation (1) will be persistent if n is sufficiently large. This means a population with any diffusion rates can persist if there is even a modest overlap between the regions where adults can reproduce and where juveniles can mature, provided that the reproductive rate of adults is sufficiently large. We characterized the asymptotic behavior as n of the principal eigenvalue of Equation (2). In the cooperative case where g c 0 , the system expressed in Equation (1) will have a unique positive equilibrium if it is persistent, and in that case we also characterized the asymptotic behavior as n of the equilibrium.
There are several directions for further research on the general topic of this paper. It would be of interest to take the approach of [4] and consider the competition between two stage structured populations described by systems such as Equation (1) that differ only in their diffusion rates. That would be somewhat challenging because it would involve systems of four equations, but at least in the cooperative case where c = g = 0 the general theory of monotone dynamical systems as in [28] and some of the ideas and methods of [10] would apply. It would also be interesting to consider models with an explicit age structure, as introduced in [11] and studied in [12,13]. Finally, it would be interesting albeit challenging to consider the case of time-periodic environments with spatial heterogeneity. Temporal variation alone is sufficient to cause faster diffusion to be favored in such environments in some cases (see [29]), but even without a stage structure, the time-dependent case is challenging, and there are many open questions.

Author Contributions

Conceptualization, R.S.C., C.C. and S.M.; mathematical analysis, R.S.C., C.C. and S.M.; writing–original draft preparation, C.C. and S.M.; writing–review and editing, R.S.C., C.C., and S.M.; funding acquisition, R.S.C., C.C. and S.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received external funding as noted in the subsequent Acknowledgments.

Acknowledgments

R.S.C. and C.C. are supported in part by NSF Awards DMS 15-14792 and 18-53478. S.M. is supported by CONICYT + PIA/Concurso de Apoyo a Centros Científicos y Tecnológicos de Excelencia con Financiamiento Basal AFB170001.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In [23], the authors considered the equilibria and dynamics of the system
u t = D L u + F ( x , u ) in Ω × ( 0 , ) , B u = 0 on Ω × ( 0 , )
where u = ( u 1 , , u n ) T is a vector of smooth functions, Ω is a bounded domain in R N with smooth boundary, u = ( u 1 ( x ) , u n ) T is a vector of smooth functions, D = d i a g ( d 1 , d n ) is a diagonal matrix of positive constants, L = d i a g ( L 1 , L n ) is a diagonal matrix of second order uniformly strongly elliptic operators of the form
L i = j , k = 1 N α j k i 2 x j x k + j = 1 N β j i x j + γ i
with smooth coefficients, and B = ( B 1 , B n ) , where for each i, B i defines a Dirichlet, Neumann, or Robin boundary condition. (They include Neumann as a case of Robin.) They also considered the associated linearized problem, which they wrote as
D L ϕ + A u ϕ = λ ϕ in Ω , B ϕ = 0 on Ω ,
where A = ( a i j ) is an n × n matrix of smooth functions with a i j 0 for i j , and ϕ = ( ϕ 1 ( x ) , ϕ n ( x ) ) T is a vector of smooth functions.
Note that, in our notation, we use the opposite sign convention to the one used in [23], such that what they denote as λ , we denote as λ .
For details on what the specific smoothness assumptions require, see [23]. Under those assumptions, by the Perron-Frobenius theorem, for each x Ω , the matrix A has a principal eigenvalue, which in our notation we denote as Λ ( x ) . The first major result of [23] is their Theorem 1.4, which can be stated as follows:
Theorem A1.
(Theorem 1.4 of [23]) The principal eigenvalue λ 1 of the system expressed in Equation (A2) with Dirichlet, Neumann, or Robin boundary conditions satisfies
lim m a x { d 1 , , d n } 0 λ 1 = max x Ω ¯ Λ ( x ) .
Except for the adjustments needed for our different notation, that theorem applies directly to our system in all cases.
The second major result of [23] gives conditions under which the system expressed in Equation (A1) with small diffusion rates has the same dynamics as the kinetic system
d U i d t = F i ( x , U 1 , , U n ) for i = 1 , , n .
The conditions can be stated as follows, noting that we have replaced the A used in [23] with L to avoid any confusion with equation numbers in this Appendix:
(L1) F i / U j 0 (i.e. the systems expressed in Equations (A1) and (A3) are cooperative).
(L2) For each x 0 Ω ¯ , the system expressed in Equation (A3) has a unique positive equilibrium α ( x 0 ) , which is globally asymptotically stable among positive solutions and is locally linearly stable, and α ( x ) depends continuously on x.
(L3) There is a δ 0 > 0 such that, for j = 1 , , n , F j ( x , U ) / U j > δ 0 for all x Ω ¯ provided 0 < U i δ 0 for i = 1 , , n .
(L4) There is a δ 0 , M > 0 such that for j = 1 , , n , F j ( x , U ) / U j < δ 0 for all x Ω ¯ provided U i M for i = 1 , , n .
The second result of [23] that we use is Proposition 5.2, which can be stated as follows:
Theorem A2.
(Proposition 5.2 of [23]) For any positive steady state w d of Equation (A1), we have that w α uniformly in Ω as m a x { d 1 , , d n } 0 , with the α given in Assumption (L2).

References

  1. Cantrell, R.S.; Cosner, C. Spatial Ecology via Reaction-Diffusion Equations; John Wiley and Sons: Chichester, UK, 2003. [Google Scholar]
  2. Cosner, C. Reaction-diffusion-advection models for the effects and evolution of dispersal. Discret. Cont. Dyn. Syst. A 2014, 34, 1701–1745. [Google Scholar] [CrossRef]
  3. He, X.; Ni, W.-M. 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]
  4. Dockery, J.; Hutson, V.; Mischaikow, K.; Pernarowski, M. The evolution of slow dispersal rates: A reaction diffusion model. J. Math. Biol. 1998, 37, 61–83. [Google Scholar] [CrossRef]
  5. Hastings, A. Can spatial variation alone lead to selection for dispersal? Theor. Popul. Biol. 1983, 24, 244–251. [Google Scholar] [CrossRef]
  6. Altenberg, L. Resolvent positive linear operators exhibit the reduction phenomenon. Proc. Natl. Acad. Sci. USA 2012, 109, 3705–3710. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Greenwood-Lee, J.; Taylor, P.D. The evolution of dispersal in spatially varying environments. Evol. Ecol. Res. 2001, 3, 649–665. [Google Scholar]
  8. Brown, K.J.; Zhang, Y. On a system of reaction–diffusion equations describing a population with two age groups. J. Math. Anal. Appl. 2003, 282, 444–452. [Google Scholar] [CrossRef] [Green Version]
  9. Cantrell, R.S.; Cosner, C.; Yu, X. Dynamics of populations with individual variation in dispersal on bounded domains. J. Biol. Dyn. 2018, 12, 288–317. [Google Scholar] [CrossRef] [Green Version]
  10. Cantrell, R.S.; Cosner, C.; Yu, X. Populations with individual variation in dispersal in heterogeneous environments: Dynamics and competition with simply diffusing populations. Sci. China Math. 2020, 63, 441–464. [Google Scholar] [CrossRef] [Green Version]
  11. Gurtin, M.E. A system of equations for age-dependent population diffusion. J. Theor. Biol. 1973, 40, 389–392. [Google Scholar] [CrossRef]
  12. Kunisch, K.; Schappacher, W.; Webb, G.F. Nonlinear age-dependent population dynamics with random diffusion. Comput. Math. Appl. 1985, 11, 155–173. [Google Scholar] [CrossRef] [Green Version]
  13. Langlais, M. A nonlinear problem in age-dependent population diffusion. SIAM J. Math. Anal. 1985, 16, 510–529. [Google Scholar] [CrossRef]
  14. Yan, S.; Guo, S. Dynamics of a Lotka-Volterra competition-diffusion model with stage structure and spatial heterogeneity. Discret. Contin. Dyn. Syst. B 2018, 23, 1559–1579. [Google Scholar] [CrossRef] [Green Version]
  15. López-Gómez, J.; Molina-Meyer, M. The maximum principle for cooperative weakly coupled elliptic systems and some applications. Differ. Int. Equ. 1994, 7, 383–398. [Google Scholar]
  16. Molina-Meyer, M. Global attractivity and singular perturbation for a class of nonlinear cooperative systems. J. Differ. Equ. 1996, 128, 347–378. [Google Scholar] [CrossRef] [Green Version]
  17. Cantrell, R.S.; Cosner, C. On the effects of spatial heterogeneity on the coexistence of competing species. J. Math. Biol. 1998, 37, 103–145. [Google Scholar] [CrossRef]
  18. Fernández-Rincón, S.; López-Gómez, J. A singular perturbation result in competition theory. J. Math. Anal. Appl. 2017, 445, 280–296. [Google Scholar] [CrossRef]
  19. Hutson, V.; López-Gómez, J.; Mischaikow, K.; Vickers, G. Limit Behavior for a Competing Speces Problem with Diffusion. In Dynamical Systems and Applications; Agarwal, R.P., Ed.; World Sci. Publishing: Hackensack, NJ, USA, 1995; pp. 343–358. [Google Scholar]
  20. Chen, S.; Shi, J. Asymptotic profiles of basic reproduction number for epidemic spreading in heterogeneous environment. arXiv 2019, arXiv:1909.10107. [Google Scholar]
  21. Magal, P.; Webb, G.F.; Wu, Y.-X. On the basic reproduction number of reaction-diffusion epidemic models. SIAM J. Appl. Math. 2019, 79, 284–304. [Google Scholar] [CrossRef]
  22. Hei, L.J.; Wu, J.H. Existence and stability of positive solutions for an elliptic cooperative system. Acta Math. Sin. Engl. Ser. 2005, 21, 1113–1120. [Google Scholar] [CrossRef]
  23. Lam, K.-Y.; Lou, Y. Asymptotic behavior of the principal eigenvalue for cooperative elliptic systems and applications. J. Dyn. Differ. Equ. 2016, 29, 29–48. [Google Scholar] [CrossRef]
  24. Cantrell, R.S.; Cosner, C.; Lou, Y. Evolution of dispersal and the ideal free distribution. Math. Biosci. Eng. 2010, 7, 17–36. [Google Scholar] [PubMed]
  25. Van den Driessche, P.; Watmough, J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2002, 108, 29–48. [Google Scholar] [CrossRef]
  26. Bhattacharyya, A. On a measure of divergence between two statistical populations defined by their probability distribution. Bull. Calcutta Math. Soc. 1943, 35, 99–110. [Google Scholar]
  27. Derpanis, K.G. The Bhattacharyya measure. Mendeley Comput. 2008, 1, 1990–1992. [Google Scholar]
  28. Smith, H. Monotone Dynamical Systems. An Introduction to the Theory of Competitive and Cooperative Systems; Mathematical Surveys and Monographs, 41; American Mathematical Society: Providence, RI, USA, 1995. [Google Scholar]
  29. Hutson, V.; Mischaikow, K.; Polácik, P. The evolution of dispersal rates in a heterogeneous time-periodic environment. J. Math. Biol. 2001, 43, 501–533. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Cantrell, R.S.; Cosner, C.; Martínez, S. Persistence for a Two-Stage Reaction-Diffusion System. Mathematics 2020, 8, 396. https://doi.org/10.3390/math8030396

AMA Style

Cantrell RS, Cosner C, Martínez S. Persistence for a Two-Stage Reaction-Diffusion System. Mathematics. 2020; 8(3):396. https://doi.org/10.3390/math8030396

Chicago/Turabian Style

Cantrell, Robert Stephen, Chris Cosner, and Salomé Martínez. 2020. "Persistence for a Two-Stage Reaction-Diffusion System" Mathematics 8, no. 3: 396. https://doi.org/10.3390/math8030396

APA Style

Cantrell, R. S., Cosner, C., & Martínez, S. (2020). Persistence for a Two-Stage Reaction-Diffusion System. Mathematics, 8(3), 396. https://doi.org/10.3390/math8030396

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