OBIC (version 4.2.1)

calc_workability: Calculate indicator for workability

Description

This function calculates the workability of soils, given as a value of relative season length between 0 and 1. A relative season length of 1 indicates that the water table is sufficiently low for the soil to be workable for the entire growing season required by the crop. The required ground water table for workability is determined by soil type and soil properties. Hydrological variables determine the groundwater table for each day of the year. The option calcyieldloss allows for calculation of yield loss based on the relative season length, differentiating in yield loss between six groups of crops Based on Huinink (2018)

Usage

calc_workability(
  A_CLAY_MI,
  A_SILT_MI,
  B_LU_BRP,
  B_SOILTYPE_AGR,
  B_GWL_GLG,
  B_GWL_GHG,
  B_GWL_ZCRIT,
  calcyieldloss = FALSE
)

Value

The workability of a soil, expressed as a numeric value representing the relative season length that the soil can be managed by agricultural activities.

Arguments

A_CLAY_MI

(numeric) The clay content of the soil (%)

A_SILT_MI

(numeric) The silt content of the soil (%)

B_LU_BRP

(numeric) The crop code from the BRP

B_SOILTYPE_AGR

(character) The agricultural type of soil

B_GWL_GLG

(numeric) The lowest groundwater level averaged over the most dry periods in 8 years in cm below ground level

B_GWL_GHG

(numeric) The highest groundwater level averaged over the most wet periods in 8 years in cm below ground level

B_GWL_ZCRIT

(numeric) The distance between ground level and groundwater level at which the groundwater can supply the soil surface with 2mm water per day (in cm)

calcyieldloss

(boolean) whether the function includes yield loss, options: TRUE or FALSE (default).

References

Huinink (2018) Bodem/perceel geschiktheidsbeoordeling voor Landbouw, Bosbouw en Recreatie. BodemConsult-Arnhem

Examples

calc_workability(A_CLAY_MI = 18,A_SILT_MI = 25,B_LU_BRP = 265,
B_SOILTYPE_AGR = 'dekzand',B_GWL_GLG = 145,B_GWL_GHG = 85,B_GWL_ZCRIT = 400,
calcyieldloss = FALSE)
calc_workability(18,25,265,'dekzand',145,85,400,FALSE)