====== Organic carbon ======
|{{:start:hype_model_description:schematic_ocmodel.png?400|Adds an ImageCaption tag}}|
|Figure 1: Schematic figure of the organic carbon (OC) model. Text in squares symbolizes earth, while text in ellipses symbolizes water.|
===== Source of organic material =====
Litter fall in forest, harvest left overs and other plant residues is a source of organic material to HYPE. The plant residues increase the immobile pool of organic carbon (OC) in the form of //fastC// in top two layers in soil. Thus the OC model does not use the fractionation of plant material (//resfast//) into fast and humus soil pools as is done for nitrogen and phosphorus in plant residues. The organic carbon addition by plant residues is defined based on crop/vegetation. Input, //resc// (//kg/ha/yr//), gives a daily supplement to the pool during the number of days determined by parameter //litterdays//.
==== Links to relevant procedures in the code ====
^Modules (file) ^Procedures ^
|[[http://hype.sourceforge.net/doxy-html/namespacenpc__soil__processes.html|npc_soil_processes (npc_soil_proc.f90)]]|soil_substance_processes|
===== Soil processes =====
====Soil pool initial values====
The initial size of organic carbon pools in the soil is dependent on land use and determined by the user. The parameters (//humusc1, humusc2, humusc3, fastc1, fastc2, fastc3//) give immobile OC content of the three soil layers. In addition the initial value for organic carbon humus soil pool may be constant with depth for some soil types (parameter //humusC0st//). If this parameter is used (>0) for some soil type, the classes of this soil type will use the humusC0st value instead of the landuse dependent parameters. The unit for these parameter values is //mg/m3//. With this information and soil layer thickness, the size of the pools in the layers are calculated. The model works with pools of the unit //kg/km2//. In addition an initial concentration of dissolved organic carbon in the soil water of different land uses may be specified (//occonc0//). The amount of DOC in the soil water is below called the //OCpool//.
====Common functions====
Many soil processes depend on temperature and soil moisture. They use the same common functions as [[start:hype_model_description:hype_np_soil#common_functions|nitrogen and phosphorus]]. The organic carbon soil transformations (both production of humusC from fastC, and turnover of fastC and turnover of humusC) use the soil moisture function with values partly given by the user instead of the given coefficients of the equation as is used used for nitrogen and phosphorus processes. The coefficient theta_low is replaced by the land-use dependent parameter //ocsoimslp// (theta_low=ocsoimslp/100), and the coefficient //satact// is replaced by land-use dependent parameter //ocsoimsat//. The coefficients theta_upp and theta_pow keep their values (i.e. 0.12 and 1.0). The percolation reduction of OC uses all the given coefficients, and is not affected by the parameters.
|{{:start:hype_model_description:organiccarboninsoil.png?400|Adds an ImageCaption tag}}|
|Figure 2: Organic carbon processes in upper soil layer. The transformation modelled are shown by hollow green arrows. The mineralization is a fraction of the organic carbon transformed between the pools, as is shown by the little perpendicular arrow. //humusC// is here denoted //slowC//|
==== Production of humusC from fastC ====
Some of the fastC of soil is converted into humusC. For HYPE this means that fastC (the pool where plant residues were added) is transformed to humusC in the uppermost soil layer. In Figure 2 the pool of humusC is denoted slowC because of its slower transformation rates.
The other soil layers (//k//) also have a transition from fastC to humusC. The loss of fastC does not all end up in the humusC pool, but a proportion (parameter //minc//) is mineralized in the process. The transformation (//fasttohumus, mg/m2/d//) depends on soil moisture (//smfcn//) and temperature (//tmpfcn//), amount of OC in the pool (//fastC//) and a vegetation dependent parameter //klh//.
fasttohumus(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 (DOC) in all soil layers (//k// = 1-3). The loss of fastC does not all go to the OCpool, 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 transfer in the opposite direction can also occur. The transformation of DOC to fastC is a decrease of DOC and a source of fastC in all soil layers (//k// = 1-3). The loss of DOC is not all ending up in fastC but a proportion (parameter //minc//) is mineralized. Turnover (//doctofast//, mg/m2/d) depends on a general parameter (//kof//) and the pool of DOC (//OCpool//). The transfer 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 DOC in all soil layers (k = 1.3). The turnover rate of humusC is lower than that of fastC, why it is sometimes called slowC (e.g. in Figure 2). The loss of humusC does not all go to the soil water OC, but 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 percolates down through the soil layers and where it is dissipated to become a regional groundwater flow. The decrease in concentration of percolating water during transport between the soil layers depends on soil moisture and temperature and a calibration parameter.
conc_perc = conc_soillayer*(1 - par*tmpfcn*smfcn)
The soil moisture function and temperature function are the general functions described for soil processes. Percolation uses the nitrogen and phosphorus coefficients for the soil moisture function, not the parameters that the OC transformations uses. The parameter, //par// in the equation, is //kcgwreg// for regional groundwater flow formation and //koc// for percolation between soil layers. Both are general parameters.
==== Delay of organic carbon in runoff ====
The OC transported by surface and groundwater runoff and tile drainage (//runoffC//) is collected in a temporary storage pool (//relpool// (kg/km2)).
relpool = relpool + runoffC
From the temporary pool organic carbon is released (//release// (kg/km2)) and then transported to the local river depending on the size of the total runoff (//runoff// (mm)). The parameters //ocfldelx// and //ocfldele// are general parameters.
release = min(relpool, relpool *(runoff {/} ocfldelx)^{ocfldele})
The released OC give the current OC concentration of runoff.
==== Links to relevant procedures in the code ====
^Modules (file) ^Procedures ^Section ^
|[[http://hype.sourceforge.net/doxy-html/namespacenpc__soil__processes.html|npc_soil_processes (npc_soil_proc.f90)]]|initiate_soil_npc|initial values|
|:::|soil_substance_processes|production of humusC from fastC, turnover |
|:::|soil_carbon_pool_transformations|:::|
|:::|doc_percolation_reduction|percolation|
|:::|carbon_runoff_delay|delay of organic carbon in runoff|
===== Riparian zone =====
Runoff from soil may 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 in OC depends on soil temperature, class altitude (//elev// (in masl)), the water table (//gwat//) 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}/{1000})*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=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}}
}{}
The activation of riparian zone processes is based on land use. It is primarily thought to act on forest runoff. The land use dependent parameter //ripz// determines the overall level of increase in concentration in the riparian zone, and if set to zero no riparian zone processes are used. In addition, two general parameters can influence the effect of the riparian processes; //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 3). It uses 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//).
|{{:start:hype_model_description:riparianzone.png?400|Adds an ImageCaption tag}}|
|Figure 3: Example of riparian zone soil moisture function, and the dependence on changes in the groundwater levels.|
==== Links to relevant procedures in the code ====
^Modules (file) ^Procedures ^
|[[http://hype.sourceforge.net/doxy-html/namespacenpc__soil__processes.html|npc_soil_processes (npc_soil_proc.f90)]]|class_riparian_zone_processes|
|:::|riparian_moisture_factor|
===== Rivers and lakes =====
The initial organic carbon concentration in rivers are assumed to be zero, while the lakes' concentration are set by the user. The parameter, //tocmean//, is lakeregion dependent, but can also be set for each lake separately.
==== 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 are calculated the same way as for nitrogen, but with its own calibration parameter (//wprodc//). The equations are repeated below. The production/mineralization depend on temperature and total phosphorus and lake area (//area//). The potential carbon transformation (//minprodCpot//, kg / day) is proportional to the potential nitrogen transformation (//minprodNpot//, see also [[start:hype_model_description:hype_np_riv_lake#primary_production_and_mineralization|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 phosphorus dependence is based on a half-saturation equation using a long-term average of total phosphorus. If phosphorus is not modelled by HYPE, a lake region dependent parameter (//tpmean//) is used instead. The long-term average concentration of phosphorus is reduced by the general parameter //limsedPP// before using it in the concentration function. The water depth (//depth//) is the lake depth, and for the river the depth calculated [[start:hype_model_description:hype_np_riv_lake#basic_assumptions|above]].
tmpfcn1 = watertemp / 20.
tmpfcn2 = (T10 - T20) / 5.
tmpfcn = tmpfcn1*tmpfcn2
TPfcn = (TPconc-limsedPP) / (TPconc-limsedPP + halfsatTPwater)
minprodNpot = wprodc * TPfcn * tmpfcn * area * depth
minprodCpot = minprodNpot * NCratio
==== Sedimentation ====
Sedimentation in lakes is a sink for OC and works the same way as for organic nitrogen and particulate phosphorus. Thus sedimentation (sedC_{lake}, //kg/day//) is calculated as a function of OC concentration in lake water (//conc//)) and lake area (//area//). The settling velocity parameter //sedoc// is general or can be specified for each lake.
sedC_lake = sedoc * conc * area
==== Links to relevant procedures in the code ====
^Modules (file) ^Procedures ^ Sections ^
|[[http://hype.sourceforge.net/doxy-html/namespacenpc__surfacewater__processes.html|npc_surfacewater_processes (npc_sw_proc.f90)]]|substance_processes_in_river|primary production and mineralization|
|:::|oc_production_mineralisation|:::|
|:::|calculate_lake_tpmean|:::|
|:::|calculate_river_tpmean|:::|
|:::|substance_processes_in_lake|primary production and mineralization|
|:::|:::|sedimentation|
|:::|lake_sedimentation|sedimentation|
===== Links to file reference =====
^ Section ^ Symbol ^ Parameter/Data ^ File ^
| Sources of organic material | | //resc// | [[start:hype_file_reference:cropdata.txt|CropData.txt]] |
| ::: | | //litterdays// | [[start:hype_file_reference:par.txt|par.txt]] |
| Soil processes | | //humusc1, humusc2, humusc3, fastc1, fastc2, fastc3, occonc0, koflim// | [[start:hype_file_reference:par.txt|par.txt]] |
| ::: | theta_low | //ocsoimslp///100 or 0.08 | ::: |
| ::: | //satact// | //ocsoimsat// or 0.6 | ::: |
| ::: | //minc, klh, klo, kof, kho// | //minc, klh, klo, kof, kho// | ::: |
| ::: | //par// | //kcgwreg// or //koc// | ::: |
| ::: | //ocfldelx,ocfldele// | //ocfldelx,ocfldele// | ::: |
| Riparian zone | //elev// | calculated from //mean_elev// and //dhslc_nn// | [[start:hype_file_reference:geodata.txt|GeoData.txt]] |
| ::: | //ripz//, //ripe//, //rips// | //ripz//, //ripe//, //rips// | [[start:hype_file_reference:par.txt|par.txt]] |
| ::: | //wp// | calculated from //wcwp//, //wcwp1//, //wcwp2//, //wcwp3// | ::: |
| ::: | //pw// | calculated from //wcwp//, //wcwp1-wcwp3//, //wcfc//, //wcfc1-wcfc3//, //wcep//, //wcep1-wcep3// | ::: |
| ::: | //d// | | [[start:hype_file_reference:geoclass.txt|GeoClass.txt]] |
| Rivers and lakes | //area//, //lakeregion// | | [[start:hype_file_reference:geodata.txt|GeoData.txt]] |
| ::: | //tocmean//, //wprodc//, //sedoc// | //tocmean//, //wprodc//, //sedoc// | [[start:hype_file_reference:par.txt|par.txt]] or [[start:hype_file_reference:lakedata.txt|LakeData.txt]] |
| ::: | //limsedpp// | //limsedpp// | [[start:hype_file_reference:par.txt|par.txt]] |
| ::: | | //tpmean// | ::: |
| ::: | //halfsatTPwater// | //hsatTP// | ::: |