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

In serbinsh/biocro: Simulation of Energycane and Sugarcane

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

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)

serbinsh/biocro documentation built on May 29, 2019, 6:57 p.m.