Mathematical conservation ecology: a one-predator-two-prey system as case study.

A method is presented to analyse the long-term stochastic dynamics of a biological population that is at risk of extinction. From the full ecosystem the method extracts the minimal information to describe the long-term dynamics of that population by a stochastic logistic system. The method is applied to a one-predator-two-prey model. The choice of this example is motivated by a study on the near-extinction of a porcupine population by mountain lions whose presence is facilitated by mule deer taking advantage of a change in land use. The risk of extinction is quantified by the expected time of extinction of the population.

A two-component model of host-parasitoid interactions: determination of the size of inundative releases of parasitoids in biological pest control.

A two-component differential equation model is formulated for a host-parasitoid interaction. Transient dynamics and population crashes of this system are analysed using differential inequalities. Two different cases can be distinguished: either the intrinsic growth rate of the host population is smaller than the maximum growth rate of the parasitoid or vice versa. In the latter case, the initial ratio of parasitoids to hosts should exceed a given threshold, in order to (temporarily) halt the growth of the host population. When not only oviposition but also host-feeding occurs the dynamics do not change qualitatively. In the case that the maximum growth rate of the parasitoid population is smaller than the intrinsic growth rate of the host, a threshold still exists for the number of parasitoids in an inundative release in order to limit the growth of the host population. The size of an inundative release of parasitoids, which is necessary to keep the host population below a certain level, can be determined from the two-component model. When parameter values for hosts and parasitoids are known, an effective control of pests can be found. First it is determined whether the parasitoids are able to suppress their hosts fully. Moreover, using our simple rule of thumb it can be assessed whether suppression is also possible when the relative growth rate of the host population exceeds that of the parasitoid population. With a numerical investigation of our simple system the design of parasitoid release strategies for specific situations can be computed.

Stochastic epidemics: the probability of extinction of an infectious disease at the end of a major outbreak.

The aim of this study is to derive an asymptotic expression for the probability that an infectious disease will disappear from a population at the end of a major outbreak ('fade-out'). The study deals with a stochastic SIR-model. Local asymptotic expansions are constructed for the deterministic trajectories of the corresponding deterministic system, in particular for the deterministic trajectory starting in the saddle point. The analytical expression for the probability of extinction is derived by asymptotically solving a boundary value problem based on the Fokker-Planck equation for the stochastic system. The asymptotic results are compared with results obtained by random walk simulations.

Stochastic epidemics: major outbreaks and the duration of the endemic period.

A study is made of a two-dimensional stochastic system that models the spread of an infectious disease in a population. An asymptotic expression is derived for the probability that a major outbreak of the disease will occur in case the number of infectives is small. For the case that a major outbreak has occurred, an asymptotic approximation is derived for the expected time that the disease is in the population. The analytical expressions are obtained by asymptotically solving Dirichlet problems based on the Fokker-Planck equation for the stochastic system. Results of numerical calculations for the analytical expressions are compared with simulation results.