Periodic Traveling Waves in an Integro-Difference Equation With Non-Monotonic Growth and Strong Allee Effect
Abstract
We derive sufficient conditions for the existence of a periodic traveling wave solution to an integro-difference equation with a piecewise constant growth function exhibiting a stable period-2 cycle and strong Allee effect. The mean traveling wave speed is shown to be the asymptotic spreading speed of solutions with compactly supported initial data under appropriate conditions. We then conduct case studies for the Laplace kernel and uniform kernel.
Key words: Integro-difference equation, period two cycle, Allee effect, periodic traveling wave.
AMS Subject Classification: 92D40, 92D25.
1 Introduction
Integro-difference equations in the form
(1.1) |
are of great interest in the studies of invasions of populations with discrete generations and separate growth and dispersal stages. They have been used to predict changes in gene frequency [8, 9, 10, 14, 17], and applied to ecological problems [2, 3, 4, 5, 7, 11, 12, 13]. Previous rigorous studies on integro-difference equations have assumed that the growth function is nondecreasing [17, 18], or is non-monotone without Allee effect [10, 16]. The results show existence of constant spreading speeds and travelng waves with fixed shapes and speeds. Sullivan et al. [15] demonstrated numerically that an integro-difference equation with a non-monotone growth function exhibiting a strong Allee effect can generate traveling waves with fluctuating speeds, and Otto [13] showed that such an equation can have robust non-spreading solutions. In this paper we give conditions for the existence of a periodic traveling wave with two intermediate speeds for an integro-difference equation with a piecewise constant growth function that exhibits a strong Allee effect and a period-2 cycle.
Piecewise constant growth functions have been used in the studies of integro-difference equations; see for example [6, 11, 13, 15]. Such an equation is analytically tractable and it can provide specific insights into the dynamics of solutions. For the piecewise constant growth function
(1.2) |
with , which exhibits a strong Allee effect and monotonicity, Kot et al. [6] and Lutscher [11] (Section 6.1) investigated the range expansion of the population and the spreading speed, and Sullivan et al. [15] studied oscillations in spreading speeds when the dispersal kernel is density dependent. More discussions about integro-difference equations with piecewise growth functions can be found in Lutscher [11] (Chapter 15).
In this paper, we consider the integro-difference equation (1.1) with
(1.3) |
with . is a piecewise constant non-monotone growth function exhibiting a strong Allee effect [1]. Specifically, it has a stable fixed point at zero and a stable period two cycle with the Allee threshold value. This is the function considered in Otto [13] where non-spreading solitions are studied. It may be viewed as an extension of (1.2). A graphical demonstration of defined by (1.3) is given by Figure 1.

We rigorously construct periodic traveling waves with two intermediate speed for (1.1) with given by (1.3) and a proper dispersal kernel . To the best of our knowledge, this is the first time that traveling waves with oscillating speeds have been analytically established for scalar spatiotemporal equations with constant parameters. We also show that the mean traveling wave speed is the spreading speed of solutions with compactly supported initial data under appropriate conditions. We finally conduct case studies for the Laplace kernel and uniform kernel.
2 Results
In this section, we present our main results about our model. Throughout, we assume all population density functions are bounded elements of , the space of real-valued continuous functions. It is straightforward to verify this space is closed under iteration of . We denote the norm of a function by .
We make the following assumptions about the dispersal kernel.
Hypothesis 2.1.
is a piecewise differentiable function satisfying
-
(H1)
for all , and ,
-
(H2)
for all ,
-
(H3)
for all for some ,
-
(H4)
for all , for all with , the expression
has sign changes, and
-
(H5)
for all , for all with , the expression
has sign changes for sufficiently large .
2.1 Periodic traveling waves
Let be bounded, continuous functions defined by
(2.1) |
and
(2.2) |
where is a constant defined by
(2.3) |
The following theorem shows that and generate a traveling wave-like solution to (1.1). We term this solution a periodic traveling wave since it alternates between two distinct wave profiles (with different limits at ) while traveling with mean speed .
Theorem 2.2.
Let be a sequence of bounded, continuous functions defined by
(2.4) |
If , then for all .
Proof.
For , define the right inverse of by
(2.5) |
Note that is finite for , where , and by the assumption of this theorem.
By the translation invariance property of , it suffices to show
(2.6) |
and
(2.7) |
for all . The first equation follows immediately from definition (2.2). To prove the second, observe that is monotonically decreasing and satisfies and (see for example Figure 2). Thus,
Applying the integro-difference operator yields
We can then differentiate with respect to :
From (H4), has at most 1 sign change; thus, has at most 1 turning point. We also have for all , hence .
We claim
(2.8) |
where is defined by (2.3). This can be shown by arguing in cases depending on the number of turning points of .
-
Case 1.
If has no turning points, then is monotone decreasing with for all . It follows by the intermediate value theorem (IVT) that satisfies (2.8).
-
Case 2.
If has one turning point, say , then must be increasing for and decreasing for . It follows that is a global maximum, and is increasing on and decreasing on . We have for . By the IVT, is the unique solution to (2.8).
Taking the convolution of equation (2.8) with yields
This completes the proof. ∎


Theorem 2.2 describes a rightward periodic traveling wave, but we can define the corresponding leftward periodic traveling wave via . This is due to the fact that is symmetric and translation invariant.
We investigated the intermediate wave speed generated by the periodic traveling wave solution (2.4). To do this, we define a sequence of continuous functions by
(2.9) |
This expression yields the rightmost intersection of the curve with the horizontal line . By computing the first difference, we obtain a new function sequence (with the same domain and range) defined by
(2.10) |
From Theorem 2.2, we have that converges to a period- cycle for every .
Our results about the intermediate wave speed sequence are summarized in the following corollary.
Corollary 2.3.
Assume the hypothesis of Theorem 2.2 holds. Then for all , the sequence satisfies
(2.11) |
Proof.
Thus, for each fixed , we can define intermediate traveling wave speeds and by
(2.12) |
and
(2.13) |
It follows that has a periodic traveling wave with wave profiles and , intermediate wave speeds and , and mean wave speed . Note that the intermediate wave speeds depend on , but their sum is identically equal to .
2.2 Periodic spreading solutions
By applying the results of Theorem 2.2, we were able to prove the asymptotic spreading speed of solutions with compact initial data. In particular, we showed that is the asymptotic spreading speed of compactly supported initial data with sufficient weight above the Allee threshold but below the overcompensation threshold.
Theorem 2.4.
Assume the hypothesis of Theorem 2.2 holds. Suppose is bounded, non-negative, and has compact support. If
-
i.
,
-
ii.
the set is connected and has sufficiently large length, and
-
iii.
,
then the sequence defined by and for spreads asymptotically with mean speed .
Proof.
We define two continuous functions by
(2.14) |
and
(2.15) |
Observe that is symmetric with respect to , and is symmetric with respect to . This notation can be justified by observing that and as for all . We also set
(2.16) |
and are the right-inverse of and , respectively. For a fixed , the value of is finite if , with the latter bound converging to as . They satisfy
We now apply the growth function to , assuming is sufficiently large enough so that ; this is guaranteed by taking . We have
Note that the above expression converges to as for all . We then apply the convolution operator:
Computing the derivative of this expression yields
By (H5), the number of turning points of is at most 3 (assuming is sufficiently large). Furthermore, is symmetric around ; thus the sign changes must come in pairs or occur exactly at . We can immediately exclude the cases of 0 and 2 turning points because is non-negative, vanishes at , and is not identically zero. As in the proof of Theorem 2.2, this leaves two cases:
-
Case 1.
If has 1 turning point, then it must occur at . We have as ; hence, it is sufficient to take large enough so that . It follows by the triangle inequality that for all .
-
Case 2.
If has 3 turning points, then again by symmetry, they must be given by , for some . We can label them in increasing order by . It follows that and are maxima with as , and is a local minima with as . Thus is monotone increasing on and monotone decreasing on . By taking sufficiently large, we can make and .
In both cases, we have on a closed interval of radius and elsewhere. Taking composition with yields
Hence,
Shifting the expression left by units yields
(2.17) |
We are now prepared to prove the theorem using an inductive argument. Let be the level set defined in the statement of the theorem. Without loss of generality, we write , for some . We have
for all . Applying the convolution operator yields
Let be a sequence of real numbers satisfying
(2.18) |
If is sufficiently large, we can show that is strictly increasing and unbounded. This is due to our assumption that . Since , it suffices to assume , i.e. is sufficiently large so that for all . Plugging this into the recurrence (2.18) yields .
This argument shows that, assuming is sufficiently large, equation (2.17) may be applied recursively. From the principle of induction, it follows that for all . Knowing that is unbounded allows us to use equation (2.18) to compute the asymptotic spreading speed as the following limit.
This completes the proof. ∎
3 Examples
In this section, we construct the periodic traveling wave solution for uniform and Laplace kernels and derive formulas for the mean spreading speed in terms of the model parameters.
3.1 Laplace kernel
We applied our main results to the Laplace dispersal kernel, defined
(3.1) |
is symmetric and has connected support, satisfying (H1), (H2), and (H3). We will now prove that satisfies (H5). Let with , let , and define . For sufficiently large, we have
is monotone on with ; thus it cannot have any sign changes there. On each of the intervals , , and , has the form . Depending on the sign of and , can have either one or zero sign changes on each interval. Summing over each interval, it follows that has at most sign changes. This proves (H5); the proof of (H4) follows a similar argument.
Assuming , the periodic traveling wave profiles are given by
(3.2) |
and
(3.3) |
with
(3.4) |
and are constants given by
(3.5) |
The formulas for and are given by
(3.6) |
and
(3.7) |
The critical spreading speed for the Laplace kernel thus falls into one of three cases:
-
Case 1.
, :
-
Case 2.
, :
-
Case 3.
, :
Figure 3 shows the results of a numerical simulation with growth parameters were , , and with the Laplace dispersal kernel, for an average spreading speed of . Figures 3(a) and 3(b) show heatmaps of the periodic traveling wave solution and a spreading solution with compactly supported initial condition. Figure 3(c) shows a numerical simulation of the periodic traveling wave solution, while Figure 3(d) shows the asymptotic convergence of the intermediate wave-speed sequence , evaluated at the Allee threshold .




3.2 Uniform kernel
The uniform dispersal kernel is given by
(3.8) |
Like the previous example, is symmetric around zero. We can verify (H4) by letting and and defining . There are two cases:
-
Case 1.
If , then
-
Case 2.
If , then
In both cases, the number of sign changes is exactly 1. To extend this argument to (H5), observe that since has compact support, we have
for sufficiently large .
The alternating wave profiles are given by
(3.9) |
and
(3.10) |
with
(3.11) |
has a global maximum at so that . Thus, a sufficient condition for the existence of the periodic traveling wave is .
The intermediate wave speeds are given by
(3.12) |
and
(3.13) |
for . By evaluating the intermediate wave speeds at , we obtain if , and if . So for , the traveling wave is periodic with two different intermediate wave speeds. Furthermore, the difference between these two intermediate speeds is increasing with respect to , the Allee threshold, and decreasing with respect to , the overcompensation threshold.
Figure 4(a) depicts the periodic traveling wave solution (2.4) with the uniform kernel with growth parameters are , , and . The mean spreading speed can be analytically calculated as . The intermediate speeds are and . Figure 4(b) demonstrates the spreading phenomena with initial domain size . Since our uniform kernel is compact with support , a sufficient lower bound on the initial domain size is .


References
- [1] W. C. Allee. 1931. Animal Aggregations. A Study on General Sociology. University of Chicago Press, Chicago, IL.
- [2] A. Hastings and K. Higgins. 1994. Persistence of transients in spatially structured ecological models. Science 263: 1133-1136.
- [3] M. Kot and W. M. Schaffer. 1986. Discrete-time growth-dispersal models. Math. Biosci. 80: 109-136.
- [4] M. Kot. 1989. Diffusion-driven period doubling bifurcations. Biosystems 22: 279-287.
- [5] M. Kot. 1992. Discrete-time traveling waves: Ecological examples. J. Math. Biol. 30: 413-436.
- [6] M. Kot, M. A. Lewis, and P. van den Driessche. 1996. Dispersal data and the spread of invading organisms. Ecology 77: 2027-2042.
- [7] M. Kot. 2001. Elements of Mathematical Ecology. Cambridge University Press. Cambridge, United Kingdom.
- [8] R. Lui. 1982. A nonlinear integral operator arising from a model in population genetics. I. Monotone initial data. SIAM. J. Math. Anal. 13: 913-937.
- [9] R. Lui. 1982. A nonlinear integral operator arising from a model in population genetics. II. Initial data with compact support. SIAM. J. Math. Anal. 13: 938-953.
- [10] R. Lui. 1983. Existence and stability of traveling wave solutions of a nonlinear integral operator. J. Math. Biol. 16:199-220.
- [11] F. Lutscher. 2019. Integrodifference Equations in Spatial Ecology. Springer.
- [12] M. Neubert, M. Kot, and M. A. Lewis. 1995. Dispersal and pattern formation in a discrete-time predator-prey model. Theor. Pop. Biol. 48: 7-43.
- [13] G. Otto. 2017. Non-spreading Solutiona in a Integro-Difference Model Incorporating Allee and Overcompensation Effects. Ph. D thesis, University of Louisville.
- [14] M. Slatkin. 1973. Gene flow and selection in a cline. Genetice 75: 733-756.
- [15] L. L. Sullivan, B. Li, T. E. X. Miller, M. G. Neubert, and A. K. Shaw. 2017. Density dependence in demography and dispersal generates fluctuating invasion speeds. Proc. Natl. Acad. Sci. USA 114: 5053-5058.
- [16] M. H. Wang, M. Kot, and M. G. Neubert. 2002. Integrodifference equations, Allee effects, and invasions. J. Math. Biol. 44: 150-168.
- [17] H. F. Weinberger. 1978. Asymptotic behavior of a model in population genetics. Nonlinear Partial Differential Equations and Applications, ed. J. M. Chadam. Lecture Notes in Mathematics 648: 47-96. Springer-Verlag, Berlin.
- [18] H. F. Weinberger. 1982. Long-time beahvior of a class of biological models. SIAM. J. Math. Anal. 13: 353-396.