VFSAPLE: Link the VFS and APLE models.

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

View source: R/VFSAPLE.R

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)

Example output

 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 

VFS documentation built on May 2, 2019, 8:58 a.m.