NPZ model
An NPZ model is the most basic abstract representation, expressed as a mathematical model, of a pelagic ecosystem which examines the interrelationships between quantities of nutrients, phytoplankton and zooplankton as time-varying states which depend only on the relative concentrations of the various states at the given time.[1]
One goal in pelagic ecology is to understand the interactions among available nutrients (i.e. the essential resource base), phytoplankton and zooplankton. The most basic models to shed light on this goal are called nutrient-phytoplankton-zooplankton (NPZ) models. These models are a subset of Ecosystem models.[2]
Example
[edit]An unrealistic but instructive example of an NPZ model is provided in Franks et al. (1986)[3] (FWF-NPZ model). It is a system of ordinary differential equations that examines the time evolution of dissolved and assimilated nutrients in an ideal upper water column consisting of three state variables corresponding to amounts of nutrients (N), phytoplankton (P) and zooplankton (Z). This closed system model is shown in the figure to the right which also shows the "flow" directions of each state quantity.
These interactions, assumed to be spatial homogeneous (and thus is termed a "zero-dimensional" model) are described in general terms as follows
Rate of Change in P | = | (Nutrient uptake from N to P) minus (Mortality of P) minus (Grazing by Z) |
Rate of Change in Z | = | (Growth Efficiency) times (Grazing by Z) minus (Mortality of Z) |
Rate of Change in N | = | minus (Nutrient uptake from N to P) plus (1-Growth Efficiency) x (Grazing by Z) plus (Mortality of P) plus (Mortality of Z). |
This NPZ model can now be cast as a system of first order differential equations:
where the parameters and variables are defined in the table below along with nominal values for a "standard environment"
Variable | Definition | Value Used in Example |
---|---|---|
maximum phytoplankton growth rate | ||
nutrient concentration | Initially Mole N liter | |
half saturation constant for nutrients | Mole N liter | |
phytoplankton stock size | Initially Mole N liter | |
phytoplankton mortality rate | ||
zooplankton stock size | Initially Mole N liter | |
zooplankton growth efficiency | ||
maximum zooplankton ration | ||
Ivlev constant | Mole N liter | |
zooplankton mortality rate |
An example of a 60 day sequence for the values shown is depicted in the figure to the right. Each state is color coded (Nutrient – black, Phytoplankton – green and Zooplankton – blue). Note that the initial nutrient concentration is rapidly consumed resulting in a phytoplankton bloom until the zooplankton begin aggressive grazing around day 10. Eventually both populations drop to a very low level and a high nutrient concentration remains. In the next section more sophistication is applied to the model in order to increase realism.
More Sophisticated NPZ Models
[edit]The Franks et al. (1986) [3] work has inspired significant analysis from other researchers[4][5] but is overly simplistic to capture the complexity of actual pelagic communities. A more realistic NPZ model would simulate control of primary production by incorporating mechanisms to simulate seasonally varying sunlight and decreasing illumination with depth. Evans and Parslow (1985) developed an NPZ model which includes these mechanisms[6] and forms the basis of the following example (see also Denman and Pena (1999)[7]).
A 200 day sequence resulting from this configuration of the FWF-NPZ model is shown in the figure to the right. Each state is color coded (Nutrient – black, Phytoplankton – green and Zooplankton – blue). Several interesting features in the model output are easily observed. First, a spring bloom occurs in the first 20 days or so, where the high nutrient concentrations are consumed by the phytoplankton causing an inverse relationship which is halted by a rise in zooplankton concentration eventually settling into a sustained steady-state solution for the remainder of the summer. Another bloom, not as pronounced as in the spring, occurs in the fall with a remixing of nutrients into the water column.
References
[edit]- ^ C.B. Miller and P.A. Wheeler, Biological Oceanography, Second Edition, Wiley-Blackwell, 2012.
- ^ Marina Lévy (16 February 2015). "Exploration of the critical depth hypothesis with a simple NPZ model". ICES Journal of Marine Science. 72 (6). Oxford University Press: 1916–1925. doi:10.1093/icesjms/fsv016.
- ^ a b c P.J.S. Franks, J.S. Wroblewski and G.R. Flierl, Behavior of a simple plankton model with food-level acclimation by herbivores, Marine Biology, 91, pp. 121-129, 1986.
- ^ S. Busenberg, S.K. Kumar, P. Austin and G. Wake, The dynamics of a model of plankton-nutrient interaction, Bulletin of Mathematical Biology, 52, pp. 677-696, 1990.
- ^ C.A. Edwards, H.P. Batchelder and T.M. Powell, Modeling microzooplankton and macrozooplankton dynamics within a coastal upwelling system, Journal of Plankton Research, 22, pp. 1619-1648, 2000.
- ^ G.T. Evans and J.S. Parslow, A model of annual plankton cycles, Biological Oceanography, 3, pp. 327-347, 1985.
- ^ K.L. Denman and M.A. Pena, A coupled 1-D biological/physical model of the northeast subarctic Pacific Ocean with iron limitation, Deep Sea Research II, 46, pp. 2877-2908, 1999.