User Tools

Site Tools


start:hype_model_description:hype_orgc

Organic carbon

Adds an ImageCaption tag
Figure 1: Schematic figure of the TOC-model.

Source of organic material

Litter production leads to the production of organic material, which increases levels fastC in top two layers in soil. The organic carbon litter fall is determined by crop resc (kg/ha/yr) gives a daily supplement to the pool during the number of days determined by parameter litterdays.

Modules (file) Procedures
npc_soil_processes (npc_soil_proc.f90)soil_carbon_processes

Soil processes

Soil pools – initial values

The initial pool size is dependent on land use and is determined by the user in the file par.txt. The parameters (humusc1, humusc2, humusc3, fastc1, fastc2, fastc3) give OC content of the three soil layers. The unit for these parameter values is mg/m3. With this information, the pools the size in the different layers are calculated. The model transforms pools into the unit kg/km2 by taking into account the thickness in the layers.

Common functions

Many soil processes depend on temperature and soil moisture. These use the same functions as for nitrogen and phosphorus (see above). Organic carbon soil transformations (production of humusC from fastC, turnover of fastC, and turnover of humusC) use the soil moisture function with parameters given by the user instead of the coefficients described for nutrients. The coefficient thetalow is replaced by the land-use dependent parameter ocsoilslp, and the coefficient satact is replaced by land-use dependent parameter ocsoilsat.

Adds an ImageCaption tag
Figure 21: Organic carbon processes in soil.

Production of humusC from fastC

Some of the litter fall (added to fastC) is converted into humus. This means that the model has a transformation of fastC to humusC in the uppermost soil layer.

In the other layers (k) is also a transition from fastC to humusC. The loss of fastC does not all go to humusC but a proportion (parameter minc) is mineralized. The transformation depends on soil moisture and temperature, amount of fastC and vegetation dependent parameter klh.

littertohumus(k) = klh * tmpfcn(k) * smfcn(k) * fastC(k)

Turnover of fastC

Turnover of fastC is a sink for fastC and a source of dissolved OC in soil water in all soil layers (k = 1-3). The loss of fastC does not all go to the OC, but a proportion (parameter minc) is mineralized. Turnover (transfC, mg/m2/d) depends on a general parameter (klo), the temperature function (tempfcn), humidity function (smfcn) and the pool of fastC (fastC).

transfC(k) = klo * tempfcn(k) * smfcn(k) * fastC(k)

In dry conditions a flow in the opposite direction can also occur. The transformation of OC to fastC is a decrease of OC and a source of fastC in all soil layers (k = 1-3). The loss of OC is not all to fastC but a proportion (parameter minc) is mineralized. Turnover (doctofast, mg/m2/d) depends on a general parameter (kof) and the pool of OC (OCpool). The flow is limited that the soil layer temperature must be less than 5 °C, the soil moisture (sm) must be less than field capacity and moisture function (smfcn) must be less than the parameter koflim.

doctofast(k) = kof * OCpool(k)

Turnover of humusC

Turnover of humusC is a sink for humusC and a source of OC in all soil layers (k = 1.3). The loss of humusC does not all go to the TOC, ut a proportion (parameter minc) is mineralized. Turnover (transhC, mg/m2/d) depends on a general parameter (kho), temperature function (tempfcn), humidity function (smfcn) and the pool of humusC (humusC).

transhC(k) = kho * tempfcn(k) * smfcn(k) * humusC(k)

Percolation

Organic carbon is lost from the soil water as it flows down through the soil layers and where it is dissipated to become a regional groundwater flow. The decrease in concentration depends on soil moisture and temperature and a calibration parameter.

conc = conc*(1. - par*tmpfcn*smfcn)

The soil moisture function and temperature function are the general functions described for soil processes. Percolation uses the coefficients for soil moisture function, not the parameters as the transformations. The parameter par in the equation is called kcgwreg for regional groundwater flow formation and koc for percolation between soil layers. Both are general parameters.

Modules (file) Procedures Section
npc_soil_processes (npc_soil_proc.f90)initiate_soil_npcinitial values
soil_carbon_processesproduction of humusC from fastC, turnover
soil_carbon_pool_transformations
doc_percolation_reductionpercolation

Riparian zone

Runoff from soil layers flows through a riparian zone before it reaches the local river. Surface runoff and drainage water from drainage pipes reaches the local river without passing through the riparian zone. In the riparian zone the levels of OC are affected, while flows remain unchanged. The change depends on soil temperature, class altitude (elev (in masl)), the water table (qwat) and its recent change, season and soil moisture (sm). The runoff concentration (conc(i)) of each soillayer (k) increases with the factor:

f(k)=1+ripz×tmpfcn(k)×({elev}/{100})×f_{grw}×f_{season}×f_{sm}

conc(k)=f×conc(k), ~~   k=1..3

The temperature function (tmpfcn) is the usual of soil processes (see above). The following equations describe the other process functions:

f_{grw}=e^{ripe×gwat}

f_{season}=delim{lbrace}{
 matrix{2}{2}{
   rips autumn
   1 otherwise}
   }{}

f_{sm}=delim{lbrace}{
 matrix{3}{2}{
   0 sm≤wp
  {f_2 (sm)} {wp<sm<wp} 
  satact sm≥pw}
   }{}

f_2(sm)=delim{lbrace}{
 matrix{2}{2}{
  {min(1,satact+(1-satact)×{pw-sm}/{d×\Theta_{upp}},{sm-wp}/{d×\Theta_{low} })} {rising grw}
  {min(1,satact+(1-satact)×{pw-sm}/{d×\Theta_{upp}},satact×{sm-wp}/{d×\Theta_{low} })} {sinking grw}}
   }{}

Three general parameters can be changed; ripz which determines the overall level of increase in concentration in the riparian zone, ripe which determines the groundwater level dependence, and rips which determines the seasonal influence. Season division is determined by ten-day and twenty-day averages of air temperature (T10, T20). Autumn is defined as the period when T10 is less than T20. The soil moisture function is different for an increasing (rising) and sinking ground water table (figure 2). It contains coefficients \Theta_{upp} = 0.12, \Theta_{low} = 0.08 and saturation (satact = 0.6). It depends on the soil moisture of all layers together (sm) and the water-holding capacity parameters; wp - wilting border and pw - total pore volume, in fractions of total soil layer thickness (d).

Adds an ImageCaption tag
Figure 2: Example of riparian zone soil moisture function, and the dependence on changes in the groundwater levels.
Modules (file) Procedures
npc_soil_processes (npc_soil_proc.f90)class_riparian_zone_processes
riparian_moisture_factor

Rivers and lakes

Primary production and mineralization

Primary production is a source of organic carbon in rivers and lakes, while mineralization is a sink. Primary production and mineralization is calculated the same way as for nitrogen, but with its own calibration parameter (wprodc). The potential carbon transformation (minprodCpot, kg / day) is proportional to the potential nitrogen transformation (minprodNpot, see NP section) with a transformation rate that depends on the carbon-nitrogen ratio (NCratio = 5.7). The calculated mineralization of organic carbon is limited to a maximum of 50% of the available OC pool. The parameters that affect the process of organic carbon are the rate parameter wprod and the lake region dependent parameter TPmean. If phosphorus is not modelled a long-term average total phosphorus concentration as a lake region dependent parameter (tpmean) is used. If set, the long-term average concentration is reduced by the general parameter limsedPP before using it in the concentration function.

tmpfcn1 = watertemp / 20.

tmpfcn2 = (T10 - T20) / 5.

tmpfcn = tmpfcn1*tmpfcn2

TPfcn = (TPconc-limsedPP) / (TPconc-limsedPP + halfsatTPwater)

minprodNpot = wprodc * TPfcn * tmpfcn * area

minprodCpot = minprodNpot * NCratio

Sedimentation

Sedimentation in lakes is a sink for OC and works the same way as for organic nitrogen and particulate phosphorus. Sedimentation (sedOC, kg/day) is calculated as a function of water concentration and lake area (area). The parameter sedoc is general or can be specified for each lake.

sedOC = sedoc * waterconcOC * area

Modules (file) Procedures Sections
npc_surfacewater_processes (npc_sw_proc.f90)oc_processes_in_riverprimary production and mineralization
oc_production_mineralisation
calculate_lake_tpmean
calculate_river_tpmean
oc_processes_in_lakeprimary production and mineralization
sedimentation
oc_sedimentationsedimentation
start/hype_model_description/hype_orgc.txt · Last modified: 2017/01/11 15:33 by cpers