222: Implementation of the Yasso07 model

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

Description

This function creates a model for five pools as described in Tuomi et al. (2009)

Usage

1
2
3
4
5
Yasso07Model(t, ks = c(kA = 0.66, kW = 4.3, kE = 0.35, kN = 0.22, 
    kH = 0.0033), p = c(p1 = 0.32, p2 = 0.01, p3 = 0.93, p4 = 0.34, 
    p5 = 0, p6 = 0, p7 = 0, p8 = 0, p9 = 0.01, p10 = 0, p11 = 0, 
    p12 = 0.92, pH = 0.04), C0, In, xi = 1, solver = deSolve.lsoda.wrapper, 
    pass = FALSE)

Arguments

t

A vector containing the points in time where the solution is sought.

ks

A vector of lenght 5 containing the values of the decomposition rates for each pool.

p

A vector of length 13 containing transfer coefficients among different pools.

C0

A vector containing the initial amount of carbon for the 5 pools. The length of this vector must be 5.

In

A single scalar or data.frame object specifying the amount of litter inputs by time

xi

A scalar or data.frame object specifying the external (environmental and/or edaphic) effects on decomposition rates.

solver

A function that solves the system of ODEs. This can be euler or ode or any other user provided function with the same interface.

pass

if TRUE forces the constructor to create the model even if it is invalid

Value

A Model Object that can be further queried

Author(s)

Carlos A. Sierra, Markus Mueller

References

Tuomi, M., Thum, T., Jarvinen, H., Fronzek, S., Berg, B., Harmon, M., Trofymow, J., Sevanto, S., and Liski, J. (2009). Leaf litter decomposition-estimates of global variability based on Yasso07 model. Ecological Modelling, 220:3362 - 3371.

See Also

ICBMModel, RothCModel

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
years=seq(0,50,0.1) 
C0=rep(100,5)
In=0

Ex1=Yasso07Model(t=years,C0=C0,In=In)
Ct=getC(Ex1)
Rt=getReleaseFlux(Ex1)

plotCPool(years,Ct,col=1:5,xlab="years",ylab="C pool",
          ylim=c(0,max(Ct)))
legend("topright",c("xA","xW","xE","xN","xH"),lty=1,col=1:5,bty="n")

plotCPool(years,Rt,col=1:5,xlab="years",ylab="Respiration",ylim=c(0,50))
legend("topright",c("xA","xW","xE","xN","xH"),lty=1,col=1:5,bty="n")

SoilR documentation built on May 4, 2017, 9:08 p.m.