Difference between revisions of "Phosphorus model details"

From JModels
Jump to: navigation, search
(tidy up image)
(tweaks)
Line 1: Line 1:
 
[[Image:pcycle model.jpg|right|thumb|300px|A schematic of the modelled phosphorus cycle.]]
 
[[Image:pcycle model.jpg|right|thumb|300px|A schematic of the modelled phosphorus cycle.]]
  
'''More Detailed Description of the Phosphorus Model.'''
+
This article provides information on '''phosphorus model details''', including the ordinary differential equations that describe the model.
  
 +
==Overview==
 
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.
 
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'').
 
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'').
  
 +
==State variables and equations==
 
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 (''P<sub>S</sub>)'', (2) the concentration of phosphate in the lower (deep) box (''P<sub>D</sub>)'', and (3) the concentration of phytoplankton in the surface box only (''A)''. The ordinary differential equations for the three state variables are:
 
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 (''P<sub>S</sub>)'', (2) the concentration of phosphate in the lower (deep) box (''P<sub>D</sub>)'', and (3) the concentration of phytoplankton in the surface box only (''A)''. The ordinary differential equations for the three state variables are:
  

Revision as of 09:07, 27 March 2008

A schematic of the modelled phosphorus cycle.

This article provides information on phosphorus model details, including the ordinary differential equations that describe the model.

Overview

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

State variables and equations

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:

Peqns.png

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.

Further information


Further Reading: