Intrinsic dynamics and stability properties of size-structured pelagic ecosystem models

Ivan D. Lima, Donald Olson, Scott C. Doney

Research output: Contribution to journalArticle

47 Citations (Scopus)

Abstract

Although pelagic ecosystem models, coupled with ocean circulation models, are being widely used to quantify fluxes of nutrients and carbon in the ocean at regional and basin scales, relatively little work has been done on understanding their intrinsic dynamics independently of physical forcing. In this study, the dynamics of three common formulations for the NPZD class of models (nutrient, phytoplankton, zooplankton and detritus) using two different types of predation functional response are analysed and compared. Our goal is to characterize the stability properties of this class of models with respect to variations in light and total nutrient concentrations. Despite important structural differences, the different model formulations all show asymptotic stable equilibria at low total nutrient concentrations and high to moderate light intensities, and limit cycle oscillations at low light intensities and high total nitrogen concentrations. Limit cycles are formed through a Hopf bifurcation as a phase lag develops between predator (Zi) and prey (Pi or Zi-1) due to a relatively sharp increase in the growth rate of the prey in relation to that of the predator. The use of variable preferences in the functional response provides a density-dependent mechanism that allows the system to self-regulate, increasing system stability considerably, but does not eliminate the instabilities completely. Instabilities occur at light and nutrient levels that correspond to those observed on the bottom of the euphotic zone near the nutricline, where the deep chlorophyll maximum is usually located. This suggests that the deep chlorophyll maximum might be dynamically unstable. This dynamic disequilibrium in species composition and the characteristically long transient times would allow species persistence in the presence of seasonal and mesoscale variations and provide a mechanism for species coexistence in the relatively homogeneous open ocean environment, thereby providing a potential solution for the 'paradox of plankton'. In the multi-species models, the higher diversity of species (wider range of values for the biological parameters) allows biological activity (photosynthesis, grazing and predation) to occur under a wider range of light and nutrient conditions, resulting in higher primary and secondary production and lower nutrient concentrations at light intensities equivalent to those in the upper part of the euphotic zone, than in the single-species model.

Original languageEnglish (US)
Pages (from-to)533-556
Number of pages24
JournalJournal of Plankton Research
Volume24
Issue number6
StatePublished - Jul 23 2002

Fingerprint

pelagic ecosystem
ecosystems
nutrient
light intensity
nutrient content
functional response
oceans
euphotic zone
nutrients
chlorophyll
predation
predator
predators
species diversity
secondary productivity
secondary production
bifurcation
open ocean
disequilibrium
detritus

ASJC Scopus subject areas

  • Ecology, Evolution, Behavior and Systematics
  • Ecology
  • Aquatic Science

Cite this

Intrinsic dynamics and stability properties of size-structured pelagic ecosystem models. / Lima, Ivan D.; Olson, Donald; Doney, Scott C.

In: Journal of Plankton Research, Vol. 24, No. 6, 23.07.2002, p. 533-556.

Research output: Contribution to journalArticle

@article{af0ea4fbd5f14f2e9e35110273d19b53,
title = "Intrinsic dynamics and stability properties of size-structured pelagic ecosystem models",
abstract = "Although pelagic ecosystem models, coupled with ocean circulation models, are being widely used to quantify fluxes of nutrients and carbon in the ocean at regional and basin scales, relatively little work has been done on understanding their intrinsic dynamics independently of physical forcing. In this study, the dynamics of three common formulations for the NPZD class of models (nutrient, phytoplankton, zooplankton and detritus) using two different types of predation functional response are analysed and compared. Our goal is to characterize the stability properties of this class of models with respect to variations in light and total nutrient concentrations. Despite important structural differences, the different model formulations all show asymptotic stable equilibria at low total nutrient concentrations and high to moderate light intensities, and limit cycle oscillations at low light intensities and high total nitrogen concentrations. Limit cycles are formed through a Hopf bifurcation as a phase lag develops between predator (Zi) and prey (Pi or Zi-1) due to a relatively sharp increase in the growth rate of the prey in relation to that of the predator. The use of variable preferences in the functional response provides a density-dependent mechanism that allows the system to self-regulate, increasing system stability considerably, but does not eliminate the instabilities completely. Instabilities occur at light and nutrient levels that correspond to those observed on the bottom of the euphotic zone near the nutricline, where the deep chlorophyll maximum is usually located. This suggests that the deep chlorophyll maximum might be dynamically unstable. This dynamic disequilibrium in species composition and the characteristically long transient times would allow species persistence in the presence of seasonal and mesoscale variations and provide a mechanism for species coexistence in the relatively homogeneous open ocean environment, thereby providing a potential solution for the 'paradox of plankton'. In the multi-species models, the higher diversity of species (wider range of values for the biological parameters) allows biological activity (photosynthesis, grazing and predation) to occur under a wider range of light and nutrient conditions, resulting in higher primary and secondary production and lower nutrient concentrations at light intensities equivalent to those in the upper part of the euphotic zone, than in the single-species model.",
author = "Lima, {Ivan D.} and Donald Olson and Doney, {Scott C.}",
year = "2002",
month = "7",
day = "23",
language = "English (US)",
volume = "24",
pages = "533--556",
journal = "Journal of Plankton Research",
issn = "0142-7873",
publisher = "Oxford University Press",
number = "6",

}

TY - JOUR

T1 - Intrinsic dynamics and stability properties of size-structured pelagic ecosystem models

AU - Lima, Ivan D.

AU - Olson, Donald

AU - Doney, Scott C.

PY - 2002/7/23

Y1 - 2002/7/23

N2 - Although pelagic ecosystem models, coupled with ocean circulation models, are being widely used to quantify fluxes of nutrients and carbon in the ocean at regional and basin scales, relatively little work has been done on understanding their intrinsic dynamics independently of physical forcing. In this study, the dynamics of three common formulations for the NPZD class of models (nutrient, phytoplankton, zooplankton and detritus) using two different types of predation functional response are analysed and compared. Our goal is to characterize the stability properties of this class of models with respect to variations in light and total nutrient concentrations. Despite important structural differences, the different model formulations all show asymptotic stable equilibria at low total nutrient concentrations and high to moderate light intensities, and limit cycle oscillations at low light intensities and high total nitrogen concentrations. Limit cycles are formed through a Hopf bifurcation as a phase lag develops between predator (Zi) and prey (Pi or Zi-1) due to a relatively sharp increase in the growth rate of the prey in relation to that of the predator. The use of variable preferences in the functional response provides a density-dependent mechanism that allows the system to self-regulate, increasing system stability considerably, but does not eliminate the instabilities completely. Instabilities occur at light and nutrient levels that correspond to those observed on the bottom of the euphotic zone near the nutricline, where the deep chlorophyll maximum is usually located. This suggests that the deep chlorophyll maximum might be dynamically unstable. This dynamic disequilibrium in species composition and the characteristically long transient times would allow species persistence in the presence of seasonal and mesoscale variations and provide a mechanism for species coexistence in the relatively homogeneous open ocean environment, thereby providing a potential solution for the 'paradox of plankton'. In the multi-species models, the higher diversity of species (wider range of values for the biological parameters) allows biological activity (photosynthesis, grazing and predation) to occur under a wider range of light and nutrient conditions, resulting in higher primary and secondary production and lower nutrient concentrations at light intensities equivalent to those in the upper part of the euphotic zone, than in the single-species model.

AB - Although pelagic ecosystem models, coupled with ocean circulation models, are being widely used to quantify fluxes of nutrients and carbon in the ocean at regional and basin scales, relatively little work has been done on understanding their intrinsic dynamics independently of physical forcing. In this study, the dynamics of three common formulations for the NPZD class of models (nutrient, phytoplankton, zooplankton and detritus) using two different types of predation functional response are analysed and compared. Our goal is to characterize the stability properties of this class of models with respect to variations in light and total nutrient concentrations. Despite important structural differences, the different model formulations all show asymptotic stable equilibria at low total nutrient concentrations and high to moderate light intensities, and limit cycle oscillations at low light intensities and high total nitrogen concentrations. Limit cycles are formed through a Hopf bifurcation as a phase lag develops between predator (Zi) and prey (Pi or Zi-1) due to a relatively sharp increase in the growth rate of the prey in relation to that of the predator. The use of variable preferences in the functional response provides a density-dependent mechanism that allows the system to self-regulate, increasing system stability considerably, but does not eliminate the instabilities completely. Instabilities occur at light and nutrient levels that correspond to those observed on the bottom of the euphotic zone near the nutricline, where the deep chlorophyll maximum is usually located. This suggests that the deep chlorophyll maximum might be dynamically unstable. This dynamic disequilibrium in species composition and the characteristically long transient times would allow species persistence in the presence of seasonal and mesoscale variations and provide a mechanism for species coexistence in the relatively homogeneous open ocean environment, thereby providing a potential solution for the 'paradox of plankton'. In the multi-species models, the higher diversity of species (wider range of values for the biological parameters) allows biological activity (photosynthesis, grazing and predation) to occur under a wider range of light and nutrient conditions, resulting in higher primary and secondary production and lower nutrient concentrations at light intensities equivalent to those in the upper part of the euphotic zone, than in the single-species model.

UR - http://www.scopus.com/inward/record.url?scp=0036315225&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0036315225&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0036315225

VL - 24

SP - 533

EP - 556

JO - Journal of Plankton Research

JF - Journal of Plankton Research

SN - 0142-7873

IS - 6

ER -