Nitrogen-phosphorus model pros
Strengths of the Nitrogen-Phosphorus Model
In computer modeling there is a continuum between (1) conceptual-type models that aim to simplify a system down to include only its most important elements, in particular those that govern its behaviour, and (2) simulation-type models which are highly detailed and try as much as possible to avoid making assumptions about what is important; instead they try to include as much detail as possible in the expectation that nothing important will be omitted. The nitrogen-phosphorus model (NP model) is very much a conceptual-type model.
Advantages of such an approach include fast run times. The most detailed (fine spatial resolution) simulation models can, at the extreme, take many months of real time to simulate only a few decades. In contrast, this conceptual model can simulate hundreds of thousands of years in only a few seconds.
Leading on from this, a second advantage is the ease of experimentation and sensitivity testing with the quicker models such as this one. If each model run takes many months then the ability to carry out many different runs in series is limited.
Another advantage of simpler models is the ability to understand why the model produces a particular behaviour. This is much harder, and sometimes impossible, with very complex models. It is not always possible to disentangle the large network of interacting effects in a complex model in order to work out which one is key, whereas it is easier in a simple model.
This model of the nitrogen and phosphorus cycles includes the quantitatively largest fluxes in the phosphorus cycle in the same way as the P-only model:
- uptake of P by phytoplankton
- release of P from phytoplankton organic matter upon decay
- downwards transport of P within sinking organic particles
- net upwards transport of P due to physical exchange of (nutrient-scarce) surface and (nutrient-rich) deep waters
- input of P to the ocean as a whole via river input
- loss of P from the ocean as a whole when organic particles sink to the seafloor, eventually get covered by further layers of sinking particles, and finally get compressed into marine sedimentary rocks
It also includes all of the largest fluxes in the nitrogen cycle:
- uptake of fixed nitrogen (nitrate, nitrite, ammonium) by phytoplankton
- uptake of dinitrogen (N2) by nitrogen-fixing phytoplankton
- release of N (initially in the form of ammonium) from phytoplankton organic matter upon decay
- bacterial conversion of ammonium to nitrite and then nitrate (nitrification)
- downwards transport of N as part of sinking organic particles
- net upwards transport of fixed nitrogen due to physical exchange of (nutrient-scarce) surface and (nutrient-rich) deep waters
- input of fixed nitrogen to the ocean as a whole via river input
- input of fixed nitrogen to the ocean as a whole via atmospheric deposition
- loss of fixed nitrogen from the ocean as a whole when organic particles sink to the seafloor, eventually get covered by further layers of sinking particles, and finally get compressed into marine sedimentary rocks
- loss of fixed nitrogen from the ocean as a whole through the processes of denitrification and anammox.
The model tracks the population dynamics of two separate groups of phytoplankton: the N2-fixers and the non-fixers (those which are obliged to use fixed nitrogen as their nitrogen source). The model includes the growth rate dependencies on phosphate and fixed nitrogen availability (the N2-fixers are not affected by availability of fixed nitrogen). By including the competition between these two groups, and its relation to nutrient concentrations, the model is able to include a feedback that links the nitrogen and phosphorus cycles together and provides a potentially stabilising negative feedback loop.
Some models are "closed-system" in the sense of forcing all substances to remain within the model ocean. Sinking material that reaches the seafloor, for instance, will be remineralised at the bottom (or top) of the ocean, and rivers will not be modelled. The boundaries of these models are absolute. This model is more realistic in the sense of being open-boundary (open-system). That is to say, it doesn't just simulate internal fluxes (moving N and P around the ocean) but also the fluxes into and out of the ocean as a whole.
Mass-balance checks are included in the model (behind the scenes) to make sure that the model equations do not lead to any improper loss or gain of P or N from the system as a whole.