African Journal of
Agricultural Research

  • Abbreviation: Afr. J. Agric. Res.
  • Language: English
  • ISSN: 1991-637X
  • DOI: 10.5897/AJAR
  • Start Year: 2006
  • Published Articles: 6914

Full Length Research Paper

Sensitive parameters for EPIC model evaluation and validity under soil water and nutrients limited conditions with NERICA cropping in West Africa

Omonlola Nadine Worou
  • Omonlola Nadine Worou
  • West African Science Service for Climate Change and Adapted Land Use.
  • Google Scholar
Thomas Gaiser
  • Thomas Gaiser
  • University of Bonn, Institute of Crop Science and Resource Conservation, Katzenburgweg 5, 53115 Bonn, Germany.
  • Google Scholar
Sylvester Oikeh
  • Sylvester Oikeh
  • African Agricultural Technology Foundation (AATF), P. O. Box 30709-00100, Nairobi, Kenya.
  • Google Scholar


  •  Received: 09 December 2014
  •  Accepted: 26 April 2015
  •  Published: 28 May 2015

 ABSTRACT

Crop models known to be based on the theory of crop physiology for describing the dynamic process of crop growth are recently explored for their uncertainties in model application under resource limited conditions. The aim of this study was to test Environmental Policy Integrated Climate (EPIC), on upland land rice production by taking into account seasonal variability in Guinean and Guinean-Sudanian zones in Benin and Nigeria (West Africa). A range of data available under farmer or experimental conditions in rainfed agriculture were measured or used from literature. The results show the accuracy of the model to simulate LAI, total above ground biomass and grain yield of upland rice for 2 NERICA rice cultivars. After calibration, the model showed average mean relative error between 0.06 and 0.15 with the model efficiency up to 0.98% in the case of LAI. The assessment of the model performances about sensitivity to N or P fertilizer application is also discussed under Ultisols. Large root mean square (RMSE) in calibration and the validation (>100) process suggested that robustness of the model became restrictive under severe environmental conditions such as in drought or flooding condition. Performance of the model at large scale should be executed of with land marginality classification.

 

Key words: Environmental Policy Integrated Climate (EPIC), modeling, upland rice, West Africa.


 INTRODUCTION

The operation of crop growth models is of interest for filling gap between information needed and that created by traditional experimental trials in soil and agronomic research or for extrapolating results gained on experimental stations leading to better integration of knowledge.
 
Beside, simulation modeling represents a research tool for assessing climatic change patterns and  their  impacts on crop growth and yield. Modeling a cropping system requires to understand the complex crop-water-soil interaction and to suggest some empirical parameters which can be applicable to diverse conditions and environments. However, the attempt to use crop growth   models under extremely unfavourable growth conditions that is, water scarcity combined with low soil fertility or with    indigenous   management   practices   remains    achallenge in tropical cropping systems such as in Africa or in Latin America (de Barros et al., 2004; Gaiser et al., 2010).
 
Indeed, for the single crop rice that has a relatively long history in modeling, model development is now geared to the issue of resources limitation due to expansion of rainfed rice systems. For instance the water and nitrogen modules in the latest version of ORYZA2000 formerly developed for estimating potential rice production suggest repeated model simulations with real-world data in order to increase the confidence in the suitability of the model for a certain purpose (Bouman and Van Laar, 2006). Even the agroecological system models such as Environmental Policy Integrated Climate (EPIC), which addresses crop simulation in response to weather, and nutrient cycling, is still not widely used to explore management strategies (Probert, 2004). As result, in rainfed low-input systems such as smallholder farms in West Africa, models developed for optimal management conditions fail to meet the needs of researchers and extension workers. That can be a key issue in Africa where about 80% of the rice production depends on rainfed conditions.
 
Although the basic use of crop models was to calculate crop growth and development for a single field, there is increasing interest in studies that concern multiple fields (Leenhardt et al., 2007; Hartkamp et al., 1999). This depends on the assumption that field scale model can be useful for evaluating management strategies at a broader scale. In rainfed upland systems in West Africa, rice yield is seldom above 2 Mgha-1. Constraints are shown to be in addition to rainfall uncertainty in West Africa, weeds and soil nitrogen availability. Indeed soil nutrient availability for upland rice culture has also been described to be related to land use and ecology (Becker and Johnson, 2001). It is then important for crop modeling targeted on upland rice to be tested on various environmental and management conditions to provide more confidence for further upscaling exercises. The adopted upland rice is NERICA cultivars that have been showed relatively high yieldsyield which vary in a controlled environment from 4.0 to 7.0 Mgha-1 (Akintayo et al., 2008). NERICAs are developed from interspecific crosses and tested for their ability to overcome drought (Asch et al., 2005) or to tolerate temporary inundation via flash flooding (Kawano et al., 2009) or for low nitrogen environment  (Saito and Futakuchi, 2009; Oikeh et al., 2008). New cultivars including NERICA 1 and NERICA 4 are assessed through the participatory varietal selection (PVS) in order to identify genotypes that perform well across or within a specific target environment.
 
The objective of this study was a multisite calibration and validation of the EPIC model for NERICA cultivars in contrasting ecological condition and the identification of site-specific model sensitive parameters. For this purpose we analyzed the sensitivity of the crop model to fertilizer and water input with data from experimental and on farm fields in Guinean and Sudanian zones in Benin and Nigeria (West Africa). The EPIC model was chosen due to its capacity to consider the effect of limiting water stress and of nutrients such as nitrogen and phosphorus on rice production (Gaiser et al., 2010). 


 MATERIALS AND METHODS

Study area
 
The model evaluation followed a calibration and validation process. Experimental data were collected from 8 experiments carried out in 2004-2010 in Benin and in Nigeria, West Africa (Table 1). The locations are listed from South to North: (The International Institute for Tropical Agriculture (IITA), Cotonou) (4), Tohoué (8), Ikenne (1, Nigeria), Niaouli (3), Bohicon (2), Kpakpazoumé (6), Tchankpéhoun (7) and Pingou (5). The calibration dataset is obtained from sites 1, 4, 6, 7 and 8. Validation plots concerned sites 2, 3 and 5.
 
The experiments ranged from the Guinean to the Sudanian Zone (Table 2) toward the North. The Guinean zone has two rainy seasons and the Sudanian Zone has a monomodal rainfall distribution. The annual precipitation is about 1100 mm with slightly declining rainfall in northward direction; there is regionally higher rainfall close to the Atacora mountain range for the case of Pingou and Tchankpéhoun (Röhrig, 2008). The length of growing season decreases also from South to North, from 250 to 130 days. In general, the rainfall distribution allows cultivation of two crops per year in the southern areas (Igué, 2000). 
 
Model data input and source
 
Crop management dates are summarized in Table 3. For the experiments in site 1 through 4 experimental layouts and number of replications are described and published in previous works (Oikeh et al., 2008; Saito and Futakuchi, 2009; Sokei et al., 2010; Koné et al., 2008). For 5 to 8 sites, the experimental design varies according to the area. Plots of 3 m x 15 m were used in Kpakpazoumé and Pingou. The farmland in Tohoué occupied 1250 m² and 5000 m² in Tchankpéhoun. In all experiments cultivar  ‘NERICA1’ or ‘NERICA4’ were used. 
 
Soil information was provided from soil profiles dug during the fallow period in 2009 for sites 8, 3, 2 and 4. Atchade (2006) reported chemical and physical characteristics of soil profiles in IITA, Niaouli and Bohicon (Cana Sud). The top soil properties (0 to 15 cm) were adapted according to Saito and Futakuchi (2009) at IITA. Two fields were used at IITA: one with low soil fertility (IITAlow) and the other with high soil fertility (IITAhigh). Soil data in Ikenne were obtained from Heuberger (1998).  The profile was described during the fallow period in Kpakpazoumé, Pingou and Tchankpéhoun in 2009.  Topsoil sample were randomly collected from the fields at 5 points of 0 to 20 cm depth along a profile down to the root-table. The samples were sieved to pass through a 2 mm mesh before analysis. The pH was determined using a soil-water ratio of 1:2. The organic carbon and organic N were analysed using the elemental analysis for Kpakpazoumé and Pingou (Fujine, 2014). The dichromate oxidation method of Walkley and Black (1934) was used for Kobli and Tanguiéta. Exchangeable bases (Mg, K, Ca, and Na) were extracted with 1 mol L-1 NH4 Acetate; Ca and Mg in the extract were measured using the atomic absorption spectrophotometer (AAS) while Na and K were determined by flame photometry. The potential cation exchangeable capacity was determined by extraction with 1 mol L-1 BaCl2. Ikenne, Niaouli have sandy textured topsoil. However, except Tohoué, all the sites have loamy to clayey subsoils (Alfisols and Ultisols). 
 
 
 
 
 
 
The soils were usually acid with low nitrogen content except in Bohicon and IITAhigh. The soil type was sandy to sandy loam (Table 2). Down to 15 cm, soil organic carbon level was classified in the order:
 
Bohicon>IITAhigh>Kpakpazoumé>Pingou>Ikenne>Tchankpéhoun>Niaouli>IITAlow> Tohoué. 
 
Daily meteorological data (maximum and minimum air temperature and global solar radiation) were collected from the synoptic weather station nearest to the fields (Table 2). For synoptic data in Ikenne, the model weather generator was used from FAO climate database LocClim for monthly mean for temperature. Solar radiation at Ikenne was derived from Apkadio and Etuk (2002) and the Hargreaves (Hargreaves and Samani, 1985) method was used for ETP estimation. For all other sites, Penman Monteith Method was applied. All these methods have been successfully tested for ETo estimation in (Rahimi et al., 2015; Valipour, 2014a, b, c;d; e; f;g;h, Valipour, 2015). Daily rainfall was retrieved from the closest rainfall gauge. 
 
Modeling with EPIC
 
The version 3060 of the EPIC model (Williams et al, 1990) was used to simulate rice productivity. EPIC is a field-based model designed to simulate crop production based on information about soil, crop rotation and management system. A full description is presented in the model documentation by Izaurralde et al (2006). Among various subroutines, the model considers N and P cycling by flows between inorganic and organic stocks. For N mineralization, EPIC couples C and N cycling in the soil. Simulated C and N compounds in EPIC are stored in either biomass, slow, or passive soil organic matter pools. Direct interaction is simulated between these pools as the function of soil moisture, temperature, nutrient content and clay content functions (Izaurralde et al., 2006; Gaiser et al., 2010). For P mineralization, the model uses the approach given in Jones et al. (1984). Two sources of mineralization are considered: the fresh organic P pool, associated with crop residue and microbial biomass, and the stable organic P pool, associated with the soil humus. The mineral P is then transferred among three pools: labile (which comprises fertilizer), active mineral, and stable mineral. Flow between the labile and active mineral pools is governed by the equilibrium equation that implies the mineral P flow, the amount in the active mineral P pool and P sorption coefficient defined as the fraction of fertilizer P remaining in the labile pool after the initial rapid phase of P sorption is completed.
 
Model evaluation
 
The evaluation of the model was firstly done by graphical presentation of the agreement between measured and calculated values for crop model by producing linear regressions between measured and simulated variables and calculating the coefficient of determination (R²) derived from 1:1 regression line where data are considered to meet independence assumption. The different comparison methods in Table 4 that highlight the feature of data and the model response were also used. The mean error (ME), the mean relative error (MRE), the mean absolute error (MAE) and the root mean square were presented where n is the sample number, x is the observed, y is the simulated value. The MRE is positive when the model in total overestimates compared to observation. The negative sign is related to underestimation. The root mean square error (RMSE) estimates the precision and reliability of the prediction for single yield estimation points. Model efficiency is used to assess the predictive power of the model taking into account the variability inside the observation data set. All these calculations were done in Microsoft Excel sheet.
 
 
 


 RESULTS AND DISCUSSION

Calibration of crop parameters
 
The calibration and validation  runs started with  a  warm up period of 8-9 years in order to stabilize the soil organic carbon pools in the model. The approach used for the calibration was to modify some initial values of the model parameters in order to iteratively fit simulation values as close as possible to the observed yield values. Therefore, we adjusted the default crop parameters for rice to the NERICA cultivars. However no cultivar distinction was taken into account in the crop file. The NERICA 1 and 4 passport data published by the Africa Rice Centre represented no feature for distinguishing the two cultivars in the crop file of the EPIC model such as the number of days to maturity which determine the Potential Heat Unit (PHU) or flowering age.
 
In the process of LAI calibration, the parameters DLAP1 and the DLAP2 were used to control the crop growth. The DLAP1 was changed from 30.01 to 30.20 and the DLAP2 from 70.95 to 60.95 for the two cultivar cultivars.  The plant population density has been also modified from 12600 to 50.600 in PPC1 and 250.600 to 250.900 in PPC2. Félix (2006) considered that the sub-model of EPIC for LAI development is based on a strong amount of empiricism, as the mechanism that controls the rate of development of LAI is not yet well understood. 
 
As result, the model outputs and the observations with regard to the LAI before and after the calibration were graphically compared. Figure 1 shows that the model first underestimated the values of LAI with a negative mean error of -0.22 (Table 6). After calibration the average relative difference between the observed values and simulated LAI was approximately 6% with a model efficiency of 98%. The LAI development was rather satisfactory calibrated similar to Yoshida et al (2007) using a complex and detailed phenological model as a function of relative crop growth rate, leaf nitrogen content and air temperature.
 
 
The LAI was estimated under full irrigation at relatively high soil fertility level compared to farmers field as described in Table 2 (Corg =19.6 gkg-1 and total nitrogen up to 2.2 gkg-1). The observed value is average of 5 cultivars including NERICA1 grown under high soil fertility conditions (Saito and Futakuchi, 2009). The authors detected no difference in rice cultivars in LAI at 42 and 56 DAS and no traits from the early vegetative stage were observed related to grain yield. The relative increase of measured LAI at mid stage (DLAP1) compared to the default value in the model is in line with the high weed competitiveness feature demonstrated for NERICA lines (Ekeleme et al., 2009).
 
Model sensitivity for soil parameter
 
Before calibration, the model showed low sensitivity to the supply of inorganic N and P on a highly weathered and strongly acid low-activity clay soils at Ikenne (Figure 2), as the experimental layout was made to test the effect of fertilizer application in the humid forest zone on Ultisols (Table 3). Leenhardt et al. (2006) suggested the use of pedotransfer functions to estimate soil properties during the simulation process as solution for unavailable data. However, Gaiser et al. (2010) using sensitivity analysis estimated the fraction of microbial biomass across some different soil types under cropland in West Africa. The fraction of biomass in the soil organic matter pool (FBM) triggers the mineralisation of soil nitrogen, which is the main growth constraint in low-input small-holder systems in West Africa. The authors set a value of FBM to 0.01 instead of 0.04 that is more representative for soils with high organic matter content (Niu et al., 2009). The recommended value of 0.01 was then used for all sites. 
 
In addition, the fraction of humus in the passive pool expresses the proportion of carbon (and nitrogen) in the soil organic matter pool that has a low turnover rate. It was set to 0.99 making less nitrogen available to the plant, thus generating more response of the crop to additional nitrogen supply. 
 
More sensitivity of yield to fertilizer P application in the model was found when initial labile phosphorus concentration in the first layer for the acid Ultisol was set to a value of 0.05 ppm. Labile phosphorus (CSP) is considered to have contributions to be correlated to P uptake (Sharpley, 1985). The labile P concentration factor allows optimum uptake rates when CSP  is  above 20 ppm which was the default value used as critical labile P concentrations for a range of crops and soils. The soils in Ikenne are classified by USDA as Ultisols: They are considered to be low in CEC and bases due to the translocation of the clay in the subsoil and high leaching. They present a high P sorption to Fe and Al-hydroxides in the subsoil (Mokwunye, 1979) or kaolinite in the clay fraction (Wisawapipat et al., 2009). Daroub et al. (2004) in developing a soil-plant P model for highly weathered soils recorded for maize an overestimation of the P uptake by the model. Apparently, their model was not able to reproduce P fixation which is much higher than in less-acid soils developing in temperate climates.
 
 
The analysis of rainfed upland system refers also to evaluation of the water availability with depends on soil texture. The coarse fragments influence soil physical hydraulic properties. In EPIC, the role of this parameter addresses directly to the water erosion engine but it has soil functioning oriented for estimation of water storage capacity at the same stand as the bulk density. In fact Chow et al. (1997) observed by incorporating 10 to 30% coarse fragments into the plow layer of the Northen American Podsol, it increased significantly the soil bulk density and this increase reduced the porosity and soil water-holding capacity. The sensitivity analysis of coarse fragments content was done in 4 sites where substantial coarse fragments were identified in soil profile to show influence of this parameter the grain yield. 
 
Figure 3 shows at 2 to 3 soil layers over 4 sites (Bohicon, Kpakpazoumé, Pingou and Tchankpéhoun) variation in coarse fragments content from 0 to 80%. It appeared that   strong   influence   of   coarse   fragment content was obtained when all the layers are concerned by the limitation and the upper layers showed the higher sensitivity for grain yield.
 
 
Simulation of total aboveground biomass and grain yield
 
Figure 2 shows that the model reflects after calibration the effect of N and P application on NERICA yield on Ultisols when P and N are limiting. This is in accordance with Nigerian forest agroecosystems where high split application of 90 to 120 kg N ha–1 has been recommended for the rice cultivars to optimize yields (Enwezor et al., 1989). The model results in 7 over 8 treatments were not significantly different from the observed yield of NERICA. The P stress has been simulated adequately to allow the expression of nitrogen stress among the treatments with application and without application of P. By simulating adequately the processes in P deficient soils, the model agrees with the results of Sahrawat et al. (1995) suggesting that P fertilization of acid-tolerant upland rice can significantly improve its productivity under Ultisols. Before the calibration and over all five sites, grain yield and total aboveground  plant biomass had a RMSE of more than 30% (Table 6). Farmland fields contributed the most to overestimation of model by 88% for 2 years on average (Figure 4).  For the calibration in Tchankpéhoun, the plant density was reduced from the theoretical plant population to the measured plant density at maturity. During the two seasons, heterogeneous planting density translated by the missing of plant hills led to the reduction of the total yield observed.  Affholder (2001) pointed out that a model developed in high input environment such as the US where the planting density is very homogeneous need numerous modifications to be applied under the conditions of West Africa where the variability of the densities of sowing is a big factor of the variability of the productivities. Oikeh et al. (2009) did not show relationship between the grain yield and NERICA density whereas density effects appeared only for tiller and panicle densities. This study was undertaken with seasonal differences in rainfall distribution and moisture availability that may reduce simple effects of N and spacing (plant density). Tchankpéhoun got adequate monomodal rainfall supply for the two years. After calibration, the goodness of fit of the model was improved for both total aboveground biomass and the grain yield (Figure 4). Lower RMSE after calibration indicated that a higher fraction of the measured variations were accounted by the model (Table 6).
 
 
 

Model validation

 

The calibration of the EPIC model for upland rice was focused on nitrogen and phosphorus as main constraints to crop growth. The validation was carried out on three sites (Niaouli, Bohicon and Pingou) over two seasons. At Niaouli the experiment tested different levels of N and P input, at Bohicon NPK application was tested and Pingou was an on-farm field experiment (reference in Table 3).

 

The validation of the model showed that a relatively high gap between averages simulated and observed yield (Table 7). The mean error was 1.2 Mg ha-1 whereas the mean relative error was 3.0 Mg ha-1, which showed a very large overestimation of the simulated yields at plot level. The variation of the individual plots was also quite high resulting in root mean square error  of  more than 100%. The observed grain mean grain yield was lower than the average in the calibration suggesting various stress effects. Indeed some causes of rice failure attributed to floods and drought were reported b for NERICA evaluation on 5 locations with similar pedoclimatic conditions to these experiments in Benin republic (JAICAF, 2007). Therefore, before the use of the model in the assessment of impacts of and adaptations to climate variability and climate change in spatial studies, there is still a need for improvement in the amount and quality of available data collection.

 

Figure 5 showed a scatter plot of the observed and estimated value of sites used for model validation. The average yield in plots used for validation was relatively low, this is due to crop failure in 2006 in Niaouli where the average yield was below 1 Mgha-1 leading to the model overestimation. In fact, the experimental design was originally set up to evaluate the tolerance to drought with nutrients application for NERICA cultivar. Niaouli is located in the sub humid zone with  bimodal  rainfall   pattern.   The   mid-season rainfall pattern associated with the sandy topsoil texture induced severe drought stress. The soil type “terre de barre” was described by Azontonde (1991) as soil with good physical hydraulic characteristics but with low water storage and their structure can be rapidly destroyed when there is no proper technique for maintaining organic matter.

 

The sensitivity of NERICA cultivar to water stress is well documented. Akinbile et al. (2007) showed that with NERICAs yield decreased under optimal satisfactory conditions almost linearly with evapotranspiration thus indicating that water application remained dominant factor at all the stages of production. In EPIC model, the potential harvest index is adjusted daily according to water stress suffered by the crop (Williams, 1995). During the calibration, the sensitivity of model was increased by setting the water stress impact (WSYF parameter), which allowed harvest index to drop to 0.01 in case of severe drought. The effect of water stress could be in fact limited to HI reduction. Fuji  et  al.  (2004)  reported  that some Nerica lines showed high dry matter production under drought condition among other rice cultivars, and this have been correlated with stomata conductance (r=0.63**).

 

 

However, intensive rains of short duration followed by long dry spells that occurred during the flowering period which lead to increased sterility and decrease in grain weight (O´Toole and Moya, 1981).  De Barros et al. (2005) observed slight overestimation of grain yield by the EPIC simulations was attributed to high rates of floral abortion caused by dry spells during the flowering periods since this factor is not considered in the model. 

 

 

After removing the plots with crop failure induced by drought in 2006 in Niaouli, the goodness of the fit of the model was improved from 0.01 to 0.51. Table 8 lists the simulation results of the remaining plots. Pingou siteo held in 2010 the lowest yield below 1 Mgha-1. In this year it was observed a shallow groundwater during the wet season at sowing (end of July) and was followed by transplanting. Therefore the first possibility for the model overestimation is that the model could not consider transplanting shock that caused a delay in phenological development resulting in reduced vegetation period in the field. No reported analyses on negative impact of flooding on upland NERICA were available. In controversy Fofana (2008) highlighted the recovery and the improvement possibility of NERICA1 production after short and intense moisture stress at the seedling emergence stage. The presence of ferric cuirasses in Pingou may result in low saturated  conductivity  at  the  middle   soil   depth,   thus increasing the submergence and runoff risk. The relatively high soil moisture should have been the explanation of low yield due to direct seeding in 2009. Indeed Ogunremi et al. (1986) demonstrated direct-seeded rice was adversely affected by the transient flooding conditions during the seedling stage on Ultilisol in Southern Nigeria. The obtained grain yield decreased indeed with increasing penetrometer resistance.
 
The tendancy of overestimation remained of yield response to fertilizer was observed in Bohicon and Niaouli. Even at Niaouli in 2005, where experiment also applied relatively high amount of NP (100 kg ha-1) that were less efficient in observation than in modeling.  Under limiting water conditions there could be less capacity of crop to continue water uptake that can probably reduce transport to the roots through mass flow. Indeed some traits of upland rice (japonica type) related to less adventious roots per hill result in relatively week ability in N uptake (Zhang, 2008). In addition, a severe drought that occurred just after the application of the first split of N could have induced urea volatilization resulting from the lack of N dissolution Oikeh et al. (2008). 
 
From data drawn from the upland experiments, the EPIC model was calibrated and parameterized for a multisite evaluation, which is particularly important for rice production because of its high dependency to nutrients and water. In the model validation the variations of individual plots also were higher than the observed. The model overestimated the yield under drought condition in the site of Niaouli. The model indeed estimated effective lower water retention in the first soil depths by the use of texture. Figure 6 shows the relatively good agreement of the model estimation for water storage capacity in the Niaouli site that allowed condition for water stress experiment. However, it has been also reported in the sensitivity analysis of the model that coarse fragment content had more or less high influence on water storage capacity in the soil layers. In four sites out of eight sites, the model was parameterized with coarse fragment content limiting water storage capacity in different layers.  
 
Adapted value was needed in case of Bohicon while the model still overestimated the water retention capacity with the measured variables. The results of final calibration are shown in Figure 7. Indeed, the model estimation of water storage required modification in Bohicon and reduced the yield gap between observed and simulated value from 3 to 2 Mgha-1 on average. As consequence, when simulating rainfed rice it appear to be a prerequisite to provide detailed site-specific soil input parameters including water retention among soil physical characteristics. 
 
At the level of discussion of plant nutrition causes on the overestimation of the model, if the model was adequately tested to respond to NP fertilizer from calibration there is still issue of the simulation of micronutrients. Several experiments conducted on highly weathered soils in Africa showed that as sufficient N and P is applied in maize, micro nutrient deficiencies may appear (Gaiser et al., 1999). Voortman et al. (2000) estimated micronutrients deficiencies on about 60% of the cropland in sub-Saharan Africa. This confirms the requirement to improve micronutrient effects in crop models.  Furthermore, Koné et al. (2009) proved that there was also a significant (P = 0.004) decreasing effect of Zn (28%), N (34%) and K (36%) exclusion on the mean grain yield in the Ferralsol soils Benin. These results attested the existence of Zn and K deficiencies which may reduce the sustainability of upland rice production.  
 
In the study  through  the  modeling  tool,  interspecific genotypes were evaluated under with local farmer’s conditions which include in addition to low inherent soil fertility, the occurrence of drought or flood. This multisite evaluation did not show any type of interaction between variety and environment that should be seen in the simulation outputs. Several kind of stress may limit varietal selection progress under unfavourable environments (Banziger and Cooper, 2001). In fact, Mandel et al. (2010) found smaller genetic variance for grain yield under low-input conditions in India. 
 
The higher effect of soil texture input on grain yield confirmed the importance of water capacity retention estimation as a key  soil parameter indicating that management of water supply (bunds building or reduction drainage) are likely to be useful for improving rice production.
 

 

 

 

 

 
 
 
 
 
 
 
 
 
 
 
 


 CONCLUSION

The EPIC model was tested with mix of collected and secondary data from an experimental station, on farm research field and farmland in the tropical humid with P fixing soils to semi-humid zones. The study showed the relevant crop and soil parameters for calibration in upland rice cropping conditions, in particular different pools of nitrogen and phosphorus in the soil. In calibration, the model presents a good response to the application of N and P mineral fertilizers. A multiple year calibration for multiple variables such as plant biomass, leaf area index and yield improve our confidence in the model calibration. Although the EPIC shows the sensitivity of rice to seasonal rainfall, its robustness under severe water stress become limited. With a multiple year calibration for multiple variables such as plant biomass, leaf area index and yield, the uncertainty in the model prediction in validation is related at first place, to the lack or quality on input data (estimation of impact drought spells on grain yield) and secondly on the model weakness (reduction of HI linked to crop phenology and fertilizer responsiveness under low input environments). 


 CONFLICT OF INTEREST

The authors have not declared any conflict of interest.



 REFERENCES

Akinbile CO, Sangodoyin AY, Nwilene FE (2007). Growth and Yield Responses of Upland Rice (NERICA) under Different Water Regimes in Ibadan, Nigeria. J. Appl. Irrigat. Sci. 42:199-206.
 
Akintayo I, Cisse B, Zadji L (2008). Guide pratique de la culture des NERICA de plateau Centre du riz pour l'Afrique (ADRAO). Juillet 2008; P. 26.
 
Affholder F (2001). Modélisation de culture et diagnostic agronomique régional Mise en point d’une méthode au cas du mais chez les petits producteurs du Brésil central Thèse de doctorat s l’institut national agronomique Paris-Grigon
 
Asch F, Dingkuhn M, Sow A, Audebert A (2005). Drought-induced changes in rooting patterns and assimilate partitioning between root and shoot in upland rice. Field Crops Res. 93:223–236.
Crossref
 
Atchade NS (2005). Caractérisation des stations pour la riziculture pluviale au sud-Bénin. Cas de IITA (Godomey), Niaouli et Bohicon. Master thesis. University of Abomey Clavi. P. 88.
 
Azontonde HA (1991). Propriétés physiques et hydrauliques des sols au Bénin. Soil Water Balance in the SudanoSahelian Zone. Proceedings of the Niamey Workshop, February 1991. IAHS Publ. P. 199.
 
Banziger M, Cooper M (2001). Breeding for low-input conditions and consequences for participatory plant breeding: examples from tropical maize and wheat. Euphytica 122:503–519.
Crossref
 
de Barros I, Williams JR, Gaiser T (2005). Modeling soil nutrient limitations to crop production in semiarid NE of Brazil with a modified EPIC version II: Field test of the model. Ecol. Model. 181:567–580.
 
Becker M, Johnson DE (2001). Cropping intensity effects on upland rice yield and sustainability in West Africa. Nutr. Cycl. Agrosyst. 59:107-117.
 
Bouman BAM,Van Laar HH (2006). Description and evaluation of the rice growth model ORYZA2000 under nitrogen-limited conditions. Agric. Syst. 87:249-273. 
Crossref
 
Chow TL, Rees HW, Monteith JO, Toner P, Lavoie J (2007). Effects of coarse fragment content on soil physical properties, soil erosion and potato production. Can. J. Soil Sci. 87:565–577.
Crossref
 
Daroub SH, Gerakis A, Ritchie JT, Friesen DK, Ryan J (2004). Development of a soil-plant Phosphorus simulation model for calcareous and weathered tropical soils. Agric. Syst. 76:1157–1181.
Crossref
 
Ekeleme F, Kamara AY, Oikeh SO, Omoigui LO, Amaza P, Chikoye D (2009). Response of upland rice cultivars to weed competition in the savannas of West Africa. Crop Prot. 28:90–96.
Crossref
 
Enwezor WO, Udo EJ, Usoroh NJ, Ayotade KA, Adepeju JA, Chude VO, Udegbe CL (1989). Fertilizer use and management practices for crops in Nigeria (Series No. 2). Fertilizer Procurement and Distribution Division. Federal Ministry of Agriculture, Water Resources and Rural Development, Lagos, Nigeria.
 
FAO (2002) Local climate estimator, LocClim 1.0 CD-ROM.
 
Félix R (2006). Optimisation de l'irrigation et de la fertilisation azotée et la recherche d'une stratégie limitant la pollution des nappes souterraines par les nitrates. PhD thesis. Faculty of Agronomy, Gembloux
 
Fofana M (2008). Effets de la sécheresse sur les paramètres agronomiques et la qualité des grains du riz (Oryza sp.) Présentée en vue de l'Obtention du grade de Docteur de l'Université de Lomé. P. 161.
 
Fuji M, Andoh C, Ishihara S (2004) Drought resistance of Nerica (New Rice for Africa) compared with Oryza sativa L. and millet evaluates by stomatal conductance and soil water content. New direction on a diverse planet: Proceedings of the 4th International Crop Science Congress. Brisbane, Australia. ISBN 1 920842209
 
Fujine K (2014). Elemental Analyzer (CHNS): User Guide Manual Information, International Ocean Discovery Program FlashEA 1112, P. 12.
 
Gaiser T, de Barros I, Sereke F, Lange FM (2010). Validation and reliability of the EPIC model to simulate maize production in small-holder farming systems in tropical sub-humid West Africa and semi-arid Brazil. Agriculture, Ecosyst. Environ.135:318–327
Crossref
 
Hargreaves GH, Samani ZA (1985). Reference crop evapotranspiration from temperature. Appl. Eng. Agric. 1:96–99.
Crossref
 
Hartkamp AD, White JW, Hoogenboom G (1999). Interfacing geographic information system with agronomic modeling: a review. Agron. J. 91:761-77
Crossref
 
Heuberger H (1998). Nitrogen efficiency in tropical maize (Zea mays L.)- indirect selection criteria with special emphasis on morphological root characteristics. PhD thesis university of Hannover. P. 182.
 
Igué AM (2000). The use of a soil and terrain database for land evaluation procedures: Case study of Central Benin. Hohenheimer Bodenkundliche Hefte 58. University of Hohenheim. P. 235.
 
Izaurralde RC, Williams JR, McGill WB, Rosenberg NJ, Quiroga Jakas MC (2006). Simulating soil C dynamics with EPIC: Model description and testing against long-term data. Ecol. Model. 192: 362–384.
Crossref
 
Jaicaf (2007). Summary of Jaicaf –Inrab, Cooperative NERICA Trials, the Wet Season, 2007, in Benin. Cotonou Benin P. 7.
 
Jones CA (1984). Estimation of percent aluminum saturation from soil chemical data. Commun. Soil Sci. Plant Anal. 15:327-33
Crossref
 
Kawano N, Ito O, Sakagami JI (2009). Flash flooding resistance of rice genotypes of Oryza sativa L., O. glaberrima Steud., and Interspecific hybridization progeny. Environ. Exp. Bot. 63:9–18.
Crossref
 
Koné B, Ettien JB, Amadji G, Diatta S (2008). Caractérisation de la tolérance de NERICA à la sécheresse de mi-saison en riziculture pluviale. Afr. Crop Sci. J. 16:133-114.
 
Leenhardt D, Wallach D, Le Moigne P, Guérif, M, Bruand A, Casterad MA (2007). Using crop models for multiple fields. In: Wallach, D, Makowski D, Jones JW, working with dynamic crop models, evaluation, analysis, parameterization and applications. Elservier, UK.
 
Mandel NP, Sinha PK, Variar M, Shukla VD, Perraju P, Mehta A, Pathak AR, Dwivedi JL, Rathi SPS, Bhandarkar S, Singh BN, Singh DN, Panda S, Mishra NC, Singh YV, Pandya R, Singh MK, Sanger RBS, Bhatt JC, Sharma RK, Raman A, Kumar A, Atlin G (2010). Implications of genotype×input interactions in breeding superior genotypes for favorable and unfavorable rainfed upland environments. Field Crops Res. 118:135–144.
Crossref
 
Niu X, Easterling W, Hays C J, Jacobs A, Mearns L (2009.) Reliability and input-data induced uncertainty of the EPIC model to estimate climate change impact on sorghum yields in the U.S. Great Plains. Agric. Ecosyst. Environ. 129:268–276.
Crossref
 
Ogunremi LT, Lal R, Babalola O (1986). Effects of tillage and seeding methods on soil physical properties and yield of upland rice for an Ultisols in southeast Nigeria. Soil. Tillage Res. 6:305-324.
 
Oikeh SO, Nwilene F, Diatta S, Osiname O, Touré A, Okeleye KA (2008). Responses of upland NERICA rice to nitrogen and phosphorus in forest agroecosystems. Agron. J. 100:735-741.
Crossref
 
Oikeh SO, Touré A, Sidibé B, Niang A, Semon M, Sokei Y, Mariko, M (2009). Responses of upland NERICA rice cultivars to nitrogen and plant density. Arch. Agron. Soil Sci. 55:301-314.
Crossref
 
O’Toole JC, Moya TB (1981). Water deficits and yield in upland rice. Field Crops Res. 4:247-259.
Crossref
 
Probert ME (2004).  A Capability in APSIM to Model. Phosphorus Responses in Crops. CSIRO.
 
Rahimi S, Gholami Sefidkouhi MA, Raeini-Sarjaz M, Valipour M (2015). Estimation of actual evapotranspiration by using MODIS images (a case study: Tajan catchment). Arch. Agron. Soil Sci. 61(5):695–709.
Crossref
 
Röhrig J (2008). Evaluation of agricultural land resources in Benin by regionalisation of the marginality index using satellite data. PhD thesis, University of Bonn. P. 174.
 
Sahrawat KL, Jones MP, Diatta S (1995). Response of upland rice to phosphorus in an Ultisol in humid forest zone, West Africa. Fertil. Res. 41:11-17.
Crossref
 
Saito K Futakuchi K (2009). Performance of diverse upland cultivars in low and high soil fertility conditions in West Africa. Field Crops Res. 111:243-250.
Crossref
 
Sharpley AN (1985). The selective erosion of plant nutrients in runoff. Soil Sci. Soc. Am. J. 49:1527–1534.
Crossref
 
Sokei Y, Akintayo I, Doumbia Y, Gibba A, Keita S, Assigbe P (2010). Growth and yield performance of upland NERICA cultivars in West Africa . Japanese J. Crop Sci. 79:2–3.
 
Valipour M (2014a). Application of new mass transfer formulae for computation of evapotranspiration. J. Appl. Water Eng. Res. 2(1):33–46.
Crossref
 
Valipour M. (2014b). Study of different climatic conditions to assess the role of solar radiation in reference crop evapotranspiration equations. Arch. Agron. Soil Sci. 61:5.
 
Valipour M (2014c). Evaluation of radiation methods to study potential evapotranspiration of 31 provinces. Meteorol. Atmos. Phys. 127(3):289-303.
Crossref
 
Valipour M (2014d). Investigation of Valiantzas' evapotranspiration equation in Iran. Theor. Appl. Climatol.
Crossref
 
Valipour M (2014e). Temperature analysis of reference evapotranspiration models. Meteorol. Appl. doi:10.1002/met.1465
Crossref
 
Valipour M (2014f). Evaluation of radiation methods to study potential evapotranspiration of 31 provinces. Meteorol. AtmosPhys.
Crossref
 
Valipour M (2014g). Assessment of different equations to estimate potential evapotranspiration versus FAO Penman-Monteith method. Acta Advan. Agric. Sci. 2(11):14–27.
 
Valipour M (2014h). Use of average data of 181 synoptic stations for estimation of reference crop evapotranspiration by temperature based methods. Water Resour. Manage. 28(12):4237–4255.
Crossref
 
Valipour M (2015). Importance of solar radiation, temperature, relative humidity, and wind speed for calculation of reference evapotranspiration. Arch. Agron. Soil Sci. 61(2):239–255.
 
Walkley A, Black IA (1934). An examination of the Degtjareff method for determining organic carbon in soils: Effect of variations in digestion conditions and of inorganic soil constituents Soil Sci. 63:251–263.
Crossref
 
Williams JR, Jones CA, Dyke PT (1990). The EPIC Model. In: Williams, J.R. (Ed.), EPIC—Erosion Productivity Impact Calculator. 1. Model Documentation. U.S. Depart. Agric. Technol. Bull. 1768:3–86.
 
Williams, JR (1995). The EPIC model. In: V.P. Singh, Editor, Computer Models of Watershed Hydrology, Water Resources Publications, Highlands Ranch, Colorado.
 
Wisawapipat W, Kheoruenromne I, Suddhiprakarn, A, Gilkes, RJ, (2009). Phosphate sorption and desorption by Thai upland soils. Geoderma, 153:408-415.
Crossref

 




          */?>