Residue Decomposition & Mineralization

Decomposition and mineralization of the fresh organic nitrogen pool is allowed only in the first soil layer. Decomposition and mineralization are controlled by a decay rate constant that is updated daily. The decay rate constant is calculated as a function of the C:N ratio and C:P ratio of the residue, temperature and soil water content.

The C:N ratio of the residue is calculated:

εC:N=0.58rsdlyorgNfrsh,ly+NO3ly\varepsilon _{C:N}=\frac{0.58*rsd_{ly}}{orgN_{frsh,ly}+NO3_{ly}} 3:1.2.5

where εC:N\varepsilon_{C:N} is the C:N ratio of the residue in the soil layer, rsdlyrsd_{ly} is the residue in layer lyly (kg/ha), 0.58 is the fraction of residue that is carbon, orgNfrsh,lyorgN_{frsh,ly} is the nitrogen in the fresh organic pool in layer lyly (kg N/ha), and NO3lyNO3_{ly} is the amount of nitrate in layer lyly (kg N/ha).

The C:P ratio of the residue is calculated:

εC:P=0.58rsdlyorgPfrsh,ly+Psolution,ly\varepsilon_{C:P}=\frac{0.58*rsd_{ly}}{orgP_{frsh,ly}+P_{solution,ly}} 3:1.2.6

where εC:P\varepsilon_{C:P} is the C:P ratio of the residue in the soil layer, rsdlyrsd_{ly} is the residue in layer lyly (kg/ha), 0.58 is the fraction of residue that is carbon, orgPfrsh,lyorgP_{frsh,ly} is the phosphorus in the fresh organic pool in layer lyly (kg P/ha), and Psolution,lyP_{solution,ly} is the amount of phosphorus in solution in layer lyly (kg P/ha).

The decay rate constant defines the fraction of residue that is decomposed. The decay rate constant is calculated:

δntr,ly=βrsdγntr,ly(γtmp,lyγsw,ly)1/2\delta_{ntr,ly}=\beta_{rsd}*\gamma_{ntr,ly}*(\gamma_{tmp,ly}*\gamma_{sw,ly})^{1/2} 3:1.2.7

where δntr,ly\delta_{ntr,ly} is the residue decay rate constant, βrsd\beta_{rsd} is the rate coefficient for mineralization of the residue fresh organic nutrients, γntr,ly\gamma_{ntr,ly} is the nutrient cycling residue composition factor for layer lyly, γtmp,ly\gamma_{tmp,ly} is the nutrient cycling temperature factor for layer lyly, and γsw,ly\gamma_{sw,ly} is the nutrient cycling water factor for layer lyly.

The nutrient cycling residue composition factor is calculated:

γntr,ly=min[exp[0.693(εC:N25)25],exp[0.693(εC:P200)200],1.0]\gamma_{ntr,ly}=min[{exp[-0.693*\frac{(\varepsilon_{C:N}-25)}{25}] ,exp[-0.693*\frac{(\varepsilon_{C:P}-200)}{200}] , 1.0}]

3:1.2.8

where γntr,ly\gamma_{ntr,ly} is the nutrient cycling residue composition factor for layer lyly, εC:N\varepsilon_{C:N} is the C:N ratio on the residue in the soil layer, and εC:P\varepsilon_{C:P} is the C:P ratio on the residue in the soil layer.

Mineralization from the residue fresh organic N pool is then calculated:

Nminf,ly=0.8δntr,lyorgNfrsh,lyN_{minf,ly}=0.8*\delta_{ntr,ly}*orgN_{frsh,ly} 3:1.2.9

where Nminf,lyN_{minf,ly} is the nitrogen mineralized from the fresh organic N pool (kg N/ha), δntr,ly\delta_{ntr,ly} is the residue decay rate constant, and orgNfrsh,lyorgN_{frsh,ly} is the nitrogen in the fresh organic pool in layer lyly (kg N/ha). Nitrogen mineralized from the fresh organic pool is added to the nitrate pool in the layer.

Decomposition from the residue fresh organic N pool is calculated:

Ndec,ly=0.2δntr,lyorgNfrsh,lyN_{dec,ly}=0.2*\delta_{ntr,ly}*orgN_{frsh,ly} 3:1.2.9

where Ndec,lyN_{dec,ly} is the nitrogen decomposed from the fresh organic N pool (kg N/ha), δntr,ly\delta_{ntr,ly} is the residue decay rate constant, and orgNfrsh,lyorgN_{frsh,ly} is the nitrogen in the fresh organic pool in layer lyly (kg N/ha). Nitrogen decomposed from the fresh organic pool is added to the humus active organic pool in the layer.

Table 3:1-2: SWAT+ input variables that pertain to mineralization.

Variable NameDefinitionFile Name

CMN

βmin\beta_{min}: Rate coefficient for mineralization of the humus active organic nutrients

.bsn

RSDCO

βrsd\beta_{rsd}: Rate coefficient for mineralization of the residue fresh organic nutrients

.bsn

RSDCO_PL

βrsd\beta_{rsd}: Rate coefficient for mineralization of the residue fresh organic nutrients

crop.dat

Last updated