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

Description Usage Arguments Details Value References Examples

Description

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

Usage

1
2
3
4
5
  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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
## 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.