Difference between revisions of "Siliconphosphorus model worksheet"
(create page) 
(No difference)

Revision as of 15:27, 8 April 2008
The aim of this siliconphosphorus model worksheet is to analyse some of the properties of the siliconphosphorus model, as well as helping users to become familiar with running the JModels.
Worksheet
Activate the siliconphosphorus model by clicking on this link.
Yool & Tyrrell (2003) state that the siliconphosphorus model: (1) regulates Si concentrations homeostatically; (2) is insensitive to most parameter values and structural/process assumptions; (3) but is dependent on two main assumptions, that μ_{D} > μ_{O} and μ_{O} > M (where μ denotes a maximum growth rate, subscript ‘D’ denotes diatoms, subscript ‘O’ denotes other (nonsiliceous) algae, and M is the constant mortality). Test these statements as follows:
 Try starting the model with: (a) deep P 10fold lower, all else the same; (b) deep Si 10fold lower, all else the same; (c) deep Si 10fold higher, all else the same; and see whether/how the model converges to steady state.
 Change the following parameters by +/ 25% and see if the model still converges to a similar steadystate: R_{org}, the ideal Si:P ratio in diatoms; SR, the surface remineralisation fraction of organic phosphorus; SR_{s}, the surface remineralisation fraction of biogenic silica; SF_{s}, the fraction of biogenic silica that is buried; RS, the riverine input of silicic acid.
 Change the following parameters by +/ 25% and see if the model still converges to a similar steadystate: μ_{D}, diatom maximum growth rate; μ_{O}, other algae maximum growth rate; and M, phytoplankton mortality rate.
 How quickly will the ocean rid itself of excess silicic acid? Test this timescale by initialising the model with 1200 mmol Si m^{3} (this is the threshold above which inorganic silica spontaneously precipitates out of seawater) and examine the model results to see how long it takes for the increased silica burial rate to return the ocean to near steadystate.
 Make a tenfold reduction in the initial concentration of deep silicate. How long does it take the deep silicate concentration to return to equilibrium? (calculate this as the time until 95% of the original perturbation is removed)
 First double and then halve the deep silicate concentration, while leaving all other variables at their default values. Save the data and then examine the burial fluxes of solid silica shortly after the alterations. How do they compare?
Once you have completed these questions, follow this link for the answers.
 Siliconphosphorus model overview
 Siliconphosphorus model details
 Siliconphosphorus model pros
 Siliconphosphorus model cons
 Silica burp hypothesis
References
 Yool, A. & Tyrrell, T. (2003). Role of diatoms in regulating the ocean's silicon cycle. Global Biogeochemical Cycles 17, 1103, doi:10.1029/2002GB002018.