Rmiscanmod: RUE-based model to calculate miscanthus growth and yield.

In BioCro: Simulation of Biomass Crops

Description

Simple model to calculate crop growth and yield based on MISCANMOD (see references).

Usage

Rmiscanmod(data, RUE = 2.4, LER = 0.01, Tb = 10, k = 0.67, LAIdrd = 0.8, 
    LAIStop = 1.8, RUEdrd = 1.3, RUEStop = 2.5, SMDdrd = -30, 
    SMDStop = -120, FieldC = 45, iWatCont = 45, a = 6682.2, b = -0.33, 
    soildepth = 0.6) 

Arguments

data	data.frame or matrix described in details.
RUE	Radiation use efficiency (g/MJ).
LER	Leaf expansion rate LAI/GDD.
Tb	Base Temperature (Celsius).
k	extinction coefficient of light in the canopy.
LAIdrd	Leaf Area Index 'down regulation decline'.
LAIStop	Leaf Area Index 'down regulation decline' threshold .
RUEdrd	Radiation Use Efficieny 'down regulation decline'.
RUEStop	Radiation Use Efficieny 'down regulation decline' threshold.
SMDdrd	Soil Moisture Deficit 'down regulation decline'.
SMDStop	Soil Moisture Deficit 'down regulation decline' threshold.
FieldC	Soil field capacity.
iWatCont	Initial water content.
a	Soil parameter.
b	Soil parameter.
soildepth	Soil depth.

Details

The data.frame or matrix should contain

column 1: year column 2: month column 3: day column 4: JD column 5: max T (Celsius) column 6: min T (Celsius) column 7: PPFD or solar radiation divided by 2 (MJ/m2) column 8: Potential evaporation column 9: precip (mm)

Value

returns a list

PotEvp	Potential Evaporation.
Deficitp	Deficitp
SMDp	Soil Moisture Deficit (potential)
AE.PE	Actual Evaporation / Potential Evaporation
Deficitp2	Deficitp2
SMDa	Soil Moisture Deficit (actual)
diffRainPE	difference between Rainfall and potential evaporation.
H2oper	H2O percent.
SoilMoist	Soil Moisture.
SoilMatPot	Soil Matric Potential.
WL.LER	Water limited Leaf Expansion Rate.
WL.RUE	Water limited Radiation Use Efficiency.
DDaTb	Degree Days above base Temperature.
DDcum	Degree Days (cumulative).
adjSumDD	adjusted Sum of Degree Days.
LAI	Leaf Area Index.
pLI	proportion of light intercepted.
Yield	Yield (dry biomass) (g/m2) to convert to Mg/ha divide by 100.

References

Clifton-Brown, J. C.; Neilson, B.; Lewandowski, I. and Jones, M. B. The modelled productivity of Miscanthus x giganteus (GREEF et DEU) in Ireland. Industrial Crops and Products, 2000, 12, 97-109.

Clifton-brown, J. C.; Stampfl, P. F. and Jones, M. B. Miscanthus biomass production for energy in Europe and its potential contribution to decreasing fossil fuel carbon emissions. Global Change Biology, 2004, 10, 509-518.

Examples

## Need to get an example data set and then run it
## Not run: 
data(WD1979)

res <- Rmiscanmod(WD1979)

## convert to Mg/ha

Yld <- res$Yield / 100

xyplot(Yld ~ 1:365 ,
       xlab="doy",
       ylab="Dry biomass (Mg/ha)")

## although the default value for Field Capacity is 45
## a more reasonable value is closer to 27



## End(Not run)

BioCro documentation built on May 2, 2019, 6:15 p.m.