Modeling Population Dynamics: Applications and Recent Developments
Org:
Kunquan Lan (Ryerson),
Gunog Seo (Colgate University) and
Gail S.K. Wolkowicz (McMaster)
[
PDF]
 JACQUES BELAIR, Université de Montréal
Dynamics of a quail population with seasondependent reproduction rates [PDF]

We present a model, inspired by field data, for the semiannual evolution of a population of quails with weatherdependent reproduction rates for each of the intermediate seasons, the density of cover making each of the rates monotone decreasing and bounded away from zero. The ensuing model takes the form of a system of two nonlinear finitedifference equations which was derived previously (Milton and Bélair, Theoretical Population Biology [37]) but considered then under considerably simplifying assumptions. A geometric method (generalized cobwebbing) is presented to explicitly determine the asymptotic behaviour of the solutions, and numerical as well as analytical methods are used to establish stability, and bifurcations, of equilibria.
 YUMING CHEN, Wilfrid Laurier University
Dirichlet problem for a diffusive logistic population model with two delays [PDF]

We investigate a diffusive logistic equation with nonzero boundary Dirichlet condition and two delays. We first exclude the existence of positive heterogeneous steady states, which implies the uniqueness of constant positive steady states. Then, we analyze the local stability and local Hopf bifurcation at the unique constant positive steady state.
We show that multiple delays can induce multiple stability switches. Furthermore, we prove global stability of the positive steady state under certain conditions and obtain global Hopf bifurcation results. Numerical simulations have been carried out to illustrated the obtained theoretical results. This is a joint work with Xuejun Pan and Hongying Shu.
 HARRY GAEBLER, University of Guelph
Modeling Uranium Bioreduction by Multispecies Biofilms [PDF]

In oxidized groundwater, soluble uranium (U[VI]) has the potential to rapidly disperse and lead to widespread contamination. In this work, we develop and discuss two models to investigate the use of insitu microorganisms, in the form of bacterial biofilms, as an effective treatment technology for uranium contaminated soil and groundwater.
The first model is a reactive transport model that describes uranium bioreduction in porous media. This model is a multiscale model developed first by describing microbial processes at the mesoscale, then upscaling these processes to the macroscale (reactor scale). Processes included in the mesoscale model are hydrodynamics and transport of substrates in the reactor, two biofilm and suspended bacteria species and their respective growth kinetics in the pore space through consumption of appropriate substrates coupled with the bioreduction of uranium, attachment of suspended cells to each of the biofilms, detachment of biofilm cells, and cell lysis.
The second model is a traditional chemostat model of the system described above, which is used to investigate the importance of specific reaction kinetics on the overall effectiveness of uranium bioreduction.
 TINGHAO HSU, University of Miami
SlowFast Systems and Regime Shifts [PDF]

In ecology, regime shifts are continual rapid change between different longlasting dynamics. An example is disease outbreak, where a system exhibits qualitative changes after long periods of apparent quiescence. Another example is rapid ecoevolutionary dynamics, which have been observed widely both in predators and in prey. For systems modeling those examples, we derive conditions under which regime shifts occur. Our approach was based on extending the socalled entryexit function to multidimensional slowfast systems using geometric singular perturbation theory. This is joint work with Shigui Ruan and Gail Wolkowicz.
 HAO KANG, University of Miami
Nonlinear agestructured population models with nonlocal diffusion and nonlocal boundary conditions [PDF]

In this paper, we develop the basic theory for agestructured population models with nonlocal diffusion and nonlocal boundary conditions. We first apply the theory of integrated semigroups and nondensely defined operators to the linear equations, study the spectrum, and analyze the asymptotic behavior via asynchronous exponential growth. Then we consider a semilinear equation and a nonlinear equation and establish the existence and stability of the nontrivial equilibria.
 KUNQUAN LAN, Ryerson University
Parabolic boundary value problems and applications to population models with spatially varying growth rate and external forcing [PDF]

This presentation is based on recent work on parabolic boundary value problems with nonlinearity changing signs and applications to population models with spatially varying growth rate and external forcing. The fixed point index theory for $r$nowhere normaloutward compact maps obtained by Yang and Lan in 2016 will be presented.This index theory is a generalization of the classical fixed point index theory for compact selfmaps defined in cones in Banach spaces, and is a good tool to deal with the above boundary value problems and population models.
This is joint work with Professor Wei Lin at Fudan University
 ZHENYANG LI, Honghe University
DYNAMICS OF MAPS WITH MEMORY [PDF]

We consider a map of the unit square which is not 11, defined as follows: $x_{n+1}=M_{\alpha }(x_{n1},x_{n})=\tau (\alpha \cdot x_{n}+(1\alpha )\cdot x_{n1}),$ where $\tau$ is a onedimensional map on $I=[0,1]$ and $0<\alpha <1$ determines how much memory is being used. We let $\tau $ to be the symmetric tent map. To study the dynamics of $M_\alpha$, we consider the twodimensional map $$ G_{\alpha }:[x_{n1},x_{n}]\mapsto [x_{n},\tau (\alpha \cdot x_{n}+(1\alpha )\cdot x_{n1})]\, .$$ The map $G_\alpha$ for $\alpha\in(0,3/4]$ studied in our paper is shown to have various properties depending on the $\alpha$ applying Tsujii's result. For $\alpha\in(3/4,1)$ the map $G_\alpha$ admits a singular SinaiRuelleBowen measure in another paper, we do this by applying Rychlik's results for the Lozi map. Unlike the Lozi map, the maps $G_\alpha$ are not invertible.
 CHIUJU LIN
A Novel Delayed Logistic Equations with DecayConsistent Delay in Growth [PDF]

We derive an alternative expression for a delayed logistic equation, assuming that the rate of change of the population depends on three components: growth, death, and intraspecific competition, with the delay in the growth component. Following Arino et al. (J Theor Biol 241(1):109–119, 2006), in our formulation we incorporate the delay in the growth term in a manner consistent with the rate of instantaneous decline in the population. For our new model, we provide a complete global analysis showing that no sustained oscillations are possible. We further extend our formulation to a two species competition and use adaptive dynamics to conclude that, assuming certain tradeoffs between the growth rate and the delay, the evolutionary trend is to make the delay approaching an optimal value.
 WEI LIN, Fudan, China
 XINZHI LIU, University of Waterloo
BioInspired Coordination and Consensus of MultiAgent Systems [PDF]

Many biological systems, such as swarms of birds, flocks of beasts, schools of fish, armies of ants, and colonies of bees, exhibit fascinating collective behavior. It is observed that in such systems, each individual acts as an autonomous agent and interacts only with its nearby neighbors, while the entire group displays coordinated behavior and can accomplish very complex tasks. Inspired by such collective intelligence of animal groups in nature, there has been an increased research interests in multiagent systems around the world in recent years. This talk discusses the distributed consensus problems of multiagent systems with both fixed and switching topologies. Hybrid consensus protocols are proposed to take into consideration of continuous communications among agents and intermittent information exchanges on a sequence of discrete times. Based on the proposed algorithms, the multiagent systems with the hybrid consensus protocols achieve consensus by employing results from matrix theory and algebraic graph theory. Our results show that the hybrid consensus protocols can solve the consensus problem if the union of continuoustime and discretetime interaction digraphs contains a spanning tree frequently enough. Simulations are provided to demonstrate the effectiveness of the proposed consensus protocols.
 ROBERT MCCANN, University of Toronto
When do interacting organisms gravitate to the vertices of a regular simplex? [PDF]

Flocking and swarming models often assume that organisms interact through a force which is attractive over large distances yet
repulsive at short distances. Assume this force is given as a difference of power laws
and normalized so that its unique minimum occurs at unit separation. For a range of exponents corresponding to mild repulsion and strong attraction, we show
that the minimum energy configuration is uniquely attained  apart from translations and rotations  by equidistributing
the organisms over the vertices of a regular topdimensional simplex (i.e. an equilateral triangle in two dimensions and regular tetrahedron in three).
If the attraction is not assumed to be strong, we show these configurations are at least local energy minimizers in the relevant
$d_\infty$ metric from optimal transportation, as are all of the other
uncountably many unbalanced configurations with the same support.
We infer the existence of phase transitions.
One ingredient from the proof will be described in another session, namely the establishment of
a simple isodiametric variance bound which characterizes regular simplices:
it shows that among probability measures on ${\bf R}^n$ whose supports
have at most unit diameter, the variance around the mean is maximized precisely by
those measures which assign mass $1/(n+1)$ to each vertex of a (unitdiameter) regular simplex.
Based on preprint with Tongseok Lim at {\tt https://arxiv.org/abs/1907.13593}
 SHIGUI RUAN, University of Miami
On a Network Model of Two Competitors with Applications to the Invasion and Competition of Aedes Albopictus and Aedes Aegypti Mo [PDF]

The two prominent mosquito species, Aedes aegypti and Ae. albopictus, are the primary vectors that transmit several arboviral diseases, including chikungunya, dengue fever, yellow fever, and Zika. The world is presently experiencing a series of major outbreaks of these vectorborne diseases, so it is very important and necessary to understand the current distributions and movements of these mosquito vectors for successful surveillance and control programs. Based on the invasion of the Ae. albopictus mosquitoes and the competition between Ae. Albopictus and Ae. aegypti mosquitoes in the US, we study a twospecies competition model in a network, that is with discrete Laplacian diffusion. In the case of strongweak competition where the invasive competitor is stronger than the local one, it is shown that the invasive species (Ae. albopictus) wins over the local species (Ae. aegypti) and the solutions converge uniformly to the semipositive equilibrium such that the invasive species survives while the local species becomes extinct, and vice versa. In the case of weakweak competition, the solutions converge uniformly to the positive equilibrium so that both invasive and local species coexist. By using numerical simulations, we apply the twospecies competition model in a network to explain the invasion and competition of Ae. Albopictus and Ae. Aegypti in the US. We also show that discrete Laplacian diffusion induces different spreading speeds in different invasive directions.
 GUNOG SEO, Colgate University
A PredatorPreyModel with a Holling Type I Functional Response Including a Predator Mutual Interference [PDF]

The most widely used functional response in describing predatorprey relationships is the Holling type II functional response, where per capita predation is a smooth, increasing, and saturating function of prey density. Beddington and DeAngelis modified the Holling type II response to include interference of predators that increases with predator density. In this talk, I introduce a predatorinterference term into a Holling type I functional response. I explain the ecological rationale for the response and note that the phase plane configuration of the predator and prey isoclines differs greatly from that of the BeddingtonDeAngelis response; for example, in having three possible interior equilibria rather than one. In fact, this new functional response seems to be quite unique. I use analytical and numerical methods to show that the resulting system shows a much richer dynamical behavior than the BeddingtonDeAngelis response, or other typically used functional responses. For example, cyclicfold, saddlefold, homoclinic saddle connection, and multiple crossing bifurcations can all occur. I then use a smooth approximation to the Holling type I functional response with predator mutual interference to show that these dynamical properties do not result from the lack of smoothness, but rather from subtle differences in the functional responses. This is joint work with Donald L. DeAngelis at the University of Miami.
 ZHONGWEI SHEN, University of Alberta
Transient dynamics in stochastic population models [PDF]

Transient dynamics, often observed in multiscale systems, are roughly defined to be the interesting dynamical behaviours that display over finite time periods. For a class of randomly perturbed dynamical systems that arise in chemical reactions and population dynamics, and that exhibit persistence dynamics over finite time periods and extinction dynamics in the long run, we use quasistationary distributions (QSDs) to rigorously capture the transient states governing the transient dynamics. We study the noisevanishing concentration of the QSDs to gain information about the transient states.
 JUNPING SHI, College of William & Mary
Spectral Monotonicity of Perturbed Quasipositive Matrices with Applications in Population Dynamics [PDF]

The reduction principle for structured population model states that fast dispersal will reduce the growth rate for a spatially distributed population. This reduction principle is proved in two different ways: one using graph theoretical method, and the one using analytic method. Applications of the reduction principle are shown for scalar model, competition model and predatorprey model. This is a joint work with Shanshan Chen, Zhisheng Shuai and Yixiang Wu.
 ROBERT SMITH?, The University of Ottawa
Using nonsmooth models to determine thresholds for microbial pest management [PDF]

Releasing infectious pests could successfully control and eventually maintain the number of pests below a threshold level. To address this from a mathematical point of view, two nonsmooth microbial pestmanagement models with threshold policy are pro posed and investigated in the present paper. First, we establish an impulsive model with statedependent control to describe the cultural control strategies, including releasing infectious pests and spraying chemical pesticide. We examine the existence and stability of an order1 periodic solution, the existence of orderk periodic solutions and chaotic phenomena of this model by analyzing the properties of the Poincaré map. Secondly, we establish and analyze a Filippov model. By examining the sliding dynamics, we investigate the global stability of both the pseudoequilibria and regular equilibria. The findings suggest that we can choose appropriate threshold levels and control intensity to maintain the number of pests at or below the economic threshold. The modelling and control outcomes presented here extend the results for the system with impulsive interventions at fixed moments.
 LINDA WAHL, University of Western Ontario
The lysislysogeny decision in variable environments [PDF]

Bacteriophage, viruses that infect bacteria, are the most abundant life form on the planet, and are critical to both ecosystem function and human health. After infecting a host cell, temperate bacteriophage have two lifehistory strategies: they can immediately kill the cell and release progeny viral particles (lysis); alternatively, they can integrate their genetic material in the host genome (lysogeny), persisting as an integral part of the host for many generations. Thus, in lysogeny, the predator can choose to become genetically part of the prey. Similar to a bethedging strategy, recent experimental evidence suggests that lysogeny may be favoured in variable environments. We examine the evolutionary stability of the lysislysogeny decision by analysing invasion reproductive numbers in this intriguing predatorprey system, evolving in the context of a periodic environment.
 HAO WANG, University of Alberta
Multiscale and qualitative analysis of a stoichiometric algae model in the absence or presence of bacteria [PDF]

Algal blooms are becoming a global concern due to the increasing prevalence of eutrophication. Often, algae and bacteria interact, within the wellmixed epilimnion, in a loose commensalistic way. Here we analyze stoichiometric models for algal dynamics and for bacteriaalgae dynamics. The algaeonly dynamics exhibit rich transient behavior, and the driving biological mechanisms are studied and understood via a multiple timescale analysis. We further perform global qualitative analysis of both models. There are three dynamical scenarios determined by the basic reproductive numbers of algae and bacteria. We use these models to make specific predictions about how the relative balance of algae and bacteria should change in response to varied nutrient and light availability. The bacteriaalgae model successfully reproduces empirical respiration data.
 LIN WANG, University of New Brunswick
Multitype bistability and long transients in a delayed budworm model [PDF]

For a delayed budworm model, we show that the delay induce four types of bistability: (i) a nodenode bistability (one stable positive equilibrium coexists with another stable positive equilibrium); (ii) a nodecycle bistability (one stable equilibrium coexists with a stable periodic solution); (iii) a cyclecycle bistability (one stable periodic solution coexists with another stable periodic solution); and (iv) a nodechaotic attractor bistability (one stable positive equilibrium coexists with a chaotic attractor). Moreover, we also show that delay can induce long transients.
 SHIROU WANG, University of Alberta
(De)synchronization for Markov random networks arising from Markov perturbations [PDF]

A physical network is naturally subject to noise influences from both external (extrinsic) and internal (intrinsic) sources. The extrinsic noises are usually environmentally related, while the intrinsic ones are typically due to internal uncertainties. A discretetime, discretestate (dtds) network with only extrinsic noises is commonly modeled by a discrete random dynamical system (RDS), but the one with only intrinsic noises is often modeled by a Markov chain. In this talk, we will consider a dtds network with both extrinsic and intrinsic noises under the framework of the socalled Markov random network (MRN). In particular, we will discuss the phenomenon and mechanism of (de)synchronizations for MRNs which arise as Markovperturbations of a synchronized discrete RDS. Characterization of (de)synchronizations will be given from the view points of both probability distributions and dynamical systems.
 JAMES WATMOUGH, University of New Brunswick
A general settlement distribution for dispersal across patchy landscapes [PDF]

Many species have a life stage that disperses larger distances than its other stages. For many marine organisms this is a pelagic larval stage. For many familiar insect species, it is an adult stage. Dispersal is across landscapes that are often patchy, with distinct interfaces between suitable and unsuitable habitat for the organism. Some of the central questions in mathematical ecology concern on the effect of patch spacing on persistence and spread of the species. In many analyses, dispersal is modelled by an integral kernel describing the distribution of individuals that survive the dispersal stage. These kernels are complicated by the fact that while dispersing across patchy landscapes, individuals are subject to differential mortality and settlement probabilities in different habitats. These differential rates introduce asymmetries into the kernel. In addition to differential rates across the landscape, individuals may change their behaviour at the patch interfaces. This interface bias may increase or decrease retention within a patch. We extend previous formulations of dispersal kernels with bias at patch boundaries and derive a general multiexponential dispersal kernel for a patchy landscape with biased movement across patch interfaces. The kernel is piecewise exponential, with coefficients that can be determined by a straightforward matrix product of the patch parameters. We speculate this matrix formulation of the kernel will simplify numerical computations of integrodifference equation models for population dynamics on patchy landscapes.
 GAIL S. K. WOLKOWICZ, McMaster University
Analysis of a Symplified Model of Anaerobic Digestion [PDF]

Anaerobic digestion is a complex naturally occurring process used for waste and wastewater treatment to produce biogas as a renewable source of energy. The detailed Anaerobic Digestion Model No. 1 (ADM1) includes 32 state variables and is not mathematically tractable.
Bornhoft, HankeRauschenbach, and Sundmacher [Nonlinear Dyn., 73 (2013)] introduced a qualitative simplification of the this model. The global dynamics of this simplified model will be obtained by first analyzing the limiting system, a model of single species growth in the chemostat in which the response function is nonmonotone and the species decay rate is included to complete the theory for the chemostat in this case. Using a Lyapunov function argument and the theory of asymptotically autonomous systems, it will be shown that even in the parameter regime where there is bistability, no periodic orbits exist and every solution converges to one of the equilibrium points. Then two algorithms for stochastically perturbing the parameters of the model will be described. Simulations done with these two algorithms will be compared with simulations done using the Gillespie and tauleaping algorithms. Finally the dynamics will be compared with the dynamics predicted by the ADM1 model and other simplifications of the ADM1 model.
This is joint work with Tyler Meadows and Marion Weedermann.
 YINGFEI YI, University of Alberta
A Mesoscopic Ergodic Theorem [PDF]

The talk will discuss ergodic properties of mesoscopic systems described by stochastic ordinary differential equations. An ergodic theorem will be presented for systems with rough coefficients which do not necessarily generate flows or semiflows. Applications to mechanical and biological
systems will be discussed.
 PEI YU, Western University
Study of SlowFast Motions Using Dynamical Systems Approach [PDF]

In this talk, we present a method to analyze certain slowfast motions in dynamical systems. For singular perturbed dynamical systems,
the wellknown Geometric Singular Perturbation Method (GSPM) is usually applied to find the special limit cycles  slowfast periodic
solutions. However, many practical problems might be not able or very difficult to be put in the form of singular perturbed equations,
but they still exhibit slowfast motions. For such cases, based on dynamical system theory, we developed a method to identify and analyze slowfast motions. We will use several biological examples to illustrate our method, and give a comparison between the GSPM and our method.
 ZHITAO ZHANG, Academy of Mathematics and Systems Science, the Chinese Academy of Sciences
HénonLaneEmden conjecture and related Schrödinger systems [PDF]

We have proved HenonLaneEmden conjecture is true for space dimension $N=3$ by scaling invariant of the solutions and Sobolev embedding on $S^{N1}$. Then we obtained new Liouvilletype theorems and showed HenonLaneEmden conjecture for polyharmoic system holds in a new region, and also proved the generalized H\'{e}nonLaneEmden conjecture in $\mathbb{R}^2$ and $\mathbb{R}^3$. Moreover, we prove some new results for existence of solutions to related Schr\"{o}dinger equations and systems.
 HUAIPING ZHU, York University
Local and global dynamics of models for mosquitoborne diseases [PDF]

Mosquitoes and mosquitoborne diseases (MBDs) have become a severe
burden of public health. Among many concerns, public health concerns
most about the triggers of an outbreak and conditions of repeated
infestations. In this talk, I will present some of our dynamical studies
to answer the two questions mathematically. I will introduce the local
dynamics including local stability and lower codimension bifurcations
and their association with the triggering conditions of an outbreak.
Besides, I will mainly focus on explaining the global stability,
existence and nonexistence of periodic solutions, and one related open
problem for MBDs.
 XINGFU ZOU, University of Western Ontario
On some RD models for the sterile insect release method [PDF]

Sterile insect release method (SIRM) is a biological control method for troublesome insects. As the fast advancement in genetic techniques and other biological techniques, its range of applications has been greatly enlarged. In this talk, I will begin with a basic ordinary differential equation model for SIRM, and then, move forward to a some reaction diffusion models on a bounded domain. Two releasing strategies for the sterile strain will be be discussed, one is domainwise release and the other is release only on the boundary. By analyzing the models, we obtain some conditions that ensure the success of the sterile release method in eradicating the insects.
© Canadian Mathematical Society