Asymptotic solutions of differential equations with time delay in population dynamics

Oscillations governed by generalized Volterra-Gause-Witt equations to include retardation effects in population dynamics are considered. Models containing either small or significant time delay are discussed. The Krylov-Bogoliubov-Mitropolskii perturbation method and its extension for differential equations with retarded argument is used.     

State-dependent life-history equations

Matrix population models provide a natural tool to analyse state-dependent life-history strategies. Reproductive value and the intrinsic rate of natural increase under a strategy, and the optimal life-history strategy can all be easily characterised using projection matrices. The resultant formulae, however, are not directly comparable with the corresponding formulae for age structured populations such as Lotka's equations and Fisher's formula for reproductive value. This is because formulae involving projection matrices lose track of what happens to an individual over its lifetime and are only concerned with expected numbers of descendants one time step in the future. In contrast the usual age-dependent formulae explicitly followed a single individual through from birth to death.In this paper I show how the state-dependent formulae can be rewritten to be directly comparable with the standard age-structured formulae. Although the formulae are intuitively obvious the decomposition into current and future reproductive success differs from that previously given and is, I suggest, a more natural definition. The derivation of appropriate equations for optimal life-histories relies on results from dynamic programming theory; and is much more general and easier than previous derivations.The value of rewriting projection matrix results in terms of the lifetime of an individual organism is illustrated by an example in which the optimal plastic response to an environment is derived.     

Research of population with impulsive perturbations based on dynamics of a neutral delay equation and ecological quality system

This paper studies the global behaviors of a nonlinear autonomous neutral delay differential population model with impulsive perturbation. This model may be suitable for describing the dynamics of population with long larval and short adult phases. It is shown that the system may have global stability of the extinction and positive equilibria, or grow without being bounded under some conditions.     

Stochastic hybrid delay population dynamics: well-posed models and extinction

Nonlinear differential equations have been used for decades for studying fluctuations in the populations of species, interactions of species with the environment, and competition and symbiosis between species. Over the years, the original non-linear models have been embellished with delay terms, stochastic terms and more recently discrete dynamics. In this paper, we investigate stochastic hybrid delay population dynamics (SHDPD), a very general class of population dynamics that comprises all of these phenomena. For this class of systems, we provide sufficient conditions to ensure that SHDPD have global positive, ultimately bounded solutions, a minimum requirement for a realistic, well-posed model. We then study the question of extinction and establish conditions under which an ecosystem modelled by SHDPD is doomed.     

Neural rate equations for bursting dynamics derived from conductance-based equations

A method of obtaining rate equations from conductance-based equations is developed and applied to fast-spiking and bursting neocortical neurons. It involves splitting systems of conductance-based equations into fast and slow subsystems, and averaging the effects of fast terms that drive the slowly varying quantities by showing that their average is closely proportional to the firing rate. The dependence of the firing rate on the injected current is then approximated in the analysis. The resulting behavior of the slow variables is then substituted back into the fast equations, with the further approximation of replacing the fast voltages in these terms by effective values. For bursting neurons the method yields two coupled limit-cycle oscillators: a self-exciting oscillator for the slow variables that commences limit-cycle oscillations at a critical current and modulates a fast spike-generating oscillator, thereby leading to slowly modulated bursts with a group of spikes in each burst. The dynamics of these coupled oscillators are then verified against those of the conductance-based equations. Finally, it is shown how to place the results in a form suitable for use in mean-field equations for neural population dynamics.     

Modelling multi-pulse population dynamics from ultrafast spectroscopy

Current advanced laser, optics and electronics technology allows sensitive recording of molecular dynamics, from single resonance to multi-colour and multi-pulse experiments. Extracting the occurring (bio-) physical relevant pathways via global analysis of experimental data requires a systematic investigation of connectivity schemes. Here we present a Matlab-based toolbox for this purpose. The toolbox has a graphical user interface which facilitates the application of different reaction models to the data to generate the coupled differential equations. Any time-dependent dataset can be analysed to extract time-independent correlations of the observables by using gradient or direct search methods. Specific capabilities (i.e. chirp and instrument response function) for the analysis of ultrafast pump-probe spectroscopic data are included. The inclusion of an extra pulse that interacts with a transient phase can help to disentangle complex interdependent pathways. The modelling of pathways is therefore extended by new theory (which is included in the toolbox) that describes the finite bleach (orientation) effect of single and multiple intense polarised femtosecond pulses on an ensemble of randomly oriented particles in the presence of population decay. For instance, the generally assumed flat-top multimode beam profile is adapted to a more realistic Gaussian shape, exposing the need for several corrections for accurate anisotropy measurements. In addition, the (selective) excitation (photoselection) and anisotropy of populations that interact with single or multiple intense polarised laser pulses is demonstrated as function of power density and beam profile. Using example values of real world experiments it is calculated to what extent this effectively orients the ensemble of particles. Finally, the implementation includes the interaction with multiple pulses in addition to depth averaging in optically dense samples. In summary, we show that mathematical modelling is essential to model and resolve the details of physical behaviour of populations in ultrafast spectroscopy such as pump-probe, pump-dump-probe and pump-repump-probe experiments.     

Numerical bifurcation analysis of delay differential equations arising from physiological modeling

This paper has a dual purpose. First, we describe numerical methods for continuation and bifurcation analysis of steady state solutions and periodic solutions of systems of delay differential equations with an arbitrary number of fixed, discrete delays. Second, we demonstrate how these methods can be used to obtain insight into complex biological regulatory systems in which interactions occur with time delays: for this, we consider a system of two equations for the plasma glucose and insulin concentrations in a diabetic patient subject to a system of external assistance. The model has two delays: the technological delay of the external system, and the physiological delay of the patient's liver. We compute stability of the steady state solution as a function of two parameters, compare with analytical results and compute several branches of periodic solutions and their stability. These numerical results allow to infer two categories of diabetic patients for which the external system has different efficiency.     

Complex discrete dynamics from simple continuous population models

Nonoverlapping generations have been classically modelled as difference equations in order to account for the discrete nature of reproductive events. However, other events such as resource consumption or mortality are continuous and take place in the within-generation time. We have realistically assumed a hybrid ODE bidimensional model of resources and consumers with discrete events for reproduction. Numerical and analytical approaches showed that the resulting dynamics resembles a Ricker map, including the doubling route to chaos. Stochastic simulations with a handling-time parameter for indirect competition of juveniles may affect the qualitative behaviour of the model.     

逻辑斯谛增长导致的复杂种群动态及其实验验证

教科书通常用“S”型曲线表示逻辑斯谛增长模型导致的种群动态,事实上,这一模型可以产生包括稳定平衡、周期性振荡、混沌等多种多样的种群动态模式。介绍了如何用Kicker图示的方法分析逻辑斯谛增长可能导致的种群动态模式。同时介绍了相关的实验研究案例,这些实验工作大都是以生活史周期较短的昆虫为材料在相对简单的实验室环境中完成的。     

Steady states and outbreaks of two-phase nonlinear age-structured model of population dynamics with discrete time delay

In this paper we study the nonlinear age-structured model of a polycyclic two-phase population dynamics including delayed effect of population density growth on the mortality. Both phases are modelled as a system of initial boundary values problem for semi-linear transport equation with delay and initial problem for nonlinear delay ODE. The obtained system is studied both theoretically and numerically. Three different regimes of population dynamics for asymptotically stable states of autonomous systems are obtained in numerical experiments for the different initial values of population density. The quasi-periodical travelling wave solutions are studied numerically for the autonomous system with the different values of time delays and for the system with oscillating death rate and birth modulus. In both cases it is observed three types of travelling wave solutions: harmonic oscillations, pulse sequence and single pulse.     

Dose-structured population dynamics

Applied population dynamics modeling is relied upon with increasing frequency to quantify how human activities affect human and non-human populations. Current techniques include variously the population's spatial transport, age, size, and physiology, but typically not the life-histories of exposure to other important things occurring in the ambient environment, such as chemicals, heat, or radiation. Consequently, the effects of such 'abiotic' aspects of an ecosystem on populations are only currently addressed through individual-based modeling approaches that despite broad utility are limited in their applicability to realistic ecosystems [V. Grimm, Ten years of individual-based modeling in ecology: what have we learned and what could we learn in the future? Ecol. Model. 115 (1999) 129-148][1]. We describe a new category of population dynamics modeling, wherein population dynamical states of the biotic phases are structured on dose, and apply this framework to demonstrate how chemical species or other ambient aspects can be included in population dynamics in three separate examples involving growth suppression in fish, inactivation of microorganisms with ultraviolet irradiation, and metabolic lag in population growth. Dose-structuring is based on a kinematic approach that is a simple generalization of age-structuring, views the ecosystem as a multi-component mixture with reacting biotic/abiotic components. The resulting model framework accommodates (a) different memories of exposure as in recovery from toxic ambient conditions, (b) differentiation between exogenous and endogenous sources of variation in population response, and (c) quantification of acute or sub-acute effects on populations arising from life-history exposures to abiotic species. Classical models do not easily address the very important fact that organisms differ and have different experiences over their life cycle. The dose structuring is one approach to incorporate some of these elements into the existing structures of the classical models, while retaining many of the features (and other limitations) of classical models.     

Modeling reflex asymmetries with implicit delay differential equations

Neuromuscular reflexes with time-delayed negative feedback, such as the pupil light reflex, have different rates depending on the direction of movement. This asymmetry is modeled by an implicit first-order delay differential equation in which the value of the rate constant depends on the direction of movement. Stability analyses are presented for the cases when the rate is: (1) an increasing and (2) a decreasing function of the direction of movement. It is shown that the stability of equilibria in these dynamical systems depends on whether the rate constant is a decreasing or increasing function. In particular, when the asymmetry has the shape of an increasing step function, it is possible to have stability which is independent of the value of the time delay or the steepness (i.e., gain) of the negative feedback.     

Small population effects in stochastic population dynamics

The discreteness of units of small populations can produce large fluctuations from a classical continuous representation, especially when null populations play a crucial role. These belong to what are here referred to as emergent and evanescent species. A few model biological systems are introduced in which this is the case, as well as a toy model that suggests a path to avoid the associated mathematical complexities. The corresponding division into null and non-null population sectors is carried out to leading order for the model systems, with promising results. Supported in part by DOE Office of Basic Science, Chemical Division. Reported at SMB03, the August 2003 meeting of the Society for Mathematical Biology.     

The equations of population growth

A solution in quadratures is found for a special class of differential equations suggested by studies in the growth of competing populations.     

The effect of habitat fragmentation on cyclic population dynamics: a reduction to ordinary differential equations

Habitat fragmentation is known to be a key factor affecting population dynamics. In a previous study by Strohm and Tyson (Bull Math Biol 71:1323?C1348, 2009), the effect of habitat fragmentation on cyclic population dynamics was studied using spatially explicit predator?Cprey models with four different sets of reaction terms. The difficulty with spatially explicit models is that often analytical tractability is lost and the mechanisms behind the behaviour of the models are difficult to analyse. In this study, we employ a simplification procedure based on a Fourier series first-term truncation of the spatially explicit models Strohm and Tyson (Bull Math Biol 71:1323?C1348, 2009) to obtain spatially implicit models. These simpler models capture the main features of the spatially explicit models and can be used to explain the dynamics observed by Strohm and Tyson. We find that the spatially implicit models and the spatially explicit models produce similar responses to habitat fragmentation for larger high-quality patch sizes. Additionally, we find that the critical patch size of the spatially implicit models provides an upper bound on the critical patch size of the spatially explicit models. Finally, we derive an approximation of the multi-patch habitat by a single-patch habitat with partial flux boundary conditions which allows for a lower bound on the critical patch size to be calculated.