INTRODUCTION
Several studies have shown that soils in West Africa are structurally degraded and that the nutrients are depleted (Pichot et al. Reference Pichot, Sedogo, Poulain and Arrivets1981; Bado et al. Reference Bado, Sedogo, Cescas, Lompo and Bationo1997; de Ridder et al. Reference de Ridder, Breman, Van Keulen and Stomph2004; Schlecht et al. Reference Schlecht, Buerkert, Tielkes, Bationo, Bationo, Waswa, Kihara and Kimetu2007). Furthermore, agricultural systems and farming techniques contribute to the degradation of environmental resources (Bandre & Batta Reference Bandre and Batta1998). In recent years, interest in the sequestration and turnover of organic matter (OM) in soils has increased because soil organic matter (SOM) affects the stabilization of soil structure and the cycling of organically held plant nutrients (Oades & Waters Reference Oades and Waters1991; Tiessen et al. Reference Tiessen, Cuevas and Chacon1994; Lal et al. Reference Lal, Follett, Stewart and Kimble2007). In addition, carbon (C) held in or released from SOM could change atmospheric carbon dioxide (CO2) levels (Schlesinger Reference Schlesinger1991; Reeves Reference Reeves1997; Tiessen et al. Reference Tiessen, Feller, Sampaio and Garin1998). The exchange of C between the land surface and the atmosphere is responsible for a large fraction of the global variation in atmospheric CO2, both within and between years (Tucker et al. Reference Tucker, Fung, Keeling and Gammon1986). Changes in vegetation productivity and C storage due to climate change and change in land use could have an important impact on the global C cycle (Mooney & Koch Reference Mooney and Koch1994).
Sorghum (Sorghum bicolor) is the main cereal crop in West Africa. Cotton (Gossypium spp.) is considered to be more profitable economically and plays an important role in Burkina Faso's economy, but can only be grown in areas where there is sufficient rainfall. Peanut (Arachis hypogoea) and cowpea (Vigna unguiculata) are the most commonly grown legumes (Bado Reference Bado2002). In many cases, food crops are grown as dual purpose crops in which the seeds are used for human consumption, whereas the stover is used for feeding animals. Bado (Reference Bado2002) reported that rotations with peanut or cowpea increased the amount of available nitrogen (N) for the next crop, resulting in a considerable increase in sorghum yield. Short-term use of chemical fertilizer, followed by maintenance fertilization using chemical and animal manure, are required to improve soil fertility and increase productivity, but in some cases socio-economic constraints have largely prevented the improvement of most semi-arid agricultural systems to date (Tiessen et al. Reference Tiessen, Feller, Sampaio and Garin1998; Bationo et al. Reference Bationo, Kihara, Vanlauwe, Waswa and Kimetu2007; Doraiswamy et al. Reference Doraiswamy, McCarty, Hunt, Yost, Doumbia and Franzluebbers2007).
Bado et al. (Reference Bado, Bationo, Lompo, Segda, Cescas, Sedogo, Zapata and Nguyen2006b) emphasized the importance of identifying a cost-effective way of improving soil fertility and productivity in Burkina Faso and other countries in Sub-Saharan Africa. In a crop rotation study conducted in Burkina Faso, Bado et al. (Reference Bado, Bationo, Lompo, Segda, Cescas, Sedogo, Zapata and Nguyen2006b) found that high yields were achieved by simultaneous application of organic and mineral fertilizers or when mineral fertilizers were associated with dolomite. Nitrogen fixing legumes, such as peanut and cowpea, are important components of the cropping system of resource-poor farmers and serve as food as well as cash crops (Bado et al. Reference Bado, Bationo and Cescas2006a). Crop rotations of continuous sorghum and cotton–peanut–sorghum resulted in a decrease in soil organic carbon (SOC). Low quantities of manure, e.g. 5 tonnes(t)/ha every 2 years, did not significantly increase SOC after 6 years of cultivation. Ouédraogo et al. (Reference Ouédraogo, Mando and Stroosnijder2006) reported that combining recalcitrant organic amendments and N fertilizer is the best option in sustaining crop production in semi-arid West Africa due to a positive change in soil particulate OM and a reduction in the decline of soil C, especially for the more stable fraction that is responsible for the soil structure maintenance. Hien et al. (Reference Hien, Gandry and Oliver2006) recommended selection of crops with a good rooting system, the application of organic manures, which are resistant to rapid biodegradation, reduction of tillage to avoid creating conditions for rapid biodegradation and regular liming in order to reduce the loss of soil C. A long-term rotation study was conducted from 1960 to 1978 in a degraded ferruginous tropical soil of Burkina Faso: soil C was measured in 1969 and 1978 and it was found that the soil C in the 0–200 mm layer decreased from 2·9 to 2·5 g/kg for the control treatment without fertilizers, from 2·9 to 2·4 g/kg for the treatment with applications of inorganic fertilizers and from 3·1 to 2·5 g/kg for the treatment with inorganic fertilizers and crop residues (sorghum) incorporated into the soil. The soil C increased from 3·1 to 3·5 g/kg for the treatment with inorganic fertilizers and 5 t/ha/yr of manure and increased from 5·3 to 6·6 g/kg for the treatment with inorganic fertilizers and 40 t/ha/yr of manure (Pichot et al. Reference Pichot, Sedogo, Poulain and Arrivets1981). Lal (Reference Lal, Lal, Kimble and Stewart2000) observed annual rates of soil C increase under no-till management ranging from 363 kg/ha/yr to more than 1000 kg/ha/yr (for one severely depleted soil) over a 3-year experiment to restore soil C in western Nigeria.
In other regions, such as India, it was found that the total system productivity was greater in groundnut–groundnut and groundnut–chickpea systems in comparison with other systems, but the SOC declined (Ghosh et al. Reference Ghosh, Manna, Dayal and Wanjari2006). In a study conducted in Spain, Herencia et al. (Reference Herencia, Ruiz, Melero, Garcia Galavis and Maqueda2008) reported that in general, treatments with manure resulted in increased soil nutrient values, which had an effect on several later crop cycles, and that the use of organic composts resulted in an increase in SOC and the storage of nutrients, which can provide long-term fertility benefits.
However, soils, climate and management practices vary over space and time, creating an almost infinite combination of factors that interact and influence how much C is stored in soils. Quantifying soil C sequestration under widely varying conditions is complicated. Thus, simulation models can be used to complement information gained from experiments to help understand and predict SOC and yield responses to soil, climate and crop management. Typically, long-term effects can be assessed by the models to predict environmental risks and to evaluate alternative management practices for alleviating these risks (Muchow et al. Reference Muchow, Hammer, Carberry, Muchow and Bellamy1991; Thornton et al. Reference Thornton, Hoogenboom, Wilkens and Bowen1995; Chipanshi et al. Reference Chipanshi, Ripley and Lawford1997; Nijbroek et al. Reference Nijbroek, Hoogenboom and Jones2003). The Decision Support System for Agrotechnology Transfer (DSSAT) is a comprehensive decision support system (DSS) for assessing management options (Tsuji et al. Reference Tsuji, Jones, Hoogenboom, Hunt, Thornton, Tsuji, Uehara and Balas1994; Jones et al. Reference Jones, Hoogenboom, Porter, Boote, Batchelor, Hunt, Wilkens, Singh, Gijsman and Ritchie2003; Hoogenboom et al. Reference Hoogenboom, Jones, Wilkens, Porter, Batchelor, Hunt, Boote, Singh, Uryasev, Bowen, Gijsman, Du Toit, White and Tsuji2004). The cropping system simulation model (CSM) included in DSSAT v4.0 (Jones et al. Reference Jones, Hoogenboom, Porter, Boote, Batchelor, Hunt, Wilkens, Singh, Gijsman and Ritchie2003; Hoogenboom et al. Reference Hoogenboom, Jones, Wilkens, Porter, Batchelor, Hunt, Boote, Singh, Uryasev, Bowen, Gijsman, Du Toit, White and Tsuji2004), is process-oriented, dynamic and simulates growth, development and yield for more than 25 different crops. The Century model (Parton et al. Reference Parton, Scurlock, Ojima, Gilmanov, Scholes, Schimel, Kirchner, Menaut, Seastedt, Garcia Moya, Kamnalrut and Kinyamario1993) simulates the long-term dynamics of C for different plant–soil systems. It was developed to analyse a wide range of crop rotations and tillage practices and to determine the impact of crop management on the productivity and sustainability of agroecosystems. The Century model has been used successfully to simulate SOC across a variety of land uses and climate types (Kelly et al. Reference Kelly, Parton, Crocker, Grace, Klir, Korschens, Poulton and Richter1997; Diagana et al. Reference Diagana, Antle, Stoorvogel and Gray2007; Ogle et al. Reference Ogle, Breidt, Easter, Williams and Paustian2007). Gijsman et al. (Reference Gijsman, Hoogenboom, Parton and Kerridge2002) incorporated a Century-based module into CSM to allow for more flexibility in handling different agricultural systems for long-term simulations of crop rotations. This capability is particularly important to enable CSM to be used for predicting yields in low-input cropping systems where soils tend to be deficient in OM and nutrients. However, even under high-input management systems where large amounts of soluble nutrients are supplied, hence, minimizing the role of OM in supplying nutrients, its role in enhancing soil productivity via the modification of soil processes cannot be overlooked (Porter et al. Reference Porter, Jones, Adiku, Gijsman, Gargiulo and Naab2010).
There are a number of uncertainties associated with simulation models and their use, and if one does not adequately address these uncertainties, the simulated results will be meaningless (Jones et al. Reference Jones, Graham, Wallach, Bostick and Koo2004). These uncertainties are due to the fact that models are simplifications of reality; there are uncertainties in model parameters and in the data inputs. Thus, research is needed to ensure that models can reproduce the responses that are observed in field experiments. However, soil C measurements are also uncertain and costly, and the errors could be much larger than the annual changes in SOC. By combining measurements with model predictions, a more accurate estimate of SOC can be obtained (Koo et al. Reference Koo, Bostick, Jones, Gijsman and Naab2003). The goal of the present study was to determine the yield variability and SOC dynamics for different rotations and fertilizer levels of cropping systems in a semi-arid region of West Africa. Specific objectives included the simulation of eight crop rotations and three fertilizer treatments, and the analysis of observed and simulated yield trends and SOC dynamics for the different crop rotations and management treatments.
MATERIALS AND METHODS
Crop management
A crop rotation study was initiated in Farakô-ba (4°20′W, 11°6′N), Burkina Faso in 1993 (Bado Reference Bado2002). The rotations consisted of various cropping sequences that included sorghum, cotton, peanut and maize, which are important food and cash crops in the region, and native grassland as a fallow. The experiment was initiated on a 6-year native grassland fallow and was an incomplete-factorial, split-plot design arranged in randomized blocks with eight crop rotations as whole plot factors, nine input levels as split-plot factors and four replications. The input levels were distinguished by the type of soil amendment that was added (Bado Reference Bado2002; Bostick et al. Reference Bostick, Bado, Bationo, Soler, Hoogenboom and Jones2007).
The eight rotations plus three management treatments were analysed from 1993 to 2004 (Table 1). The crop rotations were: (1) continuous fallow, (2) continuous sorghum, (3) continuous peanut, (4) continuous cotton, (5) peanut–cotton–sorghum, (6) peanut–sorghum–cotton, (7) maize–sorghum–cotton and (8) sorghum–fallow. The three management treatments were: (1) no supplemental N fertilizer (control treatment), (2) the application of inorganic N and (3) the application of manure. The recommended application rates of inorganic N, P and K prior to planting were 37, 10 and 11 kg/ha, respectively, for sorghum; 44, 15 and 17 kg/ha, respectively, for cotton; 14, 10 and 11 kg/ha, respectively, for peanut and 60, 12 and 7 kg/ha, respectively, for maize. For the manure treatments, air-dried cow manure (184 g C/kg, 18 g N/kg, 3·1 g P/kg and 1·6 g K/kg) at a rate of 3000 kg/ha was applied each year at planting.
Crops planted by year for the different crop rotations
Climate
The climate of Burkina Faso is similar to other countries in the region and is characterized by a dry season, which lasts, on average, from mid-November to mid-April, and a wet season, which lasts from May to mid-September. Most of the country (centre and south) lies in the Sudanese climatic zone, while the north is influenced by the climate of the Sahel. Rainfall varies from 400 mm in the north to 1000 mm in the southwest. However, the most important aspects of the rainfall are its temporal and spatial variability (Bandre & Batta Reference Bandre and Batta1998).
According to the FAO eco-climatic classification, Farako-bâ is located in the Sudanian climatic zone, defined by one rainy season per year (FAO 2005). The average seasonal cumulative precipitation and s.d. from May to October from 1993 to 2004 at the study site were 927 and 171 mm, respectively. The total precipitation for the cropping season (May–October) was variable for each individual year (Fig. 1). The maximum total seasonal rainfall was 1174 mm for 1998 and the minimum was 639 mm for 2002. The cropping season is also characterized by high temperatures. For example, the average maximum and minimum temperatures for June were 32 and 22·3°C, respectively, while for July they were 30 and 21·5°C, respectively. However, the highest temperatures usually occur during the months of March and April, with average maximum temperature of 38°C.
Total precipitation for the cropping season (May–October) in Farakô-ba, Burkina Faso.
Soil characteristics and measurements
The soil at the experimental site was a weathered ferruginous tropical soil (Bado Reference Bado2002), which is an Alfisol or Ultisol in the USDA soil taxonomy classification (Soil Survey Staff 1999). This soil covers 0·40 of the land area in the sub-humid zone of West Africa (6–12°N; Adeoye & Mohamed-Saleem Reference Adeoye and Mohamed-Saleem1990). These ferruginous tropical soils are characterized by a sandy surface horizon that is low in OM and base exchange capacity (ILCA 1979). The soil texture at the experimental site was loamy sandy, with sand, silt and clay contents of 74, 19 and 7% by volume, respectively (Bado Reference Bado2002; Bostick et al. Reference Bostick, Bado, Bationo, Soler, Hoogenboom and Jones2007). The soil profile inputs were complemented with data from a ferruginous tropical soil from northern Ghana that was obtained from DSSAT V4.0 database (Gijsman et al. Reference Gijsman, Thornton and Hoogenboom2007). The next soil horizon used as input for the model was similar to the upper horizon in terms of texture, and was characterized by sand, silt and clay contents of 66, 21 and 13% by volume, respectively (0·2–0·40 m depth). The next soil horizon, 0·4–1·2 m depth, contained 55, 22 and 23% by volume of sand, silt and clay, respectively. The bulk density in the top 0·2 m of the soil was estimated as 1·46 g/cm3 from soil texture components (Rawls Reference Rawls1983). This top layer of the soil had low values for total OM (11 g/kg), a weak cation exchange capacity (CEC), was slightly acidic (pH water=6·5) and was N and P poor (Bado Reference Bado2002).
Three composite samples were collected from the top 0·2 m of soil before planting in 1993 to characterize the initial SOC of the native grassland. The average SOC was 5·5 g/kg. SOC was also measured after harvest in 1998 and 2003 for some crop rotations only (continuous fallow, continuous sorghum, sorghum–fallow, continuous cotton and maize–sorghum–cotton) and input levels (control, N-fertilizer and manure applications). The Walkley–Black method (Walkley & Black Reference Walkley and Black1934) was used for all SOC determinations. For the 1998 and 2003 SOC measurements, four composite soil samples were taken and analysed for each treatment that was measured.
Simulation of crop rotations
The crop rotation or sequence tool of CSM-DSSAT V4.0 (Thornton et al. Reference Thornton, Hoogenboom, Wilkens and Bowen1995) was used to simulate the eight crop rotations and three fertilizer treatments. The CSM-DSSAT sequence analysis tool allows the user to conduct simulations of crop rotations or crop sequences and to analyse the results (Bowen et al. Reference Bowen, Thornton, Hoogenboom, Tsuji, Hoogenboom and Thornton1998). The main aspect of the sequence analysis is the consideration of experiments that are conducted across multiple cropping seasons. Therefore, the carry-over of the soil water and nutrient status is affected from one cropping season or crop to the subsequent one (Thornton et al. Reference Thornton, Wilkens, Hoogenboom, Jones, Tsuji, Uehara and Balas1994). The Century soil OM module is integrated into the CSM-DSSAT V4.0 in order to better model the dynamics of soil organic nutrient processes (Gijsman et al. Reference Gijsman, Hoogenboom, Parton and Kerridge2002; Porter et al. Reference Porter, Jones, Adiku, Gijsman, Gargiulo and Naab2010). The current version of CSM does not consider the P and K balance in the soil (Jones et al. Reference Jones, Hoogenboom, Porter, Boote, Batchelor, Hunt, Wilkens, Singh, Gijsman and Ritchie2003); thus, these elements were not defined as inputs.
Daily rainfall data were obtained from a conventional weather station located at the research station of the Environmental and Agricultural Research Institute (INERA), located adjacent to the experimental site, in Farakô-ba. The daily maximum and minimum temperature and solar radiation were obtained from the nearest weather station (10 km), located in Bobo Dioulasso. Using weather data from these two stations was the best way to obtain accurate daily temperatures and also rainfall at the experimental site.
The soil parameters needed to characterize the profile were obtained from Bado (Reference Bado2002). The parameters related to the different SOC decomposition rates that were required to run the Century component of the CSM were set according to previous reports on similar soils from Burkina Faso (Hien et al. Reference Hien, Gandry and Oliver2006; Barthes et al. Reference Barthes, Brunet, Hien, Enjalric, Conche, Freschet, D'Annunzio and Toucet-Louri2008). Therefore, the SOC corresponding with the fraction of stable organic C (SOM3, also called passive) was set equal to 0·26. Hien et al. (Reference Hien, Gandry and Oliver2006) reported 2·6 g/kg of SOC for a soil in the Savannah region of the Southwest of Burkina Faso after 13 years of continuous cropping, with removal of the crop residues. The intermediate SOC (SOM2) was set equal to 0·72. The model was set up to initiate the simulations 6 years prior to the beginning of the experiment. During this period, the growth of a native grassland fallow was simulated, mimicking the same crop that was grown prior to the start of the rotation experiment. Thereafter, the eight rotations and the three management treatments were simulated from 1993 to 2004 (Table 1).
For model evaluation, the simulated dates of flowering and maturity as well as yield were compared with the observed values. The phenology was evaluated using growth and development data for 2004 (data not shown). The cultivar coefficients are specific for each crop, and were obtained sequentially, starting with the phenological development parameters related to flowering and maturity dates, followed by the crop growth parameters related with kernel filling rate and kernels number per plant (Hunt & Boote Reference Hunt, Boote, Tsuji, Hoogenboom and Thornton1998). An iterative procedure (Hunt et al. Reference Hunt, Pararajasingham, Jones, Hoogenboom, Imamura and Ogoshi1993) was used to select the most appropriate value for each coefficient. A detailed description of the cultivar coefficients used by the CSM-DSSAT V4.0 can be found in Hoogenboom et al. (Reference Hoogenboom, Jones, Wilkens, Batchelor, Bowen, Hunt, Pickering, Singh, Godwin, Baer, Boote, Ritchie, White, Tsuji, Uehara and Balas1994) and Hunt & Boote (Reference Hunt, Boote, Tsuji, Hoogenboom and Thornton1998). The combination of cultivar coefficients that resulted in an average simulated yield that was the most similar to the average observed yield was selected. The selected varieties used in the present study were: Sariaso-02 for sorghum, RMP12 for peanut, FK37 for cotton and SR22 for maize.
Statistical analysis
An analysis of yield variability, total biomass incorporated into the soil and SOC was conducted for the different crop rotations and fertilizer treatment combinations. The observed and simulated yield variations were examined by plotting the normalized yield for each crop in each rotation and treatment according to Meinke & Hammer (Reference Meinke and Hammer1995).
where Y n is the normalized yield, Yi is the yield for each individual growing season, 
is the average yield for each crop across all years of the rotation and s.d. is the standard deviation of the yield for each crop across all years of the rotation.
For the simulated yield, an analysis of the trends was conducted for the different rotations. The results were examined by plotting the progression with time of the cumulative deviations from the mean for the different crops of the different rotations and treatments. The plots start and end with zero, and the year to year fluctuations correspond to periods above or below the average yield for different crops. In the present study, the plots show long-term yield trends for different crops without removing variability (Russell Reference Russell1981; Hammer et al. Reference Hammer, Woodruff and Robinson1987; Meinke & Hammer Reference Meinke and Hammer1995; Garcia et al. Reference Garcia, Hoogenboom, Guerra, Paz and Fraisse2006). For the simulated SOC dynamic, an analysis of the different rotations and treatments was conducted, which was presented in a graphical format depicting also the observed average and s.d. of the SOC values for the top soil layer (0–0·2 m) for 1993, 1998 and 2003.
RESULTS
Yield
The crop rotation or sequence tool of DSSAT V4.0 was able to simulate the observed yields for different crop rotations (Table 2). The average observed yields varied between 720 kg/ha for the continuous sorghum rotation and 2005 kg/ha for maize crop in the cotton–maize–sorghum rotation. The simulated average yields varied between 801 kg/ha for the continuous cotton rotation and 2672 kg/ha for the maize crop in the cotton–maize–sorghum rotation. The Pearson coefficient between observed and simulated yield values was 0·84 and the root mean square error (RMSE, expressed as a proportion) was 0·26 indicating a fairly good simulation. The variation in the simulated yield was similar to the variation of the observed yield for all crop rotations and the three treatments (Figs 2–4). There were a few anomalies in which the simulated variation in yield differed from the observed, e.g. in 1993 the models in general under-predicted the yield. The differences between simulated and observed values were partially due to the fact that the models do not consider all biotic and abiotic factors that impact crop growth and development. Differences between observed and simulated values were probably due to pest and diseases, but there were no detailed experimental records to support this on particular years and crops. In general, the observed yields were higher than the average only at the beginning of the cropping system rotations. This can be explained by the better soil properties for the first year of the rotations, which were favoured by the 6 previous years with native grassland. However, for the remainder of the years for the different rotations, the variation in the simulated yield was, in general, similar to the observed yield. Lack of agreement between observed and simulated accumulated deviations occurred in the control treatments for the continuous sorghum and sorghum–fallow rotations in 2001, in which the observed yields were below the average, but the simulated yield values were above average (Fig. 2). For the peanut–cotton–sorghum rotation without N fertilizer (control), the simulated peanut yield for 1996 and 2002 were below the average simulated, while the observed yields were above the average observed. In 1999, the opposite occurred; the simulated yield was above the average and the observed yield was below average. These examples indicate that the model was unable to simulate yield trends accurately for a few specific years for the control rotation (without N). Similar to these specific circumstances, there were also a few predictions in simulated yield trends that were different from observed for specific cases for crop rotations with N fertilizer and with manure applications. This included the peanut–cotton–sorghum rotation with N fertilizer and also the rotation with manure applications for 1999 where the observed peanut yield was below the average, while the simulated yield was above the average (Figs 3 and 4).
Simulated and observed yield expressed as s.d. from the mean for the different rotations without N fertilizer (control).
Simulated and observed yield expressed as s.d. from the mean for the different rotations with N fertilizer.
Simulated and observed yield expressed as s.d. from the mean for the different rotations with manure applications.
Observed and simulated average yields for different crop rotations
The progression across years of the cumulative deviations from the mean for the different crops of the various rotations and treatments were plotted for simulated and observed yield (Figs 5–7). The plots start and end with zero, and the year to year fluctuations correspond to periods above or below the average yield. Thus, the yield trends expressed as the cumulative deviation from the mean yield for each crop of a rotation and treatment had, for most of the years, crops and treatments, a good agreement between the simulated and observed values (Figs 5–7). There was a decrease in observed yield for the final years of the crop rotations (2002–4), with a below-average yield for the different crops. In general, the decrease in yield for the final years of the rotation was also simulated. The yield decrease can be explained by the reduction in soil fertility after years of continuous cropping systems.
Simulated and observed yield expressed as cumulative deviation from the mean (cumulative normalized yield) for the different rotations without N fertilizer (control).
Simulated and observed yield expressed as cumulative deviation from the mean (cumulative normalized yield) for the different rotations with N fertilizer.
Simulated and observed yield expressed as cumulative deviation from the mean (cumulative normalized yield) for the different rotations with manure applications.
Biomass
For the crop rotations involving N fertilizer and manure applications, the simulated average biomass and s.d. bars varied among rotations and crops. However, this variation was similar between N fertilized treatments and manure applied treatments, as represented in the average and s.d. bars of Fig. 8. When comparing average biomass values it was found that peanut was the crop with the largest simulated biomass values. The growing conditions were better for peanut because it is a legume and can fix N from the rhizobia that are formed in the root nodules. Sorghum was affected by N stress even when N was applied, because the amount was insufficient to supply the N needed by the crop. In the maize–sorghum–cotton rotation, the sorghum biomass with N fertilizers was lower than sorghum in the same sequence with manure. The manure applied increased the N available in the soil and this explains the high average values of biomass in comparison to the N fertilizer treatment. The aboveground biomass for the crop rotations without fertilizer (control) was removed from the experimental plot at harvest and not incorporated into the soil, causing a reduction in SOC across years.
Simulated average aboveground biomass incorporated into the soil for the different rotations with N fertilizer and manure applications.
SOC
Both the observed and simulated SOC content for all crop rotations decreased across the years of the experiment, except for the continuous fallow (native grassland) in which the SOC remained stable around 5·0–5·5 g/kg (Fig. 9) and the sorghum–fallow rotation with manure application. The crop rotation or sequence tool of CSM-DSSAT V4.0 (and the integrated Century SOM module) was able to simulate the SOC for the continuous fallow system accurately, with a value of 0·01 for the normalized RMSE (expressed as a proportion). In general, the simulated SOC values were within the range of the s.d. of the observed data for the other crop rotations (Fig. 9). The observed SOC rate of decline in the rotations sorghum–fallow, continuous sorghum and continuous cotton appeared to stabilize after 1999, but with an important variability in the SOC values. The model predicted a continued reduction in SOC rate of decline that eventually could impact the sequestration of SOC estimations. The simulated biomass that was incorporated into the soil for the N fertilized treatments and manure treatments resulted in a less pronounced decrease of simulated SOC across the years (Fig. 8).
Simulated and observed SOC for the different crop rotations and different managements.
The crop rotation that had the best performance with respect to simulated SOC was the sorghum–fallow rotation with manure application because this rotation had an average increase of 31 kg/ha/yr for the top 0·2 m of the soil profile (Fig. 10). The same crop rotation showed smaller decreases in SOC for the N fertilizer application and for the treatment that did not have any fertilizer application in comparison with the other crop rotations. The decrease in simulated SOC was 173 kg/ha/yr for the N fertilizer application treatments and 232 kg/ha/yr for treatments that did not have any fertilizer. Decreases in SOC of more than 400 kg/ha/yr were simulated for each of the other control crop rotations (treatments without fertilizer application) (Fig. 10). When N fertilizer was applied, the model simulated a decrease in SOC that ranged from 250 to 380 kg/ha/yr. When manure was applied, the model simulated a decrease in SOC decrease that ranged from 160 to 325 kg/ha/yr for the top 0·2 m of the soil, for continuous peanut and continuous cotton rotations, respectively.
Simulated average SOC change per year (0–0·2 m depth) for different rotations and management practices from 1993 to 2004.
In cropping systems with limited amounts of external inputs, the SOM pool may contribute significantly to plant nutrition (Vanlauwe et al. Reference Vanlauwe, Aman, Aihou, Tossah, Adebiyi, Sanginga, Lyasse, Diels and Merckx1998). The Century model includes three SOM pools, active (microbial SOM1), slow (intermediate SOM2) and passive (SOM3), which have different potential decomposition rates (Parton et al. Reference Parton, Scurlock, Ojima, Gilmanov, Scholes, Schimel, Kirchner, Menaut, Seastedt, Garcia Moya, Kamnalrut and Kinyamario1993; Batjes Reference Batjes2001). At the end of crop rotation in 2004, the simulated average values for the different components of the SOC was high for the SOM2 fraction (intermediate) for the treatments with manure applied, accounting for 0·72 of the total (Table 3). On the other hand, the control crop rotations (without N applied) had the lowest values of SOM2 in year 2004. When analysing the SOC dynamic for each rotation, it was found that the continuous fallow (native grassland) had a high simulated SOM2 in year 2004 of 0·75. For the crop rotations without N fertilizer (control), the SOC decreased to near 0·003 after 12 years.
Simulated average SOC for the three components SOM1 (microbial), SOM2 (intermediate) and SOM3 (passive) for all treatments at a depth of 0–0·2 m for 2004 after 11 years of crop rotation simulation
DISCUSSION
The removal of aboveground biomass at harvest for the crop rotations without fertilizer (control) caused an average reduction in SOC of 46% after 12 years, as there was no incorporation of biomass and residue into the soil. The removal of crop residues has been reported to be adverse in terms of soil quality due to the depletion of the soil C pool, which negatively affects crop production (Lal Reference Lal2007). In a study conducted in Ghana, Adiku et al. (Reference Adiku, Jones, Kumaga and Tonyigah2009) found that over a 4-year period, maize rotated with bare fallow produced an average maize biomass and yield of 4·0 and 1·0 t/ha/yr, respectively; SOC declined over time, with the maize bare fallow rotation losing about 55% of the initial SOC in the 4 years. The decline in SOC was 19% for a fertilized maize–grass rotation, and other treatments lost between 33 and 44% of the initial SOC. Furthermore, results from Zimbabwe reported by Mtambanengwe & Mapfumo (Reference Mtambanengwe and Mapfumo2005) highlight that the cumulative effects of applying substantial amounts of OM to specific areas on a regular basis resulted in consistently higher levels of SOC in designated rich fields than corresponding poor fields (or field sections). The key challenge is enabling farmers to generate sufficient quantities of biomass during each cropping cycle in order to maintain the productivity of rich fields, while rehabilitating the non-productive ones (Mtambanengwe & Mapfumo Reference Mtambanengwe and Mapfumo2005). Crop residues also have been reported to be a major factor in maintaining or reducing the loss of SOC on two tropical soils of Brazil (Leite et al. Reference Leite, Doraiswamy, Causarano, Gollany, Milak and Mendonca2009). Nziguheba et al. (Reference Nziguheba, Merckx and Palm2005) found in western Kenya that organic residues can contribute to the improvement of soil nutrient cycling. However, the magnitude of their contribution depends on the biochemical properties of the residues. Regular applications of manure on rice–wheat and maize rotation also have been associated with an increase in SOC and grain yield across sites in China (Zhang et al. Reference Zhang, Xu and Zhang2009).
Zingore et al. (Reference Zingore, Manyame, Nyamugafata and Giller2005) reported that for soils on smallholdings in Zimbabwe, the organic C declined rapidly under cultivation to attain a new equilibrium within 10 years after woodland clearance. Mando et al. (Reference Mando, Ouattara, Somado, Wopereis, Stroosnijder and Breman2005) found that 10 years of continuous sorghum cultivation without organic inputs in a sandy loam Lixisol soil at Saria, Burkina Faso caused significant losses of C. However, addition of manure was effective in maintaining similar levels of C compared to fallow plots. Without manure, SOM was mainly stored as fine organic matter (FOM) with a size-fraction <0·053 mm. In plots with manure and in fallow plots, the addition of manure doubled the particulate organic matter (POM) (size-fraction 0·053–2 mm) concentration. POM was greatly affected by long-term soil management options. Similar results were also found by Nziguheba et al. (Reference Nziguheba, Merckx and Palm2005) in a soil amended with organic residues in western Kenya and by Ouattara et al. (Reference Ouattara, Ouattara, Serpantié, Mando, Sédogo and Bationo2006) in western soils of Burkina Faso. Willson et al. (Reference Willson, Paul and Harwood2001) reported that the POM amount combined with information about recently incorporated crop residues, were found to be good predictors of N mineralization potential in systems using conventional, legume-based organic and manure-based organic fertility management. In a study conducted in West Africa, Vanlauwe et al. (Reference Vanlauwe, Aman, Aihou, Tossah, Adebiyi, Sanginga, Lyasse, Diels and Merckx1998) found that the average proportion of soil N belonging to the POM pool ranged from 0·09 to 0·29 depending on the annual N inputs from maize stover and prunings. Bationo & Buerkert (Reference Bationo and Buerkert2001) recommended targeted applications of mineral fertilizers and effective recycling of organic amendments as crop residues and manure to maintain food production for a rapidly growing population in the agro-pastoral farming systems in the Sudano–Sahelian zone of West Africa. Similarly, Vågen et al. (Reference Vågen, Lal and Singh2005) and Farage et al. (Reference Farage, Ardo, Olsson, Rienzi, Ball and Pretty2007) stated that for croplands the critical issue in increasing SOC is the addition of crop residues or manure and the reduction in soil disturbance through conservation tillage or no-till systems. De Costa & Sangakkara (Reference De Costa and Sangakkara2006) highlighted the benefits of agronomic measures such as mulching, the application of green manure, the use of crop residues, fallow rotation and agroforestry as very economical and efficient methods of increasing productivity and maintaining sustainability of small farming units in Asia. Zougmore (Reference Zougmore2003) recommended combining soil and water conservation techniques with soil fertility management practices to improve the efficiency of these measures, and thus sustain soil productivity in Burkina Faso and other countries in Sub-Saharan Africa.
In the present study, the model simulated the SOC for the continuous fallow system accurately, with a normalized RMSE (expressed as a proportion) of 0·001. In general, the simulated SOC values were within the range of the s.d. of the observed data for the other crop rotations. The simulated results for SOC from the present study are similar to previous studies that were conducted with the standalone version of the Century model. For instance, the Century model was able to accurately simulate the level of total soil C for five soil types in Brazil (Cerri et al. Reference Cerri, Cerri, Paustian, Bernoux and Mellilo2004): modelling also provided a flexible and powerful way to assess how different scenarios for pasture management and land use change can affect soil C dynamics (Cerri et al. Reference Cerri, Cerri, Paustian, Bernoux and Mellilo2004). The present study highlights that the CSM is a useful tool to simulate crop growth and the development for various crop rotations and management practices.
The understanding of the SOC dynamic for different crop rotations is important in relation to soil fertility aspects, but it also complements the research focus in recent years on the global carbon cycle, as it has become clear that increased levels of CO2 and other greenhouse gasses in the atmosphere are causing changes to our climate at an increasingly rapid rate (McCarthy et al. Reference McCarthy, Canziani, Leary, Dokken and White2001). Soil C sequestration has been recognized as an effective, low-cost technology to mitigate climate change. Simulation models, alone or in combination with soil sampling and other techniques, have been used to help with monitoring changes in soil carbon levels as affected by climate, soil and crop management options (Apezteguía et al. Reference Apezteguía, Izaurralde and Sereno2009). However, as stated by Batjes (Reference Batjes2004), mitigation of climate change by increased C sequestration in soils appears particularly useful when addressed in combination with other pressing regional challenges that affect the livelihood of the people, such as combating land degradation and ensuring food security, while at the same time curtailing global anthropogenic emissions.
In summary, the CSM simulated the yield trends of various crops grown in Burkina Faso fairly well, with some inconsistencies for a few years. The simulated SOC increased slightly across the years for the sorghum–fallow rotation with manure application. However, the most accurate simulation of SOC was obtained for the continuous fallow (native grass) with no supplemental fertilizer applied, in which the SOC remained stable around 5·0–5·5 g/kg. The crop rotations that included N-fertilizer or manure applications showed an increase in the total amount of biomass that was simulated, which was incorporated back into the soil in order to reduce the loss of SOC. The present study showed that DSSAT and its associated CSM can be used for the identification of cropping systems that can enhance the agricultural production and sustainability in Burkina Faso and other countries in West Africa that have similarly degraded soils. However, it is very important to define the initial values for SOC accurately as this will affect the outcome of the simulations. Further research is needed in order to extend the adoption of CSMs as a tool for decision makers in Burkina Faso and other countries in the region aiming to recommend suitable rotations to different soils or environmental conditions.
This work was supported in part by a grant from the U.S. Agency for International Development through the Sustainable Agriculture and Natural Resources Management Collaborative Research Support Program (SANREM-CRSP) and by State and Federal Funds allocated to Georgia Agricultural Experiment Stations Hatch Project GEO01654.
Open access
