title: Three stage semelparous Leslie models.
abstract: 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.
authors.affiliation: 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
authors.forename: J M
authors.initials: JM
authors.lastname: Cushing
mesh_terms: D000367:Age Factors; D000465:Algorithms; D000818:Animals; D008018:Life Cycle Stages; D008954:Models, Biological; D017711:Nonlinear Dynamics; D010507:Periodicity; D011157:Population Dynamics; D012098:Reproduction
references.citation: Math Biosci. 1997 Nov;146(1):37-62
references.pmid: 9357293
references.citation: Biometrika. 1945 Nov;33:183-212
references.pmid: 21006835
references.citation: Trends Ecol Evol. 1994 Dec;9(12):475-80
references.pmid: 21236926
references.citation: J Math Biol. 2005 Jan;50(1):1-21
references.pmid: 15692839
references.citation: J Math Biol. 2000 May;40(5):413-31
references.pmid: 10885591
references.citation: J Math Biol. 2007 Nov;55(5-6):781-802
references.pmid: 17639397
references.citation: J Math Biol. 2004 May;48(5):515-28
references.pmid: 15133621
references.citation: J Math Biol. 2002 May;44(5):450-62
references.pmid: 12021984
references.citation: J Biol Dyn. 2007 Oct;1(4):347-62
references.pmid: 22876821
references.citation: J Math Biol. 2003 Feb;46(2):95-131
references.pmid: 12567230
references.citation: Math Biosci Eng. 2006 Jan;3(1):17-36
references.pmid: 20361805