Flemming Nielsen (FACT/Banana hill) (2009)
Main findings: Forecast of annual Jatropha dry seed yield for twenty years based on the best knowledge currently available. See table and figure at the bottom of this page for details.
As part of the development of a business plan for the Bilibiza Biofuel Centre (BBC) it is important to have yield fore cast for Jatropha. The first part of this document will assess the
Next the maximum Jatropha production per household will be assessed under the assumption that farmers continue to prioritise production of food for own consumption and only use surplus labour on Jatropha.
Next the geographical distribution of Jatropha production will be assessed using different scenarios:
An attempt was made to create a regression equation for yield.
Date for Jatropha yield based in relation to age and rainfall data was extracted from Figure 2 in Achten et al 2008.Illustration 1: From Achten et al 2008: "Fig.2–Dry seed yield in relation to average annual rainfall (mm) and age of the JCL crop. The plotted points represent a mix of provenances, site conditions and plantage or average annual rainfall. Sources:[9,11,13,52,53,56] and personal communication: Kumar, 2005 and Buisman, 2005."
The article is available in pdf format and the data is not in table form but only in two graphs in Figure 2. To extract the numeric values from the graphs a screen dump was saved as a bmp file.
The Enguage Digitizer was used to obtain numeric values from the screen dump. Comparison of values extracted from the two graphs showed a difference of max ±2%. The imprecision is likely caused by the limited screen resolution.
Matching of data from the two graphs was done manually in OpenOffice Calc and the data was saved in ASCII file named 2009-09_Yield_data.txt:
age yield rain locality
Further analysis was done in R. To import the data file:
Jatropha <- read.table("2009-09_Yield_data.txt", header=TRUE)
Misc. commands for data exploration:
Different sub-sets of the data was used but no sensible regression equation could be obtained. There are several reasons for this. First of all some of the data from areas with very low rainfall (< 600 mm) have high yields. Since Jatropha usually does not yield when the rainfall is below 600 mm those data must be from areas where Jatropha relies on water from other source, be it irrigation or high ground water tables.
Other data is problematic too. In the graph on yield by rainfall the peak yield has been used for Nicaragua and Paraguay. Other data from high rainfall area are very low but the age of the plants is not known.
Excluding all outliers leaves us with mainly the two most complete datasets namely the Nicaragua and the Paraguay data series. The rainfall is higher in Paraguay than in Nicaragua but the yield is lower. A regression equation based on these data alone will therefore predict higher yield with decreasing rainfall, i.e. grossly over-estimate the yield for the Bilibiza area.
The multiple regression method is therefore not appropriate with the current data set.
The longest data serie from Paraguay shows the S-curve shape that is characteristic for many yield curves. It can often be described with the generalised logistic curve, also called Richard's curve:
Where Y is the yield and T is the year. The five parameters are:
A, the lower asymptote;
C, the upper asymptote, i.e. the maximum yield
M, the time of maximum growth;
B, the growth rate;
T, near which asymptote maximum growth occurs.
The Paraguay data series was entered in a spreadsheet and visual curve fitting was done resulting in the following values:
The model values are compared with the Parauay data set in Illustration 2.
Illustration 2: Richards'curve fitted to imperical Jatropha yield data from Paraguay
The average yield in the project area is estimated to be peak around 800 kg/ha. This may appear low considering the data in Illustration 1. However, poor management including poor pruning practise makes it unlikely that higher yields will be achieved.
Using the above mentioned parameters for the Richards' curve but with maximum yield set to 800 kg result in the following yield figures.
Illustration 3: Jatropha yield forecast for the project area (kg/ha)
Currently the equivalent of 250 ha have been planted. Below the total harvest from the project area is forecast for two scenarios:
The two scenarios result in the following harvest forecast.
Table 1: Harvest forecast for the project area under two scenarios
Illustration 4: Harvest forecast for the project area under two scenarios
11 Jatropha Yield Forecasting
|Selection||File type icon||File name||Description||Size||Revision||Time||User|
|Achten WMJ, Verchot L, Franken YJ, Mathijs E, Singh VP, Aerts R, Muys B 2008. Jatropha bio-diesel production and use (a literature review). Biomass and Bioenergy 32(12), 1063-1084. DOI: 10.1016/j.biombioe.2008.03.003||9 Nov 2010, 13:12||Banana hill Consultants|