1. Introduction
The standard one-dimensional Ricker equation refers to the discrete population model
with parameters
r,
a and an initial condition
. This model is named after Ricker’s work fitting the parameters to reproductive curves for fish populations, including the Pacific Herring, Pink Salmon, Coho Salmon, and Sockeye Salmon, among others [
1]. From a mathematical perspective, an increase in
r leads to period doubling bifurcations and eventually chaos [
2]. In a controlled laboratory, such bifurcations were observed in Tribolium flour beetles [
3]. However, as noted in [
4], ecologists have rarely observed chaos in animal populations in nature. This suggested that, in a real-world setting, one-dimensional models such as the Ricker model are missing important factors.
In a series of papers [
5,
6,
7,
8,
9,
10,
11], Jim Cushing, along with some coauthors, established the foundations of discrete evolutionary models using the idea of Darwinian evolution strategy [
12] based on evolutionary game theory [
13]. These evolutionary game models differ in significant ways from classical theory, while classical games focus on strategies that optimize players’ payoffs, evolutionary games emphasize strategies that persist over time. Through births and deaths, players come and go, but their strategies are passed down from generation to generation [
14].
The first appearance of evolutionary models in discrete dynamical systems was in [
6], applied to the classic Beverton–Holt (discrete logistic) model. Using standard methods for modeling evolution, the one-dimensional model becomes a planar system of coupled difference equations, consisting of a Beverton–Holt type equation for the population dynamics and a difference equation for the dynamics of the mean phenotype trait. By studying the case when the trait equation uncouples from the population dynamic equation, Cushing was able to obtain criteria under which the evolutionary system has globally asymptotically stable equilibria or periodic solutions.
Recently, Ch-Chaoui and Moknni [
15] studied the coupled case of the evolutionary Beverton–Holt model and discussed the existence of the positive fixed point, exploring its asymptotic stability. Under certain parametric conditions, they discussed the Neimark–Sacker bifurcation, a phenomenon that does not occur in the classic Beverton–Holt model.
This strategy of deriving an evolutionary model was extended to more general discrete models [
7,
9,
11]. In particular, models with the Allee effect [
10], models of Ricker type [
8,
16,
17] and multi-species models [
17,
18,
19].
In all these studies, the main mathematical goals are to find and classify the stability of equilibria, identify bifurcations, and when possible, describe the global dynamics of the models. In the case of a lone locally stable equilibrium, a great deal of studying has been performed on whether this fixed point is also globally attracting. While there exist many general theorems on global stability in one-dimensional mappings [
20,
21,
22], this problem is more difficult in higher dimensions. Some successful examples of establishing global stability in higher dimensions include of those of monotone and mixed-monotone systems [
19,
23,
24,
25,
26,
27]. For non-monotone systems, singularity curves [
28,
29,
30,
31] have also been used to establish global stability if the mapping satisfies certain topological constraints [
32,
33,
34,
35]; these curves will also play an important role in the present work. However, in full generality, the question of global stability even in the planar case is far from understood. Some open problems and conjectures in two-dimensions may be found in [
36].
In a recent paper [
37], the author proposed a series of open problems and conjectures concerning the global stability of the following evolutionary discrete model of Ricker type for a single-species
x (with a single mean trait
u and individual trait
v)
where the parameters
are the sequence of growth rates (which are functions of the trait),
is the competition coefficient,
is the speed of evolution, and
is the sensitivity of the competition. A complete derivation and explanation of these parameters may be found in [
8].
Notice that if
then Model (
1) is autonomous whereas if
then Model (
1) is
-periodic.
We remark that the first Equation in (
1) is known in the literature as the inherent fertility equation while the second equation is the trait equation. Other names for the trait equation include the canonical equation for evolution [
38], Lande’s Equation [
39,
40] and Fisher’s Equation [
41].
In the autonomous case, for the parameter values
, Model (
1) is mixed-monotone, and the following result addresses the global stability of the model:
Theorem 1 ([
19])
. Let for all , such that , , and . Then, the interior equilibrium point of Model (1) is globally asymptotically stable if it is locally asymptotically stable, where . The conditions of local stability required in Theorem 1 are
and
with
,
and
as stated in Theorem 2.1 in [
19]. In particular, Theorem 1 shows that local stability implies global stability in the parameter range
. In this paper, we address the case of
, where the mapping is no longer mixed-monotone and new techniques are needed.
In this paper, we solve one of the conjectures presented in [
37]; specifically, we demonstrate the global stability of the interior fixed point of Model (
1) when
and
, thereby completing Theorem 1 in the global stability approach. Furthermore, we go beyond the concept of global stability by characterizing the invariant manifolds of the fixed points and also the bifurcation of the interior fixed point. As is seen in populations in nature, we will prove that this model does not undergo a period doubling cascade to chaos.
In order to achieve these goals, in
Section 2 we make a change in the variables in Model (
1) and write a simplified model that will help us understand both the local and global properties of the original model.
Using the concept of critical curves [
29], in
Section 3 we identify a forward-invariant and globally attracting subset under the action of the system’s mapping. Subsequently, in
Section 4 and
Section 5, we analyze the stability and bifurcation of the fixed points within this invariant set. Finally, in the last two sections, we provide the conclusion and discuss future work.
We believe that the techniques followed in this paper will aid in the study of similar models, namely, the use of critical curves and isoclines for identifying the maximal invariant set where all the dynamics take place; the characterization of the invariant (stable and unstable) manifolds that attract all orbits in the invariant set; and the classification of the Neimark–Sacker bifurcation using a geometric argument, which, to our knowledge, is the first time that a supercritical case has been classified via such an argument.
2. Notation and Preliminaries
By using a change in variables in Model (
1), it is possible to reduce the number of parameters involved, thereby simplifying the study.
Let us consider the autonomous system
equivalent to Model (
1), for
,
,
, and
.
Let
and
. Then,
and
Hence, the properties of Model (
2) are the same as the following system
Let
F be the mapping representing Model (
3) given by
where
for
.
In the following sections, we will study the dynamics of invariant sets under the mapping
F representing Model (
3), as well as some properties of the fixed points of
F.
3. Critical Curves and Invariant Sets
We start by examining the preimages of
F representing Model (
3). Initially, we follow the work of [
29]. They assume there is a mapping
F that is a noninvertible planar mapping whose preimage at any point consists of a finite number of isolated points. They then divide the plane into regions
, where
if and only if
has
k preimages under
F. Following their notation, the boundaries of these regions is denoted by
, with preimage
. For smooth mappings, they show that
can be determined by the set of points where the determinant of the Jacobian of
F vanishes, in which case
.
For the map
F representing Model (
3), we will show that
separates the first and fourth quadrant into two distinct regions
and
, as is illustrated in
Figure 1. All points in the interior of the shaded region have two distinct preimages while points in the unshaded region have no preimage, as will be shown in Theorem 2.
The Jacobian of the mapping
F is given by
Then,
and this determinant vanishes along the graph of
with
. Consequently,
and
. We now look at the geometry of
and
in more detail.
Lemma 1. Let , and let . Then the graph of is a monotone increasing and concave curve. It has a vertical asymptote at and a horizontal asymptote at .
Proof. The first claim follows by calculating derivatives of . The asymptotes follow by taking the appropriate limits. □
Lemma 2. Let , and let . Then, the curve can be decomposed into curves from to and from to where is strictly increasing and is strictly decreasing. At the intersection point , has a vertical tangent line.
Proof. The curve
has parametrization
for
. We claim that
consists of the portion of
for
and that
consists of the part with
. We calculate
and
It is clear that
for all
t. Also, for
, both
and
, so
for
. Similarly, for
, both
and
, so
for
. Since
vanishes at
, the curve has a vertical tangent line at
. □
We remark that the relative positions of the critical curves as well their properties, as described in Lemmas 1 and 2, are clearly illustrated in
Figure 2. As we can observe, the critical curves do not intersect, and there are points on the left side of
with the same
y component that have two distinct preimages symmetric about a point on the critical curve
.
Before we proceed with the analysis of our model, let us take an overall look at a concrete example:
Example 1. Let , and . Then, the critical curve and a parametrization of is written asClearly, the conclusion of the precedent theorem is satisfied since andA representation of the graphs of and is shown in Figure 2. As one can observe the curve, divides the first and fourth quadrant into two distinct regions. Now, let , , , and , all taken with the same y-coordinate, where is the intersection of the line and the curve . To calculate , we use the parametrization and solve , yielding . Then, . A forward computation shows that the first three points have two preimages as followswhile D has a single preimage on given approximately byAll of the preimages fall on the line , marked as L in Figure 2. The precedent example motivates the following result:
Theorem 2. Let , and let . Then, the curve divides the first and fourth quadrant into two regions. For the region adjacent to the y-axis, each point with positive x coordinate has two distinct preimages. The other region is not in the range of F.
Proof. The conclusion that
divides the plane into two regions follows from Lemma 2. To prove the remaining claim, we consider the set of lines with slope
, whose image we will show is a horizontal line. Let
Since
is monotonically increasing, each line intersects
exactly once, so it suffices to consider the set of lines
.
Noting that
, we have
and
Thus, the image of the line
is horizontal since the
y coordinate is fixed at
. To study what happens to the
x coordinate, we look at the multiplicative factor
When
,
, so the
x coordinate
is
, i.e., it lies on
. Since
h is increasing for
and decreasing for
. Consequently, the image of the line is always to the left of
. Moreover, it traces the line from
to
exactly once for
, and the image traces the horizontal line in the reverse direction for
. □
To summarize, every point to the left of with positive x coordinate has exactly two preimages, one to the upper left of and one to its lower right. The map’s orientation is preserving in the upper left and the orientation is reversing in the lower right, best seen by the images of the lines . From the symmetry of the image of the mapping about these line, we only need to consider the mapping F restricted to . In fact, for sufficiently large n the y coordinate of is always between 0 and 1. This motivates the final theorem of this subsection.
Theorem 3. Let U denote the closed region bounded by the coordinate axes, the curve , and the line . Let , , and . Then, U is forward invariant under F, is injective, and U is globally attracting.
Proof. First, we establish that
U is invariant. Let
. From the discussion above,
always lies to the left of
, so we must prove that
implies
. Since
x lies to the left of
, we must have
. Then, we have
establishing that
U is forward invariant.
The graph of crosses the x-axis at the point . Thus, the graph of is always to the right of U. Consequently, at most one preimage of a point in U can lie in U since the other preimage must lie to the lower right of . This shows that the restriction of F to U is injective.
Finally, we show U is globally attracting. First, observe that the first quadrant is forward invariant. Note that if , then . If never maps to the first quadrant, then the set of y coordinates of is a bounded and strictly monotone increasing sequence that must converge to a value . If , this is a contradiction since there is no invariant set of points with fixed y negative coordinate by the inequality above in g. Consequently, either the sequence has its y coordinate approach 0, or it eventually maps to a positive value.
Thus, it suffices to consider
with
. We will show there exists an integer
k such that the
y coordinate of
is less than or equal to 1. Since
(since the
x coordinate lies to the left of
), if
, then we have
If the
y value of
always exceeds 1, then similar to before, the
y coordinates of
form a bounded and strictly monotone decreasing sequence that must converge to a value
. As there is no invariant set with
, the theorem is proved. □
In
Figure 3, we depict a prototype of the closed invariant region
U mentioned in Theorem 3, where the mapping
F is injective. As is clearly shown, only a portion of the critical curve
on the boundary of this region will move strictly according to the values of the parameters, while the remaining three segments of the boundary are defined by conditions that do not depend on the parameters. This region
U is crucial for the global analysis of the model, as all the dynamics of (
3) occur within it.
4. Isoclines and Fixed Points
Let us denote the restriction of
F to
U, mentioned in the previous section, by
. The isoclines of
are given by the set of points that satisfy
and the points that satisfy
. For
with
the isoclines can be broken into three distinct sets:
and
From the isoclines, we easily obtain that there are two fixed points, i.e., points with .
Proposition 1. Let , , and . Then, the mapping F from model (3) has exactly two fixed points given by and with . Proof. To be a fixed point, one either needs or . In the first case, the fixed point satisfies with . The positive solution to this quadratic yields . The second case immediately yields the origin as the other fixed point. □
Note that the isoclines
and
each break
U into two regions, which we denote as follows:
We also define
Figure 4 illustrates the relative positions between the four regions
, the three isoclines
, and the critical curve
. This figure is crucial for understanding the proof of Theorem 4 and provides a clear insight into the dynamics of the orbits in the neighborhood of the interior fixed point
, as will be shown subsequently.
Theorem 4. Let , , and . Then, the following relations hold:
- 1.
,
- 2.
,
- 3.
,
- 4.
.
Proof. We prove (1) and (2) as the other proofs are similar with reversed inequalities.
- 1.
Let
, so that
and
. Then,
which after rearrangement yields
. That is,
.
- 2.
First, note that if , then and . Since the graph of is decreasing, it suffices to prove the claim on the left boundary of . That is, it suffices to prove the result for where and either or .
Since the derivative of
is
, it is a strictly increasing function of
x for
, and we have that
In the second case,
and
, so
□
Thus, not only do points in U proceed counterclockwise around under the mapping , but in each loop around they cannot skip any of the in their itinerary. Consequently, we have the following corollary.
Corollary 1. Let , , and . Then, has no points with minimal period 2 or 3.
Now, we turn to the stability of the two fixed points and their bifurcations.
Theorem 5. The origin is always a saddle point. The interior fixed point is locally asymptotically stable if and only if Proof. We use a Jacobian matrix. We calculate
which has eigenvalues
and
, from which it follows that the origin is a saddle point. For the interior fixed point, we have
Using
, we may also express this as
From the Jury conditions [
42], a necessary and sufficient criteria for local asymptotic stability of a fixed point is given by
We claim that the first of these inequalities is always satisfied. Indeed, for the parameter values given, the trace is always positive since
, so to verify the first inequality we simply observe that the stated inequality is equivalent to
which is equivalent to the true statement
Using the second expression for
, the second inequality requires that
Rearranging, this is equivalent to
Isolating the
term yields
as claimed. □
Example 2. Let and . Then, . From Theorem 5, is stable if and only if . Solving yields , so that the fixed point loses stability at . Indeed, at this value, we have andThe characteristic polynomial is given by which yieldsThe solutions to this quadratic yields the eigenvalues . A calculation shows that for both eigenvalues so that the fixed point loses stability as β passes through 65. Corollary 2. For and , the fixed point is locally asymptotically stable for all .
Proof. The bound
is trivially satisfied since
and the left-hand side cannot be positive. □
Corollary 3. Let and . Then, the fixed point is locally asymptotically stable if and only if Proof. The fixed point is locally asymptotically stable if and only if
Since
, the given inequality is equivalent to
Squaring both sides and isolating
gives the desired inequality. □
For fixed
and
, the mapping undergoes a bifurcation at the critical value
. We recall that for smooth maps varying with a parameter, there are three types of bifurcations at a fixed point, depending on how the fixed point loses stability as the parameter changes. A fold bifurcation occurs when one of the eigenvalues of the Jacobian matrix becomes 1, while a period doubling bifurcation occurs when one of the eigenvalues of the Jacobian matrix becomes
. A third type, Neimark–Sacker bifurcation, occurs when a pair of complex eigenvalues both have magnitude 1, such as in Example 2. For a comprehensive review of bifurcations, see [
43]. We will prove momentarily that in Model (
3), the fixed point
has a Neimark–Sacker bifurcation at
, if
and
, but does not undergo a fold (saddle node) bifurcation or flip (period doubling) bifurcation.
Proposition 2. For Model (3), let and . Then, the fixed point does not undergo a flip or fold bifurcation at any . Proof. Fix and . Let . For brevity, let and .
The characteristic equation of
is given by
. Thus, the eigenvalues
of
are given by
To have a fold bifurcation, we would need
to be real and for
or equivalently for
Likewise, for a flip bifurcation we would need
to be real and for
or equivalently for
However, in the proof of Theorem 5, we proved that
, so we cannot have a fold or flip bifurcation. □
Lemma 3. For Model (3), let and . Let and let denote the interior fixed point. Then, the eigenvalues of satisfy . Proof. Fix
and
. Let
and let
. As in the proof of Proposition 2, let
and
. The eigenvalues
of
are given by
Since
, we have
, and consequently we have
. We note
, so in a neighborhood of
the eigenvalues
are complex and are given by
with
. Thus, both eigenvalues have magnitude 1 at
. □
As passes through the bifurcation value , the fixed point loses stability, and Lemma 3 shows the possibility of a Neimark–Sacker bifurcation. One of the main features of Neimark–Sacker bifurcations is the introduction of closed invariant curves. In a supercritical Neimark–Sacker bifurcation, the fixed point loses stability and stable invariant curves appear, while in a subcritical Neimark–Sacker bifurcation, there are unstable closed invariant curves that vanish at the same time the fixed point loses stability. Whether such a bifurcation is supercritical or subcritical can be determined by converting the map to the Neimark–Sacker normal form. However, this normal form changes depending on whether two nondegeneracy conditions are satisfied. The main theorem for the generic case is given below.
Theorem 6 ([
43])
. Suppose a two-dimensional discrete-time systemwith smooth F has a fixed point , where has eigenvalueswhere and . Let the following conditions be satisfied:Then, by the introduction of a complex variable z and a parameter transformation , the mapping can be transformed into the Poincare normal formwith parameter γ. If , then the bifurcation is supercritical, and if , then the bifurcation is subcritical. Lemma 4. For Model (3), conditions (C.1) and (C.2) in Theorem 6 are satisfied for and . Proof. Fix
and
, and let
denote the fixed point and
denote the bifurcation value. Note that
and that
It then follows that
In particular,
so that condition (C.1) is satisfied.
Since
we only have
for
if
, which is not the case as
. Thus, condition (C.2) is satisfied. □
Lemmas 3 and 4 are sufficient to establish the following proposition.
Proposition 3. Let and . Then, the fixed point undergoes a Neimark–Sacker bifurcation at .
A prototype of the stability region in the parameter space is shown in
Figure 5. For
, and
within the shaded region, the interior fixed point
is locally asymptotically stable, corresponding to the region referenced in the conditions of Theorem 5. The conditions in Corollary 2 define the region on the left side of the concave portion such that
, while the conditions in Corollary 3 define the region below the concave portion such that
. The Neimark–Sacker bifurcation mentioned in Proposition 3 occurs precisely on the surface of the concave part.
While we have established the existence of a Neimark–Sacker bifurcation, the computations thus far have not yet shown if it is supercritical or subcritical. In principle, one can calculate the coefficient
from Theorem 6, but the calculation for Model (
3) is lengthy and impractical, as it does not allow one to determine the sign of
for general parameter values (see
Appendix A), due to the size of the coefficients involved. Instead, we will present a geometric argument using stable and unstable manifolds in the next section.
5. Stable and Unstable Manifolds
For a fixed point
of an invertible mapping
, we denote the stable manifold of
by
and likewise its unstable manifold by
For the mapping , the stable manifold of the origin is the portion of the y-axis that lies in U, that is, .
Some numerically computed images of the unstable manifold at the origin are shown in
Figure 6. We will now proceed to show these plots are typical.
In , the eigenvector corresponds to the unstable direction in the linearized map with eigenvalue . Since lies tangent to this eigenvector, is initially in . We also recall the cycling nature of the , the direction of the mapping in each , and that the manifold cannot self-intersect. These facts ensure the argument of is strictly increasing when viewed as centered at .
If the curve remains in , then it must monotonically approach . Otherwise, it must eventually intersect so that it enters . Again, if it remains in , then it must monotonically approach , otherwise it must intersect . Eventually, if it does not stay in some , it will intersect the isoclines an infinite number of times as it spirals counterclockwise. Viewing as the center and looking at polar coordinates , we can look at the return map that gives the radius r from generated by intersections with any fixed value. Since the manifold does not self intersect, these return maps are strictly decreasing. As they are also bounded, we conclude the following theorem.
Theorem 7. Let , , and . Then, the unstable manifold is asymptotic to either the fixed point or a closed curve that encloses .
In
Figure 6, we illustrate the concept underlying Theorem 7. As is clearly seen in the presented four examples, all the orbits behave in a counterclockwise manner, either converging to the fixed point when the absolute values of both eigenvalues are less than one or orbiting around the fixed point when both eigenvalues exceed absolute value 1. This is precisely the main idea behind the proof of Theorem 7.
Let us now fix and and see how the unstable manifold evolves as increases. First, we prove the following lemma.
Lemma 5. Let . Let denote the discriminant of the characteristic polynomial of . Then .
Proof. Recall that we have
, where
and
Note also that
.
Initially, when
, it is straightforward to prove the set
is forward invariant, so that the unstable manifold
forms a heteroclinic connection between the origin and
. In this case,
, and the eigenvalues of
are
and
, i.e.,
. The heteroclinic connection for
persists under perturbations of
and is tangent to one of the eigenvectors at
. This idea is illustrated in
Figure 7. In this example, the tangency is clearly observed. Let us examine some properties of these eigenvectors in more detail.
Lemma 6. Let . Let β be such that the eigenvalues of are real. Then, any eigenvalue λ of satisfies .
Proof. Since
is an eigenvalue of
, then
has vanishing determinant. Since
and
, the product
must be negative. Since
, if
, then both terms are positive. If
, then both terms are negative or zero. This implies
. □
If is a real eigenvalue of , then the line through that is parallel to the corresponding eigenvector has a slope of . From the previous lemma, this slope is always negative. The larger of the two eigenvalues, given by , corresponds to the line with a steeper slope. It is this line that is asymptotic to.
A calculation shows
If
, then by continuity the eigenvalues of
remain real for all
. For
, eigenvalues transition at a critical value
from real to complex. At such values, the curve
transitions from directly looping into
and instead starts an infinite spiral, intersecting the isoclines an infinite number of times.
If is asymptotic to a closed curve, then that curve is necessarily stable. Consequently, the curve cannot be an unstable curve that emerged from a subcritical Neimark–Sacker bifurcation. For , this spiral approaches and the heteroclinic orbit is maintained. At , the heteroclinic orbit breaks as a supercritical Neimark–Sacker bifurcation occurs, and a stable closed curve bifurcates from the fixed point for .
Example 3. Let and , so that . We haveandThe eigenvalues of are real providedwhich yieldswhich for positive implies that the eigenvalues are real if . Solvingfor positive β yieldsThus, . Since (see Example 2), the dynamics for and may be described as follows. For , the unstable manifold is tangent to the real eigenvector and convergence to the fixed point is guaranteed. For , the eigenvalues are now complex, and the unstable manifold infinitely spirals around the fixed point, but still converges to it. For , where the Neimark–Sacker bifurcation occurs, a new forward invariant closed curve appears and the unstable manifold spirals towards this curve instead. See also Figure 6. In all cases, since forms a barrier which points cannot cross, this unstable manifold dictates global behavior on . In particular, in the absence of a closed invariant curve, the fixed point is globally asymptotically stable.
6. Conclusions
In this paper, we studied the local properties, the global properties, and the bifurcation of a nonlinear planar system of coupled difference equations consisting of the evolution of a species of Ricker type. By studying the model when the trait equation is coupled from the population dynamic equation, we were able to conclude that, under certain restrictions on the parameters, the model is forward invariant and globally attracting in a subset of the phase space, bounded by the coordinate axes, the horizontal line , and a critical curve , which is the image under the mapping that rules the model where the Jacobian vanishes. We show that in this subset of the phase space, there are four regions where the orbits behave in a cyclic nature (counterclockwise direction) and converge towards the asymmetric unstable manifold of the origin, which will attract all the orbits.
Moreover, in this set, the model possesses two symmetric fixed points: the origin, which is always a saddle fixed point, and an interior fixed point that is locally asymptotically stable for certain values of the parameters. This non-trivial fixed point undergoes a Neimark–Sacker bifurcation for a critical parameter value, a phenomenon that does not occur in the original one-dimensional Ricker model. By using a geometric argument with the unstable manifold at the origin, we show that this bifurcation is supercritical, and thus a stable asymmetric closed curve bifurcates from the fixed point once we cross the critical value. Finally, in the absence of the referred closed invariant curve, the symmetric non-trivial fixed point is globally asymptotically stable.
We have utilized tools from both analytic and geometric approaches to achieve our goals. The analytic tools are similar to those existing in the literature for these types of models, including the relative position of critical curves using singularity theory, the dynamics of sets around the isoclines, and local stability using Jury’s conditions.
Typically, the study and classification of the Neimark–Sacker bifurcation are conducted using the normal forms of the equations. Nondegeneracy conditions are usually computed by taking certain derivatives, as explained in
Section 4. However, in Model (
3), we were unable to use
to classify the bifurcation, as is mentioned in the end of
Section 4. Therefore, we employed a geometric argument, which, to the best of our knowledge, is the first time it has been used in this context. This geometric argument allowed us to demonstrate the global stability of the positive fixed point, meaning that every orbit of a point in the interior of the first quadrant will eventually converge to the fixed point. This is a central conclusion in real-world applications, particularly in ecology.
We should mention that the tools employed in this model are extendable to other discrete evolutionary models as well, such as the evolutionary Beverton–Holt model and the logistic model.
7. Future Work
In future work, we plan to study the dynamics of the model in seasonal environments, specifically focusing on how the discrete model behaves under periodic changes. This involves analyzing how the system’s behavior is influenced by parameters that fluctuate regularly, reflecting seasonal variations in the environment. In particular, we aim to understand how these periodic parameter changes affect the long-term dynamics of the model, i.e., the study of the model
when the sequence of parameters
is
-periodic, i.e.,
for all
n and some integer
.
We believe that the orbits in the periodic model (
4) will exhibit similar behavior, with counterclockwise dynamics around the periodic cycle in phase space. Furthermore, since it will be challenging to explicitly write the periodic solution, the geometric method developed for classifying the Neimark–Sacker bifurcation will be a useful tool. As observed in simulations, if the period of the cycle equals the period of the system, the same number of cycling orbits circulate counterclockwise around the periodic points. Moreover, when the eigenvalues of the Jacobian of the composition map, evaluated at the periodic orbit, have absolute values less than one, convergence can be observed. Thus, the geometric approach will be a valuable method for classifying the nature of the invariant manifolds and the dynamics around the isoclines.
Understanding these dynamics can provide insights into the long-term behavior and global stability of the population under study. In particular, it will be interesting to determine whether the periodic evolutionary model exhibits attenuation, resonance, or neither attenuation nor resonance, a phenomenon observed by R. Sacker [
44] in the periodic one-dimensional Ricker model.
We also plan to use similar techniques to study a generalization of Model (
3) for
n competing species, resulting in a model with
equations. Our real challenge will be generalizing the geometric approach to higher dimensions, as it may become impossible to explicitly compute the fixed points after a certain number of species. Hence, a geometric approach will be crucial.
An interesting study would also investigate whether our tools are applicable to a model similar to Model (
3) but with several traits in the canonical equation for evolution.