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
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
|
<!-- Generated by pkgdown: do not edit by hand -->
<!DOCTYPE html>
<html lang="en">
<head>
<meta charset="utf-8">
<meta http-equiv="X-UA-Compatible" content="IE=edge">
<meta name="viewport" content="width=device-width, initial-scale=1.0">
<title>Calculate predicted environmental concentrations in soil — PEC_soil • pfm</title>
<!-- jquery -->
<script src="https://cdnjs.cloudflare.com/ajax/libs/jquery/3.3.1/jquery.min.js" integrity="sha256-FgpCb/KJQlLNfOu91ta32o/NMZxltwRo8QtmkMRdAu8=" crossorigin="anonymous"></script>
<!-- Bootstrap -->
<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.3.7/css/bootstrap.min.css" integrity="sha256-916EbMg70RQy9LHiGkXzG8hSg9EdNy97GazNG/aiY1w=" crossorigin="anonymous" />
<script src="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.3.7/js/bootstrap.min.js" integrity="sha256-U5ZEeKfGNOja007MMD3YBI0A3OSZOQbeG6z2f2Y0hu8=" crossorigin="anonymous"></script>
<!-- Font Awesome icons -->
<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/4.7.0/css/font-awesome.min.css" integrity="sha256-eZrrJcwDc/3uDhsdt61sL2oOBY362qM3lon1gyExkL0=" crossorigin="anonymous" />
<!-- clipboard.js -->
<script src="https://cdnjs.cloudflare.com/ajax/libs/clipboard.js/2.0.4/clipboard.min.js" integrity="sha256-FiZwavyI2V6+EXO1U+xzLG3IKldpiTFf3153ea9zikQ=" crossorigin="anonymous"></script>
<!-- sticky kit -->
<script src="https://cdnjs.cloudflare.com/ajax/libs/sticky-kit/1.1.3/sticky-kit.min.js" integrity="sha256-c4Rlo1ZozqTPE2RLuvbusY3+SU1pQaJC0TjuhygMipw=" crossorigin="anonymous"></script>
<!-- pkgdown -->
<link href="../pkgdown.css" rel="stylesheet">
<script src="../pkgdown.js"></script>
<meta property="og:title" content="Calculate predicted environmental concentrations in soil — PEC_soil" />
<meta property="og:description" content="This is a basic calculation of a contaminant concentration in bulk soil
based on complete, instantaneous mixing. If an interval is given, an
attempt is made at calculating a long term maximum concentration using
the concepts layed out in the PPR panel opinion (EFSA PPR panel 2012
and in the EFSA guidance on PEC soil calculations (EFSA, 2015, 2017)." />
<meta name="twitter:card" content="summary" />
<!-- mathjax -->
<script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js" integrity="sha256-nvJJv9wWKEm88qvoQl9ekL2J+k/RWIsaSScxxlsrv8k=" crossorigin="anonymous"></script>
<script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/config/TeX-AMS-MML_HTMLorMML.js" integrity="sha256-84DKXVJXs0/F8OTMzX4UR909+jtl4G7SPypPavF+GfA=" crossorigin="anonymous"></script>
<!--[if lt IE 9]>
<script src="https://oss.maxcdn.com/html5shiv/3.7.3/html5shiv.min.js"></script>
<script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script>
<![endif]-->
</head>
<body>
<div class="container template-reference-topic">
<header>
<div class="navbar navbar-default navbar-fixed-top" role="navigation">
<div class="container">
<div class="navbar-header">
<button type="button" class="navbar-toggle collapsed" data-toggle="collapse" data-target="#navbar" aria-expanded="false">
<span class="sr-only">Toggle navigation</span>
<span class="icon-bar"></span>
<span class="icon-bar"></span>
<span class="icon-bar"></span>
</button>
<span class="navbar-brand">
<a class="navbar-link" href="../index.html">pfm</a>
<span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">0.5.6</span>
</span>
</div>
<div id="navbar" class="navbar-collapse collapse">
<ul class="nav navbar-nav">
<li>
<a href="../index.html">
<span class="fa fa-home fa-lg"></span>
</a>
</li>
<li>
<a href="../reference/index.html">Reference</a>
</li>
</ul>
<ul class="nav navbar-nav navbar-right">
<li>
<a href="https://github.com/jranke/pfm">
<span class="fa fa-github fa-lg"></span>
</a>
</li>
</ul>
</div><!--/.nav-collapse -->
</div><!--/.container -->
</div><!--/.navbar -->
</header>
<div class="row">
<div class="col-md-9 contents">
<div class="page-header">
<h1>Calculate predicted environmental concentrations in soil</h1>
<small class="dont-index">Source: <a href='https://github.com/jranke/pfm/blob/master/R/PEC_soil.R'><code>R/PEC_soil.R</code></a></small>
<div class="hidden name"><code>PEC_soil.Rd</code></div>
</div>
<div class="ref-description">
<p>This is a basic calculation of a contaminant concentration in bulk soil
based on complete, instantaneous mixing. If an interval is given, an
attempt is made at calculating a long term maximum concentration using
the concepts layed out in the PPR panel opinion (EFSA PPR panel 2012
and in the EFSA guidance on PEC soil calculations (EFSA, 2015, 2017).</p>
</div>
<pre class="usage"><span class='fu'>PEC_soil</span>(<span class='no'>rate</span>, <span class='kw'>rate_units</span> <span class='kw'>=</span> <span class='st'>"g/ha"</span>, <span class='kw'>interception</span> <span class='kw'>=</span> <span class='fl'>0</span>,
<span class='kw'>mixing_depth</span> <span class='kw'>=</span> <span class='fl'>5</span>, <span class='kw'>PEC_units</span> <span class='kw'>=</span> <span class='st'>"mg/kg"</span>, <span class='kw'>PEC_pw_units</span> <span class='kw'>=</span> <span class='st'>"mg/L"</span>,
<span class='kw'>interval</span> <span class='kw'>=</span> <span class='fl'>NA</span>, <span class='kw'>n_periods</span> <span class='kw'>=</span> <span class='fl'>Inf</span>, <span class='kw'>tillage_depth</span> <span class='kw'>=</span> <span class='fl'>20</span>,
<span class='kw'>leaching_depth</span> <span class='kw'>=</span> <span class='no'>tillage_depth</span>, <span class='kw'>crop</span> <span class='kw'>=</span> <span class='st'>"annual"</span>,
<span class='kw'>cultivation</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>, <span class='kw'>chent</span> <span class='kw'>=</span> <span class='fl'>NA</span>, <span class='kw'>DT50</span> <span class='kw'>=</span> <span class='fl'>NA</span>, <span class='kw'>FOMC</span> <span class='kw'>=</span> <span class='fl'>NA</span>,
<span class='kw'>Koc</span> <span class='kw'>=</span> <span class='fl'>NA</span>, <span class='kw'>Kom</span> <span class='kw'>=</span> <span class='no'>Koc</span>/<span class='fl'>1.724</span>, <span class='kw'>t_avg</span> <span class='kw'>=</span> <span class='fl'>0</span>, <span class='kw'>t_act</span> <span class='kw'>=</span> <span class='kw'>NULL</span>,
<span class='kw'>scenarios</span> <span class='kw'>=</span> <span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/c'>c</a></span>(<span class='st'>"default"</span>, <span class='st'>"EFSA_2017"</span>, <span class='st'>"EFSA_2015"</span>),
<span class='kw'>leaching</span> <span class='kw'>=</span> <span class='no'>scenarios</span> <span class='kw'>==</span> <span class='st'>"EFSA_2017"</span>, <span class='kw'>porewater</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>)</pre>
<h2 class="hasAnchor" id="arguments"><a class="anchor" href="#arguments"></a>Arguments</h2>
<table class="ref-arguments">
<colgroup><col class="name" /><col class="desc" /></colgroup>
<tr>
<th>rate</th>
<td><p>Application rate in units specified below</p></td>
</tr>
<tr>
<th>rate_units</th>
<td><p>Defaults to g/ha</p></td>
</tr>
<tr>
<th>interception</th>
<td><p>The fraction of the application rate that does not reach the soil</p></td>
</tr>
<tr>
<th>mixing_depth</th>
<td><p>Mixing depth in cm</p></td>
</tr>
<tr>
<th>PEC_units</th>
<td><p>Requested units for the calculated PEC. Only mg/kg currently supported</p></td>
</tr>
<tr>
<th>PEC_pw_units</th>
<td><p>Only mg/L currently supported</p></td>
</tr>
<tr>
<th>interval</th>
<td><p>Period of the deeper mixing. The default is NA, i.e. no
deeper mixing. For annual deeper mixing, set this to 365 when degradation
units are in days</p></td>
</tr>
<tr>
<th>n_periods</th>
<td><p>Number of periods to be considered for long term PEC calculations</p></td>
</tr>
<tr>
<th>tillage_depth</th>
<td><p>Periodic (see interval) deeper mixing in cm</p></td>
</tr>
<tr>
<th>leaching_depth</th>
<td><p>EFSA (2017) uses the mixing depth (ecotoxicological
evaluation depth) to calculate leaching for annual crops where tillage
takes place. By default, losses from the layer down to the tillage
depth are taken into account in this implementation.</p></td>
</tr>
<tr>
<th>crop</th>
<td><p>Ignored for scenarios other than EFSA_2017. Only annual crops
are supported when these scenarios are used. Only crops with a single cropping
cycle per year are currently supported.</p></td>
</tr>
<tr>
<th>cultivation</th>
<td><p>Does mechanical cultivation in the sense of EFSA (2017)
take place, i.e. twice a year to a depth of 5 cm? Ignored for scenarios
other than EFSA_2017</p></td>
</tr>
<tr>
<th>chent</th>
<td><p>An optional chent object holding substance specific information. Can
also be a name for the substance as a character string</p></td>
</tr>
<tr>
<th>DT50</th>
<td><p>If specified, overrides soil DT50 endpoints from a chent object
If DT50 is not specified here and not available from the chent object, zero
degradation is assumed</p></td>
</tr>
<tr>
<th>FOMC</th>
<td><p>If specified, it should be a named numeric vector containing
the FOMC parameters alpha and beta. This overrides any other degradation
endpoints, and the degradation during the interval and after the maximum PEC
is calculated using these parameters without temperature correction</p></td>
</tr>
<tr>
<th>Koc</th>
<td><p>If specified, overrides Koc endpoints from a chent object</p></td>
</tr>
<tr>
<th>Kom</th>
<td><p>Calculated from Koc by default, but can explicitly be specified
as Kom here</p></td>
</tr>
<tr>
<th>t_avg</th>
<td><p>Averaging times for time weighted average concentrations</p></td>
</tr>
<tr>
<th>t_act</th>
<td><p>Time series for actual concentrations</p></td>
</tr>
<tr>
<th>scenarios</th>
<td><p>If this is 'default', the DT50 will be used without correction
and soil properties as specified in the REACH guidance (R.16, Table
R.16-9) are used for porewater PEC calculations. If this is "EFSA_2015",
the DT50 is taken to be a modelling half-life at 20°C and pF2 (for when
'chents' is specified, the DegT50 with destination 'PECgw' will be used),
and corrected using an Arrhenius activation energy of 65.4 kJ/mol. Also
model and scenario adjustment factors from the EFSA guidance are used.</p></td>
</tr>
<tr>
<th>leaching</th>
<td><p>Should leaching be taken into account? The default is FALSE,
except when the EFSA_2017 scenarios are used.</p></td>
</tr>
<tr>
<th>porewater</th>
<td><p>Should equilibrium porewater concentrations be estimated
based on Kom and the organic carbon fraction of the soil instead of total
soil concentrations? Based on equation (7) given in the PPR panel opinion
(EFSA 2012, p. 24) and the scenarios specified in the EFSA guidance (2015,
p. 13).</p></td>
</tr>
</table>
<h2 class="hasAnchor" id="value"><a class="anchor" href="#value"></a>Value</h2>
<p>The predicted concentration in soil</p>
<h2 class="hasAnchor" id="details"><a class="anchor" href="#details"></a>Details</h2>
<p>This assumes that the complete load to soil during the time specified by
'interval' (typically 365 days) is dosed at once. As in the PPR panel
opinion cited below (EFSA PPR panel 2012), only temperature correction using the
Arrhenius equation is performed.</p>
<p>Total soil and porewater PEC values for the scenarios as defined in the EFSA
guidance (2017, p. 14/15) can easily be calculated.</p>
<h2 class="hasAnchor" id="note"><a class="anchor" href="#note"></a>Note</h2>
<p>While time weighted average (TWA) concentrations given in the examples
from the EFSA guidance from 2015 (p. 80) are be reproduced, this is not
true for the TWA concentrations given for the same example in the EFSA guidance
from 2017 (p. 92).</p>
<p>According to the EFSA guidance (EFSA, 2017, p. 43), leaching should be
taken into account for the EFSA 2017 scenarios, using the evaluation depth
(here mixing depth) as the depth of the layer from which leaching takes
place. However, as the amount leaching below the evaluation depth
(often 5 cm) will partly be mixed back during tillage, the default in this function
is to use the tillage depth for the calculation of the leaching rate.</p>
<p>If temperature information is available in the selected scenarios, as
e.g. in the EFSA scenarios, the DT50 for groundwater modelling
(destination 'PECgw') is taken from the chent object, otherwise the DT50
with destination 'PECsoil'.</p>
<h2 class="hasAnchor" id="references"><a class="anchor" href="#references"></a>References</h2>
<p>EFSA Panel on Plant Protection Products and their Residues (2012)
Scientific Opinion on the science behind the guidance for scenario
selection and scenario parameterisation for predicting environmental
concentrations of plant protection products in soil. <em>EFSA Journal</em>
<b>10</b>(2) 2562, doi:10.2903/j.efsa.2012.2562</p>
<p>EFSA (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. <em>EFSA Journal</em> <b>15</b>(10) 4982
doi:10.2903/j.efsa.2017.4982</p>
<p>EFSA (European Food Safety Authority) (2015) EFSA guidance document for
predicting environmental concentrations of active substances of plant
protection products and transformation products of these active substances
in soil. <em>EFSA Journal</em> <b>13</b>(4) 4093
doi:10.2903/j.efsa.2015.4093</p>
<h2 class="hasAnchor" id="examples"><a class="anchor" href="#examples"></a>Examples</h2>
<pre class="examples"><div class='input'><span class='fu'>PEC_soil</span>(<span class='fl'>100</span>, <span class='kw'>interception</span> <span class='kw'>=</span> <span class='fl'>0.25</span>)</div><div class='output co'>#> scenario
#> t_avg default
#> 0 0.1</div><div class='input'>
<span class='co'># This is example 1 starting at p. 92 of the EFSA guidance (2017)</span>
<span class='co'># Note that TWA concentrations differ from the ones given in the guidance</span>
<span class='co'># for an unknown reason (the values from EFSA (2015) can be reproduced).</span>
<span class='fu'>PEC_soil</span>(<span class='fl'>1000</span>, <span class='kw'>interval</span> <span class='kw'>=</span> <span class='fl'>365</span>, <span class='kw'>DT50</span> <span class='kw'>=</span> <span class='fl'>250</span>, <span class='kw'>t_avg</span> <span class='kw'>=</span> <span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/c'>c</a></span>(<span class='fl'>0</span>, <span class='fl'>21</span>),
<span class='kw'>Kom</span> <span class='kw'>=</span> <span class='fl'>1000</span>, <span class='kw'>scenarios</span> <span class='kw'>=</span> <span class='st'>"EFSA_2017"</span>)</div><div class='output co'>#> scenario
#> t_avg CTN CTC CTS
#> 0 19.76834 13.8619 10.53795
#> 21 19.59345 13.7169 10.39882</div><div class='input'><span class='fu'>PEC_soil</span>(<span class='fl'>1000</span>, <span class='kw'>interval</span> <span class='kw'>=</span> <span class='fl'>365</span>, <span class='kw'>DT50</span> <span class='kw'>=</span> <span class='fl'>250</span>, <span class='kw'>t_av</span> <span class='kw'>=</span> <span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/c'>c</a></span>(<span class='fl'>0</span>, <span class='fl'>21</span>),
<span class='kw'>Kom</span> <span class='kw'>=</span> <span class='fl'>1000</span>, <span class='kw'>scenarios</span> <span class='kw'>=</span> <span class='st'>"EFSA_2017"</span>, <span class='kw'>porewater</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)</div><div class='output co'>#> scenario
#> t_avg CLN CLC CLS
#> 0 0.5541984 0.6779249 0.9816693
#> 21 0.5484576 0.6693125 0.9609119</div><div class='input'>
<span class='co'># This is example 1 starting at p. 79 of the EFSA guidance (2015)</span>
<span class='fu'>PEC_soil</span>(<span class='fl'>1000</span>, <span class='kw'>interval</span> <span class='kw'>=</span> <span class='fl'>365</span>, <span class='kw'>DT50</span> <span class='kw'>=</span> <span class='fl'>250</span>, <span class='kw'>t_avg</span> <span class='kw'>=</span> <span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/c'>c</a></span>(<span class='fl'>0</span>, <span class='fl'>21</span>),
<span class='kw'>scenarios</span> <span class='kw'>=</span> <span class='st'>"EFSA_2015"</span>)</div><div class='output co'>#> scenario
#> t_avg CTN CTC CTS
#> 0 21.96827 11.53750 9.145259
#> 21 21.78517 11.40701 9.017370</div><div class='input'><span class='fu'>PEC_soil</span>(<span class='fl'>1000</span>, <span class='kw'>interval</span> <span class='kw'>=</span> <span class='fl'>365</span>, <span class='kw'>DT50</span> <span class='kw'>=</span> <span class='fl'>250</span>, <span class='kw'>t_av</span> <span class='kw'>=</span> <span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/c'>c</a></span>(<span class='fl'>0</span>, <span class='fl'>21</span>),
<span class='kw'>Kom</span> <span class='kw'>=</span> <span class='fl'>1000</span>, <span class='kw'>scenarios</span> <span class='kw'>=</span> <span class='st'>"EFSA_2015"</span>, <span class='kw'>porewater</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)</div><div class='output co'>#> scenario
#> t_avg CLN CLC CLS
#> 0 0.7589401 0.6674322 0.9147861
#> 21 0.7506036 0.6590345 0.8987279</div><div class='input'>
<span class='co'># The following is from example 4 starting at p. 85 of the EFSA guidance (2015)</span>
<span class='co'># Metabolite M2</span>
<span class='co'># Calculate total and porewater soil concentrations for tier 1 scenarios</span>
<span class='co'># Relative molar mass is 100/300, formation fraction is 0.7 * 1</span>
<span class='no'>results_pfm</span> <span class='kw'><-</span> <span class='fu'>PEC_soil</span>(<span class='fl'>100</span>/<span class='fl'>300</span> * <span class='fl'>0.7</span> * <span class='fl'>1</span> * <span class='fl'>1000</span>, <span class='kw'>interval</span> <span class='kw'>=</span> <span class='fl'>365</span>, <span class='kw'>DT50</span> <span class='kw'>=</span> <span class='fl'>250</span>, <span class='kw'>t_avg</span> <span class='kw'>=</span> <span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/c'>c</a></span>(<span class='fl'>0</span>, <span class='fl'>21</span>),
<span class='kw'>scenarios</span> <span class='kw'>=</span> <span class='st'>"EFSA_2015"</span>)
<span class='no'>results_pfm_pw</span> <span class='kw'><-</span> <span class='fu'>PEC_soil</span>(<span class='fl'>100</span>/<span class='fl'>300</span> * <span class='fl'>0.7</span> * <span class='fl'>1000</span>, <span class='kw'>interval</span> <span class='kw'>=</span> <span class='fl'>365</span>, <span class='kw'>DT50</span> <span class='kw'>=</span> <span class='fl'>250</span>, <span class='kw'>t_av</span> <span class='kw'>=</span> <span class='fu'><a href='https://www.rdocumentation.org/packages/base/topics/c'>c</a></span>(<span class='fl'>0</span>, <span class='fl'>21</span>),
<span class='kw'>Kom</span> <span class='kw'>=</span> <span class='fl'>100</span>, <span class='kw'>scenarios</span> <span class='kw'>=</span> <span class='st'>"EFSA_2015"</span>, <span class='kw'>porewater</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)</div></pre>
</div>
<div class="col-md-3 hidden-xs hidden-sm" id="sidebar">
<h2>Contents</h2>
<ul class="nav nav-pills nav-stacked">
<li><a href="#arguments">Arguments</a></li>
<li><a href="#value">Value</a></li>
<li><a href="#details">Details</a></li>
<li><a href="#note">Note</a></li>
<li><a href="#references">References</a></li>
<li><a href="#examples">Examples</a></li>
</ul>
<h2>Author</h2>
<p>Johannes Ranke</p>
</div>
</div>
<footer>
<div class="copyright">
<p>Developed by Johannes Ranke.</p>
</div>
<div class="pkgdown">
<p>Site built with <a href="https://pkgdown.r-lib.org/">pkgdown</a> 1.3.0.</p>
</div>
</footer>
</div>
</body>
</html>
|
