Residue Decomposition Submodel

Model Overview

This submodel simulates the decrease in crop residue biomass due to microbial activity. The decomposition process is modeled as a first-order reaction with temperature and moisture as driving variables. Decomposition is a function of decomposition days. Under optimum temperature and moisture conditions one decomposition day per day is accumulated. Only a fraction of a decomposition day is accumulated if conditions are less than optimum. Biomass remaining after harvest is partitioned between standing, flat, buried, and root pools. Below ground biomass decomposition is calculated for each soil layer.

Residue from different crops may decompose at different rates. Since residue decomposition can require a long period of time, crop residue biomass from sequential harvests is accounted for in three separate pools. Biomass from the most recently harvested crop will be in pool one, biomass from the penultimate crop in pool two, and there is a third pool for biomass from the oldest crop(s). After harvest, any residue biomass remaining from a previous crop is moved into the older age pools and residue from the just harvested crop is moved into the first residue decomposition pool.

Standing residue losses not only result from microbial activity, but also from physical forces. Physical transfer of crop residue from the standing biomass pool will reduce both the stem population and standing biomass. A daily estimate of the standing stem population is required in order to evaluate the vertical stem area that the wind encounters. This area is quantified by the stem area index, which is calculated from standing stem number, stem height and stem diameter. It affects aerodynamic resistance and, ultimately, wind erosion. Stems start to fall over after a threshold of cumulative decomposition days after harvest has been reached. Stem area index decreases proportionally with standing stem number.

Both standing and flat crop residue provide cover to the soil surface, protecting it against wind erosion. Soil cover from standing residue is estimated from stem number and stem diameter. Soil cover from flat residue is predicted from flat residue mass. Tillage may alter the amount of residue in the different pools.

For more detail on the science behind this submodel please see the WEPS technical documentation.