Description Usage Arguments Details Value References Examples
Simple model to calculate crop growth and yield based on MISCANMOD (see references).
1 2 3 4 | 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)
|
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. |
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)
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. |
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.
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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.