Personal tools

Difference between revisions of "Module SedimentQuality"

From MohidWiki

Jump to: navigation, search
(→‎Main Processes)
(→‎Main Processes)
Line 25: Line 25:
 
In the case that oxygen term is used in death rates it will become (1 - OxygenTerm) since term the value 1 represents a maximum death rate and 0 represents no death.
 
In the case that oxygen term is used in death rates it will become (1 - OxygenTerm) since term the value 1 represents a maximum death rate and 0 represents no death.
  
−
  <math>Term = min\left [\frac{Conc}{ConcOpt} ;1 \right ]</math>
+
  <math>Term = min\left [\frac{Conc}{ConcOpt} , 1 \right ]</math>
  
 
'''SubstrateTerm in Death Rates'''
 
'''SubstrateTerm in Death Rates'''

Revision as of 11:24, 2 June 2011

Overview

This is the 0D module that handles sediment biological processes made by bacteria as organic matter mineralization, nitrification, denitrification and immobilization.

Main Processes

Main processes in SedimentQuality are organic matter mineralization, nitrification, denitrification and immobilization driven by microorganisms. These processes and microorganisms death are simulated using a standard rate (s.day-1.pop-1) that is afected by temperature, pH, oxygen in case of aerobic processes and sustrate concentrations if exist. An old formulation based on RZWQM was adapted in order to be consistent in units and in the new formulation a maximum rate exists and "stress" factors (from 0 to 1) are computed accordingly to pH, aerobiose and aerobiose conditions, oxygen and susbtrate availability. Aerobiose and aerobiose conditions are computed from realtive water content.

A maximum rate is computed calculating the standard rate at the actual temperature (day-1.pop-1) times bacteria population that conducts the process. 

 MaxRate (day-1) = Standard Rate at actual Temperature (day-1.pop-1) * BacteriaPopulation 

 ActualRate (day-1) = MaxRate (day-1) * Aerobiose/AnaerobioseTerm * pHTerm * OxygenTerm * SubstrateTerm


pHTerm This term has a value of 1 in pHopt and symetric descending behaviour up to acid and basic pH's. In the case that is used in death rates the term is (1 - pHterm) since value 1 (maximum death rate) will be for the worst conditions (pH acid or basic) and 0 for best conditions. 

if pH > pHopt => pHused = 2* pHopt - pH

pHTerm=\frac{pHused}{2\times pHopt-pHused}

OxygenTerm and Sustrate in non death rates The oxygen term and substrate term represent the abundance of substance that has not a negative effect on high concentrations. As so, for higher concentrations than optimum the value will be always 1. for lower concentrations than optimum there is a linear decrease of the term up to zero concentration and zero term. In the case that oxygen term is used in death rates it will become (1 - OxygenTerm) since term the value 1 represents a maximum death rate and 0 represents no death. 

Term = min\left [\frac{Conc}{ConcOpt} , 1 \right ]

SubstrateTerm in Death Rates In the case that the substrate term is used in death rates it becomes a exponential decay with increasing substrate (higher than a mimimum) to represent the same behaviour as old formulation. For concentration lower than the minimum, the term is 1 and death rate is maximum. For concentrations higher than the minimum the term as an exponential decay showing the effect of increasing material on lowering death rate. 

Term=min\left [\frac{MinConc}{Conc} ;1 \right ]

Organic Matter Mineralization

Organic matter is mineralized to ammonia with heterotrophic bacteria and also immobilization may occur in order to maintain bacteria C:N ratios.


LPOM mineralization rate 

Insert formula here


RPOM mineralization rate

Insert formula here


Immobilization rate

Insert formula here


Nitrification

Nitrification rate

Insert formula here


Denitrification

'Denitrification rate

Insert formula here


Microorganisms Growth

Microorganisms incorporate carbon, nitrogen and phosphorus in the processes described above. However, they also have sinks as respiration and death.

Death Rates

Insert formula here

Other Features

Outputs

Output is done in terms of timeseries, HDF, and boxes following the MOHID standards.

References

Data File

Keywords

Properties Units:
dissolved                 - as from called - in MOHID Land mg/L
adsorbed/particulated     - as from called - in MOHID Land mg/kgsoil
gases:
CO2 and CH4               - same as adsorbed/particulated properties
N2                        - same as dissolved properties
exceptions:
oxygen                    - mol/L
hydrogen                  - mol/L
microorganisms population - #org/kgsoil
wind                      - km/day

OXYGEN   	            :  [0/1]       0      !Connect/Disconnect Oxygen computation 
SOL_BACTERIA 	            :  [0/1]       0	  !Connect/Disconnect Solubilizing bacteria computation 

DTSECONDS                  :  [s]       86400.   !dt to evaluate
NO3_LIMIT                  :  [mg/L]      0.     !Minimum value for denitrification or maximum value for 
                                                   !methane production in organic matter decay
NEW_RATES                  : [0/1]        0      !Connect/Disconnect new rates formulation using maximum * factors
IMOBILIZATION              : [0/1]        1      !Connect/Disconnect immobilization

EXPLICIT                   : [0/1]        1      !Explicit computation for time dicretization

Sample

New Formulation

!!!!!!!!!!!!!!!!!MAIN SEDIMENT OPTIONS!!!!!!!!!!!!!!!!
!++++++++++++++++++++++++++++++++++++++++++++++++++++++

!NITROGEN 	: 1
!CARBON   	: 1
!PHOSPHORUS  	: 1	
OXYGEN   	: 1 
SOL_BACTERIA 	: 0	 

DTSECONDS    : 86400.
NO3_LIMIT    : 0.     ![mg/L] Minimum value for denitrification or maximum value for methane production in organic matter decay
NEW_RATES    : 1      !1 - new rates using maximum * factors; 0 - old formulation
!IMOBILIZATION : 0    !1 - immobilization is connected; 0 - immobilization is disconnected   

EXPLICIT : 1          !explicit method
 
 
!!!!!!!!!!!!!!!!!!!!!!RATES!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!++++++++++++++++++++++++++++++++++++++++++++++++++++++

	!!!! Organic Matter Decay !!!!

<begin_Labil_OM_C_Rate>
AE           : 15.1         ![kcal.mole-1] Initial Activation Energy
Acoef        : 0.391E-9     ![s.day-1.pop-1] Rate Coefficient
kp 	      : 12.3	     ![L/mol] Salinity coefficient		
Temperature  : 35           ![ºC] Optimal Temperature     

OptimumpH    : 7            ![] Optimum pH for reaction
ConcOptO2    : 10           ![mg/L] Maxium conc. above rate will be maximum.
<end_Labil_OM_C_Rate>

<begin_RefractOM_C_Rate>
AE           : 15.1
Acoef        : 0.653E-11     
kp 	      : 12.3
Temperature  : 35

OptimumpH    : 7            ![] Optimum pH for reaction
ConcOptO2    : 10           ![mg/L] Maxium conc. above rate will be maximum.
<end_RefractOM_C_Rate>
 

	!!!! Microorganisms Death !!!!!

<begin_Heterotrophs_Rate>
AE           : 15.19
Acoef        : 1.01E-33    
kp 	      : 12.3
Temperature  : 35

ConcMinCarbon : 1000.       ![kg/kgsoil] - labil and refractary organic carbon. below rate is maximum
OptimumpH    : 7            ![] Optimum pH for reaction
ConcOptO2    : 10           ![mg/L] Maxium conc. above rate will be maximum.
<end_Heterotrophs_Rate>

<begin_Autotrophs_Rate>
AE           : 12.64
Acoef        : 1.35E-25    
kp 	      : 12.3
Temperature  : 35 

ConcMinAmmonia : 1.         ![mg/L] - ammonia min. below rate is maximum
OptimumpH    : 7            ![] Optimum pH for reaction
ConcOptO2    : 10           ![mg/L] Maxium conc. above rate will be maximum.
<end_Autotrophs_Rate>

<begin_Anaerobic_Rate>
AE           : 12.83
Acoef        : 4.522E-30    
kp 	      : 12.3
Temperature  : 35 

ConcMinNitrate : 1.           ![mg/L] - nitrate min. below rate is maximum
ConcMinCarbon  : 1000.        ![kg/kgsoil] - labil and refractary organic carbon. below rate is maximum
OptimumpH      : 7            ![] Optimum pH for reaction
<end_Anaerobic_Rate>
 

	!!!! Nitrification !!!!!

<begin_AmmoniaToNitrate_Rate>
AE           : 21
Acoef        : 0.     !Need for a standard value.
kp 	      : 12.3
Temperature  : 35

OptimumpH    : 7            ![] Optimum pH for reaction
ConcOptO2    : 10           ![mg/L] Maxium conc. above rate will be maximum.
<end_AmmoniaToNitrate_Rate>


	!!!! Ammonia Immobilization !!!!!

<begin_AmmoniaImobilization_Rate>
AE           : 15.1
Acoef        : 5E-07
kp 	      : 12.3
Temperature  : 35
<end_AmmoniaImobilization_Rate>


	!!!! Denitrification !!!!!

<begin_NitrateToNgas_Rate>
AE           : 12.6
Acoef        : 1.431E-10  
kp 	      : 0
Temperature  : 35 

ConcOptCarbon : 10000.      ![kg/kgsoil] - labil and refractary organic carbon. below rate is maximum
OptimumpH     : 7           ![] Optimum pH for reaction
<end_NitrateToNgas_Rate>
 

	!!!! Nitrate Immobilization !!!!!

<begin_NitrateImobilization_Rate>
AE           : 15.1
Acoef        : 5E-07
kp 	      : 12.3
Temperature  : 35
<end_NitrateImobilization_Rate>


	!!!! Urea hydrolysis !!!!!

<begin_Urea_Hydrolysis>
AE	     : 12.6					
Acoef       : 2.50E-04
kp 	     : 0			
Temperature : 35
<end_Urea_Hydrolysis>


	!!!! Phsphorus Immobilization !!!!!

<begin_PhosphorusImobilization_Rate>
AE	      : 15.1					
Acoef        : 1.955E-9    
kp 	      : 12.3			
Temperature  : 35

OptimumpH    : 7            ![] Optimum pH for reaction
ConcOptO2    : 10           ![mg/L] Maxium conc. above rate will be maximum.
<end_PhosphorusImobilization_Rate>


	!!!! Methane Production !!!!! -> testing

<begin_methane_production>
AE	     : 12.6					
Acoef       : 1.252E-12      
kp 	     : 0		
Temperature : 35 

OptimumpH   : 7            ![] Optimum pH for reaction
<end_methane_production>


	!!!! Solubilizing Bacteria  !!!!!

<begin_Sol_Rate>	               
AE           : 15.19
Acoef        : 5.00E-35
kp 	      : 12.3 
Temperature  : 35

OptimumpH    : 7            ![] Optimum pH for reaction
ConcOptO2    : 10           ![mg/L] Maxium conc. above rate will be maximum.
<end_Sol_Rate>

<begin_Labil_OM_C_sol_Rate>
AE           : 1.64E-04
Acoef        : 7.92E-10
kp 	      : 12.3
Temperature  : 35     

OptimumpH    : 7            ![] Optimum pH for reaction
ConcOptO2    : 10           ![mg/L] Maxium conc. above rate will be maximum.
<end_Labil_OM_C_sol_Rate>


<begin_RefractOM_C_sol_Rate>
AE           : 2.02E-04
Acoef        : 2.58E-10
kp 	      : 12.3
Temperature  : 35 

OptimumpH    : 7            ![] Optimum pH for reaction
ConcOptO2    : 10           ![mg/L] Maxium conc. above rate will be maximum.
<end_RefractOM_C_sol_Rate>


<begin_Solubilizing_Rate>	        
AE           : 12.6	                
Acoef        : 1.10E-13			
kp 	      : 0
Temperature  : 35 

OptimumpH    : 7            ![] Optimum pH for reaction
ConcOptO2    : 10           ![mg/L] Maxium conc. above rate will be maximum.
<end_Solubilizing_Rate>




!!!!!!!!!!!!!!!!!!MICROORGANISMS!!!!!!!!!!!!!!!!!!!!!!!
!++++++++++++++++++++++++++++++++++++++++++++++++++++++

<begin_Heterotrophs>
CN_RATIO                : 3.40        
CP_RATIO		: 12.09       
POPULATION_CARBON_RATIO : 950        ![(#org/kgsoil)/(mgC/kgsoil)]
CARBON_EFICIENCY        : 0.13        
NITROGEN_EFICIENCY      : 0           
MINIMUM_POPULATION      : 50000.     ![#org/kgsoil]
<end_Heterotrophs>

<begin_Autotrophs>
CN_RATIO                : 2.01 
CP_RATIO		: 13.82    
POPULATION_CARBON_RATIO : 9500
CARBON_EFICIENCY        : 0
NITROGEN_EFICIENCY      : 0.01
MINIMUM_POPULATION      : 500.
<end_Autotrophs>

<begin_Anaerobic>
CN_RATIO                : 3.40
CP_RATIO		: 13.82    
POPULATION_CARBON_RATIO : 9500
CARBON_EFICIENCY        : 0.1
NITROGEN_EFICIENCY      : 0.1
MINIMUM_POPULATION      : 5000.
<end_Anaerobic>

<begin_Sols>
CN_RATIO                :8
CP_RATIO		:60    
POPULATION_CARBON_RATIO :950
CARBON_EFICIENCY        :0.1
NITROGEN_EFICIENCY      :0.1
<end_Sols>

Old Formulation

!!!!!!!!!!!!!!!!!MAIN SEDIMENT OPTIONS!!!!!!!!!!!!!!!!
!++++++++++++++++++++++++++++++++++++++++++++++++++++++

!NITROGEN 	: 1
!CARBON   	: 1
!PHOSPHORUS  	: 1	
OXYGEN   	: 1 
SOL_BACTERIA 	: 0	 

DTSECONDS    : 86400.
NO3_LIMIT    : 0.     ![mg/L] Minimum value for denitrification or maximum value for methane production in organic matter decay
NEW_RATES    : 0      !1 - new rates using maximum * factors; 0 - old formulation
!IMOBILIZATION : 0    !1 - immobilization is connected; 0 - immobilization is disconnected   

EXPLICIT : 1          !explicit method
 
 
!!!!!!!!!!!!!!!!!!!!!!RATES!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!++++++++++++++++++++++++++++++++++++++++++++++++++++++

	!!!! Organic Matter Decay !!!!

<begin_Labil_OM_C_Rate>
AE           : 15.1         ![kcal.mole-1] Initial Activation Energy
Acoef        : 1E-7         ![s.day-1.pop-1] Rate Coefficient
kp 	      : 12.3	     ![L/mol] Salinity coefficient		
Temperature  : 35           ![ºC] Optimal Temperature     
<end_Labil_OM_C_Rate>

<begin_RefractOM_C_Rate>
AE           : 15.1
Acoef        : 1.67E-9     
kp 	      : 12.3
Temperature  : 35
<end_RefractOM_C_Rate>
 

	!!!! Microorganisms Death !!!!!

<begin_Heterotrophs_Rate>
AE           : 15.19
Acoef        : 5.00E-35    
kp 	      : 12.3
Temperature  : 35
<end_Heterotrophs_Rate>

<begin_Autotrophs_Rate>
AE           : 12.64
Acoef        : 6.68E-30    
kp 	      : 12.3
Temperature  : 35 
<end_Autotrophs_Rate>

<begin_Anaerobic_Rate>
AE           : 12.83
Acoef        : 4.47E-30    
kp 	      : 12.3
Temperature  : 35 
<end_Anaerobic_Rate>
 

	!!!! Nitrification !!!!!

<begin_AmmoniaToNitrate_Rate>
AE           : 21
Acoef        : 0.     !Need for a standard value.
kp 	      : 12.3
Temperature  : 35
<end_AmmoniaToNitrate_Rate>


	!!!! Ammonia Immobilization !!!!!

<begin_AmmoniaImobilization_Rate>
AE           : 15.1
Acoef        : 5E-07
kp 	      : 12.3
Temperature  : 35
<end_AmmoniaImobilization_Rate>


	!!!! Denitrification !!!!!

<begin_NitrateToNgas_Rate>
AE           : 12.6
Acoef        : 1.00E-15  
kp 	      : 0
Temperature  : 35 
<end_NitrateToNgas_Rate>
 

	!!!! Nitrate Immobilization !!!!!

<begin_NitrateImobilization_Rate>
AE           : 15.1
Acoef        : 5E-07
kp 	      : 12.3
Temperature  : 35
<end_NitrateImobilization_Rate>


	!!!! Urea hydrolysis !!!!!

<begin_Urea_Hydrolysis>
AE	     : 12.6					
Acoef       : 2.50E-04
kp 	     : 0			
Temperature : 35
<end_Urea_Hydrolysis>


	!!!! Phsphorus Immobilization !!!!!

<begin_PhosphorusImobilization_Rate>
AE	      : 15.1					
Acoef        : 5E-7    
kp 	      : 12.3			
Temperature  : 35
<end_PhosphorusImobilization_Rate>


	!!!! Methane Production !!!!! -> testing

<begin_methane_production>
AE	     : 12.6					
Acoef       : 1.0E-13      
kp 	     : 0		
Temperature : 35 
<end_methane_production>


	!!!! Solubilizing Bacteria  !!!!!

<begin_Sol_Rate>	               
AE           : 15.19
Acoef        : 5.00E-35
kp 	      : 12.3 
Temperature  : 35
<end_Sol_Rate>

<begin_Labil_OM_C_sol_Rate>
AE           : 1.64E-04
Acoef        : 7.92E-10
kp 	      : 12.3
Temperature  : 35     
<end_Labil_OM_C_sol_Rate>


<begin_RefractOM_C_sol_Rate>
AE           : 2.02E-04
Acoef        : 2.58E-10
kp 	      : 12.3
Temperature  : 35 
<end_RefractOM_C_sol_Rate>


<begin_Solubilizing_Rate>	        
AE           : 12.6	                
Acoef        : 1.10E-13			
kp 	      : 0
Temperature  : 35 
<end_Solubilizing_Rate>




!!!!!!!!!!!!!!!!!!MICROORGANISMS!!!!!!!!!!!!!!!!!!!!!!!
!++++++++++++++++++++++++++++++++++++++++++++++++++++++

<begin_Heterotrophs>
CN_RATIO                : 3.40        
CP_RATIO		: 12.09       
POPULATION_CARBON_RATIO : 950        ![(#org/kgsoil)/(mgC/kgsoil)]
CARBON_EFICIENCY        : 0.13        
NITROGEN_EFICIENCY      : 0           
MINIMUM_POPULATION      : 50000.     ![#org/kgsoil]
<end_Heterotrophs>

<begin_Autotrophs>
CN_RATIO                : 2.01 
CP_RATIO		: 13.82    
POPULATION_CARBON_RATIO : 9500
CARBON_EFICIENCY        : 0
NITROGEN_EFICIENCY      : 0.01
MINIMUM_POPULATION      : 500.
<end_Autotrophs>

<begin_Anaerobic>
CN_RATIO                : 3.40
CP_RATIO		: 13.82    
POPULATION_CARBON_RATIO : 9500
CARBON_EFICIENCY        : 0.1
NITROGEN_EFICIENCY      : 0.1
MINIMUM_POPULATION      : 5000.
<end_Anaerobic>

<begin_Sols>
CN_RATIO                :8
CP_RATIO		:60    
POPULATION_CARBON_RATIO :950
CARBON_EFICIENCY        :0.1
NITROGEN_EFICIENCY      :0.1
<end_Sols>
Retrieved from "http://wiki.mohid.com/index.php?title=Module_SedimentQuality&oldid=4761"