Recommendations of Fertilizer Formulas for the Maize Production in Northern Benin

Maize cultivation under soil conditions in Benin requires high quantity of nutrients. There is therefore a need to develop adequate fertilizer recommendations in order to achieve the level of productivity that could meet the needs of the increasing population in the rural area. The present study aims to update the mineral fertilizer formulas recommended for maize production in northern Benin. An experimental program was carried out in the year 2012 on three main soil types: ferric Luvisols, gleyic Luvisols and eutric Gleysols in two agroecolological zones of Northern Benin. The experimental design was a randomized completed bloc with four replicates, installed in farmers’ fields with the specific objective to validate five N, P, K based fertilizer formulas. The maize variety EVDT-97 STRW was used. Biophysical and economic analyses completed using the seasonal stool of the DSSAT model allowed to identify a series of efficient options. The results of variance analyses relating to the effect of different fertilizer formulas on maize grain yields showed that the rate simulated by the DSSAT model (115-30-75) produced the highest grain yields regardless of the soil types and agro-ecological zones. The ratio of observed-to-simulated values are close to 1 and the mean standard prediction error (NRMSE) between the observed and the simulated yields was comprised between 11% and 20% for gleyic Luvisols but between 21% and 30% for the other soil types. The results of the biophysical and economic analysis showed that N88P30K35 was the most efficient fertilizer formula for sustainable maize production in Northern Benin.


Introduction
Agriculture in Sub-Saharan Africa is characterized by low productivity due to a steady decline in soil fertility [1,2]. According to Douthwaite et al. one of the major constraints for agriculture in sub-Saharan Africa is the steady decline in soil fertility [3]. In developing countries, particularly in sub-Saharan Africa, the environment is not subject to excessive use of mineral fertilizers, but rather to very low or even non-use of fertilizers to compensate for crop exports. This has led to a decrease in soil fertility and therefore to a decrease in agricultural potential [4]. According to Kanté, the solution to the widespread depreciation of the "natural capital" and the decline in the production capacity of sub-Saharan Africa lands necessarily involves investments in soil fertility [5].
According to the author, to be sustainable, actions to improve soil fertility must be multifaceted, to take into account the existing diversity between agro-ecological and socio-economic situations. Serpentie and Ouattara, emphasize the notion of sustainability relating to soil fertility [6]. Thus, mineral fertilization is one of the soil improvement solutions proposed to compensate for nutrient losses and nutritional deficiencies observed in the production systems. In Benin, low crop yields are often due to unfavorable rainfall conditions, inherent soil nutrient deficiency and low use of external inputs [7].
Climate variability and land degradation are the main constraints limiting maize production in Benin [2,8,9]. The main causes of this land degradation are low organic matter content, the low use of fertilizer, poor soil fertility management practices and monocropping [2,8,10]. This is reflected in the negative nutrient balances observed on soils [11]. In Benin, fertilizer use, as in many other countries of West Africa, has been promoted to intensify crop production [10].
Indeed, excessive and inappropriate use of tillage equipment at farm level, export of crop residues and the shortening of fallow periods have created the conditions leading to the decrease in soil organic matter content and the degradation of their structure [11] and as a result to the decline of their fertility [12]. The tropical ferruginous soils (Luvisol hahlique, Luvisol gleyique, Plutthosol eutrique, Arenosol haplique) which occupy 60% of the total surface area of Benin [13] and to a lesser extent the hydromorphic soils of the depressed zones are clearly affected by this. These soils are known to be low in nitrogen and phosphorus [14].
Any strategy aiming at correcting this situation requires a sound management of agricultural lands. This involves the application of

Advances in Crop Science and Technology
mineral and organic fertilizers to restore the stock of organic matter and ensure nutrient supply to the soils. The role of organic matter in improving soil quality is widely recognized. Indeed, organic matter represents the main indicator and plays a determining role in biological activity [12].
In Benin, the fertilizer rates and formulas currently used by extension services (150 kg/ha of NPK 14-23-14 and 50 kg/ha of urea) for maize are mostly obsolete and generalized [9]. These same fertilizer rates are recommended for all agroecological zones within the country. Such practices do not take into account soil types and the specificity of farmers' cropping systems and farm ecology. These recommended standard fertilizer rates are old. Therefore, there is a need to update this fertilizer recommendation for maize production regarding each agroecological zone of Benin, soil types, and the economic profitability for the farmer.
In order to improve not only land productivity, but also maize productivity through the optimization of fertilizer use in North Benin, new site-specific fertilizer recommendations (adapted to the soil potential and optimal sowing dates related to climate potential) are necessary. Agricultural simulation models are one way to predict yield components in various agroecosystem to save time and reduce field trials [10]. Agricultural simulation models are originally developed, calibrated and validated under different agroecological conditions, and their application in other specific conditions does guarantee reliability [15]. The present research was carried out in the framework of the IFDC-Africa fertilizer research program in West Africa.
This study aims • To characterize the inherent fertility status of concretioned tropical ferruginous soils (Luvisol ferrique), modal ferruginous soils (Luvisol haplique) and hydromorphic soils (Gleysol eutrique), in the communes of Tanguiéta and Banikoara. • To determine on-farm the fertilizer rate recommended to achieve the best maize grain yields depending on soil types and agroecological zones. • To evaluate the added value of a combined application of mineral and organic fertilizers on the three types of soil studied.

Study environment
Trials for validating the recommended options using the DSSAT model were performed in two villages, one in Atacora and one in Alibori. Producers were selected from two soil units per village. These villages are located in two different agro-ecological zones: Zone 2 (Cotton Zone of North Benin) and Zone 4 (West Zone of Atacora) in Benin ( Figure 1).
The cotton zone of North Benin is a Sudanian zone with two contrasting seasons (a rainy season from June to October and a dry season from November to May). Rainfall varies between 800 and 1200 mm/year. The food crops grown are maize, sorghum, yam, cowpea; cotton is the industrial crop while perennial crops are mango and cashew nuts. Plant growth period is between 140 and 180 days. Relative humidity varies during the year according to temperatures maxima. Rains are heavier at the beginning of the season because of their stormy character and especially the absence of vegetal cover. The relief is a vast peneplain slightly developed and hardly undulating (slope between 1 to 4%) integrating mounds in a tabular form, increasingly high and in increasing numbers moving towards the Niger River. A trial was carried out in northeastern Benin, in the village of Arbonga, Banikoara district located between 10°50' and 11°30' north latitude and between 2° and 2°40' longitude east. The study area is characterized by a Sudano-Guinean climate with a long dry period and a single rainy season. The monthly averages vary between 2 mm and 280 mm (August) of rain. Rainfall varies widely from one year to another and during the vegetative period. The average temperature during the year is 27.4°C. The relative humidity varies according to the temperatures maximum (33.9°C). This study area is dominated by model tropical ferruginous soils and concretionned ferruginous soils [16]. The pedological study of the Banikoara district in the Banikoara commune allowed, using the topo sequential method, to distinguish eight soil types according to the French classification [17,18].
Agro-ecological zone 4 (zone West-Atacora) is characterized by a climatic variation of Sudano-Sahelian to Sudano-Guinean with an annual rainfall of 1000 to 1300 mm. Soils are also ferruginous, often deep, but with low water reserve. The vegetative period is between 160 and 220 days. In zone 4 the climate is very contrasted: in the west the dry season is 5 months in Natitingou and can reach seven months in Porga, in the central area the dry season also lasts seven months and the rainy season from June to September. In the eastern part, the two seasons are roughly equivalent.
The area studied is located approximately 10 km from Tanguiéta which is about 592 km from Cotonou. It is between 10°40' and 10°45' north latitude and between 1°20' and 1°22' east longitude. The climate is of the Sudano-Guinean type with a long dry season and a single rainy season. The soils of the Nanébou region in the Tanguiéta commune have very variable morphological and agronomic characteristics. This variability resulting from the heterogeneity of the parent rock, the diversity of the topographic positions and the pedological differentiation along the topo sequences, is reflected at the scale of the site mapping by the existence of combinations of soils rather than homogeneous units. The soil survey performed at Nanébou in north-west Benin allowed, using the topo sequential method, to distinguish seven types of soils according to the French classification CPCS and FAO.
These two agro-ecological zones of the study area correspond to the sub-humid agro-ecological zone of West Africa (IFDC and AFAP 2015).
The results of soil fertility evaluation showed that half of the soils in the study area are deficient in phosphorus and potassium, while one third is deficient in organic matter and nitrogen [16]. In all soils, the cation exchange capacity (CEC) is the major constraint. The use of organic matter raises this CEC in the soil. It also increases soils' nitrogen content. To correct soil phosphorus and potassium deficiencies, the use of phosphate and potash fertilizers as nutrient supplements is required.
The assessment of fertility status and class indicates that all soils require phosphate and potash fertilizers except hydromorphic leached tropical ferruginous soils (Gleyic Luvisols) which are not prevalent in the area studied. These soils have generally high to moderate levels of nitrogen and organic matter. The lower slopes and lowlands have the highest nutrient contents. It should be noted that almost all soils have severe to very severe limitations in terms of the sum of exchangeable bases and cation exchange capacity. This is probably due to the nature of the rocks on which these soils were formed. All soils belong to the low to very low fertility class except the concretionned ferruginous soils (Ferric Luvisols) which are of the medium fertility class. Table 1 presents the locations of the validation sites for the options and the types of soil on which the trials were carried out.

Plant material
The plant material used in this study is the EVDT 97 STRW, which is a 90-day open-pollinated (composite) maize variety. The ear coverage is good enough. The grains are white, half-toothed, halfstarchy and half-vitreous. Yields in the farming environment vary between 2 and 4 t/ha, while the potential yield is 6 t/ha. This maize variety is highly appreciated by producers in Benin [19].

Simulation of maize growth and development
Decision Support System for Agrotechnology Transfer (DSSAT v 4.5) was used for the simulations. This model requires a minimum of input data that can be grouped into three categories: daily climatic data (maximum temperature, minimum temperature, precipitation, insolation), site information (latitude, longitude, altitude, physicochemical properties of soils, previous cropping) and information on crop management (type of tillage, seeding rates, types of sowing, number of plants per square meter, depth of sowing, fertilizers application and genetic coefficients of cultivars determined on the basis of their physiological parameters and grain yields). Daily climate data for 33 years (1980 to 2012) were collected from ASECNA (Agence pour la Sécurité de la Navigation Aérienne en Afrique et à Madagascar) synoptic stations of Natitingou and Kandi close to the research area.
The calibration of the variety used (EVDT-97 STRW) was based on the database of soils, climate, crop characteristics and crop management practice in the study areas. The genetic coefficients were determined through the GLUE program, a utility for estimating the genetic coefficients incorporated in DSSAT [20]. To perform plant growth and development simulations, the maize CSM-CERES used six eco-physiological coefficients.
Biophysical and economic analyzes using the seasonal analysis tool of the DSSAT model were also performed in order to determine a series of efficient options. Graphic analyzes were finally carried out to evaluate the dispersion of the various formulations in order to select only those which give the best yields with a low variance (a minimum of risk).
The seasonal analysis has two components. The first is the biophysical analysis that determines the minimum and maximum yields and their variance for the different treatments. The second category is the strategic and financial analysis that requires economic data. This analysis leads to the choice of the best and efficient fertilizer option. In more detail, the mean-Gini stochastic dominance as developed by Fosu et al. [21]. The financial analysis was done by integrating as input in the model production cost and maize price collected in the study area. Maize price use was that of the market during the harvest period.

Statistical evaluation of the DSSAT model
The evaluation of the performance of the DSSAT CERES-Maize model in the prediction of plant growth and development consists in validating the values simulated by the model for the 2012 season based on data observed during on-farm experiments. To do so, a number of tools were used such as: correlation coefficient [22], actual deviations separating simulated values from values observed, mean prediction errors RMSE [23] and the mean standard prediction error NRMSE [24,25].

On-farm experimentation
The experimental design used for the trials is a four-replicate complete random block with 8 m x 5.6 m elementary plots. This design  Table 2). The vulgarized rate (T4) represents 200 kg/ha of NPK and 50 kg/ha of urea. The simulated rates (T1) represent the optimal levels of N, P and K simulated by DSSAT. The adaptability rates were determined for the validation of the optimal rates of N, P and K and their comparison with the vulgarized rate.  The method of soil preparation was flat plowing. The sowings were made on the elementary plots with spacing of 80 cm between two rows and 40 cm on the rows (a seeding density of 62500 plants/ha with two plants per pocket). Two weedings were carried out, the first between the 11th and 14th day after sowing (DAS) and the second between the 40th and 44th DAS.
Thus, the total amounts of TSP and KCl and half of the urea were applied two weeks after sowing at the first weeding, and the remainder of the urea one month later. The harvest was made at physiological maturity following the perfect drying of maize cobs on the useful area of each elementary plot after removal of the edges.
The GLM procedure of the Statistical Analysis System version 9.2 software (SAS v. 9.2) was used for statistical analyzes of data from onfarm trials. These consisted mainly of two-factor (soil type and fertilizer formula) analyses of variance by agro-ecological zone. The mean values were then compared with each other using the Student Newman Keuls test at the 5% threshold (the probability level used to refer to a significant effect).

Soils chemical properties
The results relating to the chemical parameters of the different soil types prior to the establishment of the trials are given in Tables 3 and 4.
The analysis of these tables shows that organic matter contents were lower in concretionned tropical ferruginous soils than in hydromorphic soils (Nanébou) and hydromorphic tropical ferruginous soils (Arbonga). This low level of organic matter as well as that of total nitrogen observed in hydromorphic soils and hydromorphic tropical ferruginous soils reflects repeated use of these soils, with little or no return of nutrients either by burial of harvest residues, or directly by mineral fertilization [16].   Table 3: Soils chemical properties before trial installation at Nanébou. *= P<0.05.
These results support those of Igué and Yallou et al. [19,26], which showed that the cultivation of lands decreases their organic matter contents. Igué et al. showed that the organic matter content of cultivated soils decreases according to the cropping systems [27]. In the unbalanced system (poor farmer), organic matter is a very severe limitation compared to other systems (medium and balanced) or the limitation is average. Igué also indicated that organic matter in the topsoil (0-20 cm) decreases from 0.05 to 0.08% per year depending on the type of soil [26].  Table 4: Soil chemical properties before trial installation at Arbonga.
According to Worou, the low organic matter content of hydromorphic soils on the study site can be explained by a rather dry soil climate [1]. Fikri et al. stated that organic matter has a major influence on the physical and soils chemical properties and therefore on crop yields [28]. Concretionned tropical ferruginous soils have higher levels of phosphorus than hydromorphic soils and hydromorphic tropical ferruginous soils. It was found that in the Department des Collines, the phosphorus content can increase by 10% after 10 to 25 years of continuous cultivation of maize/cotton [26]. This may be due to the regular application of phosphate fertilizers. On the other hand, the Ca/Mg and Mg/K ratios in the three soil types showed good cationic balance without any significant difference (P>0.05) in the different soil types. It remains slightly higher in concretionned tropical ferruginous soils compared to hydromorphic soils and hydromorphic tropical ferruginous soils which are slightly more acidic.

Effect of different fertilizer formulas on maize grain yield according to soil type and area
The results of the analysis of variance relating to the effect of various fertilizer formulas on maize grain yields showed that fertilizer formulas have a highly significant influence (P<0.01 to P<0.001) on maize grain yield regardless of the areas and types of soil. Figures 2 and 3 show maize grain yields by area and soil type according to fertilizer formulas. The analysis of these Figures reveals that maize grain yields increase with increasing rates of nitrogen. Nitrogen is therefore the major limiting factor to maize yield in Northern Benin. The rate simulated by the DSSAT model (115-30-75) leads to significantly higher grain yield, regardless of soil types and agro-ecological zones (Figures 2 and 3).
These results are in line with those of Saidou et al. and Balogoun et al. who showed that nitrogen is the main limiting factor to cereal crop yields [2,8]. These observations show the crucial role played by nitrogen fertilizers in improving cereal yields [8]. The high mineralization rate of the organic matter is mainly the source of lack of nitrogen in these soils [10]. Singh et al. and Brassard found that nitrogen is the most limiting nutrient for cereal production in the Sub-Saharan Africa's soils [22,29]. As mentioned also by previous studies, most of the Africa's soils have low P level due to the nature and the type of the clays that their content [30,31]. This shows the importance of the supply of N and P to improve maize production in this part of Africa knowing the complementarity of these nutrients for plant. Moreover, the application of mineral fertilizers without any organic restitution further affects soil chemical characteristics with the number of years of cultivation [32]. Ultimately, achieving good yields depends not only on the nature of the soils but also on the amount of nitrogen available for plant nutrition. Igué et al. showed that with the fertilizer formula N 42 P 30 K 35 combined with manure the highest yields were 2,940.25 ± 383.60 and 2,923.60 ± 653.26 kg/ha respectively for concretionned soils and hydromorphic soils [13]. Without external nutrient supply, the productive capacity of the plots shows drastic deficiencies in the major nutrients (NPK); yield levels are 1,246.88 ± 359.39 and 1,327.60 ± 165.05 kg/ha respectively for concretionned ferruginous soils and hydromorphic soils with the absolute control (N 0 P 0 K 0 ).
In Arbonga (Banikoara), the difference in average yields varied significantly with the different types of soils under maize cropping. Grain yields on the concretionned soils increased by 600 kg/ha compared to the hydromorphic tropical ferruginous soils. Atakora et al. showed in a study in Ghana, that the differences in maize grain yields were more related to differences in soil fertility level [33]. Igué et al. also showed that the treatment N 88 P 30 K 35 plus manure gave the highest yields on tropical ferruginous soils [13]. On the other hand, on hydromorphic ferruginous soils, the treatment N 74 P 20 K 23 plus manure gave the highest yields. These observations support the works of Balogoun et al. which showed that to achieve high maize yield in the South and Center Benin, a rate of 80.5 kg N/ha would be required [8]. Indeed, achieving good yields depends not only on the nature of the soil but also on the amount of nitrogen available for plant nutrition.
In the agroecological zone II, the soils have a good productive potential for the cultivation of the maize variety EVDT ETR 97 whose yield, without external inputs, is around 1.5 t/ha. However, soil fertility decline is a major cause of low productivity in tropical soils [2,3,5]. To redress this situation, the use of organic fertilizers was promoted all the more because mineral fertilization without any organic fertilizer negatively affects the chemical characteristics of the soils; which shows the limits of mineral fertilization. According to Viennot, acidic soils have a negative impact on maize yields and are considered to be moderately suitable for this crop [34]. On the other hand, the hydromorphic soils of the study area were subjected to heavy pressure characterized by overutilization associated with inappropriate agricultural practices. The low productivity of these types of soils without the use of fertilizers is also linked to their topographical position in the landscape, which causes the stagnation of water on the surface of the plots and, in turn, contributes to the asphyxiation of plant root system [16,35]. Evaluation of the performance of the DSSAT model Table 5 shows the comparison between observed grain yields and those simulated by DSSAT taking into account fertilizer formulas according to the area and the type of soil. Grain yields simulated by the DSSAT model are slightly underestimated in zone II whereas they are slightly over estimated in zone IV. The values of the ratio observed-to-simulated values are very close to 1. Simulated values are therefore very close to observed values. The mean standard prediction error NRMSE between observed and simulated yields is between 11% and 20% for hydromorphic ferruginous soils whereas it is between 21% and 30% for the other soil types.  Tetteh and Nurudeen reported that the mean standard prediction error (NRMSE) between simulated and observed grain yields over two years (2010 and 2011) in Ghana was 26.13% and 18.24%, respectively [36]. This supports our results. The general remark is that, the model was very sensitive to fertilizer rates as mentioned also by Tetteh and Nurudeen and Atakora et al. [33,36]. NRSME values between observed and simulated results of 21%-30% are acceptable according to Jamieson et al. and Loague and Green [24,25]. This proves that, with correct inputs of soil and varietal characteristics a decision support tool like DSSAT could perfectly be used to extrapolate fertilizer recommendation data within a large agroecological zone presenting similar climatic characteristics and soil types. The results are also consistent with study carried out by Ritchie and Alagarswamy and Soler et al. who found that the CERES-Maize was able to accurately predict the phenology and maize grain yield for a wide range of environmental conditions [37,38]. Figure 4 presents the results of the biophysical analysis of grain yields by fertilizer formula in zone II according to soil types for the period 1980 to 2012. From this figure it appears that, in general, simulated grain yields are based on fertilizer rates. Thus, the formula 115-30-75 gave the highest average grain yields during the 33 years on the two soil types. Nevertheless, the formula 88-30-35 shows acceptable grain yields with less risk during the 33 years on the two types of soil. Figure 5 presents the results of the biophysical analysis of grain yields by fertilizer formula in zone IV according to soil types for the period 1980 to 2012. It appears that, in general, simulated grain yields are based on nitrogen rates. Thus, the formula 115-30-75 gave the highest average grain yields during the 33 years on the two soil types. Nevertheless, the formula 88-30-35 shows acceptable grain yields and with less risk during the 33 years on the two types of soil.

Application of the model to the multi-year assessment of fertilizer formulas
The financial analysis of the monetary incomes from maize per hectare with the efficiency of the various fertilizer formulas during the period 1980 to 2012 by zone and soil type is presented in Table 6. The results show that the formula 115-30-75 yielded the best monetary income per hectare and the best efficiency whatever the type of soil and the agro-ecological zone.  This means that if a producer uses the 115-30-75 formula, he only earns between 20,000 and 35,000 FCFA more than the one using the 80-30-35 formula. This is not so much, given the additional expenses for the nutrients N and K. It can be concluded that the formula 80-30-35 is the best regardless of the soil types and agro-ecological zones.  Tetteh and Nurudeen showed that the formula 160-90-90 produced the highest monetary income in the Guinean savanna zone in Ghana followed by formulas 120-0-90 and 120-45-90 respectively [36]. They pointed out that this was due to the high monetary income per hectare and the Gini coefficient. However, they indicated that due to high prices of fertilizers, their availability on the market and low natural soil fertility, the 120-45-90 formula is the most economical for sustainable maize production on Lixisols in the agro-ecological zone of the Sudan Savanna zone of Ghana. Furthermore, one can also consider that the model has been rational in the economy of N utilization by suggesting a reduce quantity. This observation is in accordance with the findings of Fosu et al. who stated that a supply of high rate of N leads to N leaching and possible contamination of water and luxury consumption by the plant while reducing the net return. The same arguments justify the choice of the formula 88-30-35 against the 115-30-75 for the production of maize in Northern Benin [21].

Conclusion
Generally, the DSSAT model was used to simulate maize yields in the agro-ecological zones II and IV of Benin. The grain yields simulated by the DSSAT model are slightly underestimated in zone II whereas they are slightly over estimated in zone IV. The values of the simulated and observed values ratio are very close to 1. Simulated values are therefore very close to observed values. The mean standard prediction error NRMSE between observed yields and simulated yields is between 11% and 20% for hydromorphic ferruginous soils whereas it is between 21% and 30% for the other soil types. Formulas 115-30-75 and 88-30-35 gave the best yields on-farm. Moreover, the seasonal analyzes with the DSSAT model over 33 years showed that the same formulas 115-30-75 and 88-30-35 gave the best yields and the best monetary incomes. The study recommended the formula 88-30-35 kg/ha of NPK as the most economically and strategically efficient fertilizer formula that gave optimum yields with less risk during the 33 years in the two agro-ecological zones of Northern Benin.