SoilR (version 1.1-23)

getF14R_method__Model_14: average radiocarbon fraction weighted by carbonrelease

Description

Calculates the average radiocarbon fraction weighted by the amount of carbon release at each time step. $(14R_1(t)+...+14R_n(t)) )/(R_1(t)+...R_n(t)))$ Where $14R_i(t)$ is the time dependent release of $14C$ of pool $i$ and $R_i(t)$ the release of all carbon isotops of pool $i$. Since the result is always in Absolute Fraction Modern format wie have to convert it to Delta14C

Arguments

object

an object

Value

$FR$ for each time step.

Examples

LitterInput=700 

Ex=TwopParallelModel14(t=years,ks=c(k1=1/2.8, k2=1/35),C0=c(200,5000), 
                       F0_Delta14C=c(0,0),In=LitterInput, gam=0.7,inputFc=C14Atm_NH,lag=2)
R14m=getF14R(Ex)
C14m=getF14C(Ex)
C14t=getF14(Ex)

par(mfrow=c(2,1))
plot(C14Atm_NH,type="l",xlab="Year",ylab="Delta 14C (per mil)",xlim=c(1940,2010)) 
lines(years, C14t[,1], col=4)
lines(years, C14t[,2],col=4,lwd=2)
legend("topright",c("Delta 14C Atmosphere", "Delta 14C pool 1", "Delta 14C pool 2"),
       lty=c(1,1,1),col=c(1,4,4),lwd=c(1,1,2),bty="n")

plot(C14Atm_NH,type="l",xlab="Year",ylab="Delta 14C (per mil)",xlim=c(1940,2010)) 
lines(years,C14m,col=4)
lines(years,R14m,col=2)
legend("topright",c("Delta 14C Atmosphere","Delta 14C SOM", "Delta 14C Respired"),
       lty=c(1,1,1), col=c(1,4,2),bty="n")
par(mfrow=c(1,1))