Combining endogenous and exogenous spatial variability in analytical population models

Analytically tractable models of dynamics in continuous space rarely incorporate both endogenous and exogenous spatial heterogeneity. We use spatial moment equations in combination with simulation models to analyze the combined effects of endogenous and exogenous variability on population viability in a simple single-population model where landscape heterogeneity and local population density both affect mortality rate. The equations partition the effects of heterogeneity into an effect of local crowding and an effect of habitat association caused by differential mortality. Exogenous heterogeneity in mortality rate increases population viability through habitat association and decreases it through increased crowding; the net effect of exogenous heterogeneity is generally to improve population viability. This result is contrary to some (but not all) conclusions in the literature, which usually focus on the effects of fragmentation rather than the benefits of refuges to short-dispersing individuals.     

Analytical methods for predicting the behaviour of population models with general spatial interactions

Many biologists use population models that are spatial, stochastic and individual based. Analytical methods that describe the behaviour of these models approximately are attracting increasing interest as an alternative to expensive computer simulation. The methods can be employed for both prediction and fitting models to data. Recent work has extended existing (mean field) methods with the aim of accounting for the development of spatial correlations. A common feature is the use of closure approximations for truncating the set of evolution equations for summary statistics. We investigate an analytical approach for spatial and stochastic models where individuals interact according to a generic function of their distance; this extends previous methods for lattice models with interactions between close neighbours, such as the pair approximation. Our study also complements work by Bolker and Pacala (BP) [Theor. Pop. Biol. 52 (1997) 179; Am. Naturalist 153 (1999) 575]: it treats individuals as being spatially discrete (defined on a lattice) rather than as a continuous mass distribution; it tests the accuracy of different closure approximations over parameter space, including the additive moment closure (MC) used by BP and the Kirkwood approximation. The study is done in the context of an susceptible-infected-susceptible epidemic model with primary infection and with secondary infection represented by power-law interactions. MC is numerically unstable or inaccurate in parameter regions with low primary infection (or density-independent birth rates). A modified Kirkwood approximation gives stable and generally accurate transient and long-term solutions; we argue it can be applied to lattice and to continuous-space models as a substitute for MC. We derive a generalisation of the basic reproduction ratio, R(0), for spatial models.     

Analytic approximation of spatial epidemic models of foot and mouth disease

The effect of spatial heterogeneity in epidemic models has improved with computational advances, yet far less progress has been made in developing analytical tools for understanding such systems. Here, we develop two classes of second-order moment closure methods for approximating the dynamics of a stochastic spatial model of the spread of foot and mouth disease. We consider the performance of such ‘pseudo-spatial’ models as a function of R₀, the locality in disease transmission, farm distribution and geographically-targeted control when an arbitrary number of spatial kernels are incorporated. One advantage of mapping complex spatial models onto simpler deterministic approximations lies in the ability to potentially obtain a better analytical understanding of disease dynamics and the effects of control. We exploit this tractability by deriving analytical results in the invasion stages of an FMD outbreak, highlighting key principles underlying epidemic spread on contact networks and the effect of spatial correlations.     

On stochastic logistic population growth models with immigration and multiple births

This paper develops a stochastic logistic population growth model with immigration and multiple births. The differential equations for the low-order cumulant functions (i.e., mean, variance, and skewness) of the single birth model are reviewed, and the corresponding equations for the multiple birth model are derived. Accurate approximate solutions for the cumulant functions are obtained using moment closure methods for two families of model parameterizations, one for badger and the other for fox population growth. For both model families, the equilibrium size distribution may be approximated well using the Normal approximation, and even more accurately using the saddlepoint approximation. It is shown that in comparison with the corresponding single birth model, the multiple birth mechanism increases the skewness and the variance of the equilibrium distribution, but slightly reduces its mean. Moreover, the type of density-dependent population control is shown to influence the sign of the skewness and the size of the variance.     

Relatedness in spatially structured populations with empty sites: An approach based on spatial moment equations

Taking into account the interplay between spatial ecological dynamics and selection is a major challenge in evolutionary ecology. Although inclusive fitness theory has proven to be a very useful tool to unravel the interactions between spatial genetic structuring and selection, applications of the theory usually rely on simplifying demographic assumptions. In this paper, I attempt to bridge the gap between spatial demographic models and kin selection models by providing a method to compute approximations for relatedness coefficients in a spatial model with empty sites. Using spatial moment equations, I provide an approximation of nearest-neighbour relatedness on random regular networks, and show that this approximation performs much better than the ordinary pair approximation. I discuss the connection between the relatedness coefficients I define and those used in population genetics, and sketch some potential extensions of the theory.     

Spatially explicit neutral models for population genetics and community ecology: Extensions of the Neyman-Scott clustering process

Spatially explicit models relating to plant populations have developed little since Felsenstein (1975) pointed out that if limited seed dispersal causes clustering of individuals, such models cannot reach an equilibrium. This paper aims to resolve this issue by modifying the Neyman-Scott cluster point process. The new point processes are dynamic models with random immigration, and the continuous increase in the clustering of individuals stops at some level. Hence, an equilibrium state is achieved, and new individual-based spatially explicit neutral coalescent models are established. By fitting the spatial structure at equilibrium to individual spatial distribution data, we can indirectly estimate seed dispersal and effective population density. These estimates are improved when genetic data are available, and become even more sophisticated if spatial distribution and genetic data pertaining to the offspring are also available.     

The mechanistic basis of discrete-time population models: The role of resource partitioning and spatial aggregation

The purpose of this paper is to present a unified view to understand mechanistic basis of various discrete-time population models from the viewpoints of resource partitioning and spatial aggregation of individuals. A first-principles derivation is presented of a new population model which incorporates both scramble and contest competition using a site-based framework in which individuals are distributed over discrete resource sites. The derived model has parameters relating to the way of resource partitioning and the degree of spatial aggregation of individuals, respectively. The model becomes various population models in various limits in these parameters. This model thus provides a unified view to understand how various population models are interrelated. The dependence of the stability of the model on the parameters is also examined.