Introduction
The NPZ (Nutrients, Phytoplankton, Zooplankton) model is based on the model that is developed by Carrasco de la Cruz (2018)1 and is also published by Jiang et al. (2020)2. The model consists of 9 state variables of which 4 in the water column (DIN, Phyto, Zoo, Det) and 5 in the bottom (Mus, Coc, Oys, MPB and Bot_Det). All state variables are expressed in terms of nitrogen (N). The state variables in the water column have the unit of mmol-N m-3 and the state variables in the bottom have a unit of mmol-N m-2. The state variables in the water column can be transported between the compartments and the North Sea by tidal currents. Also they can be introduced into the Oosterschelde through the Krammersluizen. The state variables in the bottom are not exchanged between the compartments.
State variables in the water column
Dissolved inorganic nitrogen (DIN)
Dissolved inorganic nutrients (DIN) is the sum of NO2, NO3 and NH4. The unit is mmol-N m-3. DIN is used by phytoplankton and microphytobenthos for primary production (DINuptake and DINupMPB, respectively). DIN is released by mineralization of detritus in the water column (Mineralization) and in the bottom (BotMin) and respiration by zooplankton (ZooExcretion) and shellfish (DINprod).
graph TD;
DIN-- DINuptake--> Phy;
DIN-- DINupMPB--> MPB;
Det-- Mineralization--> DIN;
Bot_Det-- BotMin--> DIN;
Zoo-- ZooExcretion--> DIN;
Shellfish-- DINprod--> DIN;
Phytoplankton (Phy)
Phytoplankton is composed of microalgae in the water column that are transported by tidal currents. The unit of phytoplankton is mmol-N m-3. Phytoplankton increases by primary production (DINuptake). Phytoplankton is consumed by zooplankton (ZooGrazing) and shellfish (GrazingPhy). Also phytoplankton can sink to the bottom (SinkPhy), where it becomes part of the bottom detritus (Bot_Det).
graph TD;
DIN-- DINuptake--> Phy;
Phy-- ZooGrazing--> Zoo;
Phy-- GrazingPhy--> Shellfish;
Phy-- SinkPhy--> Bot_Det;
Zooplankton (Zoo)
Zooplankton (Zoo) are small animals (e.g. crustaceans, shellfish larvae) that live in the water column and are transported by the tidal currents. The unit is mmol-N m-3. Zooplankton consumes phytoplankton (ZooGrazing). A fraction of the nitrogen that is grazed is lost as feaces (ZooFaeces) and the rest will result in an increase of the zooplankton biomass. Zooplankton can also be consumed by shellfish (GrazingZoo). When the zooplankton dies (ZooMortality) it will become detritus (Det). Finally, zooplankton also produces dissolved inorganic nitrogen by respiration (ZooExcretion).
graph TD;
Phy-- ZooGrazing--> Zoo;
Zoo-- Zoofaeces--> Det;
Zoo-- GrazingZoo--> Shellfish;
Zoo-- ZooMortality--> Det;
Zoo-- ZooExcretion--> DIN;
Detritus (Det)
Detritus (Det) is the dead organic carbon. The unit is mmol-N m-3. Detritus is produced by mortality of zooplankton (ZooMortality) and the production of faeces by zooplankton (ZooFaeces). Detritus is mineralized in the watercolum (Mineralisation). Shellfish can also filter detritus from the water column (GrazingDet). Finally detritus can sink to the seafloor (SinkDet) where it becomes Bot_Det .
graph TD;
Zoo-- ZooMortality--> Det;
Zoo-- Zoofaeces--> Det;
Det-- Mineralization--> DIN;
Det-- GrazingDet--> Shellfish;
Det-- SinkDet--> Bot_Det;
State variables in the bottom
Microphytobenthos (MPB)
Microphytobenthos (MPD) are microalgae that live on sea floor. The unit is mmol-N m-2. MPB only occurs on the intertidal areas. The biomass of MPB increases through primary production (DINuptMPB). When the MPB dies (MortMPB) the nitrogen goes to the bottom detritus (Bot_Det).
graph TD;
DIN-- DINupMPB--> MPB;
MPB-- MortMPB--> Bot_Det;
Detritus in the bottom (Bot_Det)
Detritus in the bottom (Bot_Det, mmol-N m-2) increases by the sinking of detritus and phytoplankton (SinkDet and SinkPhy, respectively). When the microphytobenthos and shellfish (Mus, Oys and Coc) die, the nitrogen also goes to the pool of Bot_Det (MortMPB and MortShellfish, respectively). Finally, the faeces and pseudo faeces that is produced by the shellfish is also a source for the Bot_Det (PseudoFaeces). Nitrogen is released to the water column through mineralization of the Bot_Det (BotMin) where it becomes part of the DIN. A fraction of the nitrogen is lost as N2 due to denitrification (not included in the scheme below).
graph TD;
Det-- SinkDet--> Bot_Det;
Phy-- SinkPhy--> Bot_Det;
MPB-- MortMPB--> Bot_Det;
Shellfish-- MortShellfish--> Bot_Det;
Shellfish-- Faeces--> Bot_Det;
Shellfish-- PseudoFaeces--> Bot_Det;
Bot_Det-- BotMin--> DIN;
Shellfish (Mus, Oys and Coc)
Three shellfish species are included in the model, mussels (Mus), Pacific oysters (Oys) and cockles (Coc). The units are in mmol-N m-2. The mussels are only present at the mussel culture plots in the Oosterschelde, oysters are present at the oyster culture plots and at the wild oyster beds. Cockles are only present at the intertidal flats. The physiology of the shellfish is described by the function “SFG()”. This function calculates clearance rate (Clearance), filtration rate of phytoplankton, zooplankton and detritus (GrazingPhy, GrazingZoo and GrazingDet, respectively), respiration rate (DINprod), the production of faeces and pseudo faeces (PseudoFaeces) and the growth of the shellfish (dFF). The shellfish populations are controlled by a variable mortality (MortShellfish) which is a function of the total stock and the carrying capacity.
graph TD;
Phy-- GrazingPhy--> Shellfish;
Zoo-- GrazingZoo--> Shellfish;
Det-- GrazingDet--> Shellfish;
Shellfish-- DINprod--> DIN;
Shellfish-- MortShellfish--> Bot_Det;
Shellfish-- Faeces--> Bot_Det;
Shellfish-- PseudoFaeces--> Bot_Det;
-
Carrasco de la Cruz, P. M. (2018) An ecosystem model of the Oosterschelde estuary (EMOES), Ghent University. ↩
-
Jiang, L., T. Gerkema, J. C. Kromkamp, D. Van Der Wal, P. M. Carrasco de la Cruz and K. Soetaert (2020) Drivers of the spatial phytoplankton gradient in estuarine–coastal systems: generic implications of a case study in a Dutch tidal bay. Biogeosciences 17: 4135-4152. link ↩