Description Usage Arguments Details Value Author(s) See Also Examples
Uses the erosion and runoff output of VFS as input to APLE.
1 2 | VFSAPLE(x, soilP, OM, manureP = 25, manureSolids = 25, manureWEP = 50,
manureIn = 40, fertP = 10, fertIn = 40)
|
x |
VFS object. |
soilP |
soil test Mehlich 3 phosphorus (mg/kg). |
OM |
soil organic matter (%). |
manureP |
manure P applied (kg/ha). |
manureSolids |
manure solids (%). |
manureWEP |
manure water-extractable phosphorus/TP (%). |
manureIn |
manure incorporated (%). |
fertP |
fertilizer phosphorus applied (kg/ha). |
fertIn |
fertilizer incorporated (%). |
The VFSAPLE
function handles unit conversion and parameter estimation. Erosion, precipitation, runoff, and field parameters from VFS
are passed to APLE
. It also runs APLE
for both pre- and post-filter strip erosion values, so the efficacy of the filter strip at phosphorus removal can be calculated.
preVFS |
APLE object for pre-filter strip erosion values. |
postVFS |
APLE object for post-filter strip erosion values. |
pErosion |
Efficacy of the filter strip at removing erosion phosphorus (%). |
pTotal |
Efficacy of the fiter strip at removing total phosphorus (%). |
Sarah Goslee
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 | # state college GHCN data
#
# weather <- read.dly(system.file("extdata", "USC00368449.dly", package = "VFS"))
data("weather") # same object
weather.param <- wth.param(weather, method="markov")
rain.compare <- rainfall(365*2, weather.param)
temp.compare <- temperature(365*2, weather.param)
data(soildat)
data(bufferdat)
# bluegrass buffer, clay loam soil
# short simulation to cut down on time required
vfs.CL <- VFS(nyears = 2, thissoil = subset(soildat, Soil == "CL"),
rain=rain.compare, temperature=temp.compare,
thisbuffer = subset(bufferdat, Species == "bluegrass"), Duration = 2,
FieldArea = 4000, VFSwidth = 10.7, VFSslope = 0.02,
z = 1000, b = 1.5)
print(vfs.CL)
summary(vfs.CL)
aple.CL <- VFSAPLE(vfs.CL, soilP = 120, OM = 2)
print(aple.CL)
summary(aple.CL)
|
VFS model: 2 years
Mean annual load into vegetated filter strip:
b15
0.5338093
Mean annual load out of vegetated filter strip:
b15
0.3579431
Mean annual removal efficiency:
b15
33.12677
Mean annual load into vegetated filter strip (MUSLE):
x
1.930142
Mean annual load out of vegetated filter strip (MUSLE):
dat
1.239375
Mean annual removal efficiency (MUSLE):
dat
36.02485
ALR ALRsd ALRall ALRallsd APEA APEAsd
3.312677e+01 4.096398e-01 3.312677e+01 4.096398e-01 4.433833e+01 1.084648e+02
SedIn SedInsd SedOut SedOutsd TLR MUSLEIn
5.338093e-01 1.247220e-02 3.579431e-01 8.469482e-03 3.294551e+01 1.930142e+00
MUSLEInsd MUSLEOut MUSLEOutsd MUSLETLR Runoff Runoffsd
2.389202e-02 1.239375e+00 1.599706e-02 3.578840e+01 4.157537e+02 3.063517e+02
RunoffDays RunoffDayssd
4.050000e+01 2.757716e+01
Mean annual erosion P reduction by VFS: 28.467
Mean annual total P reduction by VFS: 11.867
AnnualErosionPRemoval AnnualErosionPRemovalsd
28.4669907 0.5289076
AnnualTotalPRemoval AnnualTotalPRemovalsd
11.8666069 1.0546169
AnnualLossErosionPre AnnualLossDissolvedSoilPre
2.6213403 1.2472611
AnnualLossDissolvedManurePre AnnualLossDissolvedFertPre
2.3030086 0.2670162
AnnualLossTotalPre AnnualLossErosionPost
6.4386262 1.8792016
AnnualLossDissolvedSoilPost AnnualLossDissolvedManurePost
1.2472611 2.3030086
AnnualLossDissolvedFertPost AnnualLossTotalPost
0.2670162 5.6964875
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.