Personal tools

Difference between revisions of "Module CEQUALW2"

From MohidWiki

Jump to: navigation, search
(Keywords)
(Properties)
 
(11 intermediate revisions by the same user not shown)
Line 1: Line 1:
Module CEQUALW2 is a module using the [[CE-QUAL-W2]] ecological model formulations. CE-QUAL-W2 is 2D laterally average hydrodynamic and ecological model developed at the U.S. Corps of Engineers and it's normally used to simulate water quality in reservoirs. Its ecological formulations were adapted and re-programmed into MOHID as Module CEQUALW2. This module is able to simulate 22 properties, including temperature, nutrients (nitrogen, phosphorus and silica biogeochemical cycles), oxygen  and several species of algae (microalgae). The model does not simulate macroalgae, neither the influence of zooplankton in the primary production.
+
Module CEQUALW2 is a module using the [[CE-QUAL-W2]] ecological model formulations. CE-QUAL-W2 is 2D laterally average hydrodynamic and ecological model developed at the U.S. Corps of Engineers and it's normally used to simulate water quality in reservoirs. Its ecological formulations were adapted and re-programmed into MOHID as Module CEQUALW2 which made it possible to run 3D applications with these ecological formulations. This module is able to simulate 22 properties, including temperature, nutrients (nitrogen, phosphorus and silica biogeochemical cycles), oxygen  and several species of algae (microalgae). The model does not simulate macroalgae, neither the influence of zooplankton in the primary production.
 +
 
 +
=Processes=
 +
 
 +
'''Properties'''
 +
 
 +
Temperature
 +
Microalgae (algae and epiphytes) [[File:esquema_cequal.png|thumb|center|400px|CEQUALW2 processes integrated in the MOHID system]]
 +
Cohesive sediments
 +
Detritus
 +
Oxygen
 +
Carbon dioxide
 +
Ph
 +
Alcalinity
 +
Bicarbonate
 +
Carbonate
 +
BOD
 +
Nitrate
 +
Ammonia
 +
Refractory particulate organic matter
 +
Labile particulate organic matter
 +
Refractory dissolved organic matter
 +
Labile dissolved organic matter
 +
Inorganic phosphorous
 +
Particulate silica
 +
Dissolved silica
 +
Inorganic carbon
  
 
=Keywords=
 
=Keywords=
Line 5: Line 31:
 
'''CequalW2'''
 
'''CequalW2'''
 
 
Keywords                      :Data Type    Default      Comment
+
Keywords                      :Data Type    Default      Comment
 
+
DTSECONDS                    :Real          1200        ! Time step
DTSECONDS                    :Real          1200        ! Time step
 
 
 
<begin_algae>
 
NAME                          :Char                      !Name of the algae
 
 
 
A_GROWTH                      :Real          2            !DEFAULT: 2 - Maximum algae growth
 
 
 
A_RESPIRATION                :Real          0.04        !Maximum algae respiration
 
 
 
A_EXCRETION                  :Real          0.04        !Maximum algae excretion
 
 
 
A_MORTALITY                  :Real          0.1          !Maximum algae mortality
 
 
 
A_HALFSAT_P                  :Real          0.003        !Algal half-saturation P [g m^-3]
 
 
 
A_HALFSAT_N                  :Real          0.014        !Algal half-saturation N [g m^-3]
 
 
 
A_HALFSAT_SI                  :Real          0            !Algal half-saturation N [g m^-3]
 
 
 
A_LIGHT_SAT                  :Real          75          !Algal saturating light intensity at maximum phtosynthetic rate [W m^-2]
 
 
 
A_T1                          :Real          5            !Lower temperature for algal growth (ºC)
 
 
 
A_T2                          :Real          25          ! Lower temperature for maximum algal growth (ºC)
 
 
 
A_T3                          :Real          35          ! Upper temperature for maximum algal growth
 
 
 
A_T4                          :Real          40          ! Upper temperature for algal growth
 
 
 
A_K1                          :Real          0.1          !Fraction of algal growth rate at AT1
 
 
 
A_K2                          :Real          0.99        !Fraction of maximum algal growth rate at AT2
 
 
 
A_K3                          :Real          0.99        !Fraction of maximum algal growth rate at AT3
 
 
 
A_K4                          :Real          0.1          !Fraction of algal growth rate at AT4
 
 
 
A_OK1                        :Real          13
 
 
 
A_OK2                        :Real          11          !Algal half-saturation coefficients for oxygen consumption [mgO2/l]
 
 
 
A_OK3                        :Real          2.5
 
 
 
A_OK4                        :Real          7
 
 
 
A_STOICHIOMETRY_P            :Real          0.005        !Algal stoichiometric coefficient for phosphorus
 
 
 
A_STOICHIOMETRY_N            :Real          0.08        !Algal stoichiometric coefficient for nitrogen
 
 
 
A_STOICHIOMETRY_C            :Real          0.45        !Algal stoichiometric coefficient for carbon
 
 
 
A_STOICHIOMETRY_Si            :Real          0.18        !Algal stoichiometric coefficient for silica
 
 
 
A_POM                        :Real          0.8          !Algal stoichiometric coefficient for POM
 
 
 
A_NEQUATIONNUMBER            :Integer      2            !Equation for preference factor (either 1 or 2)
 
 
 
A_AMMONIUM_PREF              :Real          0.001        !half-saturation copreference constant
 
 
 
O2_A_RESPIRATION              :Real          1.1          !Stoichiometry coefficient in respiration
 
 
 
O2_A_GROWTH                  :Real          1.4          !Stoichiometry coefficient in growth
 
 
 
<end_algae>
 
 
 
<begin_om>
 
 
 
LDOM_DECAY                    :Real          0.1          !LDOM decay rate [day^-1]
 
 
 
RDOM_DECAY                    :Real          0.001        !RDOM decay rate [day^-1
 
 
 
LRDOM_DECAY                  :Real          0.01        !Labile to refractory DOM dekay rate [day^-1]
 
 
 
LPOM_DECAY                    :Real          0.08        !LPOM decay rate [day^-1]
 
 
 
RPOM_DECAY                    :Real          0.001        !RPOM decay rate [day^-1]
 
 
 
LRPOM_DECAY                  :Real          0.01        !LRPOM decay rate [day^-1]
 
 
 
OM_STOICHIOMETRY_P            :Real          0.005        !Stoichiometric coef. for phosphorus
 
 
 
OM_STOICHIOMETRY_N            :Real          0.08        !Stoichiometric coef. for nitrogen
 
 
 
OM_STOICHIOMETRY_C            :Real          0.45        !Stoichiometric coef. for carbon
 
 
 
OM_STOICHIOMETRY_SI          :Real          0.18        !Stoichiometric coef. for silica
 
 
 
OM_T1                        :Real          4            !Organic Matter Temperature Rate Multipliers
 
 
 
OM_T2                        :Real          25
 
 
 
OM_K1                        :Real          0.1          !Organic Matter Temperature Rate Multipliers
 
 
 
OM_K2                        :Real          0.99
 
  
O2_OM                        :Real          1.4          !O2OM - Oxygen stoichiometry for organic matter decay
+
<begin_algae>
 +
NAME                          :Char                      !Name of the algae
 +
A_GROWTH                      :Real          2            !DEFAULT: 2 - Maximum algae growth
 +
A_RESPIRATION                :Real          0.04        !Maximum algae respiration
 +
A_EXCRETION                  :Real          0.04        !Maximum algae excretion
 +
A_MORTALITY                  :Real          0.1          !Maximum algae mortality
 +
A_HALFSAT_P                  :Real          0.003        !Algal half-saturation P [g m^-3]
 +
A_HALFSAT_N                  :Real          0.014        !Algal half-saturation N [g m^-3]
 +
A_HALFSAT_SI                  :Real          0            !Algal half-saturation N [g m^-3]
 +
A_LIGHT_SAT                  :Real          75          !Algal saturating light intensity at maximum phtosynthetic rate [W m^-2]
 +
A_T1                          :Real          5            !Lower temperature for algal growth (ºC)
 +
A_T2                          :Real          25          ! Lower temperature for maximum algal growth (ºC)
 +
A_T3                          :Real          35          ! Upper temperature for maximum algal growth
 +
A_T4                          :Real          40          ! Upper temperature for algal growth
 +
A_K1                          :Real          0.1          !Fraction of algal growth rate at AT1
 +
A_K2                          :Real          0.99        !Fraction of maximum algal growth rate at AT2
 +
A_K3                          :Real          0.99        !Fraction of maximum algal growth rate at AT3
 +
A_K4                          :Real          0.1          !Fraction of algal growth rate at AT4
 +
A_OK1                        :Real          13
 +
A_OK2                        :Real          11          !Algal half-saturation coefficients for oxygen consumption [mgO2/l]
 +
A_OK3                        :Real          2.5
 +
A_OK4                        :Real          7
 +
A_STOICHIOMETRY_P            :Real          0.005        !Algal stoichiometric coefficient for phosphorus
 +
A_STOICHIOMETRY_N            :Real          0.08        !Algal stoichiometric coefficient for nitrogen
 +
A_STOICHIOMETRY_C            :Real          0.45        !Algal stoichiometric coefficient for carbon
 +
A_STOICHIOMETRY_Si            :Real          0.18        !Algal stoichiometric coefficient for silica
 +
A_POM                        :Real          0.8          !Algal stoichiometric coefficient for POM
 +
A_NEQUATIONNUMBER            :Integer      2            !Equation for preference factor (either 1 or 2)
 +
A_AMMONIUM_PREF              :Real          0.001        !half-saturation copreference constant
 +
O2_A_RESPIRATION              :Real          1.1          !Stoichiometry coefficient in respiration
 +
O2_A_GROWTH                  :Real          1.4          !Stoichiometry coefficient in growth
 +
<end_algae>
 +
<begin_om>
 +
LDOM_DECAY                    :Real          0.1          !LDOM decay rate [day^-1]
 +
RDOM_DECAY                    :Real          0.001        !RDOM decay rate [day^-1
 +
LRDOM_DECAY                  :Real          0.01        !Labile to refractory DOM dekay rate [day^-1]
 +
LPOM_DECAY                    :Real          0.08        !LPOM decay rate [day^-1]
 +
RPOM_DECAY                    :Real          0.001        !RPOM decay rate [day^-1]
 +
LRPOM_DECAY                  :Real          0.01        !LRPOM decay rate [day^-1]
 +
OM_STOICHIOMETRY_P            :Real          0.005        !Stoichiometric coef. for phosphorus
 +
OM_STOICHIOMETRY_N            :Real          0.08        !Stoichiometric coef. for nitrogen
 +
OM_STOICHIOMETRY_C            :Real          0.45        !Stoichiometric coef. for carbon
 +
OM_STOICHIOMETRY_SI          :Real          0.18        !Stoichiometric coef. for silica
 +
OM_T1                        :Real          4            !Organic Matter Temperature Rate Multipliers
 +
OM_T2                        :Real          25
 +
OM_K1                        :Real          0.1          !Organic Matter Temperature Rate Multipliers
 +
OM_K2                        :Real          0.99
 +
O2_OM                        :Real          1.4          !O2OM - Oxygen stoichiometry for organic matter decay
 +
<end_om>
  
<end_om>
+
<begin_BOD>
 +
BOD_DECAY                    :Real          0.25        !CBOD Decay Rate [day^-1]
 +
BOD_T_COEF                    :Real          1.0147      !BOD Temperature Rate Multiplier
 +
BOD_RATIO                    :Real          1.85        !Ratio of CBOD5 to ultimate CBOD
 +
BOD_STOICHIOMETRY_P          :Real          0.004        !Phosphorus/CBOD stochiometric ratio
 +
BOD_STOICHIOMETRY_N          :Real          0.06        !Nitrogen/CBOD stochiometric ratio
 +
BOD_STOICHIOMETRY_C          :Real          0.32        !Carbon/CBOD stochiometric ratio
 +
<end_BOD>
  
<begin_BOD>
+
<begin_oxygen>
BOD_DECAY                     :Real          0.25        !CBOD Decay Rate [day^-1]
+
O2_METHOD                     :Integer      1            !Method to compute oxygen: 1 – Simple with O2_LIM; 2 – Monod curve; 3 – Ecotoxicity curve
BOD_T_COEF                    :Real          1.0147      !BOD Temperature Rate Multiplier
+
O2LIM                        :Real          0.1          !Dissolved oxygen concentration at which anaerobic processes begin [g m^-3]
BOD_RATIO                    :Real          1.85        !Ratio of CBOD5 to ultimate CBOD
+
O2_K1                        :Real          2.5          !Toxicity equation coefficients for oxygen consumption
 +
O2_K2                        :Real          7
 +
<end_oxygen>
  
BOD_STOICHIOMETRY_P           :Real          0.004        !Phosphorus/CBOD stochiometric ratio
+
<begin_nitro>
BOD_STOICHIOMETRY_N          :Real          0.06         !Nitrogen/CBOD stochiometric ratio
+
NH4_DECAY                    :Real          0.12        !Ammonium Decay Rate [day^-1]
BOD_STOICHIOMETRY_C           :Real          0.32         !Carbon/CBOD stochiometric ratio
+
NH4_T1                        :Real          5            !Minimum Temperature (T1)
<end_BOD>
+
NH4_T2                        :Real          25           !upper temperature (T2)
 +
NH4_K1                        :Real          0.1          !Multiplier factor for T1
 +
NH4_K2                        :Real          0.99        !Multiplier factor for T2
 +
NO3_DECAY                    :Real          0.03         !Nitrate Decay Rate [day^-1]
 +
NO3_T1                        :Real          5            !Minimum temperature (T1)
 +
NO3_T2                        :Real          25           !Optimal temperature (T2)
 +
NO3_K1                        :Real          0.1          !Multiplier factor for T1
 +
NO3_K2                        :Real          0.99        !Multiplier factor for T2
 +
O2_NH4                        :Real          4.57         !O2NH4 - Oxygen stoichiometry for nitrification
 +
<end_nitro>
  
<begin_oxygen>
+
<begin_epiphyton>
O2_METHOD                    :Integer      1           !Method to compute oxygen: 1 – Simple with O2_LIM; 2 – Monod curve; 3 – Ecotoxicity curve
+
NAME                          :Char                      !Name of the epiphyte
O2LIM                        :Real          0.1          !Dissolved oxygen concentration at which anaerobic processes begin [g m^-3]
+
E_GROWTH                      :Real          2           !Maximum epiphyte growth
O2_K1                        :Real          2.5          !Toxicity equation coefficients for oxygen consumption
+
E_RESPIRATION                :Real          0.04        !Maximum epiphyte respiration
O2_K2                        :Real          7
+
E_EXCRETION                  :Real          0.04        !Maximum epiphyte excretion
<end_oxygen>
+
E_MORTALITY                  :Real          0.1          !DEFAULT: 0.15 - Maximum epiphyte mortality
<begin_nitro>
+
E_HALFSAT_P                  :Real          0.003        ! epiphyte half-saturation P [g m^-3]
NH4_DECAY                    :Real          0.12        !Ammonium Decay Rate [day^-1]
+
E_HALFSAT_N                  :Real          0.014        ! epiphyte half-saturation N [g m^-3]
NH4_T1                        :Real          5            !Minimum Temperature (T1)
+
E_HALFSAT_SI                  :Real          0            ! epiphyte half-saturation N [g m^-3]
NH4_T2                        :Real          25          !upper temperature (T2)
+
E_LIGHT_SAT                  :Real          75          ! epiphyte saturating light intensity at maximum phtosynthetic rate [W m^-2]
NH4_K1                        :Real          0.1          !Multiplier factor for T1
+
E_T1                          :Real          5            !Lower temperature for epiphyte growth (ºC)
NH4_K2                        :Real          0.99        !Multiplier factor for T2
+
E_T2                          :Real          25          ! Lower temperature for maximum epiphyte growth (ºC)
NO3_DECAY                    :Real          0.03         !Nitrate Decay Rate [day^-1]
+
E_T3                          :Real          35          ! Upper temperature for maximum epiphyte growth
NO3_T1                        :Real          5            !Minimum temperature (T1)
+
E_T4                          :Real          40          ! Upper temperature for epiphyte growth
NO3_T2                        :Real          25          !Optimal temperature (T2)
+
E_K1                          :Real          0.1          !Fraction of epiphyte growth rate at AT1
NO3_K1                        :Real          0.1         !Multiplier factor for T1
+
E_K2                          :Real          0.99        !Fraction of maximum epiphyte growth rate at AT2
NO3_K2                        :Real          0.99        !Multiplier factor for T2
+
E_K3                          :Real          0.99         !Fraction of maximum epiphyte growth rate at AT3
O2_NH4                        :Real          4.57        !O2NH4 - Oxygen stoichiometry for nitrification
+
E_K4                          :Real          0.1         !Fraction of epiphyte growth rate at AT4
<end_nitro>
+
E_STOICHIOMETRY_P            :Real          0.005        !epiphyte stoichiometric coefficient for phosphorus
 +
E_STOICHIOMETRY_N            :Real          0.08        !epiphyte stoichiometric coefficient for nitrogen
 +
E_STOICHIOMETRY_C            :Real          0.45        !epiphyte stoichiometric coefficient for carbon
 +
E_STOICHIOMETRY_Si            :Real         0.18        !epiphyte stoichiometric coefficient for silica
 +
E_POM                        :Real          0.8          !epiphyte stoichiometric coefficient for POM
 +
E_NEQUATIONNUMBER            :Integer      2            !Equation for preference factor (either 1 or 2)
 +
E_AMMONIUM_PREF              :Real          0.001        !half-saturation copreference constant
 +
O2_E_RESPIRATION              :Real          1.1          !Stoichiometry coefficient in respiration
 +
O2_E_GROWTH                  :Real          1.4          !Stoichiometry coefficient in growth
 +
< end_epiphyton >
  
<begin_epiphyton>
+
<begin_phos>
NAME                          :Char                      !Name of the epiphyte
+
OM_STOICHIOMETRY_P            :Real          0.005        !Stoichiometric coef. for phosphorous
E_GROWTH                      :Real          2            !Maximum epiphyte growth
+
<end_phos>
E_RESPIRATION                :Real          0.04        !Maximum epiphyte respiration
 
E_EXCRETION                  :Real          0.04        !Maximum epiphyte excretion
 
E_MORTALITY                  :Real          0.1          !DEFAULT: 0.15 - Maximum epiphyte mortality
 
E_HALFSAT_P                  :Real          0.003        ! epiphyte half-saturation P [g m^-3]
 
E_HALFSAT_N                  :Real          0.014        ! epiphyte half-saturation N [g m^-3]
 
E_HALFSAT_SI                  :Real          0            ! epiphyte half-saturation N [g m^-3]
 
E_LIGHT_SAT                  :Real          75          ! epiphyte saturating light intensity at maximum phtosynthetic rate [W m^-2]
 
E_T1                          :Real          5            !Lower temperature for epiphyte growth (ºC)
 
E_T2                          :Real          25          ! Lower temperature for maximum epiphyte growth (ºC)
 
E_T3                          :Real          35          ! Upper temperature for maximum epiphyte growth
 
E_T4                          :Real          40          ! Upper temperature for epiphyte growth
 
E_K1                          :Real          0.1          !Fraction of epiphyte growth rate at AT1
 
E_K2                          :Real          0.99        !Fraction of maximum epiphyte growth rate at AT2
 
E_K3                          :Real          0.99        !Fraction of maximum epiphyte growth rate at AT3
 
E_K4                          :Real          0.1          !Fraction of epiphyte growth rate at AT4
 
E_STOICHIOMETRY_P            :Real          0.005        !epiphyte stoichiometric coefficient for phosphorus
 
E_STOICHIOMETRY_N            :Real          0.08        !epiphyte stoichiometric coefficient for nitrogen
 
E_STOICHIOMETRY_C            :Real          0.45        !epiphyte stoichiometric coefficient for carbon
 
E_STOICHIOMETRY_Si            :Real          0.18        !epiphyte stoichiometric coefficient for silica
 
E_POM                        :Real          0.8          !epiphyte stoichiometric coefficient for POM
 
E_NEQUATIONNUMBER            :Integer      2            !Equation for preference factor (either 1 or 2)
 
E_AMMONIUM_PREF              :Real          0.001        !half-saturation copreference constant
 
O2_E_RESPIRATION              :Real          1.1          !Stoichiometry coefficient in respiration
 
O2_E_GROWTH                  :Real          1.4          !Stoichiometry coefficient in growth
 
< end_epiphyton >
 
<begin_phos>
 
OM_STOICHIOMETRY_P            :Real          0.005        !Stoichiometric coef. for phosphorous
 
<end_phos>
 
  
 +
<begin_ammonia>
 +
OM_STOICHIOMETRY_N            :Real          0.08        !Stoichiometric coef. for phosphorous
 +
<end_ammonia>
  
<begin_ammonia>
+
<begin_silica>
OM_STOICHIOMETRY_N            :Real          0.08        !Stoichiometric coef. for phosphorous
+
PARTSI_DECAY                  :Real          0.3          !Particulate biogenic silica decay rate [day^-1]
<end_ammonia>
+
<end_silica>
  
<begin_silica>
+
<begin_dsi>
PARTSI_DECAY                  :Real          0.3          !Particulate biogenic silica decay rate [day^-1]
+
OM_STOICHIOMETRY_SI          :Real          0.18        !Stoichiometric coef. for silica
<end_silica>
+
<end_dsi>
  
<begin_dsi>
+
<begin_ic>
OM_STOICHIOMETRY_SI          :Real          0.18         !Stoichiometric coef. for silica
+
OM_STOICHIOMETRY_C            :Real          0.45         !Stoichiometric coef. for carbon
<end_dsi>
+
<end_ic>
  
<begin_ic>
+
'''BenthicCequalW2'''
OM_STOICHIOMETRY_C            :Real          0.45        !Stoichiometric coef. for carbon
 
<end_ic>
 
  
Keywords BenthicCequalW2
+
DTSECONDS                    :Real          1200
DTSECONDS                    :Real          1200
+
<begin_SOD>
<begin_SOD>
+
PO4R                          :Real          0.001        ! Release of PO4 by SOD rate
PO4R                          :Real          0.001        ! Release of PO4 by SOD rate
+
NH4R                          :Real          0.001        ! Release of PO4 by SOD rate
NH4R                          :Real          0.001        ! Release of PO4 by SOD rate
+
SiR                          :Real          0.1          ! Release of Silica by SOD rate
SiR                          :Real          0.1          ! Release of Silica by SOD rate
+
CO2R                          :Real          0.1          ! Release of CO2 by SOD rate
CO2R                          :Real          0.1          ! Release of CO2 by SOD rate
+
O2Consumption                :Integer      1            ! Sink of O2 by SOD rate
O2Consumption                :Integer      1            ! Sink of O2 by SOD rate
+
DefaultO2                    :Integer      1            ! If off, consumes O2 only when there is no Oxygen
DefaultO2                    :Integer      1            ! If off, consumes O2 only there is no Oxygen
+
SODT1                        :Real          4            ! T1 for temperature rate multiplier
SODT1                        :Real          4            ! T1 for temperature rate multiplier
+
SODT2                        :Real          35          ! T2 for temperature rate multiplier
SODT2                        :Real          35          ! T2 for temperature rate multiplier
+
SODK1                        :Real          0.1          ! K1 for temperature rate multiplier
SODK1                        :Real          0.1          ! K1 for temperature rate multiplier
+
SODK2                        :Real          0.99        ! K2 for temperature rate multiplier
SODK2                        :Real          0.99        ! K2 for temperature rate multiplier
+
<end_SOD>
<end_SOD>
 
  
<begin_det>
+
<begin_det>
DET_DECAY                    :Real          0.1          ! Sediment decay rate, [day^-1]
+
DET_DECAY                    :Real          0.1          ! Sediment decay rate, [day^-1]
DET_T1                        :Real          4            ! T1 for temperature rate multiplier
+
DET_T1                        :Real          4            ! T1 for temperature rate multiplier
DET_T2                        :Real          30          ! T2 for temperature rate multiplier
+
DET_T2                        :Real          30          ! T2 for temperature rate multiplier
DET_K1                        :Real          0.1          ! K1 for temperature rate multiplier
+
DET_K1                        :Real          0.1          ! K1 for temperature rate multiplier
DET_K2                        :Real          0.99        ! K2 for temperature rate multiplier
+
DET_K2                        :Real          0.99        ! K2 for temperature rate multiplier
<end_det>
+
<end_det>
  
 
[[Category:Modules]]
 
[[Category:Modules]]
 
[[Category:MOHID Base 1]]
 
[[Category:MOHID Base 1]]
 
[[Category:Biogeochemistry]]
 
[[Category:Biogeochemistry]]

Latest revision as of 12:16, 8 April 2014

Module CEQUALW2 is a module using the CE-QUAL-W2 ecological model formulations. CE-QUAL-W2 is 2D laterally average hydrodynamic and ecological model developed at the U.S. Corps of Engineers and it's normally used to simulate water quality in reservoirs. Its ecological formulations were adapted and re-programmed into MOHID as Module CEQUALW2 which made it possible to run 3D applications with these ecological formulations. This module is able to simulate 22 properties, including temperature, nutrients (nitrogen, phosphorus and silica biogeochemical cycles), oxygen and several species of algae (microalgae). The model does not simulate macroalgae, neither the influence of zooplankton in the primary production.

Processes

Properties
Temperature
Microalgae (algae and epiphytes)
CEQUALW2 processes integrated in the MOHID system
Cohesive sediments
Detritus
Oxygen
Carbon dioxide
Ph
Alcalinity
Bicarbonate
Carbonate
BOD
Nitrate
Ammonia
Refractory particulate organic matter
Labile particulate organic matter
Refractory dissolved organic matter
Labile dissolved organic matter
Inorganic phosphorous
Particulate silica
Dissolved silica
Inorganic carbon

Keywords

CequalW2

Keywords                      :Data Type     Default      Comment
DTSECONDS                     :Real          1200         ! Time step
<begin_algae>
NAME                          :Char                       !Name of the algae
A_GROWTH                      :Real          2            !DEFAULT: 2 - Maximum algae growth
A_RESPIRATION                 :Real          0.04         !Maximum algae respiration
A_EXCRETION                   :Real          0.04         !Maximum algae excretion
A_MORTALITY                   :Real          0.1          !Maximum algae mortality
A_HALFSAT_P                   :Real          0.003        !Algal half-saturation P [g m^-3]
A_HALFSAT_N                   :Real          0.014        !Algal half-saturation N [g m^-3]
A_HALFSAT_SI                  :Real          0            !Algal half-saturation N [g m^-3]
A_LIGHT_SAT                   :Real          75           !Algal saturating light intensity at maximum phtosynthetic rate [W m^-2]
A_T1                          :Real          5            !Lower temperature for algal growth (ºC)
A_T2                          :Real          25           ! Lower temperature for maximum algal growth (ºC)
A_T3                          :Real          35           ! Upper temperature for maximum algal growth
A_T4                          :Real          40           ! Upper temperature for algal growth
A_K1                          :Real          0.1          !Fraction of algal growth rate at AT1
A_K2                          :Real          0.99         !Fraction of maximum algal growth rate at AT2
A_K3                          :Real          0.99         !Fraction of maximum algal growth rate at AT3
A_K4                          :Real          0.1          !Fraction of algal growth rate at AT4
A_OK1                         :Real          13
A_OK2                         :Real          11           !Algal half-saturation coefficients for oxygen consumption [mgO2/l]
A_OK3                         :Real          2.5
A_OK4                         :Real          7
A_STOICHIOMETRY_P             :Real          0.005        !Algal stoichiometric coefficient for phosphorus
A_STOICHIOMETRY_N             :Real          0.08         !Algal stoichiometric coefficient for nitrogen
A_STOICHIOMETRY_C             :Real          0.45         !Algal stoichiometric coefficient for carbon
A_STOICHIOMETRY_Si            :Real          0.18         !Algal stoichiometric coefficient for silica
A_POM                         :Real          0.8          !Algal stoichiometric coefficient for POM
A_NEQUATIONNUMBER             :Integer       2            !Equation for preference factor (either 1 or 2)
A_AMMONIUM_PREF               :Real          0.001        !half-saturation copreference constant
O2_A_RESPIRATION              :Real          1.1          !Stoichiometry coefficient in respiration
O2_A_GROWTH                   :Real          1.4          !Stoichiometry coefficient in growth
<end_algae>
<begin_om>
LDOM_DECAY                    :Real          0.1          !LDOM decay rate [day^-1]
RDOM_DECAY                    :Real          0.001        !RDOM decay rate [day^-1
LRDOM_DECAY                   :Real          0.01         !Labile to refractory DOM dekay rate [day^-1]
LPOM_DECAY                    :Real          0.08         !LPOM decay rate [day^-1]
RPOM_DECAY                    :Real          0.001        !RPOM decay rate [day^-1]
LRPOM_DECAY                   :Real          0.01         !LRPOM decay rate [day^-1]
OM_STOICHIOMETRY_P            :Real          0.005        !Stoichiometric coef. for phosphorus
OM_STOICHIOMETRY_N            :Real          0.08         !Stoichiometric coef. for nitrogen
OM_STOICHIOMETRY_C            :Real          0.45         !Stoichiometric coef. for carbon
OM_STOICHIOMETRY_SI           :Real          0.18         !Stoichiometric coef. for silica
OM_T1                         :Real          4            !Organic Matter Temperature Rate Multipliers
OM_T2                         :Real          25
OM_K1                         :Real          0.1          !Organic Matter Temperature Rate Multipliers
OM_K2                         :Real          0.99
O2_OM                         :Real          1.4          !O2OM - Oxygen stoichiometry for organic matter decay
<end_om>
<begin_BOD>
BOD_DECAY                     :Real          0.25         !CBOD Decay Rate [day^-1]
BOD_T_COEF                    :Real          1.0147       !BOD Temperature Rate Multiplier
BOD_RATIO                     :Real          1.85         !Ratio of CBOD5 to ultimate CBOD
BOD_STOICHIOMETRY_P           :Real          0.004        !Phosphorus/CBOD stochiometric ratio
BOD_STOICHIOMETRY_N           :Real          0.06         !Nitrogen/CBOD stochiometric ratio
BOD_STOICHIOMETRY_C           :Real          0.32         !Carbon/CBOD stochiometric ratio
<end_BOD>
<begin_oxygen>
O2_METHOD                     :Integer       1            !Method to compute oxygen: 1 – Simple with O2_LIM; 2 – Monod curve; 3 – Ecotoxicity curve
O2LIM                         :Real          0.1          !Dissolved oxygen concentration at which anaerobic processes begin [g m^-3]
O2_K1                         :Real          2.5          !Toxicity equation coefficients for oxygen consumption
O2_K2                         :Real          7
<end_oxygen>
<begin_nitro>
NH4_DECAY                     :Real          0.12         !Ammonium Decay Rate [day^-1]
NH4_T1                        :Real          5            !Minimum Temperature (T1)
NH4_T2                        :Real          25           !upper temperature (T2)
NH4_K1                        :Real          0.1          !Multiplier factor for T1
NH4_K2                        :Real          0.99         !Multiplier factor for T2
NO3_DECAY                     :Real          0.03         !Nitrate Decay Rate [day^-1]
NO3_T1                        :Real          5            !Minimum temperature (T1)
NO3_T2                        :Real          25           !Optimal temperature (T2)
NO3_K1                        :Real          0.1          !Multiplier factor for T1
NO3_K2                        :Real          0.99         !Multiplier factor for T2
O2_NH4                        :Real          4.57         !O2NH4 - Oxygen stoichiometry for nitrification
<end_nitro>
<begin_epiphyton>
NAME                          :Char                       !Name of the epiphyte
E_GROWTH                      :Real          2            !Maximum epiphyte growth
E_RESPIRATION                 :Real          0.04         !Maximum epiphyte respiration
E_EXCRETION                   :Real          0.04         !Maximum epiphyte excretion
E_MORTALITY                   :Real          0.1          !DEFAULT: 0.15 - Maximum epiphyte mortality
E_HALFSAT_P                   :Real          0.003        ! epiphyte half-saturation P [g m^-3]
E_HALFSAT_N                   :Real          0.014        ! epiphyte half-saturation N [g m^-3]
E_HALFSAT_SI                  :Real          0            ! epiphyte half-saturation N [g m^-3]
E_LIGHT_SAT                   :Real          75           ! epiphyte saturating light intensity at maximum phtosynthetic rate [W m^-2]
E_T1                          :Real          5            !Lower temperature for epiphyte growth (ºC)
E_T2                          :Real          25           ! Lower temperature for maximum epiphyte growth (ºC)
E_T3                          :Real          35           ! Upper temperature for maximum epiphyte growth
E_T4                          :Real          40           ! Upper temperature for epiphyte growth
E_K1                          :Real          0.1          !Fraction of epiphyte growth rate at AT1
E_K2                          :Real          0.99         !Fraction of maximum epiphyte growth rate at AT2
E_K3                          :Real          0.99         !Fraction of maximum epiphyte growth rate at AT3
E_K4                          :Real          0.1          !Fraction of epiphyte growth rate at AT4
E_STOICHIOMETRY_P             :Real          0.005        !epiphyte stoichiometric coefficient for phosphorus
E_STOICHIOMETRY_N             :Real          0.08         !epiphyte stoichiometric coefficient for nitrogen
E_STOICHIOMETRY_C             :Real          0.45         !epiphyte stoichiometric coefficient for carbon
E_STOICHIOMETRY_Si            :Real          0.18         !epiphyte stoichiometric coefficient for silica
E_POM                         :Real          0.8          !epiphyte stoichiometric coefficient for POM
E_NEQUATIONNUMBER             :Integer       2            !Equation for preference factor (either 1 or 2)
E_AMMONIUM_PREF               :Real          0.001        !half-saturation copreference constant
O2_E_RESPIRATION              :Real          1.1          !Stoichiometry coefficient in respiration
O2_E_GROWTH                   :Real          1.4          !Stoichiometry coefficient in growth
< end_epiphyton >
<begin_phos>
OM_STOICHIOMETRY_P            :Real          0.005        !Stoichiometric coef. for phosphorous
<end_phos>
<begin_ammonia>
OM_STOICHIOMETRY_N            :Real          0.08         !Stoichiometric coef. for phosphorous
<end_ammonia>
<begin_silica>
PARTSI_DECAY                  :Real          0.3          !Particulate biogenic silica decay rate [day^-1]
<end_silica>
<begin_dsi>
OM_STOICHIOMETRY_SI           :Real          0.18         !Stoichiometric coef. for silica
<end_dsi>
<begin_ic>
OM_STOICHIOMETRY_C            :Real          0.45         !Stoichiometric coef. for carbon
<end_ic>
BenthicCequalW2
DTSECONDS                     :Real          1200
<begin_SOD>
PO4R                          :Real          0.001        ! Release of PO4 by SOD rate
NH4R                          :Real          0.001        ! Release of PO4 by SOD rate
SiR                           :Real          0.1          ! Release of Silica by SOD rate
CO2R                          :Real          0.1          ! Release of CO2 by SOD rate
O2Consumption                 :Integer       1            ! Sink of O2 by SOD rate
DefaultO2                     :Integer       1            ! If off, consumes O2 only when there is no Oxygen
SODT1                         :Real          4            ! T1 for temperature rate multiplier
SODT2                         :Real          35           ! T2 for temperature rate multiplier
SODK1                         :Real          0.1          ! K1 for temperature rate multiplier
SODK2                         :Real          0.99         ! K2 for temperature rate multiplier
<end_SOD>
<begin_det>
DET_DECAY                     :Real          0.1          ! Sediment decay rate, [day^-1]
DET_T1                        :Real          4            ! T1 for temperature rate multiplier
DET_T2                        :Real          30           ! T2 for temperature rate multiplier
DET_K1                        :Real          0.1          ! K1 for temperature rate multiplier
DET_K2                        :Real          0.99         ! K2 for temperature rate multiplier
<end_det>