1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/PEC_sw_drift.R
\name{PEC_sw_drift}
\alias{PEC_sw_drift}
\title{Calculate predicted environmental concentrations in surface water due to drift}
\usage{
PEC_sw_drift(
rate,
applications = 1,
water_depth = as_units("30 cm"),
drift_percentages = NULL,
drift_data = c("JKI", "RF"),
crop_group_JKI = "Ackerbau",
crop_group_RF = "arable",
distances = c(1, 5, 10, 20),
formula = c("Rautmann", "FOCUS"),
water_width = as_units("100 cm"),
side_angle = 90,
rate_units = "g/ha",
PEC_units = "µg/L"
)
}
\arguments{
\item{rate}{Application rate in units specified below, or with units defined via the
\code{units} package.}
\item{applications}{Number of applications for selection of drift percentile}
\item{water_depth}{Depth of the water body in cm}
\item{drift_percentages}{Percentage drift values for which to calculate PECsw.
Overrides 'drift_data', 'distances', 'applications', crop group and
formula arguments if not NULL.}
\item{drift_data}{Source of drift percentage data. If 'JKI', the \link{drift_data_JKI}
included in the package is used. If 'RF', the Rautmann drift data are calculated
either in the original form or integrated over the width of the water body, depending
on the 'formula' argument.}
\item{crop_group_JKI}{When using the 'JKI' drift data, one of the German names
as used in \link{drift_data_JKI}. Will only be used if drift_data is 'JKI'. Available
crop groups are "Ackerbau", "Obstbau frueh", "Obstbau spaet",
"Weinbau frueh", "Weinbau spaet", "Hopfenbau", "Flaechenkulturen > 900 l/ha" and
"Gleisanlagen".}
\item{crop_group_RF}{Crop group(s) as used in \link{drift_parameters_focus}, i.e.
"arable", "hops", "vines, late", "vines, early", "fruit, late", "fruit, early"
or "aerial".}
\item{distances}{The distances in m for which to get PEC values}
\item{formula}{By default, the original Rautmann formula is used. If you
specify "FOCUS", mean drift input over the width of the water body is
calculated as described in Chapter 5.4.5 of the FOCUS surface water guidance}
\item{water_width}{Width of the water body in cm}
\item{side_angle}{The angle of the side of the water relative to the bottom which
is assumed to be horizontal, in degrees. The SYNOPS model assumes 45 degrees here.}
\item{rate_units}{Defaults to g/ha. For backwards compatibility, only used
if the specified rate does not have \link[units:units]{units::units}].}
\item{PEC_units}{Requested units for the calculated PEC. Only µg/L currently supported}
}
\value{
A numeric vector with the predicted concentration in surface water.
In some cases, the vector is named with distances or drift percentages, for
backward compatibility with versions before the vectorisation of arguments
other than 'distances' was introduced in v0.6.5.
}
\description{
This is a basic, vectorised form of a simple calculation of a contaminant
concentration in surface water based on complete, instantaneous mixing
with input via spray drift.
}
\details{
It is recommened to specify the arguments \code{rate}, \code{water_depth} and
\code{water_width} using \link[units:units]{units::units} from the \code{units} package.
Since pfm version 0.6.5, the function is vectorised with respect to rates,
applications, water depth, crop groups and distances
}
\examples{
PEC_sw_drift(100)
# 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 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_RF = "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))
# The influence of assuming a 45° angle of the sides of the waterbody and the width of the
# waterbody can be illustrated
PEC_sw_drift(100)
PEC_sw_drift(100, drift_data = "RF")
PEC_sw_drift(100, drift_data = "RF", formula = "FOCUS")
PEC_sw_drift(100, drift_data = "RF", formula = "FOCUS", side_angle = 45)
PEC_sw_drift(100, drift_data = "RF", formula = "FOCUS", side_angle = 45, water_width = 200)
# The function is vectorised with respect to rates, applications, water depth,
# crop groups and distances
PEC_sw_drift(
rate = rep(100, 6),
applications = c(1, 2, rep(1, 4)),
water_depth = c(30, 30, 30, 60, 30, 30),
crop_group_JKI = c(rep("Ackerbau", 4), rep("Obstbau frueh", 2)),
distances = c(rep(5, 4), 10, 5))
# Try the same with the Rautmann formula
PEC_sw_drift(
rate = rep(100, 6),
applications = c(1, 2, rep(1, 4)),
water_depth = c(30, 30, 30, 60, 30, 30),
drift_data = "RF",
crop_group_RF = c(rep("arable", 4), rep("fruit, early", 2)),
distances = c(rep(5, 4), 10, 5))
# And with the FOCUS variant
PEC_sw_drift(
rate = rep(100, 6),
applications = c(1, 2, rep(1, 4)),
water_depth = c(30, 30, 30, 60, 30, 30),
drift_data = "RF",
formula = "FOCUS",
crop_group_RF = c(rep("arable", 4), rep("fruit, early", 2)),
distances = c(rep(5, 4), 10, 5))
}
\seealso{
\link{drift_parameters_focus}, \link{drift_percentages_rautmann}
}
\author{
Johannes Ranke
}
|
