Fix for UK drainage, some EFSA interception and washoff

- PEC_sw_drainage_UK() gave results that were a little bit too high
for the substances with the highest Koc (>4000), as I used 0.01% instead
of 0.008% for them. This must have come from an old vesion of the UK
data requirements handbook, at least I do not have another explanation

- Add EFSA interception (from 2014 DegT50 guidance) and tier 1
crop wash-off factors (from 2017 PEC soil guidance) for some major
arable crops

- Update docs

4 files changed, 104 insertions, 4 deletions

diff --git a/R/EFSA_GW_interception_2014.R b/R/EFSA_GW_interception_2014.R
new file mode 100644
index 0000000..15d7835
--- /dev/null
+++ b/R/EFSA_GW_interception_2014.R
@@ -0,0 +1,41 @@
+#' Subset of EFSA crop interception default values for groundwater modelling
+#'
+#' @name EFSA_GW_interception_2014
+#' @docType data
+#' @source European Food Safety Authority (2014) EFSA Guidance Document for
+#' evaluating laboratory and field dissipation studies to obtain DegT50 values
+#' of active substances of plant protection products and transformation
+#' products of these active substances in soil. \emph{EFSA Journal}
+#' \bold{12}(5):3662, 37 pp., doi:10.2903/j.efsa.2014.3662
+#' @format A matrix containing interception values, currently only for some selected crops
+#' @examples
+#' \dontrun{
+#'   # This is the code that was used to define the data
+#'   bbch <- paste0(0:9, "x")
+#'   crops <- c(
+#'     "Beans (field + vegetable)",
+#'     "Peas",
+#'     "Summer oilseed rape", "Winter oilseed rape",
+#'     "Tomatoes",
+#'     "Spring cereals", "Winter cereals")
+#'   EFSA_GW_interception_2014 <- matrix(NA, length(crops), length(bbch),
+#'     dimnames = list(Crop = crops, BBCH = bbch))
+#'   EFSA_GW_interception_2014["Beans (field + vegetable)", ] <-
+#'     c(0, 0.25, rep(0.4, 2), rep(0.7, 5), 0.8)
+#'   EFSA_GW_interception_2014["Peas", ] <-
+#'     c(0, 0.35, rep(0.55, 2), rep(0.85, 5), 0.85)
+#'   EFSA_GW_interception_2014["Summer oilseed rape", ] <-
+#'     c(0, 0.4, rep(0.8, 2), rep(0.8, 5), 0.9)
+#'   EFSA_GW_interception_2014["Winter oilseed rape", ] <-
+#'     c(0, 0.4, rep(0.8, 2), rep(0.8, 5), 0.9)
+#'   EFSA_GW_interception_2014["Tomatoes", ] <-
+#'     c(0, 0.5, rep(0.7, 2), rep(0.8, 5), 0.5)
+#'   EFSA_GW_interception_2014["Spring cereals", ] <-
+#'     c(0, 0, 0.2, 0.8, rep(0.9, 3), rep(0.8, 2), 0.8)
+#'   EFSA_GW_interception_2014["Winter cereals", ] <-
+#'     c(0, 0, 0.2, 0.8, rep(0.9, 3), rep(0.8, 2), 0.8)
+#'   save(EFSA_GW_interception_2014,
+#'     file = "../data/EFSA_GW_interception_2014.RData")
+#' }
+#' EFSA_GW_interception_2014
+NULL
diff --git a/R/EFSA_washoff_2017.R b/R/EFSA_washoff_2017.R
new file mode 100644
index 0000000..450c12e
--- /dev/null
+++ b/R/EFSA_washoff_2017.R
@@ -0,0 +1,42 @@
+#' Subset of EFSA crop washoff default values
+#'
+#' @name EFSA_washoff_2017
+#' @docType data
+#' @source
+#'   European Food Safety Authority (2017) EFSA guidance document for
+#'   predicting environmental concentrations of active substances of plant
+#'   protection products and transformation products of these active substances
+#'   in soil. \emph{EFSA Journal} \bold{15}(10) 4982
+#'   doi:10.2903/j.efsa.2017.4982
+#' @format A matrix containing wash-off factors, currently only for some selected crops
+#' @examples
+#' \dontrun{
+#'   # This is the code that was used to define the data
+#'   bbch <- paste0(0:9, "x")
+#'   crops <- c(
+#'     "Beans (field + vegetable)",
+#'     "Peas",
+#'     "Summer oilseed rape", "Winter oilseed rape",
+#'     "Tomatoes",
+#'     "Spring cereals", "Winter cereals")
+#'   EFSA_washoff_2017 <- matrix(NA, length(crops), length(bbch),
+#'     dimnames = list(Crop = crops, BBCH = bbch))
+#'   EFSA_washoff_2017["Beans (field + vegetable)", ] <-
+#'     c(NA, 0.6, rep(0.75, 2), rep(0.8, 5), 0.35)
+#'   EFSA_washoff_2017["Peas", ] <-
+#'     c(NA, 0.4, rep(0.6, 2), rep(0.65, 5), 0.35)
+#'   EFSA_washoff_2017["Summer oilseed rape", ] <-
+#'     c(NA, 0.4, rep(0.5, 2), rep(0.6, 5), 0.5)
+#'   EFSA_washoff_2017["Winter oilseed rape", ] <-
+#'     c(NA, 0.1, rep(0.4, 2), rep(0.55, 5), 0.3)
+#'   EFSA_washoff_2017["Tomatoes", ] <-
+#'     c(NA, 0.55, rep(0.75, 2), rep(0.7, 5), 0.35)
+#'   EFSA_washoff_2017["Spring cereals", ] <-
+#'     c(NA, 0.4, 0.5, 0.5, rep(0.65, 3), rep(0.65, 2), 0.55)
+#'   EFSA_washoff_2017["Winter cereals", ] <-
+#'     c(NA, 0.1, 0.4, 0.6, rep(0.55, 3), rep(0.6, 2), 0.4)
+#'   save(EFSA_washoff_2017,
+#'     file = "../data/EFSA_washoff_2017.RData")
+#' }
+#' EFSA_washoff_2017
+NULL
diff --git a/R/PEC_sw_drainage_UK.R b/R/PEC_sw_drainage_UK.R
index e53f179..3be4c97 100644
--- a/R/PEC_sw_drainage_UK.R
+++ b/R/PEC_sw_drainage_UK.R
@@ -29,6 +29,14 @@
 #'   "DFOP", "HS", or "IORE", or an mkinmod object
 #' @param model_parms A named numeric vector containing the model parameters
 #' @return The predicted concentration in surface water in µg/L
+#' @references HSE's Chemicals Regulation Division (CRD) Active substance
+#'   PECsw calculations (for UK specific authorisation requests)
+#'   \url{https://www.hse.gov.uk/pesticides/topics/pesticide-approvals/pesticides-registration/data-requirements-handbook/fate/active-substance-uk.htm}
+#'   accessed 2019-09-27
+#'
+#'   Drainage PECs Version 1.0 (2015) Spreadsheet published at
+#'   \url{https://www.hse.gov.uk/pesticides/topics/pesticide-approvals/pesticides-registration/data-requirements-handbook/fate/pec-tools-2015/PEC\%20sw-sed\%20(drainage).xlsx}
+#'   accessed 2019-09-27
 #' @export
 #' @author Johannes Ranke
 #' @examples
diff --git a/R/SSLRC_mobility_classification.R b/R/SSLRC_mobility_classification.R
index 2f8fb34..602f166 100644
--- a/R/SSLRC_mobility_classification.R
+++ b/R/SSLRC_mobility_classification.R
@@ -1,4 +1,4 @@
-# Copyright (C) 2015,2018  Johannes Ranke
+# Copyright (C) 2015,2018,2019  Johannes Ranke
 # Contact: jranke@uni-bremen.de
 # This file is part of the R package pfm
 
@@ -17,30 +17,39 @@
 
 #' Determine the SSLRC mobility classification for a chemical substance from its Koc
 #'
-#' This implements the method specified in the UK data requirements handbook and was 
+#' This implements the method specified in the UK data requirements handbook and was
 #' checked against the spreadsheet published on the CRC website
 #'
 #' @param Koc The sorption coefficient normalised to organic carbon in L/kg
 #' @return A list containing the classification and the percentage of the
 #'   compound transported per 10 mm drain water
+#' @references HSE's Chemicals Regulation Division (CRD) Active substance
+#'   PECsw calculations (for UK specific authorisation requests)
+#'   \url{https://www.hse.gov.uk/pesticides/topics/pesticide-approvals/pesticides-registration/data-requirements-handbook/fate/active-substance-uk.htm}
+#'   accessed 2019-09-27
+#'
+#'   Drainage PECs Version 1.0 (2015) Spreadsheet published at
+#'   \url{https://www.hse.gov.uk/pesticides/topics/pesticide-approvals/pesticides-registration/data-requirements-handbook/fate/pec-tools-2015/PEC\%20sw-sed\%20(drainage).xlsx}
+#'   accessed 2019-09-27
 #' @export
 #' @author Johannes Ranke
 #' @examples
 #' SSLRC_mobility_classification(100)
+#' SSLRC_mobility_classification(10000)
 SSLRC_mobility_classification <- function(Koc)
 {
   if (!is.numeric(Koc) | length(Koc) != 1) stop("Please give a single number")
   if (is.na(Koc)) {
     result <- list(NA, NA)
   } else {
-    result <- list("Non mobile", 0.01)
+    result <- list("Non mobile", 0.008)
     if (Koc < 4000) result <- list("Slightly mobile", 0.02)
     if (Koc < 1000) result <- list("Slightly mobile", 0.5)
     if (Koc < 500) result <- list("Moderately mobile", 0.7)
     if (Koc < 75) result <- list("Mobile", 1.9)
     if (Koc < 15) result <- list("Very mobile", 1.9)
   }
-  names(result) <- c("Mobility classification", 
+  names(result) <- c("Mobility classification",
                      "Percentage drained per mm of drain water")
   return(result)
 }