Appendix R Parameters
Following are parameters that can be specified within the WQ Module via user commands.
Tips for use:
- Use the search box to dynamically filter the table
- Sort the filtered (or unfiltered) table by clicking the arrows at the top of each column
- Use “s” on your keyboard to toggle the contents pane on and off to provide more screen width if needed
- UD = user defined
- N/A = not applicable
- For each parameter, use the links in the final column to navigate to either:
- A description of the corresponding control file command that includes the parameter (these links are all suffixed with ‘==’)
- The section describing the relevant science and/or equations (these links are all called ‘science’)
- Where numerical limits are provided for phytoplankton simulation parameters, the conversion bewteen \(\mu\)g/L Clh a and mmol/m\(^3\) or carbon is assumed to be 27.0 mg C/mg Chl a
Where parameter ranges are provided, these are not to be taken as fixed, but as indicative only. Some environmental settings may require parameterisation outside the ranges specified in this manual.
Parameter | Notation | Units | Range | Default | Description | Links |
---|---|---|---|---|---|---|
Sediment oxygen flux | \(F_{sed}^{O_2}\) |
mg O\(_2\)/m\(^2\)/d mmol O\(_2\)/m\(^2\)/d |
(-6400, 0) (-200, 0) |
0 0 |
The rate at which dissolved oxygen is exchanged between the sediments and overlying water column at 20\(^o\)C excluding the influence of overlying dissolved oxygen concentration. Can be varied spatially via use of model mesh material types. This is not the maxuimum rate of oxygen sediment flux |
oxygen flux == science |
Sediment half saturation oxygen concentration for oxygen flux | \(K_{sed-O_2}^{O_2}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(0, 4.8) (0, 150) |
4 125 |
The dissolved oxygen concentration for which the computed sediment oxygen flux is half the user specified rate, at 20\(^o\)C. Cannot be varied spatially |
oxygen benthic == science |
Sediment oxygen flux temperature coefficient | \(\theta_{sed}^{O_2}\) | [-] | (1, 1.1) | 1 | The factor to which (\(T\)-20) is raised to account for temperature effects on oxygen sediment flux. A value of 1.0 will remove the effect of temperature |
oxygen benthic == science |
Minimum dissolved oxygen concentration | \(\left[DO\right]_{min}^{O_2}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(0, 0.5) (0, 15.6) |
NA NA |
The minimum allowable dissolved oxygen concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
oxygen min max == science |
Maximum dissolved oxygen concentration | \(\left[DO\right]_{max}^{O_2}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(6, 20) (187.5, 625) |
NA NA |
The maximum allowable dissolved oxygen concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
oxygen min max == science |
Timestep | \(dt\) | seconds | (UD, UD) | 900 | The timestep at which the WQ Module is called from the hydrodynamic model |
wq dt == |
Sediment silicate flux | \(F_{sed}^{Si}\) |
mg Si/m\(^2\)/d mmol Si/m\(^2\)/d |
(0, 140) (0, 5) |
0 0 |
The rate at which silicate is exchanged between the sediments and overlying water column at 20\(^o\)C excluding the influence of overlying dissolved oxygen concentration. Can be varied spatially via use of model mesh material types. This is not the maxuimum rate of silicate sediment flux |
silicate flux == science |
Sediment half saturation oxygen concentration for silicate flux | \(K_{sed-O_2}^{Si}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(0, 4.8) (0, 150) |
4 125 |
The dissolved oxygen concentration for which the computed sediment silicate flux is half the user specified rate, at 20\(^o\)C. Cannot be varied spatially |
silicate benthic == science |
Sediment silicate flux temperature coefficient | \(\theta_{sed}^{Si}\) | [-] | (1, 1.1) | 1 | The factor to which (\(T\)-20) is raised to account for temperature effects on silicate sediment flux. A value of 1.0 will remove the effect of temperature |
silicate benthic == science |
Minimum silicate concentration | \(\left[Si\right]_{min}^{Si}\) |
mg/L Si mmol/m\(^3\) Si |
(0, 2.8) (0, 100) |
NA NA |
The minimum allowable silicate concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
silicate min max == science |
Maximum silicate concentration | \(\left[Si\right]_{max}^{Si}\) |
mg/L Si mmol/m\(^3\) Si |
(0, 140) (0, 5000) |
NA NA |
The maximum allowable silicate concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
silicate min max == science |
Sediment ammonium flux | \(F_{sed}^{NH_4}\) |
mg NH\(_4\)-N/m\(^2\)/d mmol NH\(_4\)-N/m\(^2\)/d |
(0, 140) (0, 10) |
0 0 |
The rate at which ammonium is exchanged between the sediments and overlying water column at 20\(^o\)C excluding the influence of overlying dissolved oxygen concentration. Can be varied spatially via use of model mesh material types. This is not the maxuimum rate of ammonium sediment flux |
ammonium flux == science |
Sediment half saturation oxygen concentration for ammonium flux | \(K_{sed-O_2}^{NH_4}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(0, 4.8) (0, 150) |
4 125 |
The dissolved oxygen concentration for which the computed sediment ammonium flux is half the user specified rate, at 20\(^o\)C. Cannot be varied spatially |
ammonium benthic == science |
Sediment ammonium flux temperature coefficient | \(\theta_{sed}^{NH_4}\) | [-] | (1, 1.1) | 1 | The factor to which (\(T\)-20) is raised to account for temperature effects on ammonium sediment flux. A value of 1.0 will remove the effect of temperature |
ammonium benthic == science |
Minimum ammonium concentration | \(\left[NH_4\right]_{min}^{NH_4}\) |
mg/L NH\(_4\)-N mmol/m\(^3\) NH\(_4\)-N |
(0, 0.098) (0, 7) |
NA NA |
The minimum allowable ammonium concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
ammonium min max == science |
Maximum ammonium concentration | \(\left[NH_4\right]_{max}^{NH_4}\) |
mg/L NH\(_4\)-N mmol/m\(^3\) NH\(_4\)-N |
(0, 98) (0, 7000) |
NA NA |
The maximum allowable ammonium concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
ammonium min max == science |
Sediment nitrate flux | \(F_{sed}^{NO_3}\) |
mg NO\(_3\)-N/m\(^2\)/d mmol NO\(_3\)-N/m\(^2\)/d |
(0, 140) (0, 10) |
0 0 |
The rate at which nitrate is exchanged between the sediments and overlying water column at 20\(^o\)C excluding the influence of overlying dissolved oxygen concentration. Can be varied spatially via use of model mesh material types. This is not the maxuimum rate of nitrate sediment flux |
nitrate flux == science |
Sediment half saturation oxygen concentration for nitrate flux | \(K_{sed-O_2}^{NO_3}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(0, 4.8) (0, 150) |
4 125 |
The dissolved oxygen concentration for which the computed sediment nitrate flux is half the user specified rate, at 20\(^o\)C. Cannot be varied spatially |
nitrate benthic == science |
Sediment nitrate flux temperature coefficient | \(\theta_{sed}^{NO_3}\) | [-] | (1, 1.1) | 1 | The factor to which (\(T\)-20) is raised to account for temperature effects on nitrate sediment flux. A value of 1.0 will remove the effect of temperature |
nitrate benthic == science |
Minimum nitrate concentration | \(\left[NO_3\right]_{min}^{NO_3}\) |
mg/L NO\(_3\)-N mmol/m\(^3\) NO\(_3\)-N |
(0, 0.098) (0, 7) |
NA NA |
The minimum allowable nitrate concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
nitrate min max == science |
Maximum nitrate concentration | \(\left[NO_3\right]_{max}^{NO_3}\) |
mg/L NO\(_3\)-N mmol/m\(^3\) NO\(_3\)-N |
(0, 98) (0, 7000) |
NA NA |
The maximum allowable nitrate concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
nitrate min max == science |
Nitrification rate | \(R_{nitrif}^{NH_4}\) | /d | (0, 1) |
0 0 |
The rate at which ammonium is nitrified to nitrate at 20\(^o\)C excluding the influence of ambient dissolved oxygen concentration. A value of 0.0 will switch off nitrification |
nitrification == science |
Half saturation oxygen concentration for nitrification | \(K_{nitrif-O_2}^{NH_4}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(0, 4.8) (0, 150) |
4 125 |
The dissolved oxygen concentration for which the computed nitrification rate is half the user specified rate, at 20\(^o\)C |
nitrification == science |
Nitrification temperature coefficient | \(\theta_{nitrif}^{NH_4}\) | [-] | (1, 1.1) | 1 | The factor to which (\(T\)-20) is raised to account for temperature effects on nitrification. A value of 1.0 will remove the effect of temperature |
nitrification == science |
Denitrification rate | \(R_{denit}^{NO_3}\) | /d | (0, 0.5) |
0 0 |
The rate at which nitrate is denitrified to nitrogen gas N\(_2\) at 20\(^o\)C excluding the influence of ambient dissolved oxygen and nitrate concentrations. A value of 0.0 will switch off denitrification |
denitrification == science |
Half saturation oxygen concentration for denitrification | \(K_{denit-O2-MM}^{NO_3}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(0, 4.8) (0, 150) |
4 125 |
The dissolved oxygen concentration for which the computed denitrification rate is half the user specified rate, at 20\(^o\)C and with no influence of ambient nitrate concentration. This value is used if the ‘Michaelis Menten’ denitrification model is speciifed as the first argument of the denitrification == command |
denitrification == science |
Normalising oxygen concentration for denitrification | \(K_{denit-O2-exp}^{NO_3}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(0, 4.8) (0, 150) |
4 125 |
The dissolved oxygen concentration that divides (normalises) the ambient dissolved oxygen concentration in the exponent term of the exponential denitrification model, at 20\(^o\)C and with no influence of ambient nitrate concentration. This value is used if the ‘exponential’ denitrification model is speciifed as the first argument of the denitrification == command |
denitrification == science |
Denitrification temperature coefficient | \(\theta_{denit}^{NO_3}\) | [-] | (1, 1.1) | 1 | The factor to which (\(T\)-20) is raised to account for temperature effects on denitrification. A value of 1.0 will remove the effect of temperature |
denitrification == science |
Half saturation nitrate concentration for denitrification | \(K_{denit-NO3}^{NO_3}\) |
mg/L NO\(_3\)-N mmol/m\(^3\) NO\(_3\)-N |
(0, 5) (0, 357.1) |
0.07 5 |
The nitrate concentration for which the computed denitrification rate is half the user specified rate, at 20\(^o\)C and with no influence of ambient dissolved oxygen concentration. This is not user specifiable and is set to the default value |
denitrification == science |
Denitrification model | Michaelis Menten;exponential | 0 | (, ) | 0 | The model used to simulate denitrification |
denitrification == science |
Wet atmospheric inorganic nitrogen deposition | \(\left[{TN}\right]_{rain}\) |
mg/L N mmol/m\(^3\) N |
(0, 98) (0, 7000) |
0 0 |
The concentration of total inorganic nitrogen in rainfall. This is the sum of ammonium and nitrate N and is applied as a multiplicative factor to rainfall depth to compute a mass flux |
atmospheric deposition == science |
Dry atmospheric inorganic nitrogen deposition | \(R_{atm-dry}^{TN}\) |
mg/m\(^2\)/d N mmol/m\(^2\)/d N |
(0, 98) (0, 7000) |
0 0 |
The rate of atmopsheric fallout of dry inorganic nitrogen. This is the sum of ammonium and nitrate N and is applied as a spatial and temporal constant value to compute a mass flux |
atmospheric deposition == science |
Fraction of atmospheric nitrogen deposition that is nitrate | \(f_{TN}^{NO_3}\) | [-] | (0, 1) | 0 | The fraction of N that is nitrate in atmospheric deposition. The same factor is applied to both wet and dry deposition of N. A value of 1 (0) will set wet and dry deposition to be purely nitrate (ammonium) |
atmospheric deposition == science |
Minimum FRP concentration | \(\left[FRP\right]_{min}^{FRP}\) |
mg/L PO\(_4\)-P mmol/m\(^3\) PO\(_4\)-P |
(0, 0.01) (0, 0.32) |
NA NA |
The minimum allowable FRP concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
FRP min max == science |
Maximum FRP concentration | \(\left[FRP\right]_{max}^{FRP}\) |
mg/L PO\(_4\)-P mmol/m\(^3\) PO\(_4\)-P |
(0, 10) (0, 320) |
NA NA |
The maximum allowable FRP concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
FRP min max == science |
Sediment FRP flux | \(F_{sed}^{FRP}\) |
mg FRP-P/m\(^2\)/d mmol FRP-P/m\(^2\)/d |
(0, 155) (0, 5) |
0 0 |
The rate at which FRP is exchanged between the sediments and overlying water column at 20\(^o\)C excluding the influence of overlying dissolved oxygen concentration. Can be varied spatially via use of model mesh material types. This is not the maxuimum rate of FRP sediment flux |
FRP flux == science |
Sediment half saturation oxygen concentration for FRP flux | \(K_{sed-O_2}^{FRP}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(0, 4.8) (0, 150) |
4 125 |
The dissolved oxygen concentration for which the computed sediment FRP flux is half the user specified rate, at 20\(^o\)C. Cannot be varied spatially |
FRP benthic == science |
Sediment FRP flux temperature coefficient | \(\theta_{sed}^{FRP}\) | [-] | (1, 1.1) | 1 | The factor to which (\(T\)-20) is raised to account for temperature effects on FRP sediment flux. A value of 1.0 will remove the effect of temperature |
FRP benthic == science |
Wet atmospheric FRP deposition | \(\left[{FRP}\right]_{rain}\) |
mg/L P mmol/m\(^3\) P |
(0, 10) (0, 320) |
0 0 |
The concentration of FRP in rainfall. This is only FRP P and is applied as a multiplicative factor to rainfall depth to compute a mass flux |
atmospheric deposition == science |
Dry atmospheric FRP deposition | \(R_{atm-dry}^{FRPads}\) |
mg/m\(^2\)/d P mmol/m\(^2\)/d P |
(0, 10) (0, 320) |
0 0 |
The rate of atmopsheric fallout of adsorbed FRP. This is only applied to adsorbed FRP P. If adsorbed FRP is not simulated then this value is ignored. If used, it is applied as a spatial and temporal constant value to compute a mass flux |
atmospheric deposition == science |
Minimum adsorbed FRP concentration | \(\left[FRPads\right]_{min}^{FRPads}\) |
mg/L PO\(_4\)-P mmol/m\(^3\) PO\(_4\)-P |
(0, 0.01) (0, 0.32) |
NA NA |
The minimum allowable adsorbed FRP concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
FRPads min max == science |
Maximum adsorbed FRP concentration | \(\left[FRPads\right]_{max}^{FRPads}\) |
mg/L PO\(_4\)-P mmol/m\(^3\) PO\(_4\)-P |
(0, 10) (0, 320) |
NA NA |
The maximum allowable adsorbed FRP concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
FRPads min max == science |
Linear FRP sorption partitioning coefficient | \(K_{ads-L}^{FRP}\) | L/mg | (0.01, 0.3) | 0.01 | The coefficient that determines the partitioning of FRP between free and adsorbed states. Can be estimated by dividing adsorbed FRP by free FRP and the concentration of suspended sediment, if these three quantities are known |
adsorption == science |
Quadratic FRP sorption partitioning coefficient | \(K_{ads-Q}^{FRP}\) | L/mg | (0.6, 0.8) | 0.7 | The ratio of adsorption and desorption rate coefficients for phosphorus |
adsorption == science |
Maximum adsorption capacity of phosphorus | \(Q^{FRP}_{max}\) | mg FRP-P/mg TSS | (0.005, 0.006) | 0.0051 | The maximum mass of P that can be adsorbed onto ambient suspended sediment, per unit mass of suspended sediment |
adsorption == science |
Maximum adsorption particle size | \(D50_{ads-max}^{FRP}\) | m | (, ) | 0 | Not yet used. Will be activated in future releases of the WQ Module |
|
Minimum phyto concentration | \(\left[PHY\right]_{min}^{PHY}\) |
\(\mu\)g/L Chl a mmol/m\(^3\) C |
(0.01, 0.1) (0.023, 0.225) |
0.01 0.0225 |
The minimum allowable group phytoplankton concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. It is also important the this minuimum value not be set to zero. If phytoplankton concentrations are allowed to reduce to zero then the numerical scheme will not allow the group to regrow because growth is applied as a multiplicative factor from one timestep to the next. Multiplying any growth rate by a previous timestep’s zero group concentration will therefore always result in zero group oncentration for subsequent timesteps |
min max == science |
Maximum phyto concentration | \(\left[PHY\right]_{max}^{PHY}\) |
\(\mu\)g/L Chl a mmol/m\(^3\) C |
(0.01, 50) (0.0025, 112.5) |
NA NA |
The maximum allowable group phytoplankton concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
min max == science |
Standard temperature | \(T_{std}^{phy}\) | \(^o\)C | (5, 25) | 10 | Used to compute shape of temperature limitation function between T\(_{std}\) and T\(_{max}\). Below T\(_{std}\), the phytoplankton temperature limitation function is a simple Arrhenius relationship. Must always be less than the specified T$_{max} |
temperature limitation == science |
Optimum temperature | \(T_{opt}^{phy}\) | \(^o\)C | (5, 25) | 20 | Used to compute shape of temperature limitation function between T\(_{std}\) and T\(_{max}\). Must always be less than the specified T$_{max} |
temperature limitation == science |
Maximum temperature | \(T_{max}^{phy}\) | \(^o\)C | (25, 35) | 30 | Used to compute shape of temperature limitation function between T\(_{std}\) and T\(_{max}\). Above T\(_{max}\), the phytoplankton temperature limitation function is 0. Must always be greater than the specified T\(_{std} and T\)_{opt}$ |
temperature limitation == science |
Optimum freshwater salinity | \(S_{opt-fresh}^{phy}\) | g/L | (0, 45) | 0 | Used to compute the freshwater salinity limitation function. It is the maximum salinity at which the freshwater salinity limitation function has a value of 1 |
salinity limitation == science |
Maximum freshwater salinity | \(S_{max-fresh}^{phy}\) | g/L | (0, 45) | 5 | Used to compute the freshwater salinity limitation function. It is the maximum salinity at which the freshwater salinity limitation function is greater than the vlaue of \(L_{max-fresh}^{phy}\) |
salinity limitation == science |
Limitation function at maximum freshwater salinity | \(L_{max-fresh}^{phy}\) | [-] | (0, -) | 0 | Used to compute the freshwater salinity limitation function. It is the value of the freshwater limitation function at a salinity of \(S_{max-fresh}\). If this limitation function value is set to be between zero and one (greater than one), then the computed limitation will be applied as a reductive (enhancement) factor to primary productivity (respiration). |
salinity limitation == science |
Optimum marine salinity | \(S_{opt-marine}^{phy}\) | g/L | (0, 45) | 0 | Used to compute the marine salinity limitation function. It is the minimum salinity at which the marine salinity limitation function has a value of 1 |
salinity limitation == science |
Limitation function at zero marine salinity | \(L_{zero-marine}^{phy}\) | [-] | (0, -) | 0 | Used to compute the marine salinity limitation function. It is the value of the marine limitation function at a salinity of zero. If this limitation function value is set to be between zero and one (greater than one), then the computed limitation will be applied as a reductive (enhancement) factor to primary productivity (respiration). |
salinity limitation == science |
Optimum mixed salinity | \(S_{opt-mix}^{phy}\) | g/L | (0, 45) | 0 | Used to compute the mixed salinity limitation function. It is the minimum salinity at which the mixed salinity limitation function has a value of 1 |
salinity limitation == science |
Maximum mixed salinity | \(S_{max-mix}^{phy}\) | g/L | (0, 45) | 5 | Used to compute the mixed salinity limitation function. It is the maximum salinity at which the mixed salinity limitation function has a value of 1 |
salinity limitation == science |
Limitation function at zero and (\(S_{opt-mix}^{phy} + S_{max-mix}^{phy}\)) mixed salinity | \(L_{zero-mix}^{phy}\) | [-] | (0, -) | 0 | Used to compute the mixed salinity limitation function. It is the value of the mixed limitation function at a salinity of zero and (\(S_{opt-mix}^{phy} + S_{max-mix}^{phy}\)). If this limitation function value is set to be between zero and one (greater than one), then the computed limitation will be applied as a reductive (enhancement) factor to primary productivity (respiration). |
salinity limitation == science |
Optimum estuarine salinity | \(S_{opt-est}^{phy}\) | g/L | (0, 45) | 0 | Used to compute the estuarine salinity limitation function. It is the salinity at which the estuarine salinity limitation function has a value of 1 |
salinity limitation == science |
Maximum estuarine salinity | \(S_{max-est}^{phy}\) | g/L | (0, 45) | 5 | Used to compute the estuarine salinity limitation function. It is the maximum salinity at which the estuarine salinity limitation function has a value greater than zero |
salinity limitation == science |
Power coefficient for estuarine salinity limitation | \(P_{est}^{phy}\) | [-] | (0, 5) | 0 | Used to compute the estuarine salinity limitation function. It determines the gradient of the estuarine limitation function (with salinity) above and below \(S_{opt-est}^{phy}\). Small (large) values of this parameter provide low (high) gradients |
salinity limitation == science |
Minimum concentration of N for phytoplankton uptake | \(\left[N\right]_{min}^{phy}\) |
mg/L N mmol/m\(^3\) N |
(0, 0.1) (0.00, 7.14) |
0.01 0.71 |
Used to compute the nitrogen limitation function in the basic phytoplankton model. Also used by the advanced phytoplankton model if internal nitrogen stores are exhausted. It is the minimum N concentration (as NH\(_4\)-N + NO\(_3\)-N) at which uptake of nitrogen by phytoplankton can occur |
nitrogen limitation == science |
Half saturation N concentration for phytoplankton uptake of N | \(K_{lim-N}^{phy}\) |
mg/L N mmol/m\(^3\) N |
(1, 5) (71.4, 357.14) |
1 71.43 |
Used to compute the nitrogen limitation function in the basic phytoplankton model. Also used by the advanced phytoplankton model if internal nitrogen stores are exhausted. It is the concentration of N (NH\(_4\)-N + NO\(_3\)-N) that when added to the specified minimum \(\left[N\right]_{min}^{phy}\) sets the nitrogen limitation function to a value of 0.5. A value of zero for this parameter will set the nitrogen limitation function to always be 1 |
nitrogen limitation == science |
Minimum ratio of internal phytoplankton N to biomass | \(X_{N-C-min}^{phy}\) |
mg N/mg Chl a mmol N/mmol C |
(2.97, 4.75) (0.09, 0.15) |
2.97 0.09 |
Used to compute the nitrogen limitation function in the advanced phytoplankton model if internal nitrogen stores are not exhausted. It should be based on the Redfield ratio of 106:16:1 or C:N:P (or Chla a equivalent). Note the units of mg of Chl a rather than \(\mu\)g |
nitrogen limitation == science |
Maximum ratio of internal phytoplankton N to biomass | \(X_{N-C-max}^{phy}\) |
mg N/mg Chl a mmol N/mmol C |
(4.75, 6.54) (0.15, 0.21) |
6.54 0.21 |
Used to compute the nitrogen limitation function in the advanced phytoplankton model if internal nitrogen stores are not exhausted. It should be based on the Redfield ratio of 106:16:1 or C:N:P (or Chla a equivalent). Note the units of mg of Chl a rather than \(\mu\)g |
nitrogen limitation == science |
Minimum concentration of P for phytoplankton uptake | \(\left[P\right]_{min}^{phy}\) |
mg/L P mmol/m\(^3\) P |
(0, 0.1) (0.0, 3.23) |
0.001 0.03 |
Used to compute the phosphorus limitation function in the basic phytoplankton model. Also used by the advanced phytoplankton model if internal phosphorus stores are exhausted. It is the minimum P concentration (as FRP-P) at which uptake of phosphorus by phytoplankton can occur |
phosphorus limitation == science |
Half saturation P concentration for phytoplankton uptake of P | \(K_{lim-P}^{phy}\) |
mg/L P mmol/m\(^3\) P |
(0.1, 0.5) (3.23, 16.13) |
0.1 3.23 |
Used to compute the phosphorus limitation function in the basic phytoplankton model. Also used by the advanced phytoplankton model if internal phosphorus stores are exhausted. It is the concentration of P (FRP-P) that when added to the specified minimum \(\left[P\right]_{min}^{phy}\) sets the phosphorus limitation function to a value of 0.5. A value of zero for this parameter will set the phosphorus limitation function to always be 1 |
phosphorus limitation == science |
Minimum ratio of internal phytoplankton P to biomass | \(X_{P-C-min}^{phy}\) |
mg P/mg Chl a mmol P/mmol C |
(0.33, 0.66) (0.005, 0.009) |
0.33 0.005 |
Used to compute the phosphorus limitation function in the advanced phytoplankton model if internal phosphorus stores are not exhausted. It should be based on the Redfield ratio of 106:16:1 or C:N:P (or Chla a equivalent). Note the units of mg of Chl a rather than \(\mu\)g |
phosphorus limitation == science |
Maximum ratio of internal phytoplankton P to biomass | \(X_{P-C-max}^{phy}\) |
mg P/mg Chl a mmol P/mmol C |
(0.66, 0.99) (0.009, 0.014) |
0.99 0.014 |
Used to compute the phosphorus limitation function in the advanced phytoplankton model if internal phosphorus stores are not exhausted. It should be based on the Redfield ratio of 106:16:1 or C:N:P (or Chla a equivalent). Note the units of mg of Chl a rather than \(\mu\)g |
phosphorus limitation == science |
Minimum concentration of silicate for phytoplankton uptake | \(\left[Si\right]_{min}^{phy}\) |
mg/L Si mmol/m\(^3\) Si |
(0, 0.1) (0.00, 3.57) |
0.01 0.36 |
Used to compute the silicate limitation function in the basic and advanced phytoplankton model. It is the minimum silicate concentration at which uptake of silicate by phytoplankton can occur |
silicate limitation == science |
Half saturation silicate concentration for phytoplankton uptake of silicate | \(K_{lim-Si}^{phy}\) |
mg/L Si mmol/m\(^3\) Si |
(1, 10) (35.71, 357.14) |
1 35.71 |
Used to compute the silicate limitation function in the basic and advanced phytoplankton model. It is the concentration of silicate that when added to the specified minimum \(\left[Si\right]_{min}^{phy}\) sets the silicate limitation function to a value of 0.5. A value of zero for this parameter will set the silicate limitation function to always be 1 |
silicate limitation == science |
Half saturation constant for basic model light growth limitation | \(I_{K-bas}\) | W/m\(^2\) | (1, 1500) | 20 | Used to compute the basic light limitation function for both basic and advanced phytoplankton constituent models. It divides (normalises) the photosynthetically available radiation incident on the top and bottom of a computational cell and forms the lower limit of the exponential integral |
light limitation == science |
Half saturation constant for monod model light growth limitation | \(I_{K-mon}\) | W/m\(^2\) | (1, 1500) | 20 | Used to compute the monod light limitation function for both basic and advanced phytoplankton constituent models. It divides (normalises) the photosynthetically available radiation incident at the centre of a computational cell |
light limitation == science |
Saturating light intensity for the steele model light growth limitation | \(I_{S-ste}\) | W/m\(^2\) | (1, 1500) | 20 | Used to compute the steele light limitation function for both basic and advanced phytoplankton constituent models. It divides (normalises) the photosynthetically available radiation incident at the centre of a computational cell. This limitation function includes photoinhibition |
light limitation == science |
Half saturation constant for webb model light growth limitation | \(I_{K-web}\) | W/m\(^2\) | (1, 1500) | 20 | Used to compute the webb light limitation function for both basic and advanced phytoplankton constituent models. It divides (normalises) the photosynthetically available radiation incident at the centre of a computational cell |
light limitation == science |
Half saturation constant for jassby model light growth limitation | \(I_{K-jas}\) | W/m\(^2\) | (1, 1500) | 20 | Used to compute the jassby light limitation function for both basic and advanced phytoplankton constituent models. It divides (normalises) the photosynthetically available radiation incident at the centre of a computational cell |
light limitation == science |
Half saturation constant for chalker model light growth limitation | \(I_{K-cha}\) | W/m\(^2\) | (1, 1500) | 20 | Used to compute the chalker light limitation function for both basic and advanced phytoplankton constituent models. It divides (normalises) the photosynthetically available radiation incident at the centre of a computational cell |
light limitation == science |
Saturating light intensity for the klepper model light growth limitation | \(I_{S-kle}\) | W/m\(^2\) | (1, 1500) | 20 | Used to compute the klepper light limitation function for both basic and advanced phytoplankton constituent models. It divides (normalises) the photosynthetically available radiation incident at the centre of a computational cell. This limitation function includes photoinhibition |
light limitation == science |
Saturating light intensity for the integrated model light growth limitation | \(I_{S-int}\) | W/m\(^2\) | (1, 1500) | 20 | Used to compute the integrated light limitation function for both basic and advanced phytoplankton constituent models. It divides (normalises) the photosynthetically available radiation incident on the top and bottom of a computational cell. This limitation function includes photoinhibition |
light limitation == science |
Constant ratio of internal phytoplankton N to biomass | \(X_{N-C-con}^{phy}\) |
mg N/mg Chl a mmol N/mmol C |
(2.97, 6.54) (0.09, 0.21) |
4.75 0.15 |
Used to compute the uptake of nitrogen during growth in the basic phytoplankton constituent model only. It should be based on the Redfield ratio of 106:16:1 or C:N:P (or Chla a equivalent). Note the units of mg of Chl a rather than \(\mu\)g |
uptake == science |
Constant ratio of internal phytoplankton P to biomass | \(X_{P-C-con}^{phy}\) |
mg P/mg Chl a mmol P/mmol C |
(0.33, 0.99) (0.009, 0.014) |
0.66 0.009 |
Used to compute the uptake of phosphorus during growth in the basic phytoplankton constituent model only. It should be based on the Redfield ratio of 106:16:1 or C:N:P (or Chla a equivalent). Note the units of mg of Chl a rather than \(\mu\)g |
uptake == science |
Constant ratio of internal phytoplankton silicate to biomass | \(X_{Si-C-con}^{phy}\) |
mg Si/mg Chl a mmol Si/mmol C |
(5.94, 11.89) (0.09, 0.19) |
8.92 0.14 |
Used to compute the uptake of silicate during growth in the basic and advanced phytoplankton constituent models. It should be based on the extended Redfield ratio of 106:16:1:15 or C:N:P:Si (or Chla a equivalent). Note the units of mg of Chl a rather than \(\mu\)g |
uptake == science |
Rate of inorganic nitrogen uptake by phytoplankton | \(R_{N-uptake}^{phy}\) |
mg N/mg Chl a/d mmol N/mmol C/d |
(0.080, 2.704) (0.0025, 0.086) |
0 0 |
Used to compute the uptake of nitrogen during growth in the advanced phytoplankton constituent model only. It is the nitrogen uptake rate per unit of phytoplankton (per day) without any influence of temperature or nutrient limitation. The unit of phytoplankton is either mg (not \(\mu\)g) Chl a or mmol C, depending on the units system used |
uptake == science |
Rate of inorganic phosphorus uptake by phytoplankton | \(R_{P-uptake}^{phy}\) |
mg P/mg Chl a/d mmol P/mmol C/d |
(0.011, 0.374) (0.0002, 0.005) |
0 0 |
Used to compute the uptake of phosphorus during growth in the advanced phytoplankton constituent model only. It is the phosphorus uptake rate per unit of phytoplankton (per day) without any influence of temperature or nutrient limitation. The unit of phytoplankton is either mg (not \(\mu\)g) Chl a or mmol C, depending on the units system used |
uptake == science |
Constant rate of phytoplankton settling | \(V_{settle}^{phy}\) | m/d | (-50, 0) | 0 | The constant settling velocity applied to phytoplankton using either the basic or advanced phytoplankton constituent model. It is applied as a constant (in the constant settling model) or modified by density (using the density correction model). A negative value is a downwards velocity |
settling == science |
Stokes settling light half saturation light intensity | \(I_{K-sto}\) | W/m\(^2\) | (0, 1500) | 20 | The half saturation light intensity used to compute cell density, which is in turn required for calculation of Stokes settling velocity. It can be applied to either the basic or advanced phytoplankton constituent model. This value is ignored if \(I_K\) has been specified in a light limitation model within the same phytoplankton constituent model block |
settling == science |
Motile settling velocity | \(V_{mot}^{phy}\) | m/d | (-50, 0) | 0 | The motile settling velocity applied to phytoplankton, that is then modified by local light conditions. It can be applied to only to the advanced phytoplankton constituent model. A negative value is a downwards velocity |
settling == science |
Motile settling light half saturation light intensity | \(I_{K-mot}\) | W/m\(^2\) | (0, 1500) | 20 | The half saturation light intensity used to compute direction of motility. It can be applied only to the advanced phytoplankton constituent model. This value is ignored if \(I_K\) has been specified in a light limitation model within the same phytoplankton constituent model block |
settling == science |
Carbon to chlorophyll a ratio | \(X_{cc}^{phy}\) | mg carbon / mg chlorophyll a | (20, 40) | 27 | The ratio of carbon to chlorophyll a in the cellular biomass of a phytoplankton group. The units of chlorophyll a in this ratio are milligrams (mg), not \(\mu\)g. This value is constant for a given phytoplankton group |
carbon chla ratio == science |
Nitrogen fixing rate | \(R_{nfix}^{N_2}\) |
mg N / mg Chl a / d mmol N / mmol C /d |
(0, 5) (0, 0.16) |
0 0 |
The rate at which nitrogen N is fixed from the atmosphere to complement nitrogen uptake from the water column. This rate excludes the influence of nitrogen limitation due to water column nutrient scarcity. It is the mass (or molar) rate of uptake of nitrogen N from the atmopshere, not N\(_2\). Note the units of chlorophyll a are in mg, not \(\mu\)g |
nitrogen fixing == science |
Penalty function of nitrogen fixing on primary productivity | \(f_{nfix}^{phy}\) | [-] | (0, 1) | 1 | A rate modifier that accounts for the metabolic cost of nitrogen fixing on primary productivity. This modifier can be thought of as a penalty function, with values of zero and one corresponding to no and maximum penalties, respectively. A penalty of 1 (the default) will always set primary productivity to zero |
nitrogen fixing == science |
Phytoplankton primary productivity rate | \(R_{prod}^{phy}\) | /d | (0, 3.5) | 0 | The phytoplankton primary productivity (growth) rate at 20\(^o\)C with no light, temperature, salinty, silicate or nutrient limitation. This rate is multiplied by ambient phytoplankton concentration to derive initial carbon, nutrient and silicate uptake fluxes |
primary productivity == science |
Phytoplankton primary productivity temperature coefficient | \(\theta_{prod}^{phy}\) | [-] | (1, 1.1) | 1 | The factor to which (\(T\)-20) is raised to account for temperature effects on phytoplankton growth in the standard temperature limitation function. A value of 1.0 will remove the effect of temperature |
primary productivity == science |
Phytoplankton respiration rate | \(R_{resp}^{phy}\) | /d | (0, 0.1) | 0 | The phytoplankton respiration rate at 20\(^o\)C with no temperature or salinty limitation. This rate is multiplied by ambient phytoplankton concentration to derive excretion and mortality fluxes. It is added to the parallel primary productivity and exudation rates and then multiplied by ambient phytoplankton concentration to compute biomass flux |
respiration == science |
Phytoplankton respiration temperature coefficient | \(\theta_{resp}^{phy}\) | [-] | (1, 1.1) | 1 | The factor to which (\(T\)-20) is raised to account for temperature effects on phytoplankton respiration. A value of 1.0 will remove the effect of temperature |
respiration == science |
Fraction of phytoplankton respiration that is true respiration | \(f_{true-resp}^{phy}\) | [-] | (0, 1) | 1 | The proportion of consumed carbohydrates (that were generated during photosynthesis) that are converted to energy, as opposed to waste matter. A value of one (zero) has all consumed carbohydrates converted to energy (waste), and will minimise (maximise) excretive losses |
respiration == science |
Fraction of phytoplanton loss that is excretion | \(f_{excr-loss}^{phy}\) | [-] | (0, 1) | 0 | The proportion of combined excretive and mortality losses that are excretive. If \(f_{true-resp}^{phy}\) is set to one (zero) then this parameter has no (maximum) effect. |
respiration == science |
Fraction of primary production lost to exudation | \(f_{exud}^{phy}\) | [-] | (0, 1) | 0 | The proportion of primary production that does not generate retained cellular carbon biomass, but is lost through exudation (and hence excretion). Exudation losses are summed with potentially modified respiration losses to compute excretion |
respiration == science |
Phytoplankton specific shading coefficient | \(Ke^{phy}\) |
/(\(\mu\)g Chl a/L)/m /(mmol C/m\(^3\))/m |
(0, 0.0225) (0, 0.01) |
0 0 |
The rate at which light is attenuated by phytoplankton. This coefficient is per unit concentration of phytoplankton, per metre (depth). It is multiplied by phytoplankton concentration in each computational cell to leave a value that has units of per meter, which is the same as standard extinction coefficients and therefore is able to be added to other (e.g. background) parallel extinction coefficients. This parameter is used in all light limitation == commands, but only the reference to the basic model is included in this table. Use the navigation in the commands Appendix to see linkages with other light limitation commands |
light limitation == science |
Minimum POC concentration | \(\left[POC\right]_{min}^{POC}\) |
mg/L C mmol/m\(^3\) C |
(0, 0.1) (0, 8.33) |
NA NA |
The minimum allowable labile particulate organic carbon concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
carbon min max == science |
Maximum POC concentration | \(\left[POC\right]_{max}^{POC}\) |
mg/L C mmol/m\(^3\) C |
(0, 100) (0, 8333) |
NA NA |
The maximum allowable labile particulate organic carbon concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
carbon min max == science |
Minimum DOC concentration | \(\left[DOC\right]_{min}^{DOC}\) |
mg/L C mmol/m\(^3\) C |
(0, 0.1) (0, 8.33) |
NA NA |
The minimum allowable labile dissolved organic carbon concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
carbon min max == science |
Maximum DOC concentration | \(\left[DOC\right]_{max}^{DOC}\) |
mg/L C mmol/m\(^3\) C |
(0, 100) (0, 8333) |
NA NA |
The maximum allowable labile dissolved organic carbon concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
carbon min max == science |
Minimum PON concentration | \(\left[PON\right]_{min}^{PON}\) |
mg/L N mmol/m\(^3\) N |
(0, 0.1) (0, 7.14) |
NA NA |
The minimum allowable labile particulate organic nitrogen concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
nitrogen min max == science |
Maximum PON concentration | \(\left[PON\right]_{max}^{PON}\) |
mg/L N mmol/m\(^3\) N |
(0, 100) (0, 7143) |
NA NA |
The maximum allowable labile particulate organic nitrogen concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
nitrogen min max == science |
Minimum DON concentration | \(\left[DON\right]_{min}^{DON}\) |
mg/L N mmol/m\(^3\) N |
(0, 0.1) (0, 7.14) |
NA NA |
The minimum allowable labile dissolved organic nitrogen concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
nitrogen min max == science |
Maximum DON concentration | \(\left[DON\right]_{max}^{DON}\) |
mg/L N mmol/m\(^3\) N |
(0, 100) (0, 7143) |
NA NA |
The maximum allowable labile dissolved organic nitrogen concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
nitrogen min max == science |
Minimum POP concentration | \(\left[POP\right]_{min}^{POP}\) |
mg/L P mmol/m\(^3\) P |
(0, 0.01) (0, 0.32) |
NA NA |
The minimum allowable labile particulate organic phosphorus concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
phosphorus min max == science |
Maximum POP concentration | \(\left[POP\right]_{max}^{POP}\) |
mg/L P mmol/m\(^3\) P |
(0, 10) (0, 323) |
NA NA |
The maximum allowable labile particulate organic phosphorus concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
phosphorus min max == science |
Minimum DOP concentration | \(\left[DOP\right]_{min}^{DOP}\) |
mg/L P mmol/m\(^3\) P |
(0, 0.01) (0, 0.32) |
NA NA |
The minimum allowable labile dissolved organic phosphorus concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
phosphorus min max == science |
Maximum DOP concentration | \(\left[DOP\right]_{max}^{DOP}\) |
mg/L P mmol/m\(^3\) P |
(0, 10) (0, 323) |
NA NA |
The maximum allowable labile dissolved organic phosphorus concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
phosphorus min max == science |
Sediment DOC flux | \(F_{sed}^{DOC}\) |
mg C/m\(^2\)/d mmol C/m\(^2\)/d |
(0, 120) (0, 10) |
0 0 |
The rate at which labile dissolved organic carbon is exchanged between the sediments and overlying water column at 20\(^o\)C excluding the influence of overlying dissolved oxygen concentration. Can be varied spatially via use of model mesh material types. This is not the maxuimum rate of labile DOC sediment flux |
DOC flux == science |
Sediment DON flux | \(F_{sed}^{DON}\) |
mg N/m\(^2\)/d mmol N/m\(^2\)/d |
(0, 21.126) (0, 1.509) |
0 0 |
The rate at which labile dissolved organic nitrogen is exchanged between the sediments and overlying water column at 20\(^o\)C excluding the influence of overlying dissolved oxygen concentration. Can be varied spatially via use of model mesh material types. This is not the maxuimum rate of labile DON sediment flux |
DON flux == science |
Sediment DOP flux | \(F_{sed}^{DOP}\) |
mg P/m\(^2\)/d mmol P/m\(^2\)/d |
(0, 2.914) (0, 0.094) |
0 0 |
The rate at which labile dissolved organic phosphorus is exchanged between the sediments and overlying water column at 20\(^o\)C excluding the influence of overlying dissolved oxygen concentration. Can be varied spatially via use of model mesh material types. This is not the maxuimum rate of labile DOP sediment flux |
DOP flux == science |
Sediment half saturation oxygen concentration for DOM flux | \(K_{sed-O_2}^{DOM}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(0, 4.8) (0, 150) |
4 125 |
The dissolved oxygen concentration for which the computed sediment labile DOM flux is half the user specified rate, at 20\(^o\)C. Cannot be varied spatially. This value is applied equally to labile dissolved organic carbon, nitrogen and phosphorus |
organics benthic == science |
Sediment DOM flux temperature coefficient | \(\theta_{sed}^{DOM}\) | [-] | (1, 1.1) | 1 | The factor to which (\(T\)-20) is raised to account for temperature effects on labile DOM sediment flux. A value of 1.0 will remove the effect of temperature. This value is applied equally to labile dissolved organic carbon, nitrogen and phosphorus |
organics benthic == science |
POC hydrolysis rate | \(R_{hyd}^{POC}\) | /d | (0, 0.1) |
0 0 |
The rate at which labile particulate organic carbon is hydrolysed to dissolved organic carbon at 20\(^o\)C excluding the influence of ambient dissolved oxygen concentration. This is not the maxuimum rate of POC hydrolysis |
hydrolysis == science |
PON hydrolysis rate | \(R_{hyd}^{PON}\) | /d | (0, 0.015) |
0 0 |
The rate at which labile particulate organic nitrogen is hydrolysed to dissolved organic nitrogen at 20\(^o\)C excluding the influence of ambient dissolved oxygen concentration. This is not the maxuimum rate of PON hydrolysis |
hydrolysis == science |
POP hydrolysis rate | \(R_{hyd}^{POP}\) | /d | (0, 0.0009) |
0 0 |
The rate at which labile particulate organic phosphorus is hydrolysed to dissolved organic phosphorus at 20\(^o\)C excluding the influence of ambient dissolved oxygen concentration. This is not the maxuimum rate of POP hydrolysis |
hydrolysis == science |
Half saturation oxygen concentration for POM hydrolysis | \(K_{hyd-O_2}^{POM}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(0, 4.8) (0, 150) |
4 125 |
The dissolved oxygen concentration for which the computed hydrolysis of labile POM is half the user specified rate, at 20\(^o\)C. Cannot be varied spatially. The same value is applied to all labile particulate organic constituents (carbon, nitrogen and phosphorus), and also to the breakdown of refractory particulate organic matter |
hydrolysis == science |
POM hydrolysis temperature coefficient | \(\theta_{hyd}^{POM}\) | [-] | (1, 1.1) | 1 | The factor to which (\(T\)-20) is raised to account for temperature effects on labile POM hydrolysis. A value of 1.0 will remove the effect of temperature. The same value is applied to all labile particulate organic constituents (carbon, nitrogen and phosphorus), and also to the breakdown of refractory particulate organic matter |
hydrolysis == science |
DOM mineralisation rate | \(R_{miner}^{DOM}\) | /d | (0, 0.1) |
0 0 |
The rate at which labile dissolved organic carbon, nitrogen and phosphorus are mineralised to their corresponding dissolved inorganic constituents at 20\(^o\)C excluding the influence of ambient dissolved oxygen concentration. This is not the maxuimum rate of mineralisation. The same rate is applied to all labile dissolved organic constituents (carbon, nitrogen and phosphorus) |
mineralisation == science |
Half saturation oxygen concentration for DOM mineralisation | \(K_{miner-O_2}^{DOM}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(0, 4.8) (0, 150) |
4 125 |
The dissolved oxygen concentration for which the computed aerobic and anaerobic mineralisation of labile DOM is half the user specified rate, at 20\(^o\)C, and with \(f_{an}\) equal to 1.0. If \(f_{an}\) is equal to zero then this half saturation oxygen concentration is the dissolved oxygen concentration for which the computed aerobic mineralisation of labile DOM is half the user specified rate, again at 20\(^o\)C. Cannot be varied spatially. The same value is applied to all labile dissolved organic constituents (carbon, nitrogen and phosphorus), and also to the activation of refractory dissolved organic carbon, nitrogen and phosporus |
mineralisation == science |
DOM mineralisation temperature coefficient | \(\theta_{miner}^{DOM}\) | [-] | (1, 1.1) | 1 | The factor to which (\(T\)-20) is raised to account for temperature effects on labile DOM mineralisation. A value of 1.0 will remove the effect of temperature. The same value is applied to all labile dissolved organic constituents (carbon, nitrogen and phosphorus), and also to the activation of refractory dissolved organic carbon, nitrogen and phosporus |
mineralisation == science |
Anaerobic mineralisation multiplier | \(f_{an}\) | [-] | (0, 1) | 0 |
The factor applied to the compenent of mineralisation of labile dissolved organic matter that does not consume dissolved oxygen. A value of 0.0 will entirely suppress the Michaelis-Menten weighting that non dissolved oxygen based processes have in calculating overall mineralisation of labile dissolved organic matter. Setting this parameter to zero does not suppress non dissolved oxygen mineralisation. This same factor is also applied to the activation of refractory dissolved organic carbon, nitrogen and phosphorus to their labile equivalents. If activation is to be simulated then the command mineralisation == will need to be issued to set \(f_{an}\). In this case, and if the simulation of dissolved organic matter mineralisation is not desired, then set \(R_{miner}^{DOM}\) to be zero in the same mineralisation command
|
mineralisation == science |
Half saturation nitrate concentration for nitrate consumption via anaerobic mineralisation of dissolved organic carbon | \(K_{miner-NO_3}^{NO_3}\) |
mg/L NO\(_3\)-N mmol/m\(^3\) NO\(_3\)-N |
(0, 2.1) (0, 150) |
1.75 125 |
The nitrate concentration for which the computed consumption of nitrate due to anaerobic mineralisation is half the corresponding anaerobic mineralisation of labile dissolved organic carbon. If this paramter is set to zero, then the consumption of nitrate as a source of oxygen as part of organic matter mineralisation is switched off. Cannot be varied spatially |
mineralisation == science |
POM specific shading coefficient | \(Ke^{POM}\) |
/(mg POC/L)/m /(mmol POC/m\(^3\))/m |
(0, 8.33) (0, 0.1) |
0 0 |
The rate at which light is attenuated by labile particulate organic matter. This coefficient is per unit concentration of particulate organic carbon (which is used as a proxy for particulate organic matter) per metre (depth). It is multiplied by POC concentration in each computational cell to leave a value that has units of per meter, which is the same as standard extinction coefficients and therefore is able to be added to other (e.g. background) parallel extinction coefficients |
self shading == science |
DOM specific shading coefficient | \(Ke^{DOM}\) |
/(mg DOC/L)/m /(mmol DOC/m\(^3\))/m |
(0, 8.33) (0, 0.1) |
0 0 |
The rate at which light is attenuated by labile dissolved organic matter. This coefficient is per unit concentration of dissolved organic carbon (which is used as a proxy for dissolved organic matter) per metre (depth). It is multiplied by DOC concentration in each computational cell to leave a value that has units of per meter, which is the same as standard extinction coefficients and therefore is able to be added to other (e.g. background) parallel extinction coefficients |
self shading == science |
RPOM breakdown rate | \(R_{bdn}^{RPOM}\) | /d | (0, 0.1) |
0 0 |
The rate at which refractory particulate organic matter (a single computed variable) is broken down to labile particulate organic carbon, nitrogen and phosphorus (three computed variables) at 20\(^o\)C excluding the influence of ambient dissolved oxygen concentration. This is not the maxuimum rate of breakdown |
ref breakdown == science |
Ratio of N to C in refractory particulate organic matter | \(X_N^{RPOM}\) | [-] | (0.094, 0.19) | 0.15 | The molar ratio of N to C in refractory particulate organic matter. This ratio should not differ significantly from the Redfield ratio of 16/106 |
ref breakdown == science |
Ratio of P to C in refractory particulate organic matter | \(X_P^{RPOM}\) | [-] | (0.00094, 0.019) | 0.0094 | The molar ratio of P to C in refractory particulate organic matter. This ratio should not differ significantly from the Redfield ratio of 1/106 |
ref breakdown == science |
Minimum RPOM concentration | \(\left[RPOM\right]_{min}^{RPOM}\) |
mg/L C mmol/m\(^3\) C |
(0, 0.1) (0, 8.33) |
NA NA |
The minimum allowable refractory particulate organic matter concentration. Although referred to as organic matter, it is accounted as (and in units of) carbon in the WQ Module. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
ref carbon min max == science |
Maximum RPOM concentration | \(\left[RPOM\right]_{max}^{RPOM}\) |
mg/L C mmol/m\(^3\) C |
(0, 100) (0, 8333) |
NA NA |
The maximum allowable refractory particulate organic matter concentration. Although referred to as organic matter, it is accounted as (and in units of) carbon in the WQ Module. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
ref carbon min max == science |
Minimum RDOC concentration | \(\left[RDOC\right]_{min}^{RDOC}\) |
mg/L C mmol/m\(^3\) C |
(0, 0.1) (0, 8.33) |
NA NA |
The minimum allowable refractory dissolved organic carbon concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
ref carbon min max == science |
Maximum RDOC concentration | \(\left[RDOC\right]_{max}^{RDOC}\) |
mg/L C mmol/m\(^3\) C |
(0, 100) (0, 8333) |
NA NA |
The maximum allowable refractory dissolved organic carbon concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
ref carbon min max == science |
Minimum RDON concentration | \(\left[RDON\right]_{min}^{RDON}\) |
mg/L N mmol/m\(^3\) N |
(0, 0.1) (0, 7.14) |
NA NA |
The minimum allowable refractory dissolved organic nitrogen concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
ref nitrogen min max == science |
Maximum RDON concentration | \(\left[RDON\right]_{max}^{RDON}\) |
mg/L N mmol/m\(^3\) N |
(0, 100) (0, 7143) |
NA NA |
The maximum allowable refractory dissolved organic nitrogen concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
ref nitrogen min max == science |
Minimum RDOP concentration | \(\left[RDOP\right]_{min}^{RDOP}\) |
mg/L P mmol/m\(^3\) P |
(0, 0.01) (0, 0.32) |
NA NA |
The minimum allowable refractory dissolved organic phosphorus concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
ref phosphorus min max == science |
Maximum RDOP concentration | \(\left[RDOP\right]_{max}^{RDOP}\) |
mg/L P mmol/m\(^3\) P |
(0, 10) (0, 323) |
NA NA |
The maximum allowable refractory dissolved organic phosphorus concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
ref phosphorus min max == science |
RDOM activation rate | \(R_{act}^{RDOM}\) | /d | (0, 0.1) |
0 0 |
The rate at which refractory dissolved organic carbon, nitrogen and phosphorus are activated to their corresponding labile dissolved organic constituents at 20\(^o\)C excluding the influence of ambient dissolved oxygen concentration. This is not the maxuimum rate of activation. The same rate is applied to all refractory dissolved organic constituents (carbon, nitrogen and phosphorus) |
ref activation == science |
RPOM specific shading coefficient | \(Ke^{RPOM}\) |
/(mg RPOM/L)/m /(mmol RPOM/m\(^3\))/m |
(0, 8.33) (0, 0.1) |
0 0 |
The rate at which light is attenuated by refractory particulate organic matter. This coefficient is per unit concentration of refractory particulate organic carbon (the units of RPOM) per metre (depth). It is multiplied by RPOM concentration in each computational cell to leave a value that has units of per meter, which is the same as standard extinction coefficients and therefore is able to be added to other (e.g. background) parallel extinction coefficients |
ref self shading == science |
CDOM multiplier | \(R_{CDOM}^{RDOM}\) | [-] | (0, 2) | 0 | The linear multiplicative factor applied to the calculated value of CDOM for use in organic matter self shading calculations. A value greater (less) than 1.0 increases (decreases) the CDOM estimate and hence increases (decreases) shading. . This coefficient is dimensionless but can be thought of as light attenuation per unit of CDOM, which are both in units of /m. Multiplying this parameter by CDOM in each computational cell therefore leaves a value that has units of per meter, which is the same as standard extinction coefficients and therefore is able to be added to other (e.g. background) parallel extinction coefficients |
ref self shading == science |
RDOM photolysis fraction | \(f_{photo}^{RDOM}\) | [-] | (0, 1) | 0 | The proportion of photolysed refractory dissolved organic matter that is is converted to labile dissolved organic matter. The remaining proportion is converted to dissolved inorganic carbon and nutrients. This multiplier is applied equally to refractory dissolved organic carbon, nitrogen and phosphorus |
ref photolysis == science |
Anaerobic oxidation of ammonium rate | \(k_{anmx}^{N_2}\) |
mg N/L/d mmol N/m\(^3\)/d |
(0, 0.5) (0, 0.036) |
0 0 |
The rate at which ammonium is oxidised (consuming nitrite) to nitrogen gas N\(_2\) at 20\(^o\)C excluding the influence of ambient ammonium and nitrite concentrations. This is not solely a per day rate, but is the mass (or mmol) rate of production of nitrogen in free nitrogen gas, per unit volume, per day. Anammox only operates when ambient dissolved oxygen concentrations are less than 0.1 mg/L. A value of 0.0 will switch off anammox |
anaerobic oxidation of ammonium == science |
Half saturation ammonium concentration for anammox | \(K_{anmx-NH_4}^{N_2}\) |
mg/L NH\(_4\)-N mmol/m\(^3\) NH\(_4\)-N |
(0, 10) (0, 714) |
1.75 125 |
The ammonium concentration for which the computed consumption of ammonium due to anaerobic oxidation would be half the corresponding anaerobic oxidation rate, excluding the influence of nitrite. Cannot be varied spatially |
anaerobic oxidation of ammonium == science |
Half saturation nitrite concentration for anammox | \(K_{anmx-NO_2}^{N_2}\) |
mg/L NO\(_2\)-N mmol/m\(^3\) NO\(_2\)-N |
(0, 10) (0, 714) |
1.75 125 |
The nitrite concentration for which the computed consumption of ammonium due to anaerobic oxidation would be half the corresponding anaerobic oxidation rate, excluding the influence of ammonium. Cannot be varied spatially |
anaerobic oxidation of ammonium == science |
Dissimilatory reduction of nitrate to ammonium rate | \(R_{DRNA}^{NO_3}\) | /d | (0, 0.05) |
0 0 |
The rate at which nitrate is reduced to ammonium at 20\(^o\)C excluding the influence of ambient nitrate concentrations. A value of 0.0 will switch off DRNA |
diss nitrate reduction to ammonium == science |
Half saturation oxygen concentration for DRNA | \(K_{DRNA-O_2}^{NO_3}\) |
mg/L O\(_2\) mmol/m\(^3\) O\(_2\) |
(0, 10) (0, 312.5) |
4 125 |
The oxygen concentration for which the computed reduction of nitrate to ammonium would be half the corresponding rate. Cannot be varied spatially |
diss nitrate reduction to ammonium == science |
Constant rate of adsorbed FRP settling | \(V_{settle}^{FRP}\) | m/d | (-50, 0) | 0 | The constant settling velocity applied to adsorbed FRP. Future releases of the WQ Module will dynamically link this settling velocity to calculations made by the TUFLOW FV ST Module. A negative value is a downwards velocity |
settling == science |
Labile organic matter stokes diameter | \(d_{lorg}\) | m | (1 e-09, 1 e -06) | 1 e-09 | The conceptual diameter of labile particulate organic matter applied in stokes calculations of settling velocity. Future releases of the WQ Module will dynamically link settling velocity to calculations made by the TUFLOW FV ST Module |
settling == science |
Labile organic matter stokes density | \(\rho_{lorg}\) | kg/m\(^3\) | (1000, 3000) | 2000 | The density of labile particulate organic matter applied in stokes calculations of settling velocity. Future releases of the WQ Module will dynamically link settling velocity to calculations made by the TUFLOW FV ST Module |
settling == science |
Labile organic matter settling velocity | \(V_{settle}^{lorg}\) | m/d | (-50, 0) | 0 | The settling velocity applied to labile particulate organic matter for constant or density related settling calculations. A negative value is a downwards velocity. Future releases of the WQ Module will dynamically link settling velocity to calculations made by the TUFLOW FV ST Module |
settling == science |
Refractory organic matter stokes diameter | \(d_{rorg}\) | m | (1 e-09, 1 e -06) | 1 e-09 | The conceptual diameter of refractory particulate organic matter applied in stokes calculations of settling velocity. Future releases of the WQ Module will dynamically link settling velocity to calculations made by the TUFLOW FV ST Module |
ref settling == science |
Refractory organic matter stokes density | \(\rho_{rorg}\) | kg/m\(^3\) | (1000, 3000) | 2000 | The density of refractory particulate organic matter applied in stokes calculations of settling velocity. Future releases of the WQ Module will dynamically link settling velocity to calculations made by the TUFLOW FV ST Module |
ref settling == science |
Refractory organic matter settling velocity | \(V_{settle}^{rorg}\) | m/d | (-50, 0) | 0 | The settling velocity applied to refractory particulate organic matter for constant or density related settling calculations. A negative value is a downwards velocity. Future releases of the WQ Module will dynamically link settling velocity to calculations made by the TUFLOW FV ST Module |
ref settling == science |
Minimum water quality calculation depth | \(d_{min-wq}\) | m | (0, UD) | 0.01 | The minimum depth at which water calculations will be computed. This depth is the depth of the entire computational column that immediately overlies the sediment (bed) for a given 2D cell. The depth is specified in the TUFLOW FV control file |
cell water quality depth == |
Disable water quality calculation | \(B_{disable-wq}\) | [-] | (0, 1) | 0 | A Boolean flag to disable (1) or enable (0) water quality calculations. The flag is specified in the TUFLOW FV control file |
disable water quality model == |
Non-kinetic equilibrium substep frequency | \(fess\) | Number of WQ Module timesteps | (1, UD) | 1 | The frequency (in WQ timesteps) at which the WQ Module calls non-kinetic equilibrium calculations. For example, a value of 12 will call non-kinetic calculations every twelfth time the WQ module is called |
wq equilibrium substeps == |
Pathogen dark death rate | \(k_d^{pth}\) | /d | (0, 1) | 0 | The natural mortality rate in freshwater at 20 \(^o\)C. Also referred to as the dark death rate |
mortality == science |
Salinity effect on mortality | \(C_{SM}^{pth}\) | /d /psu | (0, 0.1) | 0 | The salinity effect on mortality. Set this to zero to remove salinity effects on mortality |
mortality == science |
Parameter controlling salinity mortality dependence | \(\alpha\) | [-] | (0, 10) | 1 | An exponentially applied controller of the mortality response to salinity. Set this to 1.0 to make salinity dependence linear. Do not set this and \(C_{SM}^{pth}\) to both be zero |
mortality == science |
Parameter controlling DOC trophic mortality dependence | \(\beta\) | [-] | (0, 0) | 1 | An exponentially applied controller of the mortality response to DOC. Not implemented |
mortality == science |
Pathogen mortality temperature coefficient | \(\theta_{mor}^{pth}\) | [-] | (1, 1.2) | 1 | The factor to which (\(T\)-20) is raised to account for temperature effects on pathogen mortality. A value of 1.0 will remove the effect of temperature |
mortality == science |
Pathogen freshwater inactivation rate coefficient for exposure to visible light | \(k_{vis}^{pth}\) | m\(^2\) /MJ | (0, 1) | 0 | Inactivation effect of visible radiation. Excludes the effect of salinity and oxygen on inactivation |
visible inactivation == science |
Coefficient that enhances the inactivation effect of visible light under saline conditions | \(c_{S_{vis}}^{pth}\) | m\(^2\) /MJ /psu | (0, 0.007) | 0 | Enhancement of inactivation of visible radiation due to salinity. A value of 0.0 will remove the effect of salinity |
visible inactivation == science |
Half saturation oxygen concentration for visible light inactivation | \(K_{DO_{vis}}^{pth}\) | mg/L O\(_2\) |
(0, 8) (0, 250) |
4 125 |
The dissolved oxygen concentration for which the computed inactivation rate under visible light is half that which would occur without oxygen limitation. A value of 0.0 will remove the effect of oxygen |
visible inactivation == science |
Pathogen freshwater inactivation rate coefficient for exposure to UV-A light | \(k_{uva}^{pth}\) | m\(^2\) /MJ | (0, 20) | 0 | Inactivation effect of UV-A radiation. Excludes the effect of salinity and oxygen on inactivation |
uva inactivation == science |
Coefficient that enhances the inactivation effect of UV-A light under saline conditions | \(c_{S_{uva}}^{pth}\) | m\(^2\) /MJ /psu | (0, 0.007) | 0 | Enhancement of inactivation of UV-A radiation due to salinity. A value of 0.0 will remove the effect of salinity |
uva inactivation == science |
Half saturation oxygen concentration for UV-A light inactivation | \(K_{DO_{uva}}^{pth}\) | mg/L O\(_2\) |
(0, 8) (0, 250) |
4 125 |
The dissolved oxygen concentration for which the computed inactivation rate under UV-A light is half that which would occur without oxygen limitation. A value of 0.0 will remove the effect of oxygen |
uva inactivation == science |
Pathogen freshwater inactivation rate coefficient for exposure to UV-B light | \(k_{uvb}^{pth}\) | m\(^2\) /MJ | (0, 100) | 0 | Inactivation effect of UV-B radiation. Excludes the effect of salinity and oxygen on inactivation |
uvb inactivation == science |
Coefficient that enhances the inactivation effect of UV-B light under saline conditions | \(c_{S_{uvb}}^{pth}\) | m\(^2\) /MJ /psu | (0, 0.007) | 0 | Enhancement of inactivation of UV-B radiation due to salinity. A value of 0.0 will remove the effect of salinity |
uvb inactivation == science |
Half saturation oxygen concentration for UV-B light inactivation | \(K_{DO_{uvb}}^{pth}\) | mg/L O\(_2\) |
(0, 8) (0, 250) |
4 125 |
The dissolved oxygen concentration for which the computed inactivation rate under UV-B light is half that which would occur without oxygen limitation. A value of 0.0 will remove the effect of oxygen |
uvb inactivation == science |
Minimum alive pathogen concentration | \(\left[PTH_a\right]_{min}^{pth}\) | CFU/100mL | (0, 0) | NA | The minimum allowable alive (free) pathogen concentration. If the WQ Module detects a concentration less than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues. A value other than zero for this parameter would be unusual |
alive min max == science |
Maximum alive pathogen concentration | \(\left[PTH_a\right]_{max}^{pth}\) | CFU/100mL | (0, 1 e 12) | NA | The maximum allowable alive (free) pathogen concentration. If the WQ Module detects a concentration greater than this, then it will be reset to this value, and reported in the WQ Module log file. The root cause for requiring this reset should be addressed as this resetting is not intended to be a permanent solution and will cause mass conservation issues |
alive min max == science |
Target pathogen attachment fraction | \(f_{att}^{pth}\) | [-] | (0, 1) | 0 | The target fraction of attached pathogens. It is computed as the concentration of attached pathogens, divided by the sum of attached and alive pathogens for each computational |
target attached fraction == science |
Unattached pathogen settling velocity | \(V_{settle}^{pth}\) | m/d | (-50, 0) | 0 | The settling velocity of unattached pathogens. A negative value is a downwards velocity. Future releases of the WQ Module will dynamically link settling velocity to calculations made by the TUFLOW FV ST Module |
settling == science |
Attached pathogen settling velocity | \(V_{settle}^{pth_t}\) | m/d | (-50, 0) | 0 | The settling velocity of attached pathogens. A negative value is a downwards velocity. Future releases of the WQ Module will dynamically link settling velocity to calculations made by the TUFLOW FV ST Module |
settling == science |