Three stage semelparous Leslie models.
Nonlinear Leslie matrix models have a long history of use for modeling the dynamics of semelparous species. Semelparous models, as do nonlinear matrix models in general, undergo a transcritical equilibrium bifurcation at inherent net reproductive number R(0) = 1 where the extinction equilibrium loses stability. Semelparous models however do not fall under the purview of the general theory because this bifurcation is of higher co-dimension. This mathematical fact has biological implications that relate to a dichotomy of dynamic possibilities, namely, an equilibration with over lapping age classes as opposed to an oscillation in which age classes are periodically missing. The latter possibility makes these models of particular interest, for example, in application to the well known outbreaks of periodical insects. While the nature of the bifurcation at R(0) = 1 is known for two-dimensional semelparous Leslie models, only limited results are available for higher dimensional models. In this paper I give a thorough accounting of the bifurcation at R(0) = 1 in the three-dimensional case, under some monotonicity assumptions on the nonlinearities. In addition to the bifurcation of positive equilibria, there occurs a bifurcation of invariant loops that lie on the boundary of the positive cone. I describe the geometry of these loops, classify them into three distinct types, and show that they consist of either one or two three-cycles and heteroclinic orbits connecting (the phases of) these cycles. Furthermore, I determine stability and instability properties of these loops, in terms of model parameters, as well as those of the positive equilibria. The analysis also provides the global dynamics on the boundary of the cone. The stability and instability conditions are expressed in terms of certain measures of the strength and the symmetry/asymmetry of the inter-age class competitive interactions. Roughly speaking, strong inter-age class competitive interactions promote oscillations (not necessarily periodic) with separated life-cycle stages, while weak interactions promote stable equilibration with overlapping life-cycle stages. Methods used include the theory of planar monotone maps, average Lyapunov functions, and bifurcation theory techniques.
Cushing
J M
JM
Department of Mathematics and the Interdisciplinary Program in Applied Mathematics, University of Arizona, 617 N Santa Rita, Tucson, AZ 85721, USA. cushing@math.arizona.edu
D000367:Age Factors; D000465:Algorithms; D000818:Animals; D008018:Life Cycle Stages; D008954:Models, Biological; D017711:Nonlinear Dynamics; D010507:Periodicity; D011157:Population Dynamics; D012098:Reproduction
Math Biosci. 1997 Nov;146(1):37-62
9357293
Biometrika. 1945 Nov;33:183-212
21006835
Trends Ecol Evol. 1994 Dec;9(12):475-80
21236926
J Math Biol. 2005 Jan;50(1):1-21
15692839
J Math Biol. 2000 May;40(5):413-31
10885591
J Math Biol. 2007 Nov;55(5-6):781-802
17639397
J Math Biol. 2004 May;48(5):515-28
15133621
J Math Biol. 2002 May;44(5):450-62
12021984
J Biol Dyn. 2007 Oct;1(4):347-62
22876821
J Math Biol. 2003 Feb;46(2):95-131
12567230
Math Biosci Eng. 2006 Jan;3(1):17-36
20361805