Phosphorus model details
More Detailed Description of the Phosphorus Model.
The schematic to the right shows how the oceanic cycle of phosphorus is represented in the model. The model is a standard one-dimensional, two-box model of the global ocean, with the top layer representing the surface ocean down to the limit of the deepest wind-induced mixing during the year (down to the annual thermocline) and the bottom layer representing the deep ocean.
There is a simple representation of the physical processes that transfer phosphorus vertically in the ocean. This occurs in reality partly through mixing between the surface and deep oceans, partly through the physical movement of large masses of water between the surface and deep layers (upwelling and downwelling, ocean circulation), and partly due to diffusion. These physical processes are simulated by exchanging a slab of water between the surface and deep boxes. The rate of exchange equates to swapping a slab of thickness 3 metres every year (according to the standard value of the parameter K).
The model is stripped down to the simplest form that can still reproduce the essential aspects of how phosphorus (in the form of phosphate) cycles in the oceans. There are only three state variables in my model: (1) the concentration of phosphate in the upper (surface) box (PS), (2) the concentration of phosphate in the lower (deep) box (PD), and (3) the concentration of phytoplankton in the surface box only (A). The ordinary differential equations for the three state variables are:
Equation 1: the first term calculates uptake by phytoplankton, the second remineralisation back into the surface box, the third term net input of deep phosphate due to mixing, and the fourth term river input.
Equation 2: the first term calculates remineralisation back into the deep box, and the second term calculates net loss of phosphate due to mixing.
Equation 3: the first term is new growth of algal biomass, determined according to phosphate limitation (PH is the Michaelis-Menten half-saturation constant that defines the way in which scarcity of phosphate impacts on growth rate), and the second term is the loss of algal cells due to mortality (M).
Other symbols in the equations above signify (1) the maximum possible growth rate of the phytoplankton (i.e. when phosphate is extremely plentiful) (μmax), (2) the fraction of dead phytoplankton that get decomposed (returned to inorganic nutrients) in the surface box (SR) and the deep box (DR) with the remainder being buried, (3) the amount of phosphorus coming down rivers (RP), and (4) the depths of the surface (SD) and deep (DD) boxes.
The calculation of phytoplankton growth rates omits the effects of light, temperature and other nutrients on growth rate, all of which are important over some parts of the surface ocean at some times. To compensate, the phytoplankton are given lower maximum growth rates than laboratory and field experiments indicate should otherwise be the case.
Parameter values for the phosphorus model are the same as those used in the nitrogen-phosphorus model.