Climate Change Impact and Variability on Cereal Productivity among Smallholder Farmers under Future Production Systems in West Africa

: Agriculture in West Africa is constrained by several yield-limiting factors, such as poor soil fertility, erratic rainfall distributions and low input systems. Projected changes in climate, thus, pose a threat since crop production is mainly rain-fed. The impact of climate change and its variation on the productivity of cereals in smallholder settings under future production systems in Navrongo, Ghana and Nioro du Rip, Senegal was assessed in this study. Data on management practices obtained from household surveys and projected agricultural development pathways (through stakeholder engagements), soil data, weather data (historical: 1980–2009 and ﬁve General Circulation Models; mid-century time slice 2040–2069 for two Representative Concentration Pathways; 4.5 and 8.5) were used for the impact assessment, employing a crop simulation model. Ensemble maize yield changes under the sustainable agricultural development pathway (SDP) were − 13 and − 16%, while under the unsustainable development pathway (USDP), yield changes were − 19 and − 20% in Navrongo and Nioro du Rip, respectively. The impact on sorghum and millet were lower than that on maize. Variations in climate change impact among smallholders were high with relative standard deviations (RSD) of between 14% and 60% across the cereals with variability being higher under the USDP, except for millet. Agricultural production systems with higher intensiﬁcation but with less emphasis on soil conservation (USDP) will be more negatively impacted by climate change compared to relatively sustainable ones (SDP).


Introduction
The increasing population of West Africa, which is projected to become three-fold by 2050 over the baseline value in 2010 [1] has implications for food security in the sub Sustainability 2021, 13, 5191 2 of 22 region. Currently, food security remains a major structural challenge in the sub region with undernourishment, particularly in the rural areas. One possible solution to this challenge is improving the productivity of the farming systems. Yet, efforts to improve farm productivity is constrained by several factors, such as the sub-optimal use of external inputs, relatively fragile and poor soils, overreliance on rain-fed agriculture, and the fact that the small-scale farmers that dominate the farming systems have diverse practices that handicap the promotion of blanket recommendations. Agricultural production in the subregion has been reported in several literature as vulnerable to environmental stress factors due to the heavy reliance on natural resources to support livelihoods [2][3][4]. In particular, changes or variability in weather patterns pose a concern. An analysis of long-term weather data (1950-2010) by Owusu et al. [5] in the Volta Basin in West Africa indicated an apparent shift in the rainfall pattern in Ghana toward a longer dry season and the disappearance of short dry spells. Additionally, Sultan et al. [6] reported significant projected increases in temperature across the entire West African sub region by 2050, which could alter not only the climatological hydrology but also crop physiological behaviour. For example, increased temperatures could accelerate crop development, shorten life cycles and also increase crop respiration. These effects would result in a reduction in crop yields. Overall, the projected changes in these weather variables would pose challenges to the livelihoods of farmers. Therefore, there is the need to assess the extent of the impact of the projected climate on crop productivity to provide a basis for the development of appropriate adaptation strategies.
A number of studies on the impact of climate change on crop yields have previously been carried out in the West African region. Notable among them include that by Traore et al. [7], who reported an adverse impact of climate change on maize and millet yields. Other studies projected that by 2050, local agricultural production in West Africa will fall below the per capita of 50 kg [1]. Simulation studies have projected an increase in the West African temperature by 1 • C [8], which would result in an increased frequency of heat and rainfall extremes [9] and a reduced average yield of millet by 10-20%, and sorghum by 5-15%. In monetary terms, this would translate to losses of USD 2.33-4.02 billion for millet and USD 0.73-2.17 billion for sorghum, annually [8]. An earlier literature review [3] also noted that West Africa was experiencing rapid climate change, typified by increasing temperatures and more rampant extreme climate events, which would lead to losses in crop yields. For the Oueme River Basin (ORB) in Benin, which is representative of the middle belt of West Africa, it was indicated that occurrences of drought and crop failure were likely to increase under climate change conditions [10]. A meta-analysis by Roudier et al. [11] indicated that increasing temperatures associated with climate change will likely have the most negative impacts on crop production, even though the temperature effects could be exacerbated by changes in rainfall patterns. It was concluded in the analysis that a rise in temperature beyond 2 • C cannot be counteracted by any positive changes in rainfall amount or distribution [12].
The issue of concern here is that most of these impact assessments are carried out either on current or past agriculture production systems. This approach fails to recognize that future production systems may differ due to technological advances in agriculture, agricultural policy and development pathways, and the fact that the highly diversified smallholder farming systems may adopt varying future agricultural practices. An implication of the diversity among farmers is that climate change and variability effects may not be uniform: some farmers may be severely impacted, based on their choice of crop mix and cropping practice, among others. On the contrary, it is conceivable that some farmers may gain from climate change impacts. With regard to technological change, it is feasible that new crop varieties, such as heat and drought-tolerant varieties may offset adverse climate change impacts [13]. Thus, farmers who adopt these varieties may gain from climate change impacts. To date, very few studies have quantified the variability of the impact among diverse smallholders and the way future farming productivity would be affected by climate change. Given that adaptation strategies have to be implemented at

Model Description
The crop model used for this study was the CERES-Maize, CERES-Millet and CERES-Sorghum modules of the Decision Support System for Agro-technology Transfer (DSSAT) [15]. The model is process-based and mechanistic, running on a daily time step. It utilizes data on weather (daily minimum and maximum temperature, rainfall and solar radiation), soil parameters (profile level), crop management practices and crop genetic information to simulate soil processes, crop growth and yield. The model simulates soil-water dynamics using routines developed by Ritchie [16], soil organic carbon dynamics as described by Porter et al. [17] and nutrient balance. Potential crop growth is influenced by photosynthetic capacity, capture of solar radiation, plant development based on thermal time units and photoperiod sensitivity, while actual growth is limited by soil water and nutrient stresses [18,19]. More details on the model description can be found in MacCarthy et al. [20].
Prior to this current study, the modules of the three crops studied were calibrated in previous studies. The maize (Zea mays L.) variety simulated for Ghana is "Obatanpa" which is an open-pollinated tropically adapted intermediate maturing cultivar developed by Ghana's Crop Research Institute (CRI) in collaboration with the International Institute of Tropical Agriculture (IITA), Ibadan, Nigeria, the International Maize and Wheat Improvement Centre (CIMMYT), Mexico and the Sasakawa Global 2000 [21]. In the case of Senegal, the TZEEY-SRBC5 (released by the National Agriculture Research Systems of Senegal in collaboration with the International Institute for Tropical Agriculture (IITA)) [22] was the selected variety. The CERES maize model was earlier recalibrated and evaluated for these varieties [14]. The sorghum variety ICSV III was registered for use in Ghana in 1998, based on research by the CGIAR Centres and the national research systems [23]. The pearl millet variety used is the CIVT, which is an improved variety developed by ICRISAT together with the national systems in the semi-arid regions of West Africa [24]. The sorghum (Sorghum bicolor (L.) Moench)) and pearl millet (Pennisetum glaucum (L.) R. Br.) modules were calibrated by [25] and [2]. The genetic coefficients of the calibrated varieties of sorghum and millet were, thus, available for use in this current study.
Data on crop management practices were obtained from a farm survey conducted in 2012 and 2007 at Navrongo and Nioro du Rip, respectively. The farm survey involved 87 and 169 farmers who cultivated maize and sorghum, respectively, in Navrongo, while in Nioro du Rip, the survey involved 98 and 222 farmers who cultivated maize and millet, respectively. Data collected included time of planting or planting window, amount and time of inorganic fertiliser application, quantity of manure applied and grain yield. A subset of fifteen (15) sorghum and thirteen (13) millet farms with complete set of information was then used to evaluate the performance of the model in simulating the observed yields. Sowing time varied from May to June in Navrongo and from June to July in Nioro du Rip. The soil profile information used for the simulation was derived from a soil map of the study area for Nioro du Rip, [26] while those for the Navrongo site were from published literature [14,27]. Phenology (anthesis and maturity) data from the literature and the previous maize yield simulation by Freduah et al. [14] also guided the parametrisation of the CERES crop modules (Appendix A Tables A1 and A2).

Climate Data
Historical climate data for Ghana over the period from 1980 to 2009 were obtained from the Ghana Meteorological Agency for Navrongo. In the case of Nioro du Rip, Senegal, data were obtained from AGRHYMET, Niamey, Niger. The weather variables used in the model were the daily minimum and maximum temperature, rainfall and solar radiation. The historical data were used to project future climate scenarios based on the AgMIP protocol [28]. For this, five General Circulation Models (GCMs) were selected from 29 GCMs which adequately captured the main features of the West African climate in terms of seasonal variability. The selection of the GCMs was based on a scatterplot of the average seasonal temperature vs. the seasonal rainfall for the locations, which was used to identify the cluster of GCMs that fall into each of the two Representative Concentration Pathways (RCPs), namely RCP4.5 and 8.5, using the AgMIP protocols described by Rosenzwieg et al. [28] and Ruane and McDermid [29]. For each of the 2 RCPs, the five GCMs corresponded to five categories: hot/wet, hot/dry, cool/wet, cool/dry and middle (Table 1). It is worth noting that the GCMs selected for the RCPs for a given scenario were not always the same (Table 1). Apart from temperature and rainfall, atmospheric CO 2 concentration was another factor considered in the yield simulation. The baseline CO 2 concentration of 390 ppm (corresponding to 2010 levels) was used for the historical climate, while 499 ppm and 571 ppm were used for future climate scenarios of RCPs 4.5 and 8.5, respectively, to correspond with AgMIP protocols [28]. Given that the current and future government policies in the agriculture sector will influence crop productivity irrespective of climate change, the concept of Representa- tive Agricultural Pathways [30] was applied in this study. Two RAPs were assumed; (i) an agricultural system with sustainable low growth characterized by high sustainability and associated environmental performance and policies (SDP), and (ii) an agricultural system characterized by unsustainable high growth (USDP). To quantify the policy effects, a multi-stage stakeholder (scientists, farmers, and policy makers) engagement was carried out to identify biophysical elements that would change in the future for each RAP, and the magnitude and direction of the changes were documented. The identified biophysical elements ( Table 2) were incorporated into the crop model to assess the impacts of the two management pathways (SDP and USDP) on crop yield. The main drivers of the biophysical component considered under the SDP RAP were fertiliser application, use of improved crop genotypes (Tables A1 and A2), improved crop management and soil health ( Table 2). The simulated enhanced SDP involved the application of an additional fertilizer (average of 30 kg N ha −1 ) above the current farmer application rates for maize in both study sites (Navrongo and Nioro du Rip). The amounts of fertiliser applied per farm were varied according to the observations of the surveys ( Figure A1), and planting density was maintained at 4 plants m −2 . For sorghum and millet, the farmer survey indicated that no fertiliser was applied; hence, the simulated enhanced management involved the application of 10 kg N ha −1 across farms, and the plant population was maintained at the survey values of 4 and 2 plants m −2 for sorghum and millet, respectively.
For the USDP RAP, an average of 45 kg N ha −1 was applied to maize at both sites (Navrongo and Nioro du Rip). As with SDP RAP, the amounts of fertiliser applied per farm varied based on amounts applied in the survey ( Figure A1). For the SDP simulated enhanced management, the plant population was increased from 4 to 5.5 plants m −2 . For sorghum and millet, the survey indicated that no fertiliser was applied by farmers. Therefore, the simulated enhanced management applied 15 kg N ha −1 across farms and the plant population increased from 4 and 2 plants m −2 to 5.5 and 3 plants m −2 for sorghum and millet, respectively. Furthermore, interactions of the multi-stage stakeholders indicated, that soil health would deteriorate under the USDP RAP. To capture this effect, the topsoil depth was reduced by 5 cm to account for surface erosion due to unsustainable practices. The cultivar modifications made to the various improved crops are indicated in Tables A1 and A2.

Simulation of the Impact of Climate and Agricultural Policy on Crop Yields
To assess the impact of the projected climate change on the yields and the variability of the impact among farms, the AgMIP protocol for integrated assessment was followed. This approach involves the simultaneous use of information from multiple farms, multiple climate scenarios and multiple soil data to simulate yields for the various farms and climates. In addition, the information derived from the RAPs was superimposed to capture the policy effects. The assessment therefore involved three (3) sets of simulations. First, the survey data were used to set up the model to simulate crop yields for current management practices by farmers under historical climate data . This showed how the crop yields for each of the farms (e.g., 87 farms for maize) would have varied if the current farmer practices were operated over the historical 30-year climate conditions. Secondly, data from the farm surveys were adjusted to capture the RAP effects identified by the multistakeholders for the SDP and USDP ( Table 2). The adjusted management conditions were then implemented in the model and the modified model was used to simulate yields over the 30-year historical period. This set of simulation runs is referred to as future production systems (without climate change). Thirdly, the RAP-modified model was used to simulate the impact of projected climate change on grain yield for the near-term using projected climate (2040-2069) data, derived from the five (5) GCMs and for the two (2) RCPs (4.5 and 8.5). The simulations for the SDP were carried out using RCP 4.5, while those for the USDP, RCP 8.5 was used, as SDP is associated with a lower emission pathway (RCP 4.5) while USDP is associated with a higher emission scenario (RCP 8.5). Planting windows were specified from May 1 to June 30 for Navrongo and from June 19 to July 31 for Nioro du Rip. Simulated planting was within 2 weeks intervals in Navrongo and weekly intervals in Nioro du Rip. For each simulation run, sowing was allowed when 25 mm of rainfall was cumulated over 3 rainfall events within the defined sowing window per farm ( Figure A2). Fertiliser was applied in two splits (10 and 32 days after emergence). Soil data used for the simulations for both sites are the same as those for the evaluation (Supplementary  Materials Table S1).
In summary, taking the case of maize in Navrongo, the first set of simulations involved the 30 years average yield simulated for each of the 87 farms based on management indicated in the farmer survey, referred to as Survey X i . The second set simulated the 30 years average maize yield for each of the 87 farms based on the RAPs-mediated management under the historical climate, referred to here as RAP v X i , (i.e., historical climate data + RAPs), where v represents the type of RAPs. The overall change in crop yield under both farmers' current situations and the RAPs-impacted yields was described as P i , derived as follows: Thereafter, the grand mean (GP RAP ) of change in grain yield across farms per RAP (SDP and USDP) was estimated in Equation (2), as follows: where N is the number of farms.
To capture the effect of the farm diversity response to climate change, maize yields were simulated using the RAPs-modified model for 30 years on each farm under future climate (2040-2069) and RCPs. The yield for each farm was averaged over the 30 years and called (FP v X i ). The latter showed the interactive effect of GCMs x RCPs/RAPs. The overall impact of climate change and RAPs on yield per farm (d ijk ) where i is farm number, j is the climate scenario and k is RCP/RAP (for the Navrongo maize example, i = 87, j = 5, k = 2). This was estimated in Equation (3) as follows: The impact of climate change per farm d ijk was then ensembled across the 5 GCMs per RAP to obtain d ik and the grand ensemble mean change in grain yield across farms (GM ik ) was estimated as follows: The coefficient of variation (CV) was used to estimate the variability in grain yield among farms for each crop under current (survey-based management) and future production systems (RAPs-based management) for the current (historical) climate as well as future projected climate (2049-2069). The CV was estimated as follows: The relative standard deviation (RSD) was used to estimate the variability in climate change impact among farms. It is an indicator of the dispersion of values around the mean and defined in Equation (6) as follows:

Sensitivity of Grain Yield to Weather Parameters
The sensitivity of grain yield to weather parameters (temperature, rainfall, CO 2 ) were assessed for the two study locations. The sensitivity of grain yield to each parameter was assessed using 30 years of historical weather data. The temperature was varied from −2, 0, +2, +4, +6, +8 • C. For rainfall, the changes of 0.25, 0.5, 0.75, 1, 1.25, 1.5, 1.75 and 2 were considered, and for CO 2 effect, the considered changes were 450, 540, 630 and 720 ppm.

Model Evaluation
The performance of the model in predicting grain yield were assessed using a set of statistics, namely, the Root Mean Square Error (RMSE), the Relative Root Mean Square Error (RRMSE) and the Coefficient of Model Efficiency.
The RMSE is defined in Equation (7) as follows: The lower the RMSE, the better the model performance. The minimum value of zero implies a perfect model performance. The Relative Root Mean Square Error (RRMSE) is defined as follows: The RRMSE values ranges between 0% and 100%. Values close to 0 are preferred. The Model Efficiency (E 1 ) was determined as follows: Sustainability 2021, 13, 5191

of 22
A Model Efficiency of 1 denotes a perfect agreement between the observed and simulated data, while E 1 = 0 signifies poor agreement between the observed and simulated data [31].

Model Performance
The CERES-Sorghum and CERES-Millet models adequately reproduced the phenology (anthesis and maturity) and grain yields of millet and sorghum measured from farmer fields. Anthesis was predicted with RRMSE of 10% and 13%, while duration to physiological maturity was reproduced with RRMSE of 16% and 14% for sorghum and millet, respectively. The statistics of the evaluation of the model's performance in reproducing grain yield are presented in Figure 1. The low RMSE suggests a generally satisfactory performance of the CERES models in reproducing grain yields from farmers' fields. For the millet cultivar at Nioro du Rip, the model slightly overestimated grain yields. The RRMSE showed that the yields of both crops were satisfactorily reproduced by the CERES-Millet and CERES-Sorghum models with RRMSE values of less than 30%. The coefficient of model efficiency (EF) was 0.77 and 0.81 for sorghum and millet, respectively, which further indicated that the model adequately captured crop yields ( Figure 1).
A Model Efficiency of 1 denotes a perfect agreement between the observed and simulated data, while E1 = 0 signifies poor agreement between the observed and simulated data [31].

Model Performance
The CERES-Sorghum and CERES-Millet models adequately reproduced the phenology (anthesis and maturity) and grain yields of millet and sorghum measured from farmer fields. Anthesis was predicted with RRMSE of 10% and 13%, while duration to physiological maturity was reproduced with RRMSE of 16% and 14% for sorghum and millet, respectively. The statistics of the evaluation of the model's performance in reproducing grain yield are presented in Figure 1. The low RMSE suggests a generally satisfactory performance of the CERES models in reproducing grain yields from farmers' fields. For the millet cultivar at Nioro du Rip, the model slightly overestimated grain yields. The RRMSE showed that the yields of both crops were satisfactorily reproduced by the CERES-Millet and CERES-Sorghum models with RRMSE values of less than 30%. The coefficient of model efficiency (EF) was 0.77 and 0.81 for sorghum and millet, respectively, which further indicated that the model adequately captured crop yields (Figure 1).

Projected Changes in Temperature and Rainfall
All GCMs for each RCP agreed on the direction of change in temperature for the time slice of 2049-2069. The ensemble changes in projected temperatures for the Navrongo site were 1.50 and 2.12 °C with CVs of 22% and 21% for RCP 4.5 and RCP 8.5, respectively, while those for Nioro du Rip, Senegal were 1.44 °C and 1.88 °C with CVs of 22% and 28% under RCP 4.5 and RCP 8.5, respectively. Temperature changes under RCP 8.5 were higher than those projected for RCP 4.5 at both sites ( Figure 2). Changes in the projected total rainfall varied both in magnitude and direction among the GCMs/RCPs at both study locations. For the Navrongo site in Ghana, ensemble changes in projected rainfall amount

Projected Changes in Temperature and Rainfall
All GCMs for each RCP agreed on the direction of change in temperature for the time slice of 2049-2069. The ensemble changes in projected temperatures for the Navrongo site were 1.50 and 2.12 • C with CVs of 22% and 21% for RCP 4.5 and RCP 8.5, respectively, while those for Nioro du Rip, Senegal were 1.44 • C and 1.88 • C with CVs of 22% and 28% under RCP 4.5 and RCP 8.5, respectively. Temperature changes under RCP 8.5 were higher than those projected for RCP 4.5 at both sites ( Figure 2). Changes in the projected total rainfall varied both in magnitude and direction among the GCMs/RCPs at both study locations. For the Navrongo site in Ghana, ensemble changes in projected rainfall amount were +2.8 and +4.6% for the RCP 4.5 and RCP 8.5, respectively. For the Nioro du Rip site, the ensemble change in the rainfall amount was −11.6% for the two RCPs. Projected ensemble changes in rainfall events for the five climate scenarios were −2.6% for each of the RCPs in Navrongo with CVs of 150% and 250% (for RCP 4.5 and RCP 8.5, respectively), whereas at Nioro du Rip, the projected ensemble changes in rainfall events were −8.6% and −16.2% with CVs of 157% and 127% under RCP 4.5 and RCP 8.5, respectively. the ensemble change in the rainfall amount was −11.6% for the two RCPs. Projected ensemble changes in rainfall events for the five climate scenarios were −2.6% for each of the RCPs in Navrongo with CVs of 150% and 250% (for RCP 4.5 and RCP 8.5, respectively), whereas at Nioro du Rip, the projected ensemble changes in rainfall events were −8.6% and −16.2% with CVs of 157% and 127% under RCP 4.5 and RCP 8.5, respectively.

Simulated Grain Yield of Cereals under Current and Future Management Practices
The distribution of the simulated average ensemble yields of the cereals (maize, sorghum and millet) for the current (historical) climate and current-SDP and USDPmediated management practices are illustrated in Figure 3. Simulated yields under the SDP and USDP were generally higher than the farmer survey-based management with yield distributions shifting towards higher grain yields. Simulated average maize yields under the current historical climate and current farmer-based production systems were 1212 and 934 kg ha −1 for Navrongo and Nioro du Rip, respectively. The inclusion of improved RAPsmanagement practices, for example, modest increases in fertiliser application and planting densities, among others, resulted in increased maize yields, reaching 2421 and 2165 kg ha −1 representing yield increases of 115% and 187% at Navrongo and Nioro du Rip, respectively for the SDP pathway. Under the USDP pathway where fertiliser applications were higher, ensemble average maize yields of 2683 and 2484 kg ha −1 were simulated, representing average yield increases of 144% and 261% for Navrongo and Nioro du Rip, respectively.

Simulated Grain Yield of Cereals under Current and Future Management Practices
The distribution of the simulated average ensemble yields of the cereals (maize, sorghum and millet) for the current (historical) climate and current-SDP and USDP-mediated management practices are illustrated in Figure 3. Simulated yields under the SDP and USDP were generally higher than the farmer survey-based management with yield distributions shifting towards higher grain yields. Simulated average maize yields under the current historical climate and current farmer-based production systems were 1212 and 934 kg ha −1 for Navrongo and Nioro du Rip, respectively. The inclusion of improved RAPs-management practices, for example, modest increases in fertiliser application and planting densities, among others, resulted in increased maize yields, reaching 2421 and 2165 kg ha −1 representing yield increases of 115% and 187% at Navrongo and Nioro du Rip, respectively for the SDP pathway. Under the USDP pathway where fertiliser applications were higher, ensemble average maize yields of 2683 and 2484 kg ha −1 were simulated, representing average yield increases of 144% and 261% for Navrongo and Nioro du Rip, respectively. The variability in the simulated yields among farms under current production systems was much higher (78%) compared to those obtained under simulated future RAPsmediated production systems (36% and 33% in Navrongo, 51% and 38% in Nioro du Rip for the SDP and USDP pathways, respectively). The response to improved management practices and cultivar also varied among farms. In Navrongo, the variations were similar under the two development pathways (31% and 34% for SDP and USDP, respectively) while those of the Nioro du Rip site were 66% and 55% for the SDP and USDP pathways, respectively. The variability in the simulated yields among farms under current production systems was much higher (78%) compared to those obtained under simulated future RAPs-mediated production systems (36% and 33% in Navrongo, 51% and 38% in Nioro du Rip for the SDP and USDP pathways, respectively). The response to improved management practices and cultivar also varied among farms. In Navrongo, the variations were similar under the two development pathways (31% and 34% for SDP and USDP, respectively) while those of the Nioro du Rip site were 66% and 55% for the SDP and USDP pathways, respectively.
The simulated average yield of sorghum under the current production system was 572 kg ha −1 . For the future RAPs-mediated production system, simulated yields increased to 1075 and 1285kg ha −1 under SDP and USDP, respectively, in response to improved management practices and improved cultivar. These resulted in yield increases of 124% and 174% under SDP and USDP, respectively. The higher yield increase under the USDP is due to the higher fertiliser application as well as the high plant population associated with this scenario ( Table 2). As with maize, variability in the yields of sorghum among farms reduced from 56% to 29% and 25% (SDP and USDP, respectively). The variability in the response to management practices was similar for the two development pathways (56% and 55% for SDP and USDP respectively).
Simulated average yield of millet was 585 kg ha −1 under the current farmer-based production system and current historical climate. This increased to 1210 and 1305 kg ha −1 when RAPs effects for SDP and USDP, respectively, were simulated. This translated into ensemble average yield increases of 114% and 135% for SDP and USDP, respectively. The variation in the response to the improved practices under SDP and USDP were 16% and 25%, respectively. The variations in the changes in ensemble simulated yields reduced from 50% under the current practices to 40% and 36% under SDP and USDP, respectively. Figure 4 illustrates the simulated ensemble yield distribution of maize, sorghum and millet under both the current (historical) and future climate for the two RAPs (SDP and USDP) conditions. The spread of ensemble yield distributions for the SDP and USDP were generally similar for all crops. However, the simulated yields under the current climate were noticeably always higher than those under the future climate, especially for Navrongo, irrespective of the RAPs scenario.

Impact of Future Climate Change on Grain Yield and Variation among Smallholder Farms
to 1075 and 1285kg ha −1 under SDP and USDP, respectively, in response to improved management practices and improved cultivar. These resulted in yield increases of 124% and 174% under SDP and USDP, respectively. The higher yield increase under the USDP is due to the higher fertiliser application as well as the high plant population associated with this scenario ( Table 2). As with maize, variability in the yields of sorghum among farms reduced from 56% to 29% and 25% (SDP and USDP, respectively). The variability in the response to management practices was similar for the two development pathways (56 and 55% for SDP and USDP respectively).
Simulated average yield of millet was 585 kg ha −1 under the current farmer-based production system and current historical climate. This increased to 1210 and 1305 kg ha −1 when RAPs effects for SDP and USDP, respectively, were simulated. This translated into ensemble average yield increases of 114% and 135% for SDP and USDP, respectively. The variation in the response to the improved practices under SDP and USDP were 16% and 25%, respectively. The variations in the changes in ensemble simulated yields reduced from 50% under the current practices to 40% and 36% under SDP and USDP, respectively. Figure 4 illustrates the simulated ensemble yield distribution of maize, sorghum and millet under both the current (historical) and future climate for the two RAPs (SDP and USDP) conditions. The spread of ensemble yield distributions for the SDP and USDP were generally similar for all crops. However, the simulated yields under the current climate were noticeably always higher than those under the future climate, especially for Navrongo, irrespective of the RAPs scenario.  Simulated ensemble average maize yield among farms under future climate in Navrongo were 1954 and 1899 kg ha −1 for the SDP and USDP, respectively. In Nioro du Rip, ensemble average yields were 1659 and 1639 kg ha −1 for the SDP and USDP, respectively. These indicate a general reduction in maize grain yields across the two sites. Maize yields declined by 13% and 19% at Navrongo and 16% and 20% at Nioro du Rip under SDP and USDP, respectively. Yield reductions were higher under the USDP than the SDP at both sites, probably due to soil degradation by erosion and higher projected temperatures in RCP 8.5. This effect was included in the USDP-mediated scenario. A simulated reduction in yield in Nioro du Rip were also higher than those obtained for Navrongo. The ensemble yield reductions varied among farms ( Figure 5) with an RSD of 14% and 17% in Navrongo and 32% and 36% in Nioro du Rip under the SDP and USDP, respectively. average negative impact (−4% and −9% under SDP and USDP, respectively) of climate change on sorghum yields under both RAPs. As with maize, the yield reduction was higher under the USDP than under the SDP. The impacts of climate change among sorghum farms varied considerably with RSD, of 38% and 42% for the SDP and the USDP, respectively, which were much higher than those obtained for maize at both sites. Thus, the variability in sorghum yields among the farms was lower under SDP than the USDP. Simulated average millet yields under the future climate were 1131 and 1137 kg ha −1 for SDP and USDP, respectively. As with maize and sorghum, simulated yields under the future climate were always lower than those for the current climate (used as the reference datum). The average change in the yield of millet under future production systems (SDP and USDP) at Nioro du Rip were −5% and −10% for the SDP and USDP, respectively. Thus, the yield decline was higher under USDP by two units. The variability in the climate change impact among millet farms was moderately high, with an RSD of 60% and 36% under SDP and USDP, respectively. As with maize, simulated millet grain yields were higher under USDP than under SDP.

Impact of Future Climate Change on Grain Yield and Variation among Smallholder Farms
The factors contributing to the variation in the climate change impact among farms are shown on Figure 6. The relative importance of the three factors analysed (soil, amount of fertiliser and planting window) varied among crops and between sites. The magnitude of the contributions from the three factors to the variation in the climate change impact among farms were higher for Nioro du Rip than for Navrongo. For sorghum and millet, fertiliser did not contribute to yield variations among farms, as similar quantities were applied by all farms. The magnitude of impact for the planting window at Nioro du Rip was more than three-folds of that obtained in Navrongo for maize, irrespective of the RAPs. The simulated ensemble average yields of sorghum under the future climate were 1028 and 1173 kg ha −1 for the SDP and the USDP, respectively. These translated into an average negative impact (−4% and −9% under SDP and USDP, respectively) of climate change on sorghum yields under both RAPs. As with maize, the yield reduction was higher under the USDP than under the SDP. The impacts of climate change among sorghum farms varied considerably with RSD, of 38% and 42% for the SDP and the USDP, respectively, which were much higher than those obtained for maize at both sites. Thus, the variability in sorghum yields among the farms was lower under SDP than the USDP.
Simulated average millet yields under the future climate were 1131 and 1137 kg ha −1 for SDP and USDP, respectively. As with maize and sorghum, simulated yields under the future climate were always lower than those for the current climate (used as the reference datum). The average change in the yield of millet under future production systems (SDP and USDP) at Nioro du Rip were −5% and −10% for the SDP and USDP, respectively. Thus, the yield decline was higher under USDP by two units. The variability in the climate change impact among millet farms was moderately high, with an RSD of 60% and 36% under SDP and USDP, respectively. As with maize, simulated millet grain yields were higher under USDP than under SDP.
The factors contributing to the variation in the climate change impact among farms are shown on Figure 6. The relative importance of the three factors analysed (soil, amount of fertiliser and planting window) varied among crops and between sites. The magnitude of the contributions from the three factors to the variation in the climate change impact among farms were higher for Nioro du Rip than for Navrongo. For sorghum and millet, fertiliser did not contribute to yield variations among farms, as similar quantities were applied by all farms. The magnitude of impact for the planting window at Nioro du Rip was more than three-folds of that obtained in Navrongo for maize, irrespective of the RAPs.

Impact of Climate Change on Duration to Flowering and Maturity
The impact of climate change resulted in the decline in the durations that plants took to flower and to mature in all crops (Figure 7). The magnitude of decline was generally higher for the USDP than the SDP. Additionally, the declines were higher for the crops at the Navrongo site than for those at Nioro du Rip. The decline in duration was also higher for the maturity than for the flowering alone, an indication that the duration of the reproductive stage also declined.

Impact of Climate Change on Duration to Flowering and Maturity
The impact of climate change resulted in the decline in the durations that plants took to flower and to mature in all crops (Figure 7). The magnitude of decline was generally higher for the USDP than the SDP. Additionally, the declines were higher for the crops at the Navrongo site than for those at Nioro du Rip. The decline in duration was also higher for the maturity than for the flowering alone, an indication that the duration of the reproductive stage also declined.

Sensitivity of Grain Yield to Temperature, Rainfall and CO2 Concentration
The sensitivity of the grain yield of the crops to changes in temperature, rainfall and CO2 are presented in Figures 8-10. The sensitivities generally varied annually, irrespective of the RAPs. Increasing temperatures resulted in a decline in grain yields, and this was more prominent for maize than sorghum or millet at both sites with a higher magnitude of losses. However, the variability in the sensitivity of grain yield to temperature was higher for sorghum and millet compared to maize. Grain yields of the cereals were also sensitive to varying total rainfall amounts, annually (Figure 9). Generally, grain yields

Sensitivity of Grain Yield to Temperature, Rainfall and CO 2 Concentration
The sensitivity of the grain yield of the crops to changes in temperature, rainfall and CO 2 are presented in Figures 8-10. The sensitivities generally varied annually, irrespective of the RAPs. Increasing temperatures resulted in a decline in grain yields, and this was more prominent for maize than sorghum or millet at both sites with a higher magnitude of losses. However, the variability in the sensitivity of grain yield to temperature was higher for sorghum and millet compared to maize. Grain yields of the cereals were also sensitive to varying total rainfall amounts, annually ( Figure 9). Generally, grain yields declined under increasing rainfall amounts beyond the base rainfall amount. The variability in the sensitivities were generally higher across the three crops and management scenarios when rainfall amounts reduced below the base total rainfall. For the maize in Nioro du Rip, the magnitude of the yield decline was generally higher than that of Navrongo. On the other hand, the magnitude for grain yield reduction was higher for sorghum in Navrongo than for pearl millet in Nioro du Rip. The magnitude of sensitivity of the grain yield to CO 2 was generally much smaller than those for temperature and total rainfall amount. Cereals may not show significant positive responses to CO 2 , but the differences in grain yield sensitivity to CO 2 between SDP and USDP were obvious ( Figure 10). For maize in Nioro du Rip, the sensitivity under SDP increased marginally, while that under USDP declined.   . Sensitivity of grain yield of cereals to varying temperature at two sites (Navrongo and Nioro du Rip). The box and whiskers plots represent variation in sensitivity annually over 30 years. 'Maz' is maize, 'Sgg' is sorghum, 'Pml' is pearl millet, and SDP and USDP are sustainable and unsustainable pathways, respectively. Figure 9. Sensitivity of grain yield of cereals to varying rainfall amounts at two sites (Navrongo and Nioro du Rip). The box and whiskers plots represent variation in sensitivity annually over 30 years. 'Maz' is maize, 'Sgg' is sorghum, 'Pml' is pearl millet, and SDP and USDP are sustainable and unsustainable pathways, respectively. Figure 9. Sensitivity of grain yield of cereals to varying rainfall amounts at two sites (Navrongo and Nioro du Rip). The box and whiskers plots represent variation in sensitivity annually over 30 years. 'Maz' is maize, 'Sgg' is sorghum, 'Pml' is pearl millet, and SDP and USDP are sustainable and unsustainable pathways, respectively.
Sustainability 2021, 13, x FOR PEER REVIEW 17 of 24 Figure 10. Sensitivity of grain yield of cereals to varying CO2 concentrations at two sites (Navrongo and Nioro du Rip). The box and whiskers plots represent variation in sensitivity annually over 30 years. 'Maz' is maize, 'Sgg' is sorghum, 'Pml' is pearl millet, and SDP and USDP are sustainable and unsustainable pathways, respectively.

Projected Changes in Temperature and Rainfall
The projected increases in temperatures by the selected GCMs are in agreement with those reported by Adiku et al. [2] and Sultan et al. [6] for the sub region. Changes in the projected total rainfall varied both in magnitude and direction among the GCMs/RCPs for both study locations. The inconsistencies in the direction of change in projected rainfall amounts by GCMs used in this study have been reported by other studies [2,9,32] for rainfall in West Africa. Figure 10. Sensitivity of grain yield of cereals to varying CO 2 concentrations at two sites (Navrongo and Nioro du Rip). The box and whiskers plots represent variation in sensitivity annually over 30 years. 'Maz' is maize, 'Sgg' is sorghum, 'Pml' is pearl millet, and SDP and USDP are sustainable and unsustainable pathways, respectively.

Projected Changes in Temperature and Rainfall
The projected increases in temperatures by the selected GCMs are in agreement with those reported by Adiku et al. [2] and Sultan et al. [6] for the sub region. Changes in the projected total rainfall varied both in magnitude and direction among the GCMs/RCPs for both study locations. The inconsistencies in the direction of change in projected rainfall amounts by GCMs used in this study have been reported by other studies [2,9,32] for rainfall in West Africa.

Sensitivity of Cereal Yields to Temperature, Rainfall and CO 2 Concentration
The marginal response of maize to CO 2 fertilisation is in conformity with those reported by other studies in Sub-Saharan Africa (SSA). The simulated response of maize yields to increasing rainfall was also similar for both sites with the grain yield increase progressively declining with an increasing amount of rainfall (Figure 8), a phenomenon that can probably be attributed to increase leaching of nitrogen under high rainfall. Changes in the weather parameters also have the potential to influence mineralisation of organic matter, thereby influencing nutrient availability. The impact of flooding on crop yields has not been simulated in this study but can become a major factor in climate change studies. The differences in sensitivity of the crops to temperature can be explained by the differences in the temperature threshold as noted by Tsusaka and Otsuka [33] in their analysis of the effect of temperature and rainfall on the yields of cereals in SSA, reporting on the differences in the sensitivity of crops to these weather parameters.

Cereals Yields under Improved Management Scenarios
The percentage yield increases simulated are within the yield responses reported for fertiliser trials that are conducted with the recommended management practices. In an on-station 4-year trial in a similar environment in northern Ghana, Naab et al. [34] reported an average maize yield increase of about 314% when 60 kg ha −1 N fertiliser was applied. A study by Akponikpè et al. [24] also reported a yield increase of 104% for millet when 40 kg N ha −1 was applied in a two-year experiment in Niger. In a recent yield simulation study by Adam et al. [35], it was reported that yield increases of sorghum in response to improved management practices in Koutiala (Mali) and Navrongo could be attributed more to improved management practices and that focussing on intensification in these (suboptimised) systems is currently more important than climate change impacts. Generally, the yield changes in maize, in response to management and genotype modification, were higher in Nioro du Rip than in Navrongo, even though the cultivar of maize used in Nioro du Rip had much shorter maturity duration and, hence, had lower yield potential. The length of the growing period in Nioro du Rip is also generally short, thereby limiting farmers to short and extra short maturity duration maize varieties. Additionally, the soil in Nioro du Rip is deeper and therefore has higher water storage capacity, while those in Navrongo are shallow and more prone to the leaching of nutrients. Therefore, it appears that the soil water and nutrient stress associated with shallow soils may have limited the response of crops to management practices in Navrongo.

Climate Change Impact on Crop Yields
In general, the simulated yield loss in Nioro du Rip was higher than in Navrongo. For instance, the yield loss under SDP in Nioro du Rip was three units higher than those simulated for Navrongo under a similar management scenario (SDP), even though ensemble temperature changes were lower in Nioro du Rip. This can be explained partly by the drastic reduction in the rainfall amount and events in Nioro du Rip where the length of the growing season was already much shorter than at Navrongo. The shorter maturity duration varieties used in Nioro du Rip compared to the intermediate variety used in Navrongo may also contribute to the higher yield loss in Nioro du Rip. A comparison of the historical and projected rainfall patterns show an extension of rainfall distribution in Navrongo into drier months in the historical data. Thus, late planting in the historical climate which encountered drought at the reproductive stage now benefited from the extension of rainfall into the drier months under the projected climate enabling the reproduction stages to coincide with high rainfall periods, ensuring better crop growth.
The short length of the growing season at Nioro du Rip compared to Navrongo may explain the fluctuations in yields as observed in this study, as farmers who plant late in Nioro du Rip are likely to obtain low yields in some years. Additionally, the projected decline in rainfall amounts and events are also more severe in Nioro du Rip, thereby contributing to higher uncertainty in yields among farms. Ensemble yield reductions were generally higher under USDP than under SDP. Thus, investments into crop production are likely to yield lower benefits under USDP than under SDP, in spite of the relatively higher yield increases. The more intensified smallholder systems become, the higher the impact of climate change [9,36]. The positive yield gains expected from CO 2 fertilisation were generally offset by the higher temperature effect on plant growth, hence resulting in low grain yields. According to Chemura et al. [37], the agroclimatic suitability of maize and sorghum are projected to decline in Ghana due to climate change by 2050. They reported that the optimally suitable land area (51,323 km 2 ) for maize production will decline by 14% under RCP 8.5. Similarly, for sorghum, the optimally suitable land area (66,731 km 2 ) is projected to decline by 13% in the same time frame under the same climate scenario. Given that the population in West Africa is projected to increase by about three-fold by 2050 [1], any decline in both the optimally suitable production area and productivity poses a further threat to food security, hence requires investments in appropriate adaptation measures to offset grain yield reductions.
As with the baseline climate, the variation in grain yields under SDP was lower than that in USDP. Higher yield variability is undesirable because this denotes yield uncertainty. The better soil condition under the SDP that ensured better water and nutrient retention within the season could account for this observation. It could also be due to lower planting density coupled with lower use of fertiliser. A number of studies have attributed higher variability in yield to increased fertiliser use [19,38]. Variability in yields among farms due to climate change were similar to those under the baseline climate.
Out of the three crops simulated, maize was the most vulnerable at both study sites under the projected climate change, irrespective of the type of management scenario. This result confirms those of Faye et al. [9] and Akumaga et al. [32] who reported that millet and sorghum were less sensitive to climate change impact than maize. Policymakers' attention needs to be drawn to this observation to inform their decision making, particularly regarding the choice of crops to promote for climate change adaptation and mitigation. Until recently, the production of maize in the semi-arid regions of West Africa has been rare. There is, however, increasing adoption of maize to the detriment of the other small grains (sorghum and millet), which hitherto were dominant in this region. These efforts have been spurred further by the advances in breeding efforts of both the CGIAR Centres and the National research systems, which resulted in the breeding of shorter maturity duration maize varieties as well as other varieties that are tolerant to stresses, such as drought, striga and nutrient deficiency.

Variation in Climate Change Impact among Smallholder Farms
This study is very unique in that it assesses the variation in climate change impact among smallholders with diverse crop management practices. Normally, a single (or few) crop management practice(s) is/are used to represent smallholder systems and a lumped single value of climate change impact is produced, a phenomenon that poses challenges for targeting adaptation planning processes which have to be implemented at the farm level. Variations in climate change impact are normally limited to the inter-annual variations, which Sultan et al. [6] reported to be on the rise in West Africa. Again, Adam et al. [35], indicated in their study that the variation in inter-annual yields was less important than the variation in yields among diverse smallholder farms in two locations in West Africa. In an earlier study, Freduah et al. [14] assessed the variability in the sensitivity of smallholder farmers to climate change and illustrated a wide variation in impact among maize farmers. This study builds on that by first developing future agricultural development (RAPs) trajectories through stakeholder engagements prior to assessing climate change impact.
The two factors that contributed to the variation in climate change impact among farms are the planting date or window and the type of soil on which crops were cultivated. Given that limited access to input credit and availability of labour for cultural practices [7] have been indicated as some of the reasons for late planting, policies must aim at addressing these so as to reduce the number of farms that are planted late. Since these farms are typically characterized with lower yields in these environment [7,39], such policy interventions will ensure higher grain yield. Another reason for late planting is due to the delayed onset of rain and false start of the season, which poses challenges for farmers. To our knowledge, this study is the first to take this approach and help in envisaging potential impacts and variations in climate change impact among smallholder farmers under the two representative agricultural pathways: sustainable development practices (SDP) and unsustainable development practices (USDP) co-developed with stakeholders.

Conclusions
This study provides a quantification of the climate change impact on yields of cereals (maize, sorghum and millet) in two semi-arid communities in Ghana (Navrongo) and Senegal (Nioro du Rip), as well as the diversity in the magnitude of climate change among smallholder farmers under two future production systems scenarios. The nature of future production systems will influence the magnitude of the climate change impact among smallholder farmers in West Africa. Agricultural production systems with higher intensification (including the use of higher fertiliser amount) but with less emphasis on soil conservation (USDP) will be more negatively impacted by climate change than those that are relatively sustainable (SDP). Out of the three cereals studied, maize was the most vulnerable. Thus, current policies promoting maize cultivation in the semi-arid regions of West Africa need to consider its relative vulnerability to climate change compared to the smaller grains (sorghum and pearl millet). Climate change impacts are generally varied widely among farms irrespective of future management scenarios. Variation in planting time was the main factor contributing to the diversity of climate change impact among farms, followed by the type of soil on which the crops were cultivated. The interaction between the planting date and climate change impact provides opportunity for agricultural scientists, as well as policy and decision makers, to develop appropriate adaptation packages that target the diversity in the smallholder landscape.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data is not yet publicly available due to efforts at publishing it in a data journal.

Acknowledgments:
The authors acknowledge the research support provided by AgMIP and Columbia University and the sub-grant provided to CIWARA for this study.

Conflicts of Interest:
The authors declare no conflict of interest.  Table A2. Genetic coefficients of sorghum and millet varieties used in the study and their corresponding improved varieties used in simulation of future production systems for each site.