VFSAPLE: Link the VFS and APLE models.

Description Usage Arguments Details Value Author(s) See Also Examples

Description

Uses the erosion and runoff output of VFS as input to APLE.

Usage

1
2
VFSAPLE(x, soilP, OM, manureP = 25, manureSolids = 25, manureWEP = 50,
    manureIn = 40, fertP = 10, fertIn = 40)

Arguments

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 (%).

Details

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.

Value

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 (%).

Author(s)

Sarah Goslee

See Also

APLE, VFS

Examples

 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)

sgoslee/VFS documentation built on May 16, 2019, 11:13 p.m.