Abstract
A Kolmogorov-type competition model featuring allocation profiles, gain functions, and cost parameters is examined. For plant species that compete for sunlight according to the canopy partitioning model [R.R. Vance and A.L. Nevai, Plant population growth and competition in a light gradient: a mathematical model of canopy partitioning, J. Theor. Biol. 245 (2007), pp. 210–219] the allocation profiles describe vertical leaf placement, the gain functions represent rates of leaf photosynthesis at different heights, and the cost parameters signify the energetic expense of maintaining tall stems necessary for gaining a competitive advantage in the light gradient. The allocation profiles studied here, being supported on three alternating intervals, determine “interior” and “exterior” species. When the allocation profile of the interior species is a delta function (a big leaf) then either competitive exclusion or coexistence at a single globally attracting equilibrium point occurs. However, if the allocation profile of the interior species is piecewise continuous or a weighted sum of delta functions (multiple big leaves) then multiple coexistence states may also occur.
1. Introduction
Consider two species (i=1, 2) that obey the Kolmogorov-type competition model Citation21,Citation32
System (Equation1) describes the interactions of two clonal plant species that compete for sunlight by placing their leaves at multiple heights Citation21,Citation32. For species i, its total leaf area is denoted by x
i
, its vertical leaf area density with respect to height z (i.e. its vertical leaf profile) is indicated by s
i
(z), the rate of photosynthesis performed by its leaves when shaded from above by cumulative leaf area x is given by the energy gain function φ
i
(x), its leaf light extinction coefficient is σ
i
, its growth constant is γ
i
, and the energetic cost of maintaining its leaves, stems, and roots is represented by C
i
. The growth and cost parameters γ
i
and C
i
in turn depend on close to a dozen primary plant parameters such as tissue mass, mean stem height, and tissue metabolic rate. For a complete description of these parameters and functions and how they relate to Equation (Equation1
) see Citation32.
The application of Equation (Equation1) to light competition is not limited to terrestrial plants only. For example, Weissing and Huisman Citation34 proposed earlier a related system of equations to model phytoplankton that compete for sunlight in a lake. Both models are members of a more general class of systems
Furthermore, it is usually of great interest to establish the manner in which the allocation profiles and gain functions influence the outcome of competition in Equation (Equation1). Toward this end, it is known that
-
if two species have allocation profiles which are supported on non-overlapping intervals then either competitive exclusion or coexistence at a single globally attracting equilibrium point occurs Citation21;
-
if two species have rectangular (i.e. uniform) allocation profiles supported on overlapping but non-nested intervals, and the species share a common gain function, then either competitive exclusion or coexistence at a single globally attracting equilibrium point occurs Citation21;
-
there exist pairs of species possessing rectangular allocation profiles and distinct gain functions such that multiple non-trivial equilibrium points occur Citation19;
-
there exist pairs of species, one with a rectangular allocation profile surrounded by a bi-rectangular allocation profile belonging to the other, such that multiple non-trivial equilibrium points occur, even when the two species share the same gain function Citation19;
-
if two species possess rectangular allocation profiles, they share a common gain function, and there exists some τ>0 such that s 2(z)=s 1(z−τ) and C 2=C 1+bτ, then the competitive ability of species 2 diminishes from dominant to compatible to inferior as the allocation difference τ increases Citation22.
In summary, it has been found that when two plant species possess overlapping rectangular leaf profiles (as it occurs for many grasses, e.g. alkali seaheath (Frankenia grandifolia) and Virginia glasswort (Salicornia virginica)) then the degree to which their leaves overlap can have a marked influence on the outcome of competition. However, other kinds of leaf overlap are also possible; for example, some pairs of species interlace their leaves in an alternating fashion over multiple discrete heights, e.g. corn (Zea mays) and soybean (Glycine max). We can express the associated canopy patterns using a short-hand notation; for example, (2A–2B–A) denotes a five-leaf canopy in which the two lowest leaves belong to species A, the two middle leaves belong to species B, and the highest leaf belongs to species A. We seek the simplest such patterns for which coexistence can arise at multiple locally stable equilibrium points. As mentioned earlier, these discrete-leaf patterns cannot be of the form (kA–ℓB). The three-leaf pattern (A–B–A) is the simplest one that remains, and we remark that it also describes the simplest canopy in which the leaves of two species can interlace nontrivially. Of the simplest canopy patterns which remain are (A–2B–A) and (A–B–A–B), although we consider the former pattern to be less complex than the latter because it has fewer transitions between species.
In this article, we characterize the simplest canopy patterns that can produce coexistence at multiple equilibrium points for two theoretical species that compete for sunlight and obey model (Equation1).
We emphasize that many other factors are known to influence the outcome of plant competition, including seed dispersal and root competition for both soil nutrients and water Citation32. Here, we focus exclusively on whether competition for sunlight alone can produce multiple equilibrium points and we concentrate on examining the minimal complexity of interleaving patterns sufficient for such outcomes.
2. Asymptotic behaviour
Given any nonnegative initial condition, a nonnegative solution to Equation (Equation1) exists for all time, it is unique, and it always approaches some finite equilibrium within the closure of the first (or quadrant as
Citation13, pp. 144 and 250]). Thus, to determine the asymptotic behaviour of the dynamical system, it suffices to restrict our attention to nullclines and their with each other and with the coordinate axes. In general, the (non-trivial) nullcline
for species i is the set of all (x
1, x
2) in the closed first quadrant that satisfy the equation
So that we may state our results precisely, we must first make some preliminary definitions.
Recall from Citation32 that φ i (x) can be interpreted as the rate at which a unit area of leaves belonging to species i performs photosynthesis at a fixed height when shaded from above by cumulative leaf area x. This rate depends continuously on the overlying leaf area, it is always positive, it decreases as the overlying leaf area increases, and very large overlying leaf areas produce correspondingly very small rates of photosynthesis. Following Citation32, we define φ i (x) to be a gain function if it satisfies the following properties:
-
(P1) φ i is a continuously differentiable function of x≥0.
-
(P2) φ i (x)>0 for x≥0.
-
(P3) φ i ′(x)<0 for x≥0.
-
(P4) φ i (x)→0 as
.
To denote that φ i is a gain function, we will write φ i ∈Π. A special case arises when φ1 and φ2 both have the particular form
It will also be useful to introduce new mathematical notation so that we may later distinguish important leaf locations (such as the understory and overstory) within a canopy. Let denote the set of all allocation profiles, that is,
if and only if s(z) is nonnegative for z≥0 and
. Let Ω
s
denote the support of a given
. For
, we define the following parameterized subsets of
:
-
,
-
C(a, c)={s(z)∈M(a, c):s(z) is piecewise continuous and everywhere finite},
-
,
-
,
-
,
, and
,
-
,
, and
.
Some comments are in order. First, these definitions are not mutually exclusive. For example, . Second, if s(z)∈C(a, c) and
then
In applications of the model to plant competition for sunlight, allocation profiles in B(a) and B(a
1, a
2, …, a
m
) represent, respectively, the situation in which a plant species places its leaves exclusively at a single height or at m different heights. For allocation profiles in U(a) ⊕ O(c), the intervals (0, a) and represent the understory (U) and overstory (O), respectively.
3. Statement of the main results
Below we assume that:
-
the supports of the two species’ allocation profiles do not overlap, and
-
the support of the allocation profile belonging to species 2 surrounds the support of the allocation profile belonging to species 1 on both sides.
Our first main result concerns a situation in which s 1(z) is singular and s 2(z) is arbitrary but subject to conditions (i) and (ii).
Theorem 1 (U–B–O)
Suppose that s
1(z)∈B(a) and s
2(z)∈U(a) ⊕ O(a) with a>0, and let φ1, φ2∈Π and σ1, σ2, C
1, C
2>0 be given. Then at most one non-trivial equilibrium point exists, and when one does exist, it is always globally asymptotically stable with respect to the (open) positive quadrant. The same result also holds for s
2(z)∈M(0, a) or .
Our second main result concerns the case in which s 1(z) is piecewise continuous and finite and s 2(z) is bi-singular.
Theorem 2 (B–C–B)
Let s
1(z)∈C(a, c) and s
2(z)∈B(a, c) with and let σ1, σ2, C
1, C
2>0 be given. Then there exist φ1, φ2∈Π such that infinitely many non-trivial equilibrium points exist. Examples exist where the two nullclines coincide on an interval with the nonempty interior; there also exist examples with infinitely many locally asymptotically stable, and infinitely many unstable isolated equilibrium points. If S
2(c)C
1<C
2<C
1 then all this can even be arranged so that φ1 ≡ φ2.
The hypothesis of Theorem 2 excludes some pairs of profiles that are of type “B–M–B”. The next result partially fills in this gap by addressing situations in which s 1(z) is multi-singular.
Theorem 3 (B–mB–B)
Let s
1(z)∈B(z
1, z
2, …, z
m
) and s
2(z)∈B(a, c) with and let σ1, σ2, C
1, C
2>0 be given.
-
If S 1(z 1)<S 2(c) then there exist φ1, φ2∈Π such that infinitely many non-trivial equilibrium points exist. Examples exist where the two nullclines coincide on an interval with the nonempty interior; there also exist examples with infinitely many locally asymptotically stable, and infinitely many unstable isolated equilibrium points. If S 2(c)C 1<C 2<C 1 then all this can even be arranged so that φ1 ≡ φ2.
-
If S 1(z 1)≥S 2(c), then given any φ1, φ2∈Π at most one non-trivial equilibrium point exists, and when one does exist, it is always globally asymptotically stable with respect to the open first quadrant. This final result also holds for s 2(z)∈B(a) ⊕ O(c).
We remark that when m=1, then the assumption in part (b) of Theorem 3 is always satisfied. In this case, Theorem 3 is simply a special case of Theorem 1.
The final result concerns a special case of the situation occurring in Theorem 3 in that it restricts attention to cases in which s 1(z) and s 2(z) are both bi-singular. However, its conclusion is different because of its restriction to gain functions in parametric form.
Theorem 4 (B–2B–B)
Let and σ1, σ2, C
1>0 be given. Then there exist s
1(z)∈B(z
1, z
2), s
2(z)∈B(a, c),
and C
2>0, such that at least two non-trivial equilibrium points exist. In addition,
-
φ1 and φ2 may be chosen so that
for any given I in, J 1, J 2>0, even when C 2>0 is also given,
-
φ1 and φ2 may be chosen so that φ1≡φ2.
Note that neither part (a) nor part (b) of Theorem 4 implies the other. If , then these functions may still differ in the parameters ϕ
i
and π
i
appearing in Equation (Equation4
). Similarly, although the gain functions may be chosen to have the same parametric form, it may not be possible to specify values a priori for the parameters I
in, J
1, J
2, and C
2.
The biological implications of these results, as they pertain to plants that compete for sunlight are outlined in the Discussion.
4. Proof of Theorem 1 (U–B–O)
In this section, we prove Theorem 1. Suppose first that s
1(z)∈B(a) and s
2(z)∈U(a) ⊕ O(a) with a>0, and let φ1, φ2∈Π and σ1, σ2, C
1, C
2>0 be given. Then the species 1 nullcline satisfies
5. Proof of Theorem 2 (B–C–B)
In this section, we prove Theorem 2. Suppose that s
1(z)∈C(a, c) and s
2(z)∈B(a, c) with , and let σ1, σ2, C
1, C
2>0 be given. We will construct φ1, φ2∈Π such that Equation (Equation1
) has infinitely many equilibrium points.
Step 1: Getting started.
Define k=S 2(c)∈(0, 1).
-
(1a) Select Y 1>X 1>X 0>0 and Y 0>Y 1+(1/k−1)(X 1−X 0).
-
(1b) Define ωℓ=Y ℓ−X ℓ and Z ℓ=ωℓ+(1/k)X ℓ for ℓ∈{0, 1}.
Observe from the earlier-mentioned construction that (i) ω0>ω1, (ii) X ℓ<Y ℓ<Z ℓ for ℓ∈{0, 1}, and (iii) Z 0>Z 1. For γ∈[0, 1], define
-
(1c) X(γ)=X 0+γ(X 1−X 0) and Y(γ)=Y 0+γ(Y 1−Y 0).
-
(1d) ω(γ)=Y(γ)−X(γ) and Z(γ)=ω(γ)+(1/k)X(γ).
Observe that X(γ), Y(γ), ω(γ), and Z(γ) are all linear and they satisfy X(γ)∈[X 0, X 1], Y(γ)∈[Y 1, Y 0], ω(γ)∈[ω1, ω0], and Z(γ)∈[Z 1, Z 0] for γ∈[0, 1].
Step 2: Select φ1.
We construct φ1 as follows.
-
(2a) Define
and select
so that
-
(i) φ1(x)>0 and φ1′(x)<0 for
,
-
(ii)
for γ∈[0, 1].
(2b) Extend φ1(x) elsewhere so that φ1∈Π.
-
We remark that the operation in Equation (2a) is well-defined. To see this, first select a function ψ(x)∈C
1([X
0, X
1]) which satisfies (a) ((X
1−X
0)/(Y
1−Y
0))ψ(X
0)+C
1>0 and (b) ψ(x)>0 and ψ′(x)<0 for x∈[X
0, X
1]. Next, extend ψ to [Y
1, Y
0] by setting ψ(Y (γ))=((X
1−X
0)/(Y
1−Y
0))ψ(X(γ)) for γ∈[0, 1]. Extend ψ further to be a function in which satisfies
and ψ′(x)<0 for
. Finally, define φ1(x)=ψ(x)+C
1 for
. Observe from (1c) that
and therefore that
for γ∈[0, 1]. Since φ1(Y
0)>0, and the integral in part (ii) represents the average value of φ1 on [X(γ), Y(γ)], it is straightforward to see that φ1 satisfies the conditions in (2a).
Step 3: Select φ2.
We construct φ2 as follows. First, select λ0>λ1>μ1>μ0>0 so that kλℓ+(1−k)μℓ=C 2 for ℓ∈{0, 1} and define φ2(X(γ))=(1−γ)λ0+γλ1 and φ2(Z(γ))=(1−γ)μ0+γμ1 for γ∈[0, 1]. Next,
-
(3a) extend φ2(x) elsewhere so that φ2∈Π.
It is straightforward to show that
-
(3b) kφ2(X(γ))+(1−k)φ2(Z(γ))=C 2 for γ∈[0, 1].
5.1. An interval of equilibrium points
Observe from steps 1 and 2 with γ=0 that
5.2. Isolated equilibrium points
We now modify the earlier-mentioned construction to get infinitely many locally asymptotically stable, and infinitely many unstable isolated equilibrium points. First, select an infinite sequence γ1, γ2, … satisfying 0<γ1<γ2<⋯<1. Then
5.3. The same gain function
To complete the proof, it remains only to show that if C 2∈(k C 1, C 1) then we can modify the earlier-mentioned construction above so that φ1 ≡ φ2. Recall that φ1 was subject to a certain integral condition on [X 0, Y 0] in step 2 and that φ2 was subject to a linear condition on [X 0, X 1] and [Z 1, Z 0] in step 3. Below, we arrange for the intervals [X 0, Y 0] and [Z 1, Z 0] to be disjoint so that we will have sufficient freedom to make φ1 ≡ φ2. First, we repeat step 1 but replace (1a) with
-
(1a’) Select Y 1>X 1>X 0>0 so that X 1−X 0<((1−k)/2)X 1.
-
(1a”) Select θ∈(1/k, 2/k) and define Y 0=Y 1+θ(X 1−X 0).
These modifications allow us to directly compare Y 0 and Z 1:
6. Proof of Theorem 3 (B–mB–B)
In this section, we prove Theorem 3 using an approach that resembles the one used in the proof of Theorem 2. However, since both allocation profiles considered here are discrete, we will need notation and arguments that are very different.
Suppose that s
1(z)∈B(z
1, z
2, …, z
m
) and s
2(z)∈B(a, c) with , and let σ1, σ2, C
1, C
2>0 be given.
6.1. Proof of Theorem 3 part (a)
To prove part (a), we will again construct φ1, φ2∈Π such that Equation (Equation1) has infinitely many equilibrium points. Let [m]={1, 2, …, m}. Recall that s
1(z) has the form
The construction begins as follows. First, select , ω1>0, and
Lemma 1
There exist φ1, φ2∈Π such that
Figure 1. The relative magnitudes of certain quantities appearing in the proof of Theorem 3. Some but not all quantities can be directly compared. Quantities appearing higher are provably larger than quantities that appear lower when they are path-connected by upward lines or ellipses. For example, is greater than X* and any quantity on the bottom two diagonals but less than any other quantity on the upper m diagonals. However, the relative position of X
2(γ) with respect to the quantities in the column above
is unknown.
![Figure 1. The relative magnitudes of certain quantities appearing in the proof of Theorem 3. Some but not all quantities can be directly compared. Quantities appearing higher are provably larger than quantities that appear lower when they are path-connected by upward lines or ellipses. For example, is greater than X* and any quantity on the bottom two diagonals but less than any other quantity on the upper m diagonals. However, the relative position of X 2(γ) with respect to the quantities in the column above is unknown.](/cms/asset/2fc8293d-3205-4d81-861f-72ce90ecadba/tjbd_a_385173_o_f0001g.gif)
6.1.1. An interval of equilibrium points
Observe from Lemma 1 with γ=0 that
6.1.2. Isolated equilibrium points
The procedure to obtain infinitely many locally asymptotically stable, and infinitely many unstable isolated equilibrium points is completely analogous to the construction given in the proof of Theorem 2. First, select an infinite sequence γ1, γ2, … satisfying 0<γ1<γ2<⋯<1. Then modify φ2 on the interval so that
will cross
at (u(γ
j
), v(γ
j
)) from below to above when j is even and from above to below when j is odd. We omit the details.
6.1.3. The same gain function
To complete the proof of part (a), it remains only to show that if C 2∈(k C 1, C 1) then we can modify the earlier-mentioned construction so that φ1 ≡ φ2. The proof of the next lemma appears in Appendix A.
Lemma 2
If the following conditions are satisfied
then we can choose φ1 and φ2 in Lemma 1 so that φ1 ≡ φ2.
Observe that we can always choose ,
, and ω1
from the outset so that conditions (b) and (c) of Lemma 2 are satisfied. As in the construction presented in the proof of Theorem 2, conditions (A) and (B) in the proof of Lemma 2 give us sufficient freedom to choose φ on the intervals
and
in such a way that the two nullclines coincide along an interval with a nonempty interior, or meet at infinitely many locally asymptotically stable, and infinitely many unstable isolated equilibrium points. This completes the proof of part (a).
6.2. Proof of Theorem 3 part (b)
To prove part (b), let us assume that φ1, φ2∈Π are given. Using the same notation as in the proof of part (a), we need only to show that if h
1≥k then at most one equilibrium point exists, and that when one does exist it is always globally asymptotically stable with respect to the open first quadrant. First, we argue by contradiction to show that the slope of cannot exceed −(σ1
h
1/σ2
k). For suppose otherwise, i.e. that
. It follows from implicit differentiation of Equation (Equation13
) with respect to x
1 that
7. Proof of Theorem 4 (B–2B–B)
In this section, we prove Theorem 4. Let and σ1, σ2, C
1>0 be given.
7.1. Proof of Theorem 4 part (a)
To prove part (a), we fix C
2, I
in, J
1, J
2>0 and seek s
1(z)∈B(z
1, z
2), s
2(z)∈B(a, c), and such that Equation (Equation1
) has at least two equilibrium points. For ε>0, let
Lemma 3
There exist 0<ε≪1 and such that
Let s
1(z)=ε2δ(z−z
1)+(1−ε2)δ(z−z
2) and s
2(z)=2εδ(z−c)+(1−2ε)δ(z−a). Then s
1(z)∈B(z
1, z
2) and s
2(z)∈B(a, c), and equations for the species nullclines and
can be obtained from Equation (Equation3
). That is,
satisfies
and
satisfies
. In view of Lemma 3 and Equation (Equation16
), let u
1=(σ1)−1, u
2=(εσ1)−1, v
1=(2εσ2)−1, and v
2=(2σ2)−1. Then the two species coexist at the two equilibrium points (u
1, v
1) and (u
2, v
2) within the interior of the first quadrant. Other equilibrium points may also exist. This completes the proof of part (a).
7.2. Proof of Theorem 4 part (b)
To prove part (b), we will select s
1(z)∈B(z
1, z
2), s
2(z)∈B(a, c), , and C
2>0 such that φ1 ≡ φ2 and at least two equilibrium points exist. For ζ and ε, define
Lemma 4
There exist 0<ε≪1, C
2>0, and
such that
Let and s
2(z)=(1−ε)δ(z−c)+εδ(z−a). Then s
1(z)∈B(z
1, z
2) and s
2(z)∈B(a, c), and equations for the species nullclines
and
can be obtained from Equation (Equation3
). That is,
satisfies
8. Discussion
In analytical studies of terrestrial plant competition, it is usually assumed that model species place their leaves at a single height Citation23,Citation25,Citation27,Citation35,Citation36 or, if a coarse distinction between sun leaves and shade leaves is to be made, then at two heights (see Citation6,Citation30). Such simplified approaches involving big-leaf species are appropriate when studying multiple factors in combination that together influence plant competition (e.g. light-limitation, nutrient-limitation, and dispersal) because they can make otherwise intractable problems more amenable to mathematical methods. Yet, situations remain in which sunlight appears to be the limiting resource for plants in competition Citation10,Citation11,Citation28,Citation33,Citation34. In such cases, we consider it appropriate to include a fuller account of the competitors’ vertical leaf profiles (or VLPs) because higher leaves everywhere receive more sunlight energy than lower leaves and therefore can perform photosynthesis at higher rates Citation1,Citation9,Citation11. Furthermore, higher leaves must be supported by taller stems, which are energetically more expensive to construct and to maintain than shorter stems Citation24.
The mechanistic Canopy Partitioning Model Citation21,Citation32 can be used to predict the outcome of competition between two plant species with fixed VLPs that occupy a relatively constant environment. The defining feature of the model is that it accommodates arbitrary shapes for the competitors’ VLPs. It can therefore capture in real time subtle (and usually nonlinear) changes in leaf photosynthesis rates at different heights, which result when plant populations change in leaf area over time. Although initial applications of the model focused on species possessing VLPs that are constant functions over a height interval (rectangular profiles) or a pair of disjoint height intervals (bi-rectangular profiles) Citation19,Citation21,Citation22, it can equally well accommodate species with big-leaves positioned at multiple discrete heights in a canopy. The purpose of this study is to examine this latter situation.
In general, a central problem is to determine how many non-trivial equilibrium points the Kolmogorov-type competition model (Equation1) can have and to establish their stability or instability. Earlier work by Nevai and Vance Citation21,Citation22 delineated conditions under which at most one such equilibrium can exist and demonstrated that if this unique equilibrium exists for rectangular VLPs, it must globally attract the interior of the first quadrant when the species share the same photosynthesis gain function. In Citation19 it was shown that there exist pairs of species with rectangular VLPs and distinct gain functions for which any finite number of stable and unstable equilibrium points are possible. It was also shown that at least two equilibrium points are possible even if the species share the same gain function, provided that one of their VLPs is bi-rectangular and surrounds the other. However, no non-trivial examples of plant species pairs that obey (Equation1
) and admit infinitely many equilibrium points were previously known; neither were there known examples of such systems with multiple equilibrium points in which both gain functions are in parametric form (Equation4
) (see Figure ).
Figure 2. Allocation profiles that are arbitrary on nontrivial intervals (rounded), piecewise constant on nontrivial intervals (rectangular), or with finite support (thick horizontal lines). In the first row, the allocation profiles are assumed to be the same. If ∀ appears in the upper left corner of a cell then the number of nontrivial stable (m) and unstable (n) equilibrium points satisfies the given conditions for all pairs of model species. However, if ∃ appears in the corner of a cell then there exist pairs of model species such that m and n both satisfy the given conditions. Degenerate situations under which two species nullclines coincide everywhere are excluded. A [p] or (q) in the lower right corner indicates that the result appears in citation [p] or as Theorem q in this article, respectively.
![Figure 2. Allocation profiles that are arbitrary on nontrivial intervals (rounded), piecewise constant on nontrivial intervals (rectangular), or with finite support (thick horizontal lines). In the first row, the allocation profiles are assumed to be the same. If ∀ appears in the upper left corner of a cell then the number of nontrivial stable (m) and unstable (n) equilibrium points satisfies the given conditions for all pairs of model species. However, if ∃ appears in the corner of a cell then there exist pairs of model species such that m and n both satisfy the given conditions. Degenerate situations under which two species nullclines coincide everywhere are excluded. A [p] or (q) in the lower right corner indicates that the result appears in citation [p] or as Theorem q in this article, respectively.](/cms/asset/70e1327d-1ee7-4ba4-8458-a42666cbe9c3/tjbd_a_385173_o_f0002g.gif)
All of our theorems here concern cases in which at least one species has a VLP composed of one or more big-leaves (i.e. with finite support) and the canopy is partitioned into three height zones: an upperstory consisting only of leaves belonging to species 2, a middlestory consisting only of leaves belonging to species 1, and an understory consisting only of leaves belonging to species 2. Thus, we may refer to species 1 as the interior species and to species 2 as the exterior species. We would like to point out that this framework applies equally well to phytoplankton that compete for sunlight in lakes according to the well-known model of Weissing and Huisman Citation34.
We show in Theorem 1 that when the interior species has a single big-leaf and it is surrounded from below and above by an arbitrary VLP belonging to the exterior species, then there are only two long-term possibilities: either competitive exclusion results with the identity of the surviving species depending only on parameter values but not on initial plant population abundances, or else competitive coexistence occurs at a single globally attracting equilibrium point. Below, we will refer to either of these competitive outcomes as being of “type I.” We remark that type I outcomes do not admit the possibility of competitive exclusion owing to founder control. These same comments also apply when one species has a single big-leaf lying entirely above or below an arbitrary VLP belonging to the other, which complements a similar result obtained in Citation21 for two species with uniform but non-overlapping VLPs.
The situation changes markedly when the VLP of the interior species consists of multiple big-leaves situated at different heights.
Theorem 2 shows that there exist many pairs of model species that can experience both type I competitive outcomes and also coexistence at any one of infinitely many equilibrium points (“type II”). In particular, one may specify a priori a VLP for the interior species which is not concentrated at any single point, a VLP for the exterior species which consists of two big-leaves (one lying in the overstory and one in the understory), and arbitrary cost parameters. Then corresponding species can be found whose nullclines coincide at infinitely many points. This remarkable conclusion continues to hold even when it is required that both species share the same photosynthesis gain function. The actual construction of this gain function is performed in the proof of Theorem 2. The numerous basins of attraction that are possible in Theorem 2 constitute a form of founder control in which different initial species abundances lead to different equilibrium abundances without necessarily resulting in competitive exclusion for either species.
Theorem 3 concerns canopies in which the interior species has an arbitrary but finite number of big-leaves and the exterior species still has a single big-leaf in both the overstory and understory. Again, there exist pairs of species for which type II competitive outcomes occur. However, in these cases, the uppermost leaf belonging to the interior species must bear a smaller fraction of the plant’s total leaf area than the highest leaf of the exterior species does relative to its own total leaf area. If not, then we can deduce from Theorem 3 (b) that no combination of parameters and functions exist for which two species can coexist at multiple equilibrium points and the competitive outcome will always be of type I. The distinction that arises between the results in parts (a) and (b) warrants corroboration by empirical results.
Finally, Theorem 4 implies that if the interior species has exactly two big-leaves and they are surrounded on either side by a single big-leaf belonging to the exterior species, then parametric photosynthesis gain functions which produce coexistence at multiple equilibrium points can be chosen. The shared parameter I in represents the constant irradiance striking the upper surface of the canopy, and J 1 and J 2 are half-saturation constants Citation31 of the species. The proof of Theorem 4 uses the Intermediate Value Theorem in both parts. In the proof of part (a), we carefully choose constants and functions so that the two gain functions are in parametric form for specified values of I in, J 1, and J 2. Our construction involves two species that allocate small fractions (ε2 and 2ε) of their total leaf area to their upper leaves. When ε is very small, these two upper leaves each carry very low leaf areas. Thus, the lower leaf of the interior species experiences nearly full sunlight and is capable of casting substantial shade on the understory leaf belonging to the exterior species. In the proof of part (b), we have more freedom to obtain the parameters and functions in our construction, but we need to ensure that the two gain functions coincide. We accomplish this by letting the interior species be fairly balanced in its allocation of leaf area to its two leaves and then letting the proportion of leaf area in the overstory (which belongs to the exterior species) tend to 1.
Our results open up several new avenues for future research. In particular, it is an open question to determine which of our results have analogs for competitive systems in which both VLPs are continuous on their supports. For example, if the VLP for the interior species is rectangular and the VLP for the exterior species is bi-rectangular and surrounds the first, then can there be infinitely many equilibrium points? Can one obtain multiple equilibrium points for any non-trivial pair of VLPs, continuous on their supports, when both photosynthesis gain functions are required to be in parametric form?
Acknowledgements
The authors thank the reviewers for valuable comments and the Mathematical Biosciences Institute at The Ohio State University for support of this project. This study is based upon work supported by the National Science Foundation under Agreement No. 0112050.
Additional information
Notes on contributors
Andrew L. Nevai
Present address: Department of Mathematics, University of Central Florida, Orlando, FL 32816, USA.References
- Anten , N. P.R. and Hirose , T. 2001 . Limitations on photosynthesis of competing individuals in stands and the consequences for canopy structure . Oecologia , 129 : 186 – 196 .
- Armstrong , R. A. and McGehee , R. 1980 . Competitive exclusion . Am. Nat. , 115 : 151 – 170 .
- Butler , G. J. , Freedman , H. I. and Waltman , P. 1986 . Uniformly persistent systems . Proc. Amer. Math. Soc. , 96 : 425 – 430 .
- Cantrell , R. S. , Cosner , C. and Hutson , V. 1993 . Permanence in ecological systems with spatial heterogeneity . Proc. Roy. Soc. Edinburgh , 123A : 533 – 559 .
- Cantrell , R. S. , Cosner , C. and Hutson , V. 1996 . Ecological models, permanence and spatial heterogeneity . Rocky Mt. J. Math. , 26 : 1 – 35 .
- Pury , D. G.G. and Farquhar , G. D. 1997 . Simple scaling of photosynthesis from leaves to canopies without the errors of big-leaf models . Plant Cell Environ. , 20 : 537 – 557 .
- Freedman , H. I. and Ruan , S. G. 1995 . Uniform persistence in functional differential equations . J. Differ. Equations , 115 : 173 – 192 .
- Hardin , G. 1960 . The competitive exclusion principle . Science , 131 : 1292 – 1297 .
- Hikosaka , K. , Sudoh , S. and Hirose , T. 1999 . Light acquisition and use by individuals competing in a dense stand of an annual herb, Xanthium canadense . Oecologia , 118 : 388 – 396 .
- Hirose , T. 2005 . Development of the Monsi–Saeki theory on canopy structure and function . Ann. Bot. , 95 : 483 – 494 .
- Hirose , T. and Werger , M. J.A. 1995 . Canopy structure and photon flux partitioning among species in a herbaceous plant community . Ecology , 76 : 466 – 474 .
- Hirsch , M. W. 1988 . Systems of differential equations which are competitive or cooperative. III. Competing species . Nonlinearity , 1 : 51 – 71 .
- Hirsch , M. W. , Smale , S. and Devaney , R. L. 2004 . Differential Equations, Dynamical Systems, and an Introduction to Chaos , San Diego, CA, , USA : Academic Press .
- Hofbauer , J. and Garay , B. 2003 . Robust permanence for ecological differential equations, minimax, and discretizations . SIAM J. Math. Anal. , 34 : 1007 – 1039 .
- Hofbauer , F. , Hofbauer , J. , Raith , P. and Steinberger , T. 2004 . Intermingled basins in a two species system . J. Math. Biol. , 49 : 293 – 309 .
- Hofbauer , J. and Schreiber , S. J. 2004 . To persist or not to persist? . Nonlinearity , 17 : 1393 – 1406 .
- Hutson , V. , Lou , Y. and Mischaikow , K. 2005 . Convergence in competition models with small diffusion coefficients . J. Differ. Equations , 211 : 135 – 161 .
- Hutson , V. and Schmitt , K. 1992 . Permanence and the dynamics of biological systems . Math. Biosci. , 111 : 1 – 71 .
- Just , W. and Nevai , A. L. 2008 . A Kolmogorov-type competition model with multiple coexistence states and its applications to plant competition for sunlight . J. Math. Anal. Appl. , 348 : 620 – 636 .
- Kolmogorov , A. 1936 . Sulla teoria di Volterra della lotta per l'esistenza . Giorn. Ist. ltal. Attuari , 7 : 74 – 80 .
- Nevai , A. L. and Vance , R. R. 2007 . Plant interspecies competition for sunlight: a mathematical model of canopy partitioning . J. Math. Biol. , 55 : 105 – 145 .
- Nevai , A. L. and Vance , R. R. 2008 . The role of leaf height in plant interspecies competition for sunlight: analysis of a canopy partitioning model . Math. Biosci. Eng. , 5 : 101 – 124 .
- Perry , L. G. , Neuhauser , C. and Galatowitsch , S. M. 2003 . Founder control and coexistence in a simple model of asymmetric competition for light . J. Theor. Biol. , 222 : 425 – 436 .
- Pronk , T. E. , Schieving , F. , Anten , N. P.R. and Werger , M. J.A. 2007 . Plants that differ in height investment can coexist if they are distributing non-uniformly within an area . Ecol. Complexity , 4 : 182 – 191 .
- Rees , M. and Bergelson , J. 1997 . Asymmetric light competition and founder control in plant communities . J. Theor. Biol. , 184 : 353 – 358 .
- Rescigno , A. and Richardson , I. W. 1967 . The struggle for life. I. Two species . Bull. Math. Biophys. , 29 : 377 – 388 .
- Reynolds , H. L. and Pacala , S. W. 1993 . An analytical treatment of root-to-shoot ratio and plant competition for soil nutrient and light . Am. Nat. , 141 : 51 – 70 .
- Takada , T. and Iwasa , Y. 1986 . Size distribution dynamics of plants with interaction by shading . Ecol. Mod. , 33 : 173 – 184 .
- Tang , B. and Kuang , Y. 1996 . Permanence in Kolmogorov-type systems of nonautonomous functional-differential equations . J. Math. Anal. Appl. , 197 : 427 – 447 .
- Thornley , J. H.M. 2002 . Instantaneous canopy photosynthesis: analytical expressions for sun and shade leaves based on exponential light decay down the canopy and an acclimated non-rectangular hyperbola for leaf photosynthesis . Ann. Bot. , 89 : 451 – 458 .
- Thornley , J. H.M. and Johnson , I. 1990 . Plant and Crop Modelling – A Mathematical Approach to Plant and Crop Physiology , Oxford, , UK : Clarendon Press .
- Vance , R. R. and Nevai , A. L. 2007 . Plant population growth and competition in a light gradient: a mathematical model of canopy partitioning . J. Theor. Biol. , 245 : 210 – 219 .
- Weihe , P. E. and Neely , R. K. 1997 . The effects of shading on competition between purple loosestrife and broad-leaved cattail . Aquat. Bot. , 59 : 127 – 138 .
- Weissing , F. and Huisman , J. 1994 . Growth and competition in a light gradient . J. Theor. Biol. , 168 : 323 – 336 .
- Xiao , S. , Chen , S.-Y. and Wang , G. 2006 . An ESS for the height of a plant population, or an optimal height for an individual? – Rethinking game-theoretic models for plant height . Bull. Math. Biol. , 68 : 957 – 967 .
- Zavalaa , M. A. and Bravo de la Parra , R. 2005 . A mechanistic model of tree competition and facilitation for Mediterranean forests: scaling from leaf physiology to stand dynamics . Ecol. Mod. , 188 : 76 – 92 .
Appendix
Proof of Lemma 1
We first construct φ1 as follows.
-
(2a) Define α=h 1(h 1−h*)<0 and
.
-
(2b) Select ε>0 so that
.
-
(2c) For δ>0 define
-
(2d) Select 0<δ*≪1 so that
.
Since φ1(X*−ε;δ*)>C 1>φ1(X*+ε;δ*), we can
-
(2e) Extend φ1(x;δ*) to a function φ1∈Π defined for x≥0.
This construction implies then (with δ* now omitted) that for γ∈[0, 1]. To see that this equation holds, apply (vii) and (ix) and then (2c) and (2a) to get
We now construct φ2 as follows.
-
(3a) Select
so that
for ℓ∈{0, 1}.
-
(3b) Define φ2(X 0(γ))=(1−γ)λ0+γλ1 and
for γ∈[0, 1].
-
(3c) Extend φ2(x) elsewhere so that φ2∈Π.
It is straightforward to show that for γ∈[0, 1].
Proof of Lemma 2
It follows from condition (b) that
-
.
-
.
-
for γ∈[0, 1].
Since C 2∈(k C 1, C 1), these conditions can easily be satisfied if δ* in (2d) is chosen to be sufficiently small.
Proof of Lemma 3
Step 1: Select φ1.
To begin, let 0<ε≪1 and define X 1=1+ε2, Y 1=2ε, X 2=2, and Y 2=1/ε+ε as in the main proof. For ν>1 and γ>0, also define
Here, θ is taken to be the positive constant I in/J 1. Then
Step 2: Select φ2.
The argument to select φ2 is similar to the one used to select φ1, but it requires the use of different constants. With ε as mentioned earlier, let X
0=1, Y
0=ε, X
3=1/2ε+1, and as in the main proof. Also, for ν>1 and γ>0 (not necessarily the same as before) define
Proof of Lemma 4
Although the following argument uses similar notation as in the proof of Lemma 3, it relies on a different set of techniques. For and 0<ε≪1, consider the quantities defined in Equation (Equation17
). Since 0<ε≪1, we have (i) X
3 ≈ 4, Y
0 ≈ 0, Y
1 ≈ 2−4ζ, Y
2 ≈ 4, and Y
3=Y
2+ε3/(1−ε) ≈ 4, (ii) Y
0<X
0<X
3<Y
3, (iii) Y
1<X
1<X
2<Y
2, (iv) X
0−Y
0=1−ε2, (v)
, (vi) Y
2−X
2=2ε/1−ε+ε2 ≈ 2ε, and (vii) Y
3−X
3=ε+ε2/(1−ε) ≈ ε. For γ>0 and θ>0, let
We now show that there exists some such that
. For the moment, fix 0<δ≪1, let
, and let ε tend to zero. Then
, and
We now turn to the pair of functions M(ε, γ, ζ) and F(ε, γ, ζ, θ). Recall that
It follows from Equations (EquationA5) to (EquationA8
) that for sufficiently small ε we have both
and
. Hence, for sufficiently large θ we also have
and
. Since
depends continuously on ζ, there exists some
such that
. Thus, we have shown that there exist ε,
, and
such that
Finally, with , we observe from Equations (EquationA3
) and (EquationA9
) that