# 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 (the annual thermocline) and the bottom layer representing the deep ocean.

(4) slow mixing between the surface and deep oceans, partly due to 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 between the surface and deep boxes a slab of water of a given thickness (3 metres, according to the standard value of the parameter *K*) every year.

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 ([[Image:]]) in the upper (surface) box (*P _{S})*, (2) the concentration of phosphate ([[Image:]]) in the lower (deep) box (

*P*, and (3) the concentration of phytoplankton in the surface box only (

_{D})*A)*. The equations for the three state variables are:

where 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.

where the first term calculates remineralisation back into the deep box, and the second term calculates net loss of phosphate due to mixing.

here the first term is new growth of algal biomass, determined according to phosphate limitation (*P _{H}* 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.