Limits to prediction in a size-resolved pelagic ecosystem model

Oxford University Press
Publication Type:
Journal Article
Journal Of Plankton Research, 2010, 32 (8), pp. 1131 - 1146
Issue Date:
Full metadata record
Files in This Item:
Filename Description Size
Thumbnail2009006616OK.pdf720.29 kB
Adobe PDF
A size-resolved pelagic ecosystem model has been developed based on a continuous (with size) set of model equations and using allometric relationships to specify size-dependent physiological rates. Numerical experiments with identical model equations but different initial conditions and size-class distributions are used to investigate inherent limits to prediction of instantaneous state from an initial condition. The simulations have relatively constant physical forcings, such as solar radiation, to emphasize the dynamical properties of the size-resolved model. Initial condition experiments show that perturbations of 1, 0.1, 0.01, 0.001 and 0.0001% of the initial biomass of individual size-classes from a flat size spectrum lead to equal spread of model trajectories. The greatest divergence of trajectories occurs when a 2.7 µm equivalent spherical radius phytoplankton size-class blooms. This divergence has a finite-time Lyapunov exponent of 0.21 day-1 and a prediction time of 33 days for a precision of 10-3 mol N m-3. Large member ensembles can approximately halve the effect of growth of initial condition perturbations on prediction.
Please use this identifier to cite or link to this item: