ind_gw_recharge: Calculate groundwater recharge of a soil

In OBIC: Calculate the Open Bodem Index (OBI) Score

Calculate groundwater recharge of a soil

Description

This function calculates an index score for groundwater storage based on precipitation surplus, infiltration at saturation, sealing risk, drainage and subsoil compaction

Usage

ind_gw_recharge(
  B_LU_BRP,
  D_PSP,
  D_WRI_K,
  I_P_SE = NULL,
  I_P_CO = NULL,
  B_DRAIN,
  B_GWL_CLASS,
  D_SE = NULL,
  B_SC_WENR = NULL,
  D_P_CO = NULL
)

Arguments

B_LU_BRP	(numeric) The crop code from the BRP
D_PSP	(numeric) The precipitation surplus per crop calculated by calc_psp
D_WRI_K	(numeric) The value for top soil permeability (cm/d) as calculated by calc_permeability
I_P_SE	(numeric) The indicator value for soil sealing (deprecated)
I_P_CO	(numeric) The indicator value for occurrence of subsoil compaction (deprecated)
B_DRAIN	(boolean) Are drains installed to drain the field (options: yes or no)
B_GWL_CLASS	(character) The groundwater table class
D_SE	(numeric) The value of soil sealing calculated by calc_sealing_risk
B_SC_WENR	(character) The risk for subsoil compaction as derived from risk assessment study of Van den Akker (2006)
D_P_CO	(numeric) Compaction value based on visual soil assessment (see Details)

Details

Soil compaction risk can be assessed either with the soilcompaction risk map (B_SC_WENR) or visual soil assessment (VSA) (D_P_CO). D_P_CO itself is calculated from the VSA scoring of earthworms, compaction, rooting depth, and puddling: D_P_CO = (3 x A_EW_BCS + 3 x A_SC_BCS + 3 x A_RD_BCS - 2 x A_P_BCS - A_RT_BCS)/18 where the .*_BCS variables have the values 0, 1, or 2. When both B_SC_WENR and D_P_CO are provided, D_P_CO is used as it is considered more accurate determination of compaction.

Value

The evaluated score for the soil function to improve groundwater recharge. A numeric value between 0 and 1.

Examples

ind_gw_recharge(B_LU_BRP = 265, D_PSP = 200, D_WRI_K = 10, B_DRAIN = FALSE,
 B_GWL_CLASS = 'V', D_SE = 2, B_SC_WENR = 'Groot')
ind_gw_recharge(B_LU_BRP = 233, D_PSP = 400, D_WRI_K = 10, B_DRAIN = TRUE,
 B_GWL_CLASS = 'II', D_SE = 5, D_P_CO = 0.2)

OBIC documentation built on Feb. 11, 2026, 5:06 p.m.