Numerical equilibrium analysis for structured consumer resource models.

In this paper, we present methods for a numerical equilibrium and stability analysis for models of a size structured population competing for an unstructured resource. We concentrate on cases where two model parameters are free, and thus existence boundaries for equilibria and stability boundaries can be defined in the (two-parameter) plane. We numerically trace these implicitly defined curves using alternatingly tangent prediction and Newton correction. Evaluation of the maps defining the curves involves integration over individual size and individual survival probability (and their derivatives) as functions of individual age. Such ingredients are often defined as solutions of ODE, i.e., in general only implicitly. In our case, the right-hand sides of these ODE feature discontinuities that are caused by an abrupt change of behavior at the size where juveniles are assumed to turn adult. So, we combine the numerical solution of these ODE with curve tracing methods. We have implemented the algorithms for "Daphnia consuming algae" models in C-code. The results obtained by way of this implementation are shown in the form of graphs.

Numerical methods and parameter estimation of a structured population model with discrete events in the life history.

We consider two numerical methods for the solution of a physiologically structured population (PSP) model with multiple life stages and discrete event reproduction. The model describes the dynamic behaviour of a predator-prey system consisting of rotifers predating on algae. The nitrate limited algal prey population is modelled unstructured and described by an ordinary differential equation (ODE). The formulation of the rotifer dynamics is based on a simple physiological model for their two life stages, the egg and the adult stage. An egg is produced when an energy buffer reaches a threshold value. The governing equations are coupled partial differential equations (PDE) with initial and boundary conditions. The population models together with the equation for the dynamics of the nutrient result in a chemostat model. Experimental data are used to estimate the model parameters. The results obtained with the explicit finite difference (FD) technique compare well with those of the Escalator Boxcar Train (EBT) method. This justifies the use of the fast FD method for the parameter estimation, a procedure which involves repeated solution of the model equations.