springbarley_biocontrol: Springbarley aphid biocontrol model function

In yclough/ecodeal: Functions for economic-ecological analyses of agricultural landscapes

Description

Uses spatial data on floral and nesting quality to compute flower visitation and population dynamics for wild pollinators for a given year. The function can be used for mutliple years by running a loop and using as input the output of the previous model (see example). useraster=TRUE,cSMD,mask,parameters,landscvals,reclassin_PGL,reclassin_PAC,N00="usedefault"

Usage

springbarley_biocontrol(useraster = TRUE, cSMD, mask, parameters,
  reclassin_PGL, reclassin_PAC, N00 = "usedefault")

Arguments

useraster	whether to use a raster, if FALSE the function assumes a matrix.
cSMD	land-use data, provided as a raster or a matrix.
mask	same format as cSMD but containing binary values, 1 = produce a prediction, 0 = ignore
parameters	list of two dataframes containing parameter values
reclassin_PGL	matrix indicating for a sequence of land-use codes the grasslands (1) and the non-grassland uses (0)
reclassin_PAC	matrix indicating for a sequence of land-use codes the annual land-uses (1) and the permanent ones (0)
N00	the density of colonizing aphids, does not vary in space, if set to usedefault its value is taken from OTHER_PAR.

Value

a list containing :

mc1

mortality phase 1

mc2

mortality phase 2

ad

aphid days, calculated as N00*(1-mc1)*(exp((f-mc2)*T)-1)/(f-mc2)

cd

crop loss in kg.ha-1, calculated as 12.7*ad^.66

References

Jonsson, M., Bommarco, R., Ekbom, B., Smith, H. G., Bengtsson, J., Caballero-Lopez, B., Winqvist, C. and Olsson, O. (2014), Ecological production functions for biological control services in agricultural landscapes. Methods Ecol Evol, 5: 243-252. http://dx.doi.org/10.1111/2041-210X.12149

Examples

PRED_PAR <- read.table("C:/a1workmain/ES_synthesis/BC_springbarley/yann/PRED_PAR.txt")
OTHER_PAR <- read.table("C:/a1workmain/ES_synthesis/BC_springbarley/yann/OTHER_PAR.txt")
parameters.bcsb<-list("PRED_PAR"=PRED_PAR,"OTHER_PAR"=OTHER_PAR)
reclassin_PAC<-as.matrix(reclassin_info[,c(1,2,4)])
reclassin_PGL<-as.matrix(reclassin_info[,c(1,2,5)])
biocontrol_pred<-springbarley_biocontrol(cSMD=cSMD_Sk2020,mask=landsmask,
parameters=parameters.bcsb,reclassin_PGL=reclassin_PGL,reclassin_PAC=reclassin_PAC)

yclough/ecodeal documentation built on May 14, 2019, 2:01 a.m.