Spreading speeds and traveling waves in competitive recursion systems

This paper is concerned with the spreading speeds and traveling wave solutions of discrete time recursion systems, which describe the spatial propagation mode of two competitive invaders. We first establish the existence of traveling wave solutions when the wave speed is larger than a given threshold. Furthermore, we prove that the threshold is the spreading speed of one species while the spreading speed of the other species is distinctly slower compared to the case when the interspecific competition disappears. Our results also show that the interspecific competition does affect the spread of both species so that the eventual population densities at the coexistence domain are lower than the case when the competition vanishes.     

Global dynamics of a reaction and diffusion model for Lyme disease

This paper is devoted to the mathematical analysis of a reaction and diffusion model for Lyme disease. In the case of a bounded spatial habitat, we obtain the global stability of either disease-free or endemic steady state in terms of the basic reproduction number R?. In the case of an unbounded spatial habitat, we establish the existence of the spreading speed of the disease and its coincidence with the minimal wave speed for traveling fronts. Our analytic results show that R? is a threshold value for the global dynamics and that the spreading speed is linearly determinate.     

Spreading speed, traveling waves, and minimal domain size in impulsive reaction-diffusion models

How growth, mortality, and dispersal in a species affect the species' spread and persistence constitutes a central problem in spatial ecology. We propose impulsive reaction-diffusion equation models for species with distinct reproductive and dispersal stages. These models can describe a seasonal birth pulse plus nonlinear mortality and dispersal throughout the year. Alternatively, they can describe seasonal harvesting, plus nonlinear birth and mortality as well as dispersal throughout the year. The population dynamics in the seasonal pulse is described by a discrete map that gives the density of the population at the end of a pulse as a possibly nonmonotone function of the density of the population at the beginning of the pulse. The dynamics in the dispersal stage is governed by a nonlinear reaction-diffusion equation in a bounded or unbounded domain. We develop a spatially explicit theoretical framework that links species vital rates (mortality or fecundity) and dispersal characteristics with species' spreading speeds, traveling wave speeds, as well as minimal domain size for species persistence. We provide an explicit formula for the spreading speed in terms of model parameters, and show that the spreading speed can be characterized as the slowest speed of a class of traveling wave solutions. We also give an explicit formula for the minimal domain size using model parameters. Our results show how the diffusion coefficient, and the combination of discrete- and continuous-time growth and mortality determine the spread and persistence dynamics of the population in a wide variety of ecological scenarios. Numerical simulations are presented to demonstrate the theoretical results.     

一个具有年龄结构的单种群离散反应扩散模型的波前解

利用上下解方法研究了一个具有年龄结构的单种群离散反应扩散模型波前解的存在性,并证明了存在具有临界波速的波前解.     

Traveling wave solutions of a reaction diffusion model for competing pioneer and climax species

Presented is a reaction-diffusion model for the interaction of pioneer and climax species. For certain parameters the system exhibits bistability and traveling wave solutions. Specifically, we show that when the climax species diffuses at a slow rate there are traveling wave solutions which correspond to extinction waves of either the pioneer or climax species. A leading order analysis is used in the one-dimensional spatial case to estimate the wave speed sign that determines which species becomes extinct. Results of these analyses are then compared to numerical simulations of wave front propagation for the model on one and two-dimensional spatial domains. A simple mechanism for harvesting is also introduced.     

Repulsion Effect on Superinfecting Virions by Infected Cells

In this paper, the repulsion effect of superinfecting virion by infected cells is studied by a reaction diffusion equation model for virus infection dynamics. In this model, the diffusion of virus depends not only on its concentration gradient but also on the concentration of infected cells. The basic reproduction number, linear stability of steady states, spreading speed and existence of traveling wave solutions for the model are discussed. It is shown that viruses spread more rapidly with the repulsion effect of infected cells on superinfecting virions, than with random diffusion only. For our model, the spreading speed of free virus is not consistent with the minimal traveling wave speed. With our general model, numerical computations of the spreading speed show that the repulsion of superinfecting virion promotes the spread of virus, which confirms, not only qualitatively but also quantitatively, the experimental result of Doceul et al. (Science 327:873–876, 2010).     

Speeds of invasion in a model with strong or weak Allee effects

We study an invasion model based on a reaction-diffusion equation with an Allee effect. We use a special, piecewise-linear, population growth rate. This function allows us to obtain traveling wave solutions and to compute wave speeds for a full range of Allee effects, including weak Allee effects. Some investigators claim that linearization fails to give the correct speed of invasion if there is an Allee effect. We show that the minimum speed for a sufficiently weak Allee may, in fact, be the same as that derived by means of linearization.     

Neutral Genetic Patterns for Expanding Populations with Nonoverlapping Generations

We investigate the inside dynamics of solutions to integrodifference equations to understand the genetic consequences of a population with nonoverlapping generations undergoing range expansion. To obtain the inside dynamics, we decompose the solution into neutral genetic components. The inside dynamics are given by the spatiotemporal evolution of the neutral genetic components. We consider thin-tailed dispersal kernels and a variety of per capita growth rate functions to classify the traveling wave solutions as either pushed or pulled fronts. We find that pulled fronts are synonymous with the founder effect in population genetics. Adding overcompensation to the dynamics of these fronts has no impact on genetic diversity in the expanding population. However, growth functions with a strong Allee effect cause the traveling wave solution to be a pushed front preserving the genetic variation in the population. In this case, the contribution of each neutral fraction can be computed by a simple formula dependent on the initial distribution of the neutral fractions, the traveling wave solution, and the asymptotic spreading speed.     

Dynamics of prostate cancer stem cells with diffusion and organism response

We develop a systems based model for prostate cancer, as a sub-system of the organism. We accomplish this in two stages. We first start with a general ODE that includes organism response terms. Then, to account for normally observed spatial diffusion of cell populations, the ODE is extended to a PDE that includes spatial terms. Numerical solutions of the full PDE are provided, and are indicative of traveling wave fronts. This motivates the use of a well known transformation to derive a canonically related (non-linear) system of ODEs for traveling wave solutions. For biological feasibility, we show that the non-negative cone for the traveling wave system is time invariant. We also prove that the traveling waves have a unique global attractor. Biologically, the global attractor would be the limit for the avascular tumor growth. We conclude with comments on clinical implications of the model.     

Effects of long-range taxis and population pressure on the range expansion of invasive species in heterogeneous environments

We consider a new model for biological invasions in periodic patchy environments, in which long-range taxis and population pressure are incorporated in the framework of reaction-diffusion-advection equations. We assume that long-range taxis is induced by a weighted integral of stimuli within a certain sensing range. Population pressure is incorporated in the diffusion coefficient that linearly increases with population density. We first analyze the model in the absence of population pressure and demonstrate how the sensing length of long-range taxis influences the range expansion pattern of invasive species and its rate of spread. The effects of population pressure are examined for both homogeneous and periodic patchy environments. For the homogeneous environment, an exact and explicit traveling wave solution and the spreading speed are obtained. For the periodic patchy environment, we find numerically that a population starting from any localized distribution evolves to a traveling periodic wave if the null solution of the RDA equation is locally unstable, and that the traveling wave speed significantly increases with increasing population pressure. Furthermore, the population pressure and taxis intensity synergistically enhance the spreading speed when they are increased together.     

格上时滞Lotka-Volterra合作系统的波前解

研究了定义在格上并具有时滞的Lotka-Volterra合作系统的波前解.通过构造上下解得到了波前解的存在性,借助于比较原理和渐近传播理论得到了波前解的不存在性,进而在得到了波前解最小波速的充分条件.     

Existence of Traveling Waves for the Generalized F–KPP Equation

Variation in genotypes may be responsible for differences in dispersal rates, directional biases, and growth rates of individuals. These traits may favor certain genotypes and enhance their spatiotemporal spreading into areas occupied by the less advantageous genotypes. We study how these factors influence the speed of spreading in the case of two competing genotypes under the assumption that spatial variation of the total population is small compared to the spatial variation of the frequencies of the genotypes in the population. In that case, the dynamics of the frequency of one of the genotypes is approximately described by a generalized Fisher–Kolmogorov–Petrovskii–Piskunov (F–KPP) equation. This generalized F–KPP equation with (nonlinear) frequency-dependent diffusion and advection terms admits traveling wave solutions that characterize the invasion of the dominant genotype. Our existence results generalize the classical theory for traveling waves for the F–KPP with constant coefficients. Moreover, in the particular case of the quadratic (monostable) nonlinear growth–decay rate in the generalized F–KPP we study in detail the influence of the variance in diffusion and mean displacement rates of the two genotypes on the minimal wave propagation speed.     

On spreading speeds and traveling waves for growth and migration models in a periodic habitat

 It is shown that the methods previously used by the author [Wei82] and by R. Lui [Lui89] to obtain asymptotic spreading results and sometimes the existence of traveling waves for a discrete-time recursion with a translation invariant order preserving operator can be extended to a recursion with a periodic order preserving operator. The operator can be taken to be the time-one map of a continuous time reaction-diffusion model, or it can be a more general model of time evolution in population genetics or population ecology in a periodic habitat. Methods of estimating the speeds of spreading in various directions will also be presented. Received: 12 July 2001 / Revised version: 19 July 2002 / Published online: 17 October 2002 Mathematics Subject Classification (2000): 92D40, 92D25, 35K55, 35K57, 35B40 Keywords or phrases: Periodic – Spreading speed – Traveling wave     

Impact of Directed Movement on Invasive Spread in Periodic Patchy Environments

To address the effect of taxis of invasive animals on their spreading speed in heterogeneous environments, we deal with an advection-diffusion-reaction equation (ADR) in a periodic patchy environment. Two-types of advection that spatially vary depending on environmental heterogeneity are taken into consideration: a stepwise taxis function and a saw-like taxis function. We first analyze the ADR with the stepwise taxis advection, and derive an invasion criterion. When the invasion criterion holds, an initially localized population evolves to a traveling periodic wave (TPW). The asymptotic speed of the TPW is found to be equal to the minimal speed of the TPW analytically derived. Thus, we examine how the minimal speed is influenced by the taxis. The major results are: (1)?As the magnitude of the taxis toward favorable patches increases, invasion becomes more feasible. However, the spreading speed increases at first, and then decreases to show a one-humped curve against the magnitude of the taxis; (2)?As the scale of fragmentation in the patchy environment is increased, the spreading speed increases when the magnitude of the taxis is small, while it decreases when the magnitude of the taxis becomes sufficiently large. These characteristic features qualitatively apply to the ADR model with the saw-like taxis function.     

Integrodifference equations, Allee effects, and invasions

 Models that describe the spread of invading organisms often assume no Allee effect. In contrast, abundant observational data provide evidence for Allee effects. We study an invasion model based on an integrodifference equation with an Allee effect. We derive a general result for the sign of the speed of invasion. We then examine a special, linear–constant, Allee function and introduce a numerical scheme that allows us to estimate the speed of traveling wave solutions. Received: 14 April 2000 / Revised version: 23 December 2000 / Published online: 8 February 2002     

Propagation of HBV with spatial dependence

A mathematical model is proposed to simulate the hepatitis B virus (HBV) infection with spatial dependence. The existence of traveling waves is established via the geometric singular perturbation method. Numerical simulations show that the model admits non-monotone traveling profiles. Influences of various parameters on the minimum wave speed are also discussed.     

Refinement and Pattern Formation in Neural Circuits by the Interaction of Traveling Waves with Spike-Timing Dependent Plasticity

Traveling waves in the developing brain are a prominent source of highly correlated spiking activity that may instruct the refinement of neural circuits. A candidate mechanism for mediating such refinement is spike-timing dependent plasticity (STDP), which translates correlated activity patterns into changes in synaptic strength. To assess the potential of these phenomena to build useful structure in developing neural circuits, we examined the interaction of wave activity with STDP rules in simple, biologically plausible models of spiking neurons. We derive an expression for the synaptic strength dynamics showing that, by mapping the time dependence of STDP into spatial interactions, traveling waves can build periodic synaptic connectivity patterns into feedforward circuits with a broad class of experimentally observed STDP rules. The spatial scale of the connectivity patterns increases with wave speed and STDP time constants. We verify these results with simulations and demonstrate their robustness to likely sources of noise. We show how this pattern formation ability, which is analogous to solutions of reaction-diffusion systems that have been widely applied to biological pattern formation, can be harnessed to instruct the refinement of postsynaptic receptive fields. Our results hold for rich, complex wave patterns in two dimensions and over several orders of magnitude in wave speeds and STDP time constants, and they provide predictions that can be tested under existing experimental paradigms. Our model generalizes across brain areas and STDP rules, allowing broad application to the ubiquitous occurrence of traveling waves and to wave-like activity patterns induced by moving stimuli.     

Population spread in patchy landscapes under a strong Allee effect

Many species of invasive insects establish and spread in regions around the world, causing enormous economical and environmental damage, in particular in forests. Some of these insects are subject to an Allee effect whereby the population must surpass a certain threshold in order to establish. Recent studies have examined the possibility of exploiting an Allee effect to improve existing control strategies. Forests and most other ecosystems show natural spatial variation, and human activities frequently increase the degree of spatial heterogeneity. It is therefore imperative to understand how the interplay between this spatial variation and individual movement behavior affects the overall speed of spread of an invasion. To this end, we study an integrodifference equation model in a patchy landscape and with Allee growth dynamics. Movement behavior of individuals varies according to landscape quality. Our study focuses on how the speed of the resulting traveling periodic wave depends on the interaction between landscape fragmentation, patch-dependent dispersal, and Allee population dynamics.     

Traveling wave solutions of a nonlinear reaction–advection equation

 We establish the existence of traveling wave solutions for a nonlinear partial differential equation that models a logistically growing population whose movement is governed by an advective process. Conditions are presented for which traveling wave solutions exist and for which they are stable to small semi-finite domain perturbations. The wave is of mathematical interest because its behavior is determined by a singular differential equation and those with small speed of propagation steepen into a shock-like solutions. Finally, we indicate that the smoothing presence of diffusion allows wave persistence when an advective slow moving wave may collapse. Received: 24 November 1997 / Revised version: 13 July 1998