CAEDYM biogeochemical parameters:
! CAEDYM PARAMETER FILE: (v2.3)- Lake Benmore simulation, June 2009
! LIGHT constants *************************************************************!
!------!
!-- Extinction FractionOfIncidentSW Bandwidth
1.00000 0.410 : NIR, Near InfraRed
0.15000 0.450 : PAR, Photosynthetically Active
1.50000 0.035 : UVA, UV-A Region
2.50000 0.005 : UVB, UV-B Region
!------!
! PHYTOPLANKTON constants *****************************************************!
!------!
! Pmax (/day) : Maximum potential growth rate of phytoplankton !
0.60000 ! Phytoplankton group 1 (CAEDYM group no. 2)
1.30000 ! Phytoplankton group 2 (CAEDYM group no. 7)
! Ycc (mg C/mg chla) : Average ratio of C to chlorophyll a !
50.00000 ! Phytoplankton group 1
50.00000 ! Phytoplankton group 2
!------!
! Light limitation (2=no photoinhibition, 3=photoinhibition) !
! algt (no units) : Type of light limitation algorithm !
2 ! Phytoplankton group 1
2 ! Phytoplankton group 2
! IK (microE/m^2/s) : Parameter for initial slope of P_I curve !
70.00000 ! Phytoplankton group 1
20.00000 ! Phytoplankton group 2
! Kep (ugchlaL^-1m^-1) : Specific attenuation coefficient !
0.02000 ! Phytoplankton group 1
0.02000 ! Phytoplankton group 2
!------!
! Nutrient Parameters !
! KP (mg/L) : Half saturation constant for phosphorus !
0.00200 ! Phytoplankton group 1
0.00250 ! Phytoplankton group 2
! Po (mg/L) : Low concentrations of PO4 at which uptake ceases !
0.00000 ! Phytoplankton group 1
0.00000 ! Phytoplankton group 2
! KN (mg/L) : Half saturation constant for nitrogen !
0.01400 ! Phytoplankton group 1
0.01500 ! Phytoplankton group 2
! No (mg/L) : Low concentrations of N at which uptake ceases !
0.00000 ! Phytoplankton group 1
0.00000 ! Phytoplankton group 2
! Sicon (mg Si/mg Chla) : Constant internal Silica concentration !
0.00000 ! Phytoplankton group 1
0.00000 ! Phytoplankton group 2
! KSi (mg/L) : Half saturation constant for silica !
0.00000 ! Phytoplankton group 1
0.00000 ! Phytoplankton group 2
! Sio (mg/L) : Low concentrations of Si at which uptake ceases !
0.00000 ! Phytoplankton group 1
0.00000 ! Phytoplankton group 2
! KCa (mg/L) : Half saturation constant for carbon !
2.00000 ! Phytoplankton group 1
2.00000 ! Phytoplankton group 2
! INmin (mg N/mg Chla) : Minimum internal N concentration !
3.00000 ! Phytoplankton group 1
3.00000 ! Phytoplankton group 2
! INmax (mg N/mg Chla) : Maximum internal N concentration !
9.00000 ! Phytoplankton group 1
9.00000 ! Phytoplankton group 2
! UNmax (mg N/mg Chla/day) : Maximum rate of Phytoplankton nitrogen uptake!
3.00000E+00 ! Phytoplankton group 1
3.00000E+00 ! Phytoplankton group 2
! IPmin (mg P/mg Chla) : Minimum internal P concentration !
0.23000 ! Phytoplankton group 1
0.20000 ! Phytoplankton group 2
! IPmax (mg P/mg Chla) : Maximum internal P concentration !
2.00000 ! Phytoplankton group 1
2.00000 ! Phytoplankton group 2
! UPmax (mg P/mg Chla/day) : Maximum rate of phosphorus uptake !
0.70000E+00 ! Phytoplankton group 1
0.90000E+00 ! Phytoplankton group 2
! ICmin (mg C/mg Chla) : Minimum internal C concentration !
15.00000 ! Phytoplankton group 1
15.00000 ! Phytoplankton group 2
! ICmax (mg C/mg Chla) : Maximum internal C concentration !
80.00000 ! Phytoplankton group 1
80.00000 ! Phytoplankton group 2
! UCmax (mg C/mg Chla/day) : Maximum rate of carbon uptake !
50.00000 ! Phytoplankton group 1
50.00000 ! Phytoplankton group 2
! NFixationRate (mg N/mg Chla /day): Maximum nitrogen fixation rate !
2.00000! ! Phytoplankton group 1
2.00000 ! Phytoplankton group 2
! gthRedNFix () : Growth rate reduction under maximal N fixation !
1.00000 ! Phytoplankton group 1
1.00000 ! Phytoplankton group 2
!------!
! Temperature representation !
! vT (no units) : Temperature multiplier !
1.08000 ! Phytoplankton group 1
1.06000 ! Phytoplankton group 2
! Tsta (Deg C) : Standard temperature !
20.00000 ! Phytoplankton group 1
18.00000 ! Phytoplankton group 2
! Topt (Deg C) : Optimum temperature !
28.00000 ! Phytoplankton group 1
23.00000 ! Phytoplankton group 2
! Tmax (Deg C) : Maximum temperature !
35.00000 ! Phytoplankton group 1
30.00000 ! Phytoplankton group 2
!------!
! Respiration mortality and excretion. !
! kr (/day) : Respiration rate coefficient !
0.16000 ! Phytoplankton group 1
0.1800 ! Phytoplankton group 2
! vR (no units) : Temperature multiplier (no units) !
1.06000 ! Phytoplankton group 1
1.08000 ! Phytoplankton group 2
! Fraction of respiration relative to total metabolic loss rate !
0.70000 ! Phytoplankton group 1
0.70000 ! Phytoplankton group 2
! Fraction of metabolic loss rate that goes to DOM (rest goes to POM) !
0.50000 ! Phytoplankton group 1
0.50000 ! Phytoplankton group 2
!------!
! Vertical migration and settling (0-stokes, 1-constant, 2-motile w/o photoinhibition 3-motile w/ photoinhibition!
! phvel (no units) : Type of vertical migration algorithm !
1 ! Phytoplankton group 1
1 ! Phytoplankton group 2
! ws (ms^-1) : Constant settling velocity !
0.00000E+00 ! Phytoplankton group 1
-0.50000E-05 ! Phytoplankton group 2
!------!
! Resuspension !
! tcpy (N/m^2) : Critical shear stress !
1.00000 ! Phytoplankton group 1
1.00000 ! Phytoplankton group 2
! alpPy (mg Chla/m^2/s) : Resuspension rate constant !
0.008
! KTPy (mg Chla/m^2) : Controls rate of resuspension !
0.01000 ! Phytoplankton group 1
0.01000 ! Phytoplankton group 2
! DTphy (days) : Phytoplankton sediment survival time !
2.00000 ! Phytoplankton group 1
2.00000 ! Phytoplankton group 2
!------!
! DISSOLVED OXYGEN constants **************************************************!
!------!
1.20000 vOP (no units) : Temperature multiplier for SOD !
50.00000 PCmax (g/m^2) : Maximum limit of polychaete biomass !
2.66667 YOC (mg C/mg O) : Respiration stoichiometric ratio of C to O2!
0.10000 fox (no units) :Frac of net DO allocated to seagrass roots !
2.66667 YSG (mg seagC/mg O) : Stoichiometric factor, seagrass C : DO !
2.66667 YOJ (mg jelC/mg O) : Stoichiometric factor, jellyfish C : DO !
0.03000 oxmin (mg/L) : Minimum DO in the bottom layer (mg/L) !
1.10000 rSOs (g/m^2/day) : Static sediment exchange rate !
0.50000 KSOs (mg O/L) : 1/2 sat constant for static DO sediment flux !
0.01400 prc (no units) : Photo-respiration phytoplankton DO loss!
!------!
! C, N & P CYCLE Constants ****************************************************!
!------!
! ORGANIC PARTICLES (POM) ------! !
0.0150 ! POC1max (/day) : Max transfer of POCL->DOCL !
0.0005 ! POC2max (/day) : Max transfer of POCR->DOCR !
0.0150 ! POP1max (/day) : Max transfer of POPL->DOPL !
0.0005 ! POP2max (/day) : Max transfer of POPR->DOPR !
0.0150 ! PON1max (/day) : Max transfer of PONL->DONL !
0.0005 ! PON2max (/day) : Max transfer of PONR->DONR !
! POMDiameter (m) : Diameter of POM particles !
2.00000E-05 ! POM1 (eg. LABILE) !
5.00000E-05 !POM2 (eg. REFRACTORY) !
! POMDensity (kg/m^3) : Density of POM particles !
0.11000E+04 ! POM1 (eg.LABILE) !
0.10500E+04 !POM2 (eg. REFRACTORY) !
! tcPOM (N/m^2) : Critical shear stress for respn !
0.0200 ! POM1 (eg.LABILE) !
0.0100 ! POM2 (eg.REFRACTORY) !
! KePOC (mg/L/m) : Specific attenuation coefficient of POC !
0.02000E+00 ! POC1 (eg.LABILE) !
0.02000E+00 !POC2 (eg. REFRACTORY) !
!------! !
! DISSOLVED ORGANICS (DOM) !
0.010 ! DOC1max (/day) : Max mineralisation of DOCL->DIC !
0.0005 ! DOC2max (/day) : Max mineralisation of DOCR->DIC !
0.015 ! DOP1max (/day) : Max mineralisation of DOPL->PO4 !
0.0005 ! DOP2max (/day) : Max mineralisation of DOPR->PO4 !
0.015 ! DON1max (/day) : Max mineralisation of DONL->NH4 !
0.0005 ! DON2max (/day) : Max mineralisation of DONR->NH4 !
! KeDOC (mg/L/m) : Specific attenuation coefficient of DOC !
0.01000E+00 ! DOC1 (eg. LABILE) !
0.01000E+00 ! DOC2 (eg. REFRACTORY) !
! Salinity bounds for DOCr flocculation !
1.0000 ! Smindoc !
10.000 ! Smaxdoc !
! kfloc: flocculation rate constant !
0.50000 !
! kSWNP: Rate of DOCr Photolytic Decay !
0.30000 !
!------!
! DISSOLVED INORGANICS !
! Ionic Strength (required for DIC & pH calculation) !
0.00150 Ionic Strength (g/m3) !
! Adsorption/Desorption !
0.00000 Kadd1NH4 !
0.00000 Kadd2NH4 !
0.00000 Kadd1PO4 !
0.00000 Kadd2PO4 !
! Nitrification/Denitrification !
1.08000 vN2 (-) : Temp multiplier for denitrification !
0.02000 koN2 (/day) : Denitrification rate coefficient !
2.50000 KN2 (mg/L) : Half sat const for denitrification !
1.08000 vON (-) : Temp multiplier for nitrification !
0.16000 koNH (/day) : Nitrification rate coefficient !
2.00000 KOn (mg O/L) : Half sat constant for nitrification !
3.42857 YNH (mg N/mg O) : Ratio of O2 to N for nitrification !
!------! !
! SEDIMENT PARAMETERS !
! NUTRIENT FLUXES !
! Theta(sed) : Temp multiplier of sediment fluxes !
1.08 !
! PO4 sediment flux parameters !
0.00260 SmpPO4 (g/m2/day) : Release rate of PO4 !
0.50000 KDOS-PO4 (g/m^3) : Controls sed release of PO4 via O !
1000000.0 KpHS-PO4 (-) : Controls sed release of PO4 via pH !
! NH4 sediment flux parameters !
0.00050 !**23Apr09** was 0.001 SmpNH4 (g/m2/day) : Release rate of NH4 !
0.50000 KDOS-NH4 (g/m^3) : Controls sed release of NH4 via O !
1000000.0 KpHS-NH4 (-) : Controls sed release of NH4 via pH !
! NO3 sediment flux parameters !
-0.04 SmpNO3 (g/m2/day) : Release rate of NO3 !
0.5000 KDOS-NO3 (g/m^3) : Controls sed release of NO3 via O !
1000000.0 KpHS-NO3 (-) : Controls sed release of NO3 via pH !
! Si sediment flux parameters !
0.03850 SmpSi (g/m2/day) : Release rate of Si !
0.50000 KDOS-Si (g/m^3) : Controls sed release of Si via O !
1000000.0 KpHS-Si (-) : Controls sed release of Si via pH !
! DOC sediment flux parameters !
0.00000 !**23Apr09** was 0.092 SmpdocL (g/m2/day): Release rate of DOCL !
0.00100 SmpdocR (g/m2/day): Release rate of DOCR !
0.50000 KDOS-doc (g/m^3) : Controls sed release of DOC via O!
1000000.0 KpHS-doc (-) : Controls sed release of DOC via pH !
! DOP sediment flux parameters !
0.00000 !**23Apr09** was 0.00001 SmpdopL (g/m2/day): Release rate of DOPL!
0.00001 SmpdopR (g/m2/day): Release rate of DOPR !
0.50000 KDOS-dop (g/m^3) : Controls sed release of DOP via O !
1000000.0 KpHS-dop (-) : Controls sed release of DOP via pH !
! DON sediment flux parameters !
0.00000 !**23Apr09** was 0.00001 SmpdonL (g/m2/day): Release rate of DONL !
0.00001 SmpdonR (g/m2/day): Release rate of DONR !
0.50000 KDOS-don (g/m^3) : Controls sed release of DON via O !
1000000.0 KpHS-don (-) : Controls sed release of DON via pH !
!--- Sediment composition !
! Fraction of sediment that is organics !
0.15 sedOrganicFrac !
! Sediment Porosity (i.e. porewater fraction) !
0.10 sedPorosity !
! Composite resuspension rate (g/m2/day) !
0.005 resusRate !
! Half sat conc for resus dependence on sed mass (g) !
1.0 resusKT !
!------!
! INORGANIC PARTICLE constants ************************************************!
!------!
! Density of suspended solid particles !
0.26500E+04 deS (kg/m^3) : SSOL1 !
0.26500E+04 deS (kg/m^3) : SSOL2 !
! Diameter of suspended solids groups !
2.00000E-06 diaSS (m) : SSOL1 !
10.00000E-06 diaSS (m) : SSOL2 !
! Specific attenuation coefficient !
0.03000E+00 KeSS (mg^-1Lm^-1) : SSOL1 !
0.03000E+00 KeSS (mg^-1Lm^-1) : SSOL2 !
! Critical shear stress !
0.05000E+00 tcs (N/m2) : SSOL1 !
0.05000E+00 tcs (N/m2) : SSOL2 !
!------!
! IRON constants **************************************************************!
!------!
1.05000 vFeR (-) : Temp multiplier for reduction !
0.03000 kFeR (/day) : Maximum reduction rate !
0.50000 K_FeR (mg/L) : Half saturation for reduction !
1.05000 vFeO (-) : Temp multiplier for oxidation !
0.06000 kFeO (/day) : Maximum oxidation rate !
2.00000 K_FeO (mg/L) : Half saturation for oxidation !
0.44000 SFe (g/m^2/day) : Sediment release rate !
0.0400 KDOFe (mg/L) : Oxygen sediment half saturation!
1000000 KpHFe (-) : pH sediment half saturation !
0.04000 tcFe (N/m^2) : Critical shear stress !
-0.00001 Feset (m/s) : Settling velocity of particulate Fe !
!------!
! MANGANESE constants *********************************************************!
!------!
1.05000 vMnR (no units) : Temperature multiplier for reduction!
0.03 kMnR (/day) : Maximum reduction rate !
0.50000 K_MnR (mg/L) : Half saturation for reduction !
1.05000 vMnO (no units) : Temperature multiplier for oxidation!
0.10 kMnO (/day) : Maximum oxidation rate !
0.50000 K_MnO (mg/L) : Half saturation for oxidation!
0.08600 Smn (g/m^2/day) : Sediment release rate !
0.40000 KDOMn (mg/L) : Oxygen sediment half saturation!
1000000 KpHMn (no units) : pH sediment half saturation !
0.04000 tcMn (N/m^2) : Critical shear stress !
-0.00001 Mnset (m/s) : Settling velocity of particulate Mn !
!------!
! ALUMINIUM constants *********************************************************!
!------!
1.08000 vAlR (no units) : Temperature multiplier for reduction!
0.03000 kAlR (/day) : Maximum reduction rate !
0.50000 K_AlR (mg/L) : Half saturation for reduction !
1.08000 vAlO (no units) : Temperature multiplier for oxidation!
0.10000 kAlO (/day) : Maximum oxidation rate !
2.00000 K_AlO (mg/L) : Half saturation for oxidation !
0.00450 S_AL (g/m^2/day) : Sediment release rate !
0.50000 KDOAl (mg/L) : Oxygen sediment half saturation!
1000000 KpHAl (no units) : pH sediment half saturation!
0.04000 tcAl (N/m^2) : Critical shear stress !
-0.00001 Alset (m/s) : Settling velocity of particulate Al !
!------!
! PATHOGEN and MICROBIAL INDICATOR ORGANISM constants ************************!
!------!
!-- Mortality rate (Dark death rate) @ 20C and 0 psu !
0.03 ! Cryptosporidium oocysts !
0.71 ! Fecal Coliforms !
0.48 ! E. coli !
!-- Temperature multiplier for mortality ! !
1.145 ! Cryptosporidium oocysts !
1.060 ! Fecal Coliforms !
1.110 ! E. coli !
!-- Salinity effect on mortality ! !
0.00 ! Cryptosporidium oocysts !
0.05 ! Fecal Coliforms !
0.05 ! E. coli !
!-- Light inactivation------:!
!----
Visible UV-A UV-B ! !
0.00E-00 0.00E-00 3.50E-06 ! Cryptosporidiumoocysts !
1.48E-08 1.18E-08 6.19E-08 ! FecalColiforms !
1.48E-08 1.18E-08 6.19E-08 ! E. coli !
!-- Salinity effect on light inactivation! !
0.00 ! Cryptosporidium oocysts !
0.00 ! Fecal Coliforms !
0.00 ! E. coli !
!-- DO effect on light ! !
0.00 ! Cryptosporidium oocysts !
0.00 ! Fecal Coliforms !
0.00 ! E. coli !
!-- Loss rate due to grazing ! !
0.00 ! Cryptosporidium oocysts !
0.00 ! Fecal Coliforms !
0.00 ! E. coli !
!-- Sedimentation velocity (m/s) at 20C (-ve means down) !
-1.00E-6 ! Cryptosporidium oocysts !
-1.00E-6 ! Fecal Coliforms !
-1.00E-6 ! E. coli !
!-- Resuspension constants------: !
!----
ResusRate critStress ResusHalfSat !
0.01 0.005 1.0E8 ! Cryptosporidium oocysts !
0.01 0.010 1.0E8 ! Fecal Coliforms !
0.01 0.010 1.0E8 ! E. coli !
!------!