From 904a86fa92f1577414f2fce2bdbfe21e12381106 Mon Sep 17 00:00:00 2001 From: Ranke Johannes Date: Thu, 1 Feb 2024 10:39:53 +0100 Subject: Calculate drift percentages and PECsw using the Rautmann formula --- R/PEC_sw_drift.R | 66 +++++++++++++++++++++++++++++++++++----------- R/drift_parameters_focus.R | 6 +++++ 2 files changed, 57 insertions(+), 15 deletions(-) (limited to 'R') diff --git a/R/PEC_sw_drift.R b/R/PEC_sw_drift.R index 2dfd651..1d22e4e 100644 --- a/R/PEC_sw_drift.R +++ b/R/PEC_sw_drift.R @@ -4,15 +4,15 @@ #' concentration in surface water based on complete, instantaneous mixing #' with input via spray drift. #' +#' @inheritParams drift_percentages_rautmann_formula #' @param rate Application rate in units specified below -#' @param applications Number of applications for selection of drift percentile #' @param drift_percentages Percentage drift values for which to calculate PECsw. #' 'drift_data' and 'distances' if not NULL. #' @param drift_data Source of drift percentage data. If 'JKI', the [drift_data_JKI] #' included in the package is used. If 'RF', the Rautmann formula is used, if #' implemented for the crop type and number of applications -#' @param crop Crop name (use German names for JKI data), defaults to "Ackerbau" -#' @param distances The distances in m for which to get PEC values +#' @param crop_group_JKI When using the 'JKI' drift data, one of the German names +#' as used in [drift_parameters_focus]. #' @param water_depth Depth of the water body in cm #' @param rate_units Defaults to g/ha #' @param PEC_units Requested units for the calculated PEC. Only µg/L currently supported @@ -24,9 +24,13 @@ #' # Alternatively, we can use the formula for a single application to "Ackerbau" from the paper #' PEC_sw_drift(100, drift_data = "RF") #' -#' # This makes it possible to also use different substances +#' # This makes it possible to also use different distances #' PEC_sw_drift(100, distances = c(1, 3, 5, 6, 10, 20, 50, 100), drift_data = "RF") #' +#' # or consider aerial application +#' PEC_sw_drift(100, distances = c(1, 3, 5, 6, 10, 20, 50, 100), drift_data = "RF", +#' crop_group_focus = "aerial") +#' #' # Using custom drift percentages is also supported #' PEC_sw_drift(100, drift_percentages = c(2.77, 0.95, 0.57, 0.48, 0.29, 0.15, 0.06, 0.03)) PEC_sw_drift <- function(rate, @@ -34,7 +38,11 @@ PEC_sw_drift <- function(rate, water_depth = 30, drift_percentages = NULL, drift_data = c("JKI", "RF"), - crop = "Ackerbau", + crop_group_JKI = c("Ackerbau", + "Obstbau frueh", "Obstbau spaet", "Weinbau frueh", "Weinbau spaet", + "Hopfenbau", "Flaechenkulturen > 900 l/ha", "Gleisanlagen"), + crop_group_focus = c("arable", "hops", "vines, late", "vines, early", + "fruit, late", "fruit, early", "aerial"), distances = c(1, 5, 10, 20), rate_units = "g/ha", PEC_units = "\u00B5g/L") @@ -42,22 +50,16 @@ PEC_sw_drift <- function(rate, rate_units <- match.arg(rate_units) PEC_units <- match.arg(PEC_units) drift_data <- match.arg(drift_data) + crop_group_JKI <- match.arg(crop_group_JKI) + crop_group_focus <- match.arg(crop_group_focus) water_volume <- 100 * 100 * (water_depth/100) * 1000 # in L (for 1 ha) PEC_sw_overspray <- rate * 1e6 / water_volume # in µg/L dist_index <- as.character(distances) - RF <- list( - "1" = list("Ackerbau" = function(distance) 2.7705 * distance^-0.9787) # p. 134 - ) - if (drift_data == "RF") { - if (is.null(RF[[as.character(applications)]])) stop("Rautmann formula not included for ", applications, " applications") - if (is.null(RF[[as.character(applications)]][[crop]])) stop("Rautmann formula not included for this case") - } - if (is.null(drift_percentages)) { drift_percentages <- switch(drift_data, - JKI = pfm::drift_data_JKI[[applications]][dist_index, crop], - RF = RF[[applications]][[crop]](distances) + JKI = pfm::drift_data_JKI[[applications]][dist_index, crop_group_JKI], + RF = drift_percentages_rautmann_formula(distances, applications, crop_group_focus) ) names(drift_percentages) <- paste(dist_index, "m") } else { @@ -67,3 +69,37 @@ PEC_sw_drift <- function(rate, PEC_sw_drift <- PEC_sw_overspray * drift_percentages / 100 return(PEC_sw_drift) } + +#' Calculate the drift percentages according to the Rautmann formula +#' +#' @param distances The distances in m for which to get PEC values +#' @param applications Number of applications for selection of drift percentile +#' @param crop_group_focus One of the crop groups as used in [drift_parameters_focus] +#' @export +#' @examples +# # One application on field crops, for 1 m, 3 m and 5 m distance +#' drift_data_JKI[[1]][as.character(c(1, 3, 5)), "Ackerbau"] +#' drift_percentages_rautmann_formula(c(1, 3, 5)) +#' +#' # One application to early or late fruit crops +#' drift_data_JKI[[1]][as.character(c(3, 5, 20, 50)), "Obstbau frueh"] +#' drift_percentages_rautmann_formula(c(3, 5, 20, 50), crop_group = "fruit, early") +#' drift_data_JKI[[1]][as.character(c(3, 5, 20, 50)), "Obstbau spaet"] +#' drift_percentages_rautmann_formula(c(3, 5, 20, 50), crop_group = "fruit, late") +drift_percentages_rautmann_formula <- function(distances, applications = 1, + crop_group_focus = c("arable", "hops", "vines, late", "vines, early", "fruit, late", + "fruit, early", "aerial")) +{ + cg <- match.arg(crop_group_focus) + if (!applications %in% 1:8) stop("Only 1 to 8 applications are supported") + parms <- pfm::drift_parameters_focus[pfm::drift_parameters_focus$crop_group == cg & + pfm::drift_parameters_focus$n_apps == applications, c("A", "B", "C", "D", "hinge")] + + drift_percentages = with(as.list(parms), { + A <- ifelse(distances < hinge, A, C) + B <- ifelse(distances < hinge, B, D) + A * distances^B + } + ) + return(drift_percentages) +} diff --git a/R/drift_parameters_focus.R b/R/drift_parameters_focus.R index 6a456c6..c5e5661 100644 --- a/R/drift_parameters_focus.R +++ b/R/drift_parameters_focus.R @@ -4,6 +4,11 @@ #' using the R code given in the file `data_generation/drift_parameters_focus.R` #' installed with this package. The appendix itself is not included in the package, #' as its licence is not clear. +#' +#' For the hinge distance, `Inf` was substituted for the cases where no hinge +#' distance is given in the data, in this way parameters C and D are never +#' used for any distance if A and B are used for the case that the distance +#' is smaller than the hinge distance. #' #' @name drift_parameters_focus #' @docType data @@ -23,4 +28,5 @@ #' @keywords datasets #' @examples #' drift_parameters_focus +#' unique(drift_parameters_focus$crop_group) "drift_parameters_focus" -- cgit v1.2.1