Optimizing Traditional Cropping Systems Under Climate Change: A Case of Maize Landraces and Bambara Groundnut

Traditional crop species are reported to be drought-tolerant and nutrient-dense with potential to contribute to sustainable food and nutrition security within marginal production systems under climate change. We hypothesized that intercropping maize landraces (Zea mays L.) with bambara groundnut (Vigna subterranea (L.) Verdc.), together with optimum management strategies, can improve productivity and water use efficiency (WUE) under climate change. Using an ex-ante approach, we assessed climate change impacts and agronomic management options, such as plant ratios, and plant sequences, on yield and WUE of intercropped maize landrace and bambara groundnut. The Agricultural Production Systems sIMulator (APSIM) model was applied over four time periods; namely past (1961–1991), present (1995–2025), mid-century (2030–2060) and late-century (2065–2095), obtained from six GCMs. Across timescales, there were no significant differences with mean annual rainfall, but late century projections of mean annual temperature and reference crop evaporation (ET0) showed average increases of 3.5°C and 155mm, respectively. By late century and relative to the present, the projected changes in yield and WUE were −10 and −15% and 5 and 7% for intercropped bambara groundnut and maize landrace, respectively. Regardless of timescale, increasing plant population improved yield and WUE of intercropped bambara groundnut. Asynchronous planting increased yield and WUE for both maize landrace (5 and 14%) and bambara groundnut (35 and 47%, respectively). Most significant improvements were observed when either crop was planted 2–3 months apart. To reduce yield gaps in intercrop systems, low-cost management options like changing plant populations and sequential cropping can increase yield and WUE under projected climate change. To further increase sustainability, there is a need to expand the research to consider other management strategies such as use of other traditional crop species, fertilization, rainwater harvesting and soil conservation techniques.

Traditional crop species are reported to be drought-tolerant and nutrient-dense with potential to contribute to sustainable food and nutrition security within marginal production systems under climate change.We hypothesized that intercropping maize landraces (Zea mays L.) with bambara groundnut (Vigna subterranea (L.) Verdc.), together with optimum management strategies, can improve productivity and water use efficiency (WUE) under climate change.Using an ex-ante approach, we assessed climate change impacts and agronomic management options, such as plant ratios, and plant sequences, on yield and WUE of intercropped maize landrace and bambara groundnut.The Agricultural Production Systems sIMulator (APSIM) model was applied over four time periods; namely past (1961-1991), present (1995-2025), mid-century (2030-2060) and late-century , obtained from six GCMs.Across timescales, there were no significant differences with mean annual rainfall, but late century projections of mean annual temperature and reference crop evaporation (ET 0 ) showed average increases of 3.5 • C and 155 mm, respectively.By late century and relative to the present, the projected changes in yield and WUE were −10 and −15% and 5 and 7% for intercropped bambara groundnut and maize landrace, respectively.Regardless of timescale, increasing plant population improved yield and WUE of intercropped bambara groundnut.Asynchronous planting increased yield and WUE for both maize landrace (5 and 14%) and bambara groundnut (35 and 47%, respectively).Most significant improvements were observed when either crop was planted 2-3 months apart.To reduce yield gaps in intercrop systems, low-cost management options like changing plant populations and sequential cropping can increase yield and WUE under projected climate change.To further increase sustainability, there is a need to expand the research to consider other management strategies such as use of other traditional crop species, fertilization, rainwater harvesting and soil conservation techniques.

INTRODUCTION
Sub-Saharan Africa has a dualistic food system with the formal system taking a more national focus, and also focused on a few strategic crops while the informal system supports local food systems, which support household food and nutrition security (Tcoli, 2016;Mabhaudhi et al., 2019a).While several nations are food secure at a national level, household food insecurity remains problematic with an estimated 821 million people currently food insecure and malnourished (Abegaz, 2018;Gashu et al., 2019;Xie et al., 2019).Most of these people rely on agriculture as their mainstay; thus, the importance of agriculture within these communities provides an opportunity to improve food and nutrition security, reduce poverty, and enhance rural economic development [New Partnership for Africa's Development (NEPAD), 2014].However, current crop yields are low and challenged by worsening land degradation, especially declining soil fertility (Ukeje, 2010;Rippke et al., 2016;Badu-Apraku and Fakorede, 2017), and low water use efficiency (WUE) (Mabhaudhi et al., 2018b;O'Leary et al., 2018;Nouri et al., 2019).Furthermore, climate variability and change are adversely affecting productivity through increased incidences and intensity of droughts (Mpandeli et al., 2018;O'Leary et al., 2018;Nhamo et al., 2019).There is consensus that rural agricultural systems must increase resource use efficiencies and adopt strategies to adapt to climate risk (Isaacs et al., 2016;Matthews and McCartney, 2018).
A considerable amount of literature depicts the adoption of improved technologies such as the use of high yielding, improved crop varieties (Hammer et al., 2014;Ran et al., 2017;Mabhaudhi et al., 2019a).However, marginalized farmers have experienced several challenges when trying to adopt conventional farming practices.Chief among these include inadequate access to agrochemicals, loss in agro-biodiversity and an increase in the vulnerability of the system to climate risk (Mabhaudhi et al., 2019b;Malik and Chaudhary, 2019).The low adoption and consequent challenges have partly contributed to the widening gaps in food and nutrition security (Midega et al., 2015;Mrema et al., 2018).Within the context of marginal systems, agriculture needs to sustainably contribute to food and nutrition security and rural economic development, while reducing negative impacts on the environment or improving the environment (van Ittersum et al., 2016).Demand for more sustainable agriculture, which is less dependent on external inputs and better suited to marginal environments, has revived interest in traditional systems (Keatinge et al., 2015;Govender et al., 2016;Saharan et al., 2018).In line with this, there is a renewed focus on the inclusion of neglected and underutilized crops (NUS) as alternative crop choices in marginal cropping systems (Mabhaudhi et al., 2019a).
Neglected and underutilized crops, also referred to as underutilized indigenous and traditional crops, are defined as "plant species that are part of more substantial biodiversity, were once popular (in and out of their centers of diversity), and are neglected by users and research but remain relevant in the regions of their diversity" (Dansi et al., 2012).They are associated with high nutritional value, adaptation to marginal soils, and tolerance to drought and heat stresses (Slabbert et al., 2004;Chibarabada et al., 2015;Chimonyo et al., 2016a;Hadebe et al., 2017;Mabhaudhi et al., 2017).They often require fewer inputs such as fertilizer and agrochemicals, as they are also tolerant of several pests and diseases (Mabhaudhi et al., 2019a).Their nutritional attributes and adaptability make them suitable crops for promotion in marginal areas where poverty and food and nutrition insecurity remain high; however, their contribution to mainstay agriculture remains low (Massawe et al., 2016).As is reflected by their name, the potential of underutilized crops has not yet been fully harnessed, but most of them contribute to diversification and resilience of agroecosystems.Therefore, they have the potential for future agriculture under adverse agroclimatic conditions (Padulosi et al., 2011).Many proponents of modern agriculture and the Green Revolution have discouraged their continued production, highlighting low productivity and resource use efficiencies (Tokatlidis and Vlachostergios, 2016;Missio et al., 2018).For example, water use efficiency of bambara groundnut was reported to be 0.45 kg ha −1 mm −1 compared to 0.89 kg ha −1 mm −1 for groundnut (Chibarabada et al., 2017), while landrace sorghum varieties had 20% less WUE relative to hybrid varieties (Hadebe et al., 2019).However, the argument is to not promote them as replacement crops for high yielding major crops, but as complementary crops (Mabhaudhi et al., 2019a), especially in marginal areas where the major crops may not perform well (Massawe et al., 2016).Within these areas, NUS have potential to contribute to improving rural livelihoods and may be "better bet" technologies; however, this potential remains largely untapped due to limited information detailing their genetic, eco-physiological and agronomic performance (Chivenge et al., 2015).It is against this backdrop we hypothesize that, by optimizing resource use, yields of NUS can be sustainably increased.Intercropping involves growing of two or more crops simultaneously or overlapped on the same piece of land, which can sustainably increase WUE (Martin-Guay et al., 2018).
In this study, we hypothesize that intercropping a maize landrace (Zea mays L.) with bambara groundnut (Vigna subterranea (L.) Verdc.) is beneficial because the latter's smaller canopy offers little competition to the cereal crop (Saxena et al., 2018).As a legume, bambara groundnut also fixes atmospheric nitrogen.It contributes to soil fertility (Sprent et al., 2010), and the low cost of bambara groundnut seed makes it an exemplar crop for enhancing food and nutrition security within cereal producing households (Muhammad et al., 2016;Mayes et al., 2019).While traditional cropping systems featured multicrops (Muzari et al., 2012), intercropping maize with bambara groundnut is no longer a common practice.Little information is known about crop interaction and the impacts of climate variability and change on productivity and water productivity.While intercropping, in general, could be considered positive in terms of yield (Martin-Guay et al., 2018), the performance of each crop in an intercrop system is determined by the interaction between different crops and the availability of resources.With the impacts of climate variability and change, adapting agronomic management in response to changing resources can allow for sustainable intensification of the traditional cropping systems through improved resource use efficiency.Using an ex-ante approach in APSIM, the current study assessed the productivity and water use of a maize landrace -bambara groundnut intercrop under changing climate and in response to different management options.APSIM has been used widely to study impacts of climate change on crop growth and productivity across Africa (Beveridge et al., 2018;Duku et al., 2018;Xiao et al., 2020).However, its application for studying intercrop systems remains scanty, with no known research on its application for climate change studies.

Study Area
The study area was the University of KwaZulu-Natal's Ukulinga Research Farm (29 • 40 ′ S; 30 • 24 ′ E; 809 m a.s.l.).Ukulinga Research Farm is classified as semi-arid with 77% of the mean annual rainfall of 750 mm received mostly between October and April.The summer months are warm to hot, with an average temperature of 26.5 • C (Kunz et al., 2015).Soil textures are characterized as predominantly clay to clay loam and are moderately shallow, ranging from 0.6 to 0.8 m (Chimonyo et al., 2016a).

Brief Description of the APSIM Model
The APSIM version 7.10 is a daily time step, field-scale multiyear, a multi-crop model that provides an analytical tool for assessing the impacts of climate, soil factors and farming management on cropping system production (Holzworth et al., 2014).The model is driven by daily temperature, precipitation, and solar radiation and is capable of simulating soil carbon (C), soil water, phosphorus (P), and nitrogen (N) dynamics and their interaction (Keating et al., 2003).Management practices include sowing date, variety selection, irrigation water management, fertilizer application, crop residue management, crop rotations and conservation tillage; this makes the model ideal for assessing the impacts of various management options on resource use.APSIM also allows users to set up atmospheric CO 2 concentration (Jones et al., 2001), which is ideal for assessing climate change impacts.Furthermore, through the CANOPY module, the model can simulate resource use within intercrop systems.For detailed information on the technical workings of the APSIM model, refer to McCown et al. (1996), Dimes andRevanuru (2004), andHolzworth et al. (2006;2014).
The CANOPY module determines resources intercepted by each component of the intercrop using leaf area index (LAI), extinction coefficient and height for each crop.Arbitration for water and nitrogen uptake is done based on the module changing the order each day (on a rotational basis) in which the competing species are allowed to capture soil resources.Through the CANOPY module, the model accounts for the vertical profiles of LAI in different species in a mixture (Keating et al., 2003), and assumes a horizontally homogeneous canopy for each species (Gou et al., 2017).The CANOPY module has been published and successfully applied by Smith et al. (2016) and Snapp et al. (2018) for maize and pigeon pea; Carberry et al. (1996) for maize and bean; Chimonyo et al. (2016a) for sorghum and cowpea;and Hoffmann et al. (2020) for various maize intercrop systems.Although Nelson et al. (1998a,b) used APSIM to simulate a maize and Desmanthus virgatus intercrop system, the two crops were grown as monocultures and did use the CANOPY module.It was not clear whether Amarasingha et al. (2017) used the CANOPY module when maize and mung bean intercrop systems were simulated in APSIM.In contrast, Knörzer et al. (2011) found that APSIM was unable to simulate wheat-pea and maize-pea intercropping systems in Germany because it strongly underestimates the competitive ability of the species that was planted the first relative to the one that was planted last.In this study, we used the CANOPY module to simulate the effects of climate change on a maize landrace and bambara groundnut intercrop system.The current study, therefore, adds to the existing body of knowledge on the use of ASIM in simulating intercrop systems.It goes further to simulate different management options under different climate change impacts on the intercrop system.

Model Calibration, Testing and Application
The calibration and testing of the APSIM were carried out using observed data obtained from field experiments conducted during the 2015/16 growing season for a maize landracebambara groundnut intercrop established at the University of KwaZulu-Natal's Ukulinga Research Farm.Sub-plots comprised intercrop combinations, that is, sole maize landrace, sole bambara groundnut and maize landrace-bambara groundnut intercrop.The irrigated treatments were used for calibration, while the rainfed treatments were used to validate the model.For a detailed description of the experiment, refer to Supplementary Information 1.The simulation files were, therefore, created using observed data collected from the rainfed and irrigated treatment.

Met File
For model calibration and testing, a 10-year (2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017)(2018)(2019) weather data file that contained daily estimates of rainfall, minimum and maximum temperatures, solar radiation and reference evapotranspiration was sourced from SASRI weather site (http://sasex.sasa.org.za/irricane/tables/Ash_tables_AR.pl) using the nearest station to the location except for Ukulinga where there is a weather station on-site.With the 5-year climate file, we were able to back-calculate and estimate the initial soil water and initial soil nitrogen at planting.APSIM require an average ambient temperature (TAV) and the annual amplitude in monthly temperature (AMP).These values are calculated using long-term daily minimum and maximum temperatures by software program named "tav_amp."

Soil File
The soil file was generated using soil details at Ukulinga Rresearch Farm.Soils at the research farm have been described as being shallow clayey to clayey loam with medium fertility (Mabhaudhi et al., 2013).The soil file selected to represent this description best was Clay_Shallow_MF_101 mm.The soil module was created using information obtained from Chimonyo et al. (2016b) (Table 1), and this was matched to a pre-existing soil file available in APSIM soil module -Africa (Generic).

Crop Files
Within maize APSIM crop file, we used the maize cultivar "mwi_local" as it best described the maize landrace used in terms of days to maturity and yield potential of 3 t ha −1 .However, slight iterations to genetic coefficients were done using an iterative approach until simulated values were within 9-20% of observed values (Table 2).Since APSIM does not have a bambara groundnut crop file, the groundnut cultivar "kangwana" was modified as it closely resembled bambara groundnut in terms of physiology, growth habit and phenology (Table 2).The groundnut crop module was iterated by first adjusting the reproductive parameters within the crop life cycle (phenology, e.g., time to emergence, first leaf, reproductive stages, and maturity) to resemble what was observed from the monocropped treatment during the field experiment.After that, where simulations disagreed with observations, parameters in the groundnut module were modified in a sequential approach following the order proposed by Boote et al. (2002).The steps were: (1) leaf appearance rate, canopy height, and width, (2) specific leaf area, leaf area index, and partitioning among vegetative organs, including the rate of total biomass accumulation and lastly, (3) onset, rate, and duration of pod addition and seed growth.Besides modifications based on comparisons with the observed data, some parameter modifications were made based on a literature review.

Management File
The management file considered planting date, plant densities, fertilizer rate, irrigation and harvest rules.The plant populations used to calibrate and test the model were 2 0 and 2.2 (plants m −2 ) of the maize landrace and bambara groundnut, respectively.The plant population used represented the densities observed in the field experiment and were less than the recommended densities for dryland maize (2.6 plants m −2 ) and bambara   groundnut (4.4 plants m −2 ) production (Jensen et al., 2003).
Since the field experiment used to calibrate and test the model was conducted in one season, we used the irrigated treatments for calibration and the rainfed treatments to test model.The module "irrigate on the date" was used to apply irrigation on dates corresponding to actual irrigation dates.Observed irrigation applied per event for the field experiment was calculated to be, on average, 15 mm, which was applied thrice during the experiment.
Nitrogen fertilizer was applied automatically within 50 cm depth in the soil at a rate of 50 kg ha −1 to avoid any nitrogen stress.

Climate Scenarios
Ukulinga Research Farm is located within quinary subcatchment 4,697 of quaternary catchment U30J (Schulze et al., 2011).In addition to historical data, the study also used downscaled future climate projections for the Ukulinga quinary.
The climate projections were developed by the Council for Scientific and Industrial Research (CSIR) (Table 3) using output from six global climate models (GCMs) from the CMIP5 archive that was forced by Representative Concentration Pathway 8.5 (RCP 8.5).The climates produced under RCP 8.5 were used as they represent the most extreme scenarios.The selection of these six GCMs was based on their ability to provide a reasonable representation of the El Nino-Southern Oscillation (ENSO) phenomenon for the region.The climate projections were dynamically downscaled to improve spatial resolution to 0.5 • (∼50 km) using the CCAM regional climate model developed by the Commonwealth Scientific and Industrial Research Organization, CSIRO (McGregor andDix, 2001, 2008;McGregor, 2005).After that, a multiple-nudging strategy was followed to obtain a downscaling  Mabhaudhi et al., 2018a).Climate scenarios were then extracted for the gridded pixel that overlapped quinary sub-catchment 4,697.For application in crop modeling at a local scale, it is necessary to correct for systematic and localized biases in rainfall and temperature projections produced by the climate models.When compared to observed rainfall data from the historical quinary climate database for sub-catchment 4,697, the downscaled climate projections were found to have a substantially larger number of rain days, with many rain days having minimal rainfall depths (i.e., < 0.1 mm).Therefore, a quantile delta mapping method, as described and assessed by Cannon et al. (2015), was applied to bias correct the climate scenarios using a multiplicative factor for rainfall and an additive factor for temperature.
The bias-corrected climate data provide daily rainfall and temperature scenarios for a continuous period from 1961 to 2100.Daily reference crop evaporation (ETo) estimates were then computed as described for the historical data set (see Schulze et al., 2011).Solar radiation for each GCM for Ukulinga was then calculated as described by Schulze and Chapman (2007).The climate database, therefore, satisfied APSIM's climate file input requirements and was used to develop projections for the past (1961-1991), present (1995-2025), mid-century (2030-2060) and late-century (2065-2095) periods.Throughout the analysis, the "present" timescale was regarded as the baseline.

Management and Agronomic Scenarios
Two management scenarios were used to develop recommendations for best management practices.The scenarios were as follows.
Scenario 1: Planting Dates Maize production guidelines published by the Department of Agriculture, Forestry and Fisheries suggest that maize should be planted between October 1 and mid-December throughout South Africa [Department of Agriculture Forestry Fisheries (DAFF), 2003].As it is, South Africa exhibits a wide variation of agroecologies, both at the micro and macro level.Due to climate variability and change, this variation has increased, and there is an observed increase in the land area occupied by semi-arid arid agro-ecologies since 2000 (Cairns et al., 2013).Conversely, there is a continual need to redefine planting dates.In this study, we adopted five fixed dates between September 15 to January 15 as this approach is much easier for farmers to use.These dates were assumed to represent early to late planting.However, a significant weakness of this approach is the need to redefine the dates because of continuous shifting in agro-ecologies.

Scenario 2: Plant Populations
Model simulations were performed using plant populations that were 50% less to 50% more than the recommended values.Simulations were carried out by maintaining the recommended plant population of one component and changing the other.The total number of simulations was a 3 by 3 factorial with maize populations of 13,000, 26,000, and 39,000 plants ha −1 and bambara groundnut populations of 6,500, 13,000, and 19,500 plants ha −1 .The lower populations would reduce resource competition and improve productivity for either component crop, while higher populations assumed that there was a need to minimize unproductive resource use from the system and improve their productive use.From this, optimum plant populations were determined for both landraces.

Model Runs
For model calibration and testing, the APSIM intercrop model was run for 10 consecutive years from 2009 to 2019.The 10year run allowed for soil conditions to stabilize around what was observed in the actual experiment.During the scenario analyses management options were run independently from each other across the six climate projections to minimize the interactive effects of the scenarios.The RCPs were run continuously from 1961-2095 periods.

Data Analyses
Since APSIM does not calculate WUE directly, simulated outputs of water use (WU in mm) and yield (Y in kg ha −1 ) or biomass (B in kg ha −1 ) were used to determine water use efficiency (WUE in kg mm −1 ha −1 ) over the growing season (sowing to harvest) as follows: Within the model, WU was determined as crop water uptake from the soil profile by either maize landrace or bambara groundnut crop, i.e., maize Ep and bambara Ep and soil evaporation Es [Ep (for either maize landrace or bambara groundnut +Es)].
For model calibration and validation, model performance was evaluated by comparing simulated (S) vs. observed (O) values for phenology, leaf area index, WU, WUE B , grain yield and biomass.Model performance was evaluated using the coefficient of determination (R 2 ), root mean square error (RMSE) and normalized RMSE (nRMSE).Values of R 2 range between 0 and 1 with high values indicating less error variance.Since the interpretation of R 2 is independent, low values are only acceptable if n is large.However, R 2 values are sensitive to outliers and insensitive to additive, and proportional differences between S and O.The simulation was considered excellent when nRMSE < 10%, good if 10-20%, acceptable or fair if 20-30%, and poor if >30% of the observed mean (Jamieson et al., 1991;Granderson and Price, 2014).
Simulation outputs for yield and water use were subjected to descriptive statistics, t-test analyses and generalized linear mixed analysis (GLMM) using R statistical software (version 3.6.0).Descriptive statistics such as means, standard deviations, bubble charts and box and whisker plots were used to analyse outputs.Box and whisker plot can show stability and general distribution of data sets.The GLMM was used to identify significant factors influencing maize landrace and bambara groundnut yield.

Developing Guidelines
The Food and Agriculture Organization (FAO) suggested a list of guiding questions to review transformative elements within an intervention (Carter et al., 2018).These questions are meant to provide clarity on the adaptation planning process.In this study, we adopted selected questions to assess the implications of the research, provide actionable recommendations and provide a way forward.Key findings were summarized in Table 5 and implications outlined.

Model Performance
Comparisons of simulated and observed values for maize landrace and bambara phenology and LAI, and biomass, yield and water use (WU) and water use efficiency (WUE) are given in Figure 1 and Table 4, respectively.For phenology, the close alignment of the points to the 1:1 line indicates that the model was able to simulate the maize landrace and bambara phenology correctly.The model could explain more than 90% of the variation of either crop in phenological stages (Figure 1).During the calibration, the nRMSE for the system LAI was < 10% of the observed LAI for the maize landrace and bambara groundnut intercrop system.The nRMSE for the system LAI during model validation increased slightly to 14%; this implied good simulation for the intercrop system grown under rainfed conditions.Reasonable simulations of crop water use (WU) by the model during calibration, were also observed (RMSE = 41 mm); however, during model validation, WU was over-estimated by 48%.The output suggests that the APSIM model might not be sensitive to water.A closer look at model outputs for maize landrace and bambara groundnut simulated under irrigated (used for calibration) and rainfed (used for validation) conditions showed that transpiration was mostly unaffected by the reduction in water availability.In nature, the low availability of water results in a reduction in transpiration due to the reduction in stomatal conductivity.Field results showed no significant differences between the irrigated and rainfed treatments (see Supplementary Information 1).In this case, the model appropriately captured maize landrace and bambara groundnut physiology.
For biomass and grain yield, the model tended to overestimate the outputs for maize landrace and bambara groundnut.During model calibration, simulated yield and biomass for maize landrace as 11 and 16% higher than observed, and this implied reasonable simulation.Model simulation of maize landrace yield and biomass under rainfed conditions were satisfactory (RSME = 49 and 267 kg ha −1 ) However, simulated yield and biomass for bambara groundnut were 32 and 55% higher than those for observed yields.The performance of APSIM would suggest that, for improved model simulations, additional parameterisation may be required to simulate bambara groundnut adequately.The WUE calculated based on model simulated biomass (WUE B ) of both the maize landrace and bambara groundnut showed a good fit (1 and 4 kg mm −1 ha −1 , respectively), for simulated and observed results (Table 4).Then again, the bambara cultivar used to calibrate the crop file was a landrace selection.It could be that performance under low water availability had adverse effects on its productivity, and model did not capture this response.Considering that the model was still able to simulate low yields TABLE 4 | Calibration and validation results for observed and simulated outputs for maize landrace and bambara groundnut intercrop system for final biomass (kg ha −1 ), yield (kg ha −1 ) and intercrop system water use (mm).for the maize landrace and bambara groundnut and possible errors in the observation data (e.g., iterated cultivar parameters), the APSIM model performance was considered to be acceptable for the simulation of the intercrop system.

Change in Climate During the Growing Season
Dynamically downscaled and bias-corrected climate projections for six GCMs forced by RCP 8.5, together with an impact model APSIM 7.7, were used to simulate bambara groundnut and maize landrace yields for past, present, mid-, and late-centuries.The primary aim was to assess how climate change impacts on yield, WU and WUE of the maize landrace and bambara groundnut intercrop.Secondary to that, we assessed the impacts of various management options on mitigating the impacts of climate change.The median value of climate change projections for minimum and maximum temperatures for Ukulinga showed a consistent warming trend across all months from past to late century.Figure 2 indicates a warmer future (mid-and latecentury) with mean maximum temperature increasing by 4.5 • C relative to the baseline maximum temperature of 24 • C.This suggests an increased probability of heat stress, especially for the maize landrace.This warming trend across the selected timescales is consistent with projected trends for South Africa (Mangani et al., 2018).
The six GCMs project an increase in mean minimum temperatures in the future (mid-and late-century) that ranges from 2.0-4.8• C from a minimum baseline temperature of 13 • C. The projected increases suggest an increased probability of hot nights and longer and more frequent heatwaves.The warmer temperatures may result in a faster accumulation of heat units and a reduction in growth duration and accumulation of photosynthesis and increase in night-time respiration, all resulting in reduced crop yield (Schlenker and Roberts, 2009).Unlike bambara groundnut (C3 plant), maize (C4 plant) generally originates from warmer climates (Leff et al., 2004;Jia et al., 2016) and thus, may be more resilient to projected increases of temperature (Choudhary et al., 2019).Then again, for bambara groundnut, optimum temperatures range between 28 and 35 • C and the lethal temperature has been reported to be 50 • C (Soni et al., 2015).The wide temperature adaptation makes the crop ideally suited for building resilience to cropping systems located in areas where temperature increases have been projected.
Results across the GCMs show that mean annual rainfall (MAP) for the future is projected to remain somewhat unchanged (Figure 3).for the late-century period, data showed that ACC and CCS predict a 10.6 and 8.3% increase in MAP, respectively, while slight reductions of 3.5 and 2.5% are predicted by CNR and NOR respectively.However, the more extended box and whisper plots for ACC predict an increase in the inter-annual variability of mean rainfall (750 mm) (Figure 3).This suggests an increase in the probability of extreme weather events such as drought and floods.In all instances, projected ET o was observed to be higher (35%) than projected rainfall and is set to increase in the future (mid and late century) (Figure 3).In this regard, the rainfall: ET o ratio is projected to decrease in the near future.The increase in ET o is consistent with the projected increase in minimum and maximum temperature and suggests an increase in crop water stress (Zhao et al., 2017).Then again, intercrop systems with cereals and legumes are advantageous as the cereal over-story can lower canopy temperature and minimize evaporative losses (Eskandari, 2011;Chimonyo et al., 2016b).The modification of microclimate within intercrop systems makes it an ideal system to mitigate against projected temperature and ET 0 increases.Our results suggest that, while the tolerances of traditional crops to high temperatures may vary, intercropping crop species with different physiological and morphological traits can be a strategy to increase the resilience of marginalized production systems to projected temperature ET 0 increases.

Yield, Water Use and Water Use Efficiency
Across the GCMs, yield trends for intercropped bambara groundnut showed a gradual reduction toward the late century by 24% when compared to the baseline yield of 365 kg ha −1 .The observed trend for simulated bambara groundnut yield was late-century (285 ± 57) < mid-century (323 ± 62) < present (365 ± 67) < past (450 ± 65 kg ha −1 ) (Figure 4).Across the GCMs, the magnitude of change in simulated bambara groundnut yield during the mid-and late-century periods was consistent with corresponding projected increases in ET o and temperature.On the other hand, the mean yield trends for intercropped maize landrace across the GCMs and time scales were inconsistent with projected increases in ET o and temperature.The observed trend for simulated maize landrace yield was past (845) < late-century (855) < present (923) < mid-century (967 kg ha −1 ).For mid-century, although the model predicted a slight increase in maize landrace yield, results also showed larger yield variations relative to past and present.Standard deviations for the intercropped maize landrace yield were past (288) < present (363) < mid-century (351) < late-century (436 kg ha −1 ) (Figure 4).These results are in line with the increased probability of extreme weather events such as drought and floods (Schulze, 2011).Although bambara groundnut yield decreased across the time scales, the magnitude of yield variations within each timescale and GCM was somewhat consistent with an average standard deviation of 63 kg ha −1 .Within each timescale, our results would suggest that yields of bambara groundnut are more stable to climate fluctuation; however, it could be more sensitive to significant climate changes.
The trend for water use in the intercropped maize landrace and bambara groundnut was inconsistent across the GCM and  timescales.Overall, CCS predicted the highest water use (265 and 253 mm), while the lowest was under NOR (242 and 235 mm) for maize landrace and bambara groundnut, respectively (Figure 5).Differences in simulated WU across the GCMs could be that each climate model has been developed based on its assumptions and unique mathematical representations of physical climate system processes, providing different climate projections (Confalonieri et al., 2016).There were slight reductions in WU across the timescale; however, based on the pairwise t-test analysis, the reductions were not significant (P > 0.05).On the other hand, simulated results for crop water use efficiency (WUE) for intercropped bambara groundnut showed a reduction across time scales.The trend was such that past (1.78 ± 0.45) > present (1.52 ± 0.41) > mid-century (1.37 ± 0.38) > late century (1.16 ± 0.42 kg ha −1 mm 1 ) (Figure 6).The observed trend was consistent with the observed reduction in future yield.Large inconsistencies were observed for maize landrace WUE across GCMs and time scales.For example, CNR predicted the highest water use (4.01 ± 1.98 kg ha −1 mm −1 ), while the lowest was under NOR (3.25 ± 1.05 kg ha −1 mm −1 ).The trend for maize landrace WUE across the timescale was such that present (3.58 ± 1.25) < past (3.62 ± 0.81) < late century (4.37 ± 1.38) < mid-century (4.56 ± 1.82 kg ha −1 mm −1 ) (Figure 6).The observed trend was consistent with the simulated improvements of maize yield.

Impacts of Planting Density on Yield and Water Use Efficiency of a Maize Landrace and Bambara Groundnut Intercrop System
Simulation results of yield and WU for intercropped maize landrace and bambara groundnut across the six GCM were not significantly (P > 0.05) different; therefore, the results presented in this section are average values across the six GCMs.Across timescales, the trend for maize landrace yield was past (850 ± 288) < late-century (853 ± 443) < present (893 ± 359) < mid-century (959 ± 362 kg ha −1 ) (Figure 7).Increasing maize landrace plant population resulted in a significant increase in mean yield but did not affect WUE (Figure 8).Regardless of bambara groundnut plant population, increasing maize landrace plant population resulted in a 12% reduction in its mean yield while reducing the population resulted in an 8% improvement in its mean yield (Figure 8).On the other hand, simulated yield and WUE for intercropped bambara groundnut was significantly (P < 0.05) affected by timescales and by the interaction between maize landrace and bambara groundnut planting date.
There was a reduction in simulated mean yield for bambara groundnut with the increase in maize landrace plant population.The trend for bambara groundnut yield was 841 ± 305 > 720 ± 379 > 657 ± 422 (kg ha −1 ) when intercropped with maize landrace at plant populations of 1.3, 2.6 and 3.9 plants m −2 , respectively (Figure 7).Similar to the simulated yield trend of intercropped bambara, increasing maize plant population resulted in a reduction of calculated bambara WUE and an increase in its variability (standard deviation) (Figure 7).The reduction of simulated yield and WUE maxima and minima for bambara groundnut and the increase in yield variability under high maize landrace plant populations could be attributed to increased competition for resources with the maize landrace.Peake et al. (2008) observed that increasing maize plant populations beyond a specific limit could increase the risk of crop failure due to an increase in competition for water and FIGURE 8 | Simulated yield response of maize landrace and bambara groundnut to plant population (plant m −2 ) for climate scenarios obtained from six GCM (ACC, CSS, CNR, GFD, NOR, and MPI).Different colors of boxplots represent the bambara groundnut plant density (plant m −2 ).The effect of the interaction between maize landrace plant density and bambara groundnut plant density on maize landrace grain yield -F-statistic: 62.47 on 8 and 891 DF, P = 0.000.The effect of the interaction between maize landrace plant density and bambara groundnut plant density on bambara groundnut grain yield -F-statistic: 38.93 on 24 and 875 DF, P = 0.000.FIGURE 9 | Calculated water use efficiency maize landrace and bambara groundnut plant population (plant m −2 ) for climate scenarios obtained from six GCM (ACC, CSS, CNR, GFD, NOR, and MPI).The colored boxplots represent the bambara groundnut plant density (plant m −2 ).The effect of maize landrace plant density on maize landrace WUE -F-statistic: 6.78 on 2 and 891 DF, P = 0.001.The effect of the interaction between maize landrace plant density and bambara groundnut plant density on bambara groundnut WUE -F-statistic: 38.93 on 24 and 875 DF, P = 0.000.solar radiation.In cases were both yields of the maize landrace and bambara groundnut are desired by a farmer, it might be worthwhile to reduce maize landrace plant populations to maximize yield for bambara groundnut.Alternatively, there is a need to improve water availability through rainwater harvesting and conservation techniques to reduce competition for water within the intercrop.

Impacts of Planting Dates on Yield and Water Use Efficiency of Maize Landrace and Bambara Groundnut Intercrop System
Simulation results for yield and WU for maize landrace and bambara groundnut across the six GCM were also not significantly different (P > 0.05); therefore, the results presented in this section were average values across the six GCMs.Simulated yield for maize landrace and bambara groundnut was significantly (P < 0.05) affected by the interaction of their planting dates (Figures 10, 11).Overall, early planting (September) of maize landrace or bambara groundnut resulted in higher simulated yields relative to late planting (January).Across the planting dates, mean yield trends for intercropped bambara groundnut was September (992 ± 296) > October (889 ± 357) > November (681 ± 383) > December (548 ± 301) >January (486 ± 283 kg ha −1 ).On the other hand, calculated WUE trend for intercropped bambara groundnut was September (2.81 ± 0.78) > October (2.53 ± 0.96) > November (2.46 ± 1.15) > January (1.99 ± 1.12) > December (1.94 ± 1.01 kg mm −1 ).For maize landrace, mean yield trends was September (1,052 ± 116) > October (902 ± 197) > November (835 ± 213) > December (648 ± 261) >January (596 ± 283 kg ha −1 ).On the other hand, calculated WUE trends for intercropped bambara groundnut was September (3.01 ± 0.88) > October (2.73 ± 0.76) > November (2.33 ± 0.95) > January (2.09 ± 1.00) > December (1.84 ± 0.66 kg mm −1 ).The calculated WUE was lower (24-107%) than the calculated WUE baseline of 4.12 kg mm 1 .This could be attributed to a reduction in the lower quartile values (Figure 10), which would suggest an increase in maize landrace yield gap with later planting.According to several research outputs, climate change is expected to reduce the length of the growing season and increase the occurrence of dry spells (Mitchell et al., 2015;Ajetomobi, 2016;Paff and Asseng, 2018).Despite the loss of growing days, our result suggests that, when planting on the same day, early planting (September) may ensure stable yields and WUE are obtained.Rezvani Moghaddam et al. (2014) found that early planting could be used as an adaptation strategy for maize under future climate in arid regions of Iran.Hussain et al. (2018) also highlight that, regardless of planting date, yield responses are highly dependent on resource availability and distribution, in this case, rainfall.
Intercropping bambara ground at different planting dates to maize landrace improved the mean yield (Figure 10), and WUE (Figure 11) provided it was done before November 15.Planting bambara groundnut a month earlier than maize landracefor instance, planting in the former September and maize landrace in October, resulted in a 176 and 57% increase in its mean yield and WUE, respectively, relative to the baselines.Planting bambara groundnut 2 and 3 months earlier than maize landrace resulted in a 184% increase in yield (Figure 10) and improved WUE by 61% increase in WUE (Figure 11).Planting maize landrace a month earlier than bambara groundnut -for instance, planting in September and bambara groundnut in October, resulted in the most significant mean yield increase (56%) relative to the baseline (Figure 11).The asynchronous or sequential planting did not result in the overlap in critical phenological stages for both the maize landrace and bambara groundnut.This minimized the competition for water and other resources and maximized resource use through extending canopy duration, therefore improving yield and WUE for maize landrace and bambara groundnut within the intercrop.When critical periods overlap, Yu et al. (2016)suggested that the competitive balance in cereal-legume intercrops can be maintained by planting the legumes earlier than the cereals.This can be viewed as a strategy to minimize the risk of yield loss in the event of intermittent dry spells within the season.However, sequential cropping in rainfed systems is constrained by the length of growing period (Inthavong et al., 2011;Kotir, 2011;Vadez et al., 2012;Duku et al., 2018;Minda et al., 2018).In this study, we did not assess the impacts of climate change on changes in the length and shifts of the growing season, nor the probability of dry spell occurrence and duration.

WAY FORWARD AND RECOMMENDATIONS
Overall, crop simulation models (CSM) and climate scenarios provided a monitoring and surveillance system to identify climate trends and associated impacts on intercropped maize landrace and bambara groundnut yield and WUE.In this regard, the use of a CSM driven by climate projections from six GCMs provided an opportunity to assess the suitability and sustainability of intercropping traditional crops as a potential climate adaptation strategy under low input-low output production systems.Our study demonstrated that the availability of a range of GCM outputs provided useful indications of the potential magnitude of yield and WUE changes and the temporal variation that could occur for the intercrop system.This type of analysis was, therefore, helpful in improving our understanding of the type of climate risk on the maize landrace and bambara groundnut intercrop system (Table 5).We recommend that the use of a CSM with GCM output should be considered when assessing the applicability of agricultural adaptation strategy.
Our results further showed that, at present, functional crop diversity could enhance crop productivity, stability, and thus food security, through efficient water utilization.Also, the adoption of asynchronous or sequential planting and moderating plant populations of either maize landrace or bambara groundnut can be viewed as a low-cost option to improve productivity and WUE under increasing temperature.This allows for the identification of short, medium and long-term strategies to aid in mitigating the impacts of climate change on the productivity and WUE of maize landrace and bambara groundnut intercrop system (Table 5).However, these approaches do not represent the diversity and breadth of adaptation strategies that can be adopted by marginal farmers.
To better represent adaptation, there is a need to expand the research to consider other management strategies (e.g., other traditional crop species, different cropping sequences, fertilization, rainwater harvesting and soil conservation techniques) (Seyoum et al., 2017).In addition, more system (agroecosystem) and place-based approaches that can represent local context, knowledge and aspects of food and nutrition security other than availability (e.g., nutrition, access, utilization and stability) may be required (Beveridge et al., 2018).To increase the contribution of agriculture to improving food and nutrition security, poverty reduction, and enhance rural economic development, climate impact modeling studies should be coupled with social, economic and environmental system models.This will ensure that traditional crops and associated cropping systems are assessed in a holistic manner that informs their sustainable integration into existing cropping systems.However, the adoption of traditional crops and intercropping should not be viewed as a panacea to solve all climate adaptation challenges, nor is it the only adaptation strategy.The inclusion of traditional crops into cropping systems should be considered as a complementary strategy to increasing climate resilience in marginal cropping systems.
FIGURE 10 | Simulated yield for maize landrace and bambara groundnut across different planting date combinations under rainfed conditions obtained from six GCMs (ACC, CSS, CNR, GFD, NOR, and MPI).The x-axis represents the maize landrace planting dates and the colored boxplots represent the planting date for bambara groundnut.The effect of the interaction between maize landrace planting date and bambara groundnut planting date on maize landrace grain yield -F-statistic: 49.93 on 24 and 875 DF, P = 0.000; The effect of the interaction between maize landrace planting date and bambara groundnut planting date on bambara groundnut grain yield -F-statistic: 75.37 on 24 and 875 DF, P = 0.000.FIGURE 11 | Calculated water use efficiency (kg ha −1 mm −1 ) for maize and bambara groundnut across different planting date combinations under rainfed conditions obtained from six GCMs (ACC, CSS, CNR, GFD, NOR, and MPI).The x-axis represents the maize landrace planting dates and the colored boxplots represent the planting date for bambara groundnut.The effect of the interaction between maize landrace planting date and bambara groundnut planting date on maize landrace WUE -F-statistic: 38.93 on 24 and 875 DF, P = 0.000 and The effect of the interaction between maize landrace planting date and bambara groundnut planting date on maize landrace WUE -F-statistic: 35.16 on 24 and 875 DF, P = 0.000.
A gap between the potential and practical realization of adaptation exists, and the evidence from our study supports the view that adaptation strategies need to be both climateinformed and context-specific to be viable (Beveridge et al., 2018;Carter et al., 2018).The cultivation of traditional crops has been done for millennia; however, to our knowledge, no study has quantified the yield and WUE responses in an intercrop system and under the impacts of climate change.Further to this, the FAO guidelines and key questions provided a useful framework to contextualize the observed results in an informative manner (Table 5) and less prescriptive.With the impacts of climate variability and change, our results provide evidence that adapting agronomic management could allow for sustainable intensification of the traditional systems through improved resource use efficiencies.However, we acknowledge that this type of study should be repeated across other agroecologies different from that of Ukulinga, allowing for more robust crop management practices and adaptation strategies to be identified.The calibration and validation process concerning the APSIM maize landrace and bambara groundnut intercrop study was the first attempt to evaluate the impacts of climate change on growth and water use.In this study, data to calibrate and validate the model were obtained from irrigated and rainfed experiments, respectively, in the same growing season.This may not fully meet strict requirements for using an independent data set for model validation.Therefore, future studies should repeat the experiment across various agro-ecologies and time scales different from that of the calibration data set; this will allow for better validation of model performance and robustness.

CONCLUSIONS
There is a high probability that yield and WUE for intercropped bambara groundnut will decrease in the near to distant future if current management options are maintained.Assuming future rainfall remains mostly unchanged, the primary limitations to intercropped bambara groundnut yield and WUE will be temperature and ET o under minimal rainfall changes.However, projected changes in temperature and ET o will increase yield and WUE variability for a maize landrace and bambara groundnut intercrop system.Improving WU, through increased plant density or asynchronous planting of the maize landrace and bambara groundnut mitigated the negative impacts of climate change on yield and WUE.In this regard, optimum plant management can optimize traditional production systems.Thus, intercrop systems of maize landrace and bambara groundnut should be promoted as a potential future system for climate change adaptation in rainfed production systems.
While the results of these simulations are limited to one agroecology and a single intercrop system, the findings confirm the views that several traditional crops are drought tolerant and thus, are suitable for cultivation in marginal agricultural production areas.Furthermore, intercropping them can increase system resilience under climate change.The concept of WUE, among other parameters, has been suggested in selecting management options that can sustainably increase productivity under climate changes, heat and water stress, and interactions among them.
Intercropping maize landraces and bambara groundnut with the appropriate place-based management practices can be used as an adaptation strategy in environments that are projected to face increasing water scarcity.Reduced land and water demand from intercropping maize landraces and bambara groundnut and improved water use efficiencies mitigate the risks associated with increasing climate variability and extreme events such as drought.For resource-poor farmers that are inherently risk-averse, the production of traditional crops such as maize landraces and bambara groundnut, and their optimisation through inexpensive management strategies present an opportunity to build resilient cropping systems.Our results have important implications on how traditional crops and cropping systems should be viewed, in that their incorporation into marginal production systems can be an alternative adaptation strategy that may lead to sustainable intensification outcomes under increasing climate risk.

FIGURE 1 |
FIGURE 1 | Comparison of observed and simulated (A) phenology (days after sowing) and (B) leaf area index for maize landrace and bambara groundnut during model calibration and validation.Red triangles represent bambara groundnut, and associated numbers 1, 2, and 3 represent phenological stages; emergence, the onset of flowering and start of grain filling respectively.Blue triangles represent maize landrace, and associated numbers 1, 2, 3, 4, and 5 represent phenological stages as emergence, floral initiation, flag leaf formation, the onset of flowering, and start of grain filling respectively.

FIGURE 2 |
FIGURE 2 | Distribution of average monthly minimum (A) and maximum (B) temperature data for the different timescales (past, present, mid-, and late-century) as simulated by the six GCMs (ACC, CSS, CNR, GFD, NOR, and MPI).

FIGURE 3 |
FIGURE 3 | Rainfall data representative of four different timescales (past, present, mid-, and late-century) as simulated by the six GCMs (ACC, CSS, CNR, GFD, NOR, and MPI).The average yearly rainfall calculated from observed rainfall data received between 2004 and 2019 was used as the mean annual rainfall.

TABLE 1 |
Soil water properties at different depths for soil at the experimental site.
a Bulk density; b Hygroscopic moisture content; c Permanent wilting point; d Field capacity; e Total available water; f Saturation; g Hydraulic conductivity.

TABLE 2 |
Modification of groundnut crop coefficients based on experimental data and data obtained from the literature.

TABLE 3 |
Global climate models used in this study.

TABLE 5 |
FAO guidelines and key questions for assessing the impacts of adaptation strategy.