The SimComponent SlimRoots estimates the daily increase in the biomass of seminal and lateral roots in a variable number of soil layers based and converts it into root length per layer. The calculations depend on the input of assimilates from the shoot which must be provided by another SimComponent. Before emergency, assimilates can also be supplied by the seeds (SimComponent SeedstoSprout). The assimilates provided by the shoot or the seeds are in the first placed used for the growth of seminals which determine the vertical penetration into the soil. In addition to the supply of assimilates by the shoot, vertical penetration of seminals depends on soil temperature, soil dryness and soil strength which reduce the crop specific maximum daily elongation rate. The remaining assimilates are used for the development of lateral roots.
The growth of seminals starts in the soil layer where seeds have been placed which is determined by the sowing depth (SEEDD in m). The assimilates supplied by the shoot (RWRT in g m-2) or the seeds are in the first place used for the growth of seminals. When daily mean soil temperature is above 0 °C, it is assumed that the daily elongation of seminals (in m) is proportional to the mean soil temperature using the RTFAC coefficient of 0.12 m K-1. Below a soil temperature of 0 °C no seminal growth occurs. The temperature dependent daily elongation is then converted into assimilate demand (ASROOT) by applying the specific root weight of seminals (WSROOT in g m-1), which is defined by the user in a separate input file for each crop (SlimFile or SlimProperties).The daily demand of assimilates is limited to the crop-specific maximum elongation rate per day (MaximalSeminalRootLengthPerDay), which can be defined by the user in the same input file. Finally, the demand for assimilates can be further reduced by the effect of soil strength. To account for the soil strength two approaches have been implemented in SlimRoots to calculate a BulkDensitySeminalRootReductionFactor in each soil layer:
The user can chose either of the two methods by specifying the required parameters in the input file SlimProperties or switching them off by supplying negative values for these parameters. The BulkDensitySeminalRootReductionFactor reduces the demand for assimilates by the seminals when soil bulk density increases thus restricting the daily vertical root elongation rate. However, the presence of macropores in a soil layer can counteract the effect of a compacted soil layer. Therefore, SlimRoots offers a routine to account for the effect of soil macropores on the BulkDensitySeminalRootReductionFactor as proposed by Gaiser et al. (2013).
If the demand for assimilates by the seminal roots (ASROOT) is lower than the supply by the shoot (RWRT) the remaining assimilates (ALROOT) are used to produce lateral roots. ALROOT is distributed stepwise to the soil layers which contain seminal roots from the top layer to the deepest soil layer with seminals present. For a soil layer thickness of 5 cm, the first soil layer receives 16.33% of ALROOT. Then the remaining assimilates for lateral roots (REMLT) are recalculated and the subsequent soil layer receives again the same portion (PLRASS) of 16.33% until all assimilates are consumed (REMLT =0). PLRASS`can be reduced when soil layers are dry by the dryness factor `DRYFAC. The dryness factor DRYFAC depends on the soil moisture in each soil layer and is provided by any soil water balance SimComponent. After assimilates have been distributed to seminals and larterals in each soil layer, the root length of seminals and laterals are calculated based on the specific root weight of laterals (WLROOT) and seminals (WSROOT). Examples of specific weight of lateral roots for wheat are given in Jamieson & Ewert (1999) and Van Noordwijk and Brouwer (2004) (0.006134 to 0.004065041 g/m). The default values are presently set to 0.015 g m-1 for WSROOT and 0.0035 g m-1 for WLROOT which can be used for wheat roots (Jamieson & Ewert, 1999). Finally, the root length density in each layer i (RootLengthDensityPerLayer(i)) can be calculated.
Root length density and root age are important factor for the magnitude of crop water uptake. Therefore, SlimRoots calculates a root related restriction factor for soil water uptake (FRR) as
\[ \begin{eqnarray} FRR(i) &=& 1-(e^{-0.3 \cdot RDA4AGE(i)}) & & for & & FRR <=1 \end{eqnarray} \]where RDA4AGE is depending on the age of the roots and the root length density in layer i with
\[ \begin{eqnarray} RDA4AGE(i) &=& MAX(RD, RD \cdot \frac{YOUTH}{RLAGE(i)}) \end{eqnarray} \]RD is the root length density (cm cm-3) and RLAGE is the age of the roots in days in layer i. YOUTH is the period in days when roots are considered to be active (the default for YOUTH is set to 50 days).
Root start to decay after a user-defined developments stage (cDeadRootsDevStage) has occurred. In cereal crops, this is usually after anthesis (DVS=1.0). However, drought stress can delay the decay. The accumulated assimilates in the laterals in each layer (ASSLR(i)) are multiplied with a user-defined daily death rate (DailyDeadRootsRate) and then the root length density of the decayed roots is subtracted from the total root length density in each layer.
References:
Content Type | Name | Description | Data Type | Unit | Min Value | Max Value | Default Value |
---|---|---|---|---|---|---|---|
constant | cBiomassLossRateAtPartitioning | Rate of root biomass, that is not distributed to lateral roots | DOUBLE | 1 | - | - | 0.005 |
constant | cBulkDensity | Bulk density of soil layer i | DOUBLEARRAY | g/cm3 | 0.0 | 20.0 | - |
constant | cBulkDensityThresholdLower | Below this soil bulk density the reduction factor for the elongation of seminal roots is 1 (no impedance for root elongation) | DOUBLE | g/cm3 | 0.0 | 20.0 | 0.05 |
constant | cBulkDensityThresholdMiddle | Between this soil bulk density and the upper threshold the reduction factor for the elongation of seminal roots decreases linearly (increasing impedance for root elongation) | DOUBLE | g/cm3 | 0.0 | 20.0 | 0.05 |
constant | cBulkDensityThresholdUpper | Above this soil bulk density the reduction factor for the elongation of seminal roots is at its minimum (maximum impedance for root elongation) | DOUBLE | g/cm3 | 0.0 | 20.0 | 0.05 |
constant | cDailyDeadRootsRate | Daily fraction of dying roots | DOUBLE | 1 | 0.001 | 20.0 | 0.02 |
constant | cDeadRootsDevStage | Crop development stage in DVS when roots start to die | DOUBLE | 1 | 0.001 | 20.0 | 0.02 |
constant | cMacroporeRootGrowthFactor | Above this volumetric fraction of macropores in the soil seminal root elongation is not restricted by soil bulk density | DOUBLE | m3/m3 | - | - | 0.0 |
constant | cMacroporeVolume | Volumetric fraction of macropore volume (pores > 2mm) in layer i | DOUBLEARRAY | m3/m3 | 0.0 | 100.0 | - |
constant | cMaximalSeminalRootLengthPerDay | Crop-specific maximum root elongation rate per day | DOUBLE | m | 0.0 | 20.0 | 0.033 |
constant | cMinimalBDSeminalRootReductionFactor | Bulk density reduction factor for the elongation of seminal roots when soil bulk density is above the upper threshold (Set to -1 to suppress calculation of the reduction factor by method 1) | DOUBLE | 1 | -1.0 | 1.0 | -1.0 |
constant | cNSRPP | number of seminal roots per m2 | DOUBLE | 1 | - | - | 5.0 |
constant | cRINPOP | number of emerged plants per m2 | DOUBLE | 1 | - | - | 200.0 |
constant | cRTFAC | temperature factor driving penetration of seminal roots | DOUBLE | 1 | - | - | 0.12 |
constant | cSEEDD | Soil depth where seeds have been placed | DOUBLE | m | 0.001 | 20.0 | 0.02 |
constant | cSandPercentage | Proportion of sand fraction in the fine earth in soil layer i | DOUBLEARRAY | g/hg | 0.0 | 100.0 | - |
constant | cSoilStrengthParameter | Parameter to increase sensitivity of the calculation of the bulk density reduction factor for seminar root elongation in Method 2 (Between 1.2 and 2. Set to -1 to suppress calculation of the reduction factor by method 2) | DOUBLE | 1 | -1.0 | 2.0 | -1.0 |
constant | cWLROOT | specific weight of lateral roots (in g m-1) | DOUBLE | g/m | - | - | 0.0035 |
constant | cWSROOT | specific weight of seminal roots in g m-1 | DOUBLE | g/m | - | - | 0.015 |
input | cSoilLayerDepth | Depth of the bottom of each soil layer after soil layer transformation | DOUBLEARRAY | m | - | - | - |
input | iDevStage | Crop development stage (DVS) | DOUBLE | 1 | 0.0 | 3.0 | 0.0 |
input | iDryFactor | Dryness factor in each soil layer (input from soil water balance SimComponent) | DOUBLEARRAY | 1 | 0.0 | 1.0 | - |
input | iIsHarvest | Flag to indicate day of harvest | BOOLEAN | 1 | - | - | false |
input | iRWRT | Daily rate of biomass partitioned to the roots | DOUBLE | g/m2 | 0.0 | 20.0 | 0.0 |
input | iSoilTemperature | Daily mean temperature in soil layer i | DOUBLEARRAY | °C | -50.0 | 50.0 | - |
input | iTMAX | Daily maximum air temperature (used if soil temperature is not available) | DOUBLE | °C | -50.0 | 50.0 | - |
input | iTMIN | Daily minimum air temperature (used if soil temperature is not available) | DOUBLE | °C | -50.0 | 50.0 | - |
state | sLateralRootLength | Length of lateral roots in layer i (in m) | DOUBLEARRAY | m | 0.0 | 500.0 | - |
state | sRootAge | Age of seminal roots in layer i in days | DOUBLEARRAY | d | 0.0 | 500.0 | - |
state | sSEEDDAdjusted | Depth of sowing adjusted to layer thickness (Sowing must always occur in the first soil layer) | DOUBLE | m | 0.0 | 20.0 | 0.02 |
state | sSeminalRootDepth | Actual depth of seminal roots (deepest soil layer with roots) | DOUBLE | m | 0.0 | 50.0 | 0.0 |
state | sSeminalRootLength | Length if seminal roots in layer i (in m) | DOUBLEARRAY | m | 0.0 | 500.0 | - |
rate | rSRERR | Daily increment in seminal root depth [m] | DOUBLE | m/d | 0.0 | 2.0 | 0.0 |
out | BiomassBalance | Daily mass balance of biomass partitioned to the roots | DOUBLE | g/m2 | -100000.0 | 10000.0 | 0.0 |
out | DeadRootRatePerLayer | Daily increase of dead root biomass per layer | DOUBLEARRAY | g/m2 | 0.0 | 2000.0 | - |
out | MD95 | Actual depth of seminal roots (depth of soil with crop roots) expressed as number of layers | INT | 1 | 0 | 1000 | 0 |
out | RelativeRootLengthPerLayer | Relative rootlength in layer i (sum of lat.roots and sem.roots), in relation to the total rootlength in all layers | DOUBLEARRAY | m | 0.0 | 20.0 | - |
out | RootLengthDensityPerLayer | Rootlength per soil volume in layer i (sum of lat.roots and sem.roots divided by soil volume) | DOUBLEARRAY | centimetre_per_cubic_centimetre | 0.0 | 20.0 | - |
out | RootRestrictionFactor | Factor restricting uptake of water and mineral nitrogen in soil layer i due to root density and root age | DOUBLEARRAY | 1 | 0.0 | 1.0 | - |
out | THICKL | Layer thickness in m | DOUBLE | m | 0.001 | 20.0 | 0.0 |
out | TotalRootBiomass | Total Root Biomass over the soil profile | DOUBLE | g/m2 | -10000.0 | 10000.0 | 0.0 |
out | TotalRootBiomassPerLayer | Total root biomass per soil layer | DOUBLEARRAY | g/m2 | 0.0 | 200.0 | - |
out | TotalRootLengthPerLayer | Total root length per layer | DOUBLEARRAY | m | 0.0 | 20.0 | - |
out | UnusedRootBiomass | Daily remnant of biomass partitioned to the roots which can not be converted into root growth due to root growth restrictions (is returned to shoot on day i+1) | DOUBLE | g/m2 | 0.0 | 200.0 | 0.0 |