Skip to contents

Function to calculate and evaluate the N leaching risks in soils in view of water purification for groundwater quality

Usage

bln_wat_nrisk_gw(
  ID,
  B_LU_BRP,
  B_SOILTYPE_AGR,
  B_AER_CBS,
  B_GWL_CLASS,
  B_SC_WENR,
  B_FERT_NORM_FR = 1,
  A_CLAY_MI,
  A_SAND_MI,
  A_SILT_MI,
  A_SOM_LOI,
  A_P_AL,
  A_P_WA,
  A_P_CC,
  A_PH_CC,
  A_CEC_CO,
  A_K_CO_PO,
  A_K_CC,
  M_GREEN
)

Arguments

ID

(character) A field id

B_LU_BRP

(numeric) The crop code

B_SOILTYPE_AGR

(character) The agricultural type of soil

B_AER_CBS

(character) The agricultural economic region in the Netherlands (CBS, 2016)

B_GWL_CLASS

(character) The groundwater table class

B_SC_WENR

(character) The risk for subsoil compaction as derived from risk assessment study of Van den Akker (2006).

B_FERT_NORM_FR

(numeric) The fraction of the application norm utilized

A_CLAY_MI

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

A_SAND_MI

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

A_SILT_MI

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

A_SOM_LOI

(numeric) The organic matter content of the soil (%)

A_P_AL

(numeric) The P-AL content of the soil

A_P_WA

(numeric) The P-content of the soil extracted with water (mg P2O5 / 100 ml soil)

A_P_CC

(numeric) The plant available P content, extracted with 0.01M CaCl2 (mg / kg)

A_PH_CC

(numeric) The acidity of the soil, measured in 0.01M CaCl2 (-)

A_CEC_CO

(numeric) The cation exchange capacity of the soil (mmol+ / kg), analyzed via Cobalt-hexamine extraction

A_K_CO_PO

(numeric) The occupation of the CEC with K (%)

A_K_CC

(numeric) The plant available K content, extracted with 0.01M CaCl2 (mg / kg)

M_GREEN

(boolean) A soil measure. Are catch crops sown after main crop (optional, option: yes or no)

Examples

bln_wat_nrisk_gw(
ID = 15,
B_LU_BRP = c(233,259,2014,308),
B_SOILTYPE_AGR = rep('rivierklei',4),
B_GWL_CLASS = rep('GtIII',4),
B_AER_CBS = rep('LG06',4),
B_SC_WENR = rep(902,4),
B_FERT_NORM_FR = rep(1,4),
A_SOM_LOI = rep(3,4),
A_CLAY_MI = rep(20,4),
A_SAND_MI = rep(20,4),
A_SILT_MI = rep(60,4),
A_P_AL = c(25, 20,10,15),
A_P_CC = c(1.5, 3,4,6),
A_P_WA = c(24,20,15,10),
A_PH_CC = c(5.6,7.04,5.65,6.22),
A_CEC_CO = c(122,191,157,142),
A_K_CC = c(95,107,73,89),
A_K_CO_PO = c(1.6,2.9,3.1,2.4),
M_GREEN = FALSE
)
#> [1] 0.9866600 0.9743001 0.9554349 0.9922947