Next Article in Journal
Winter Malting Barley Growth, Yield, and Quality following Leguminous Cover Crops in the Northeast United States
Previous Article in Journal
Effect of Organic and Inorganic Sources of Nitrogen on Growth, Yield, and Quality of Beetroot Varieties in Nepal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluating APSIM-and-DSSAT-CERES-Maize Models under Rainfed Conditions Using Zambian Rainfed Maize Cultivars

1
School of Natural Resources, Plant and Environmental Sciences, Copperbelt University, P. O. Box 21692 Kitwe, Zambia
2
School of Agricultural Sciences, Department of Soil Science, University of Zambia, P. O. Box 32379 Lusaka, Zambia
*
Author to whom correspondence should be addressed.
Submission received: 25 December 2020 / Revised: 11 September 2021 / Accepted: 13 September 2021 / Published: 23 September 2021

Abstract

:
Crop model calibration and validation is vital for establishing their credibility and ability in simulating crop growth and yield. A split–split plot design field experiment was carried out with sowing dates (SD1, SD2 and SD3); maize cultivars (ZMS606, PHB30G19 and PHB30B50) and nitrogen fertilizer rates (N1, N2 and N3) as the main plot, subplot and sub-subplot with three replicates, respectively. The experiment was carried out at Mount Makulu Central Research Station, Chilanga, Zambia in the 2016/2017 season. The study objective was to calibrate and validate APSIM-Maize and DSSAT-CERES-Maize models in simulating phenology, mLAI, soil water content, aboveground biomass and grain yield under rainfed and irrigated conditions. Days after planting to anthesis (APSIM-Maize, anthesis (DAP) RMSE = 1.91 days; DSSAT-CERES-Maize, anthesis (DAP) RMSE = 2.89 days) and maturity (APSIM-Maize, maturity (DAP) RMSE = 3.35 days; DSSAT-CERES-Maize, maturity (DAP) RMSE = 3.13 days) were adequately simulated, with RMSEn being <5%. The grain yield RMSE was 1.38 t ha−1 (APSIM-Maize) and 0.84 t ha−1 (DSSAT-CERES-Maize). The APSIM- and-DSSAT-CERES-Maize models accurately simulated the grain yield, grain number m−2, soil water content (soil layers 1–8, RMSEn ≤ 20%), biomass and grain yield, with RMSEn ≤ 30% under rainfed condition. Model validation showed acceptable performances under the irrigated condition. The models can be used in identifying management options provided climate and soil physiochemical properties are available.

1. Introduction

Maize (Zea mays L.) is the third important cereal crop after wheat and rice in the world [1]. It is grown both in the warmer temperate regions and humid subtropics [2]. In tropical Africa, nearly all maize grain is used for human food and prepared and consumed in many ways [1,2]. The generation of new research data using field experiments and the publication of such findings are not adequate to meet the increasing needs for the transference of new agro-technologies [3]. Furthermore, traditional fertilizer field trials’ findings are seasonal, site-specific, laborious and costly. The experimental data for one site may not be replicated at another due to spatial variations in the soil type and climate. In contrast, most maize agronomic recommendations are based on specific sets of field trials that are rarely repeatable due to environmental and seasonal variations [4].
The literature revealed that there is a need to transfer agrotechnology that enables the solving of agricultural problems through less time-consuming and expensive system-based research approaches. The use of crop simulation models (CSMs) is an efficient complement to agronomic field trails, and these are used to characterize, develop and assess crop production systems [5]. The use of system-based experiments to transfer agro-technology to farmers, policymakers and planners is feasible by using decision support tools (DSTs), such as the Agricultural Production Systems Simulator (APSIM) [6,7] and Decision Support Systems for Agrotechnology Transfer (DSSAT) [3,8,9]. Policymakers and land managers can use crop models to identify management options, provided climatic data, soil physiochemical properties, socioeconomic information and management are available [10]. The interactions between weather, management and soil physiochemical properties that affect crop growth and yield can be assessed using crop models [11]. APSIM and DSSAT allow users to conduct onscreen experiments in minutes on a workstation, server, laptop or desktop computer [12].
The APSIM has been used to simulate wheat evapotranspiration, yield, nitrogen uptake and soil evaporation [13]. APSIM has been evaluated in Asia under different environments, crops and management scenarios. The model performed well in simulating different cropping systems and multi-crop sequences [14]. However, the model was unable to accurately simulate aspects related to high temperatures, salinity, conservation agriculture and soil micronutrients deficiencies. APSIM has been used in predicting soil water deficits and drainage in New Zealand [15]. Additionally, the APSIM model more sensitive to management practices and environmental conditions.
The APSIM-Maize and CERES-Maize models were assessed [16]. The APSIM and DSSAT models are different in the way they model incoming solar radiation and photosynthetic active radiation (PAR). APSIM uses the Beer’s law approach [17], while DSSAT uses the approach by [18]. The evaporation algorithm developed [19] is used in APSIM-Maize and CERES-Maize models. The APSIM-maize and DSSAT-CERES-Maize both require many inputs, with many model parameters not readily available at the farm level. APSIM parameterization and calibration requires obtaining better growth and development descriptions of the low-yielding cultivars found in African farming systems [20]. APSIM and DSSAT-CERES-Maize models have been used in long-term cropping systems trials, on-farm decision-making, farming systems design for production, seasonal climate forecasting, risk assessment and education in tropical and semi-arid areas [6,21,22].
The DSSAT-CERES-Maize model has been evaluated under a wide range of field practices and environmental conditions [3,8,23]. It has also been evaluated in Zimbabwe, Malawi, Kenya, South Africa and Zambia [2,4,24]. The ZMS606 had previously been calibrated for the DSSAT-CERES-Maize model under rainfed conditions in the 2013/2014 season [2]. ZMS606 is a medium (125–130 days) maturing cultivar commonly grown by small-scale farmers throughout Zambia. The DSSAT-CERES-Maize model accurately predicts the crop response to nitrogen, nitrogen uptake and crop yield variability, and it has been used to explore the potential of new cultivars for new areas and optimal sowing dates without establishing costly field trials [25].
The use of APSIM-Maize and DSSAT-CERES-Maize models is limited in Sub-Saharan Africa due to the scarcity of suitable model input data for model parameterization, calibration, testing, validation and application. The parameterization and calibration of the APSIM-Maize and DSSAT-CERES-Maize models is necessary for new cultivars, environments and subsequent model evaluations to ensure the credibility of the model performances [26]. The performances of the APSIM and DSSAT-CERES-Maize models have been tested in different locations, environments and under long-term cropping systems and climate change [6,9]. Although APSIM plant and soil modules have been extensively calibrated and evaluated under Australian and Eastern and Southern African conditions [6,22], they have not been comprehensively calibrated and tested under Zambian conditions. Therefore, the study objectives were: (i) to determine the CSPs of three maize cultivars for parameterization of the APSIM-Maize and DSSAT-CERES-Maize models under rainfed conditions; (ii) to evaluate the APSIM-Maize and DSSAT-CERES-Maize models in predicting the grain, biomass yield, leaf area index and soil water content by maize under three sowing dates (SDs) and nitrogen fertilizer rates under rainfed conditions and (iii) to validate the APSIM-Maize and DSSAT-CERES-Maize models using N fertilizer levels under irrigated conditions.

2. Materials and Methods

2.1. Study Site and Soil Data

The study was at the Mount Makulu Central Research Station in Chilanga, Zambia (latitude: 15° 550′ S, longitude: 28° 250′ E, altitude: 1213 m). The daily weather data (solar radiation, maximum and minimum temperature and rainfall) were collected from the automatic meteorological stations for the Zambia Meteorological Department (ZMD) near the experimental site. The crop models required accurate solar radiation to drive their simulation of photosynthesis and the carbon balances that govern plant growth. The soil type at Mount Makulu is classified in Soil Taxonomy as Ustic Paleustalf [27].

2.2. Weather Data under Rainfed Condition

The daily precipitation, solar radiation, maximum and minimum temperature data are presented in Figure 1 for the 2016/2017 season. The total precipitation, mean, maximum and minimum temperature during the field experiment period were 930.17 mm, 21.83 °C, 28.29 °C and 15.36 °C, respectively. As the season progressed, the precipitation, solar radiation, maximum and minimum temperatures reduced. SD1 (precip: 850.37 mm, tmax: 27.47 °C, tmin: 17.14 °C, tmean: 21.72 °C and Srad: 17.38 MJ m−2 day−1) recorded higher meteorological parameters compared to SD2 (precip: 763.27 mm, tmax: 27.02 °C, tmin: 16.32 °C, tmean: 21.67 °C and Srad: 17.22 MJ m−2 day−1) and SD3 (precip: 515.27 mm, tmax: 26.88 °C, tmin: 17.16 °C, tmean: 21.24 °C and Srad: 17.38 MJ m−2 day−1). All the three SDs experienced mean maximum temperatures below 38 °C during the growing period from emergence to physiological maturity. The critical temperature that can have detrimental effects on the maize yield is 32 °C, as reported [28]. Maize during its growing period requires about 450–600 mm of water season−1, which is taken from the soil moisture reserves. SD3 received little water during the season, and this most likely affected the plant growth and yield.

2.3. Weather Data under Irrigation Condition

The weather data (rainfall, maximum, minimum temperature and solar radiation) used in the validation was obtained from the Zambia Meteorological Department (ZMD), as shown in Figure 2. The weather data was taken from May 2016 to November 2016. The Tmin, Tmax, Tmean and precipitation were 13.76 °C, 28.64 °C, 21.21 °C and 71.70 mm for the whole period, respectively. The mean temperature ranged from 13.20 °C to 31.05 °C under the irrigated conditions. It has been observed that maize grows well in between mean temperatures of 15 and 35 °C [29].

2.4. Experimental Design, Treatments, Soil Physical and Chemical Analysis

A rainfed split–split plot design was set up at the Mount Makulu Central Research Station in Chilanga, Zambia with three sowing dates (SDs); three maize cultivars (ZMS 606, PHB 30G19 and PHB 30B50); three nitrogen (67.20, 134.40 and 201.60 kg N ha−1) fertilizer rates and three replicates, respectively. The main plot was the SD (3 sowing dates [SD]) at 14-day intervals, while the subplots and sub-subplots were the maize cultivars and N rates. Individual plot sizes were 6 m (7 rows) by 5 m (30 m2). The plots were separated from each other by a 2-m distance to prevent the cross-contamination of treatments. Three seeds were sown by hand at a 5-cm depth in a flat seedbed in 0.75-m row spacing and 0.50-m spacing between plants per station and later thinned to two plants. The initial soil conditions were sampled using a soil auger at 20-cm intervals until a 100-cm depth, weighed and oven-dried at 105 °C.
The soil parameters presented in Table 1 were analyzed at the Mount Makulu Central Research Station in Chilanga, Zambia. The soil pH in calcium chloride ranged from 6.2 to 6.80, considered to be within the neutral to optimal range for crop growth. The pH was considered to be favorable for availability of most plant nutrients. The total nitrogen percentage values ranged from 0.031 to 0.061 and, therefore, were considered to be very low. The soil organic carbon percent (OC) was analyzed using the Walkley and Black method [30], and the values ranged from 0.35 to 0.80 and were considered to be low. The critical OC value is 1.58%. The phosphorus (P) and potassium (K) ranged from 10 to 18 mg kg−1 and 0.59 to 1.12 mg kg−1, respectively. The P ranged from very low to low, while the K was very low. The soil texture and organic matter were used as inputs into the Soil–Plant–Air–Water (SPAW) pedo-transfer functions [31,32] to determine the lower limit (LL), drained upper limit (DUL), saturation (SAT), hydraulic conductivity and bulk density (Table 1). Higher organic matter increased the water-holding capacity and conductivity, as a result of pore space distribution and soil aggregation [32]. In addition, the organic matter and soil texture are two variables that affect the soil water content.

2.5. Soil Water Content Measurement

The Diviner 2000 series II probe (Senteck Sensor Technologies, Stepney, SA, Australia; https://sentektechnologies.com/product-range/soil-data-probes/diviner-2000/) is a portable and robust device measuring soil water at multiple depths (at 10-cm intervals) in the profile [34]. The Diviner 2000 series II probe consists of a probe and handheld data logging display unit. It measures the relative changes in the volumetric soil moisture content, which is used to show the most important soil water trends with time. Access tubes were installed in each experimental plot and soil water content measurement, taken two to three times per week using a Diviner 2000 series II probe. The series of readings taken by the Diviner 2000 series II assisted in showing trends of crop water use in the soil profile. The volumetric soil water content collected using the Diviner 200 series II probe were used as inputs and simulated with the APSIM-Maize and DSSAT-CERES-Maize models.

2.6. Plant Materials

The maize cultivars PHB 30G19, PHB 30B50 and ZMS 606 have medium maturity with the comparative relative maturity of 120–130 days. These cultivars were selected as the major cultivars planted by small-scale farmers and their long commercial lives. PHB30B50 is recommended to be grown under irrigation, while PHB30G19 and ZMS 606 can be grown under irrigated and rainfed conditions in Zambia. PHB30G19 (white) and PHB30B50 (yellow) are produced by DuPont Pioneer (Lusaka, Zambia). ZMS 606 is an exceptionally good drought-tolerant maize cultivar produced by ZamSeed (Lusaka, Zambia). The selected cultivars can be grown in all the three agro-ecological regions of Zambia.

2.7. Plant Growth Analysis

A plant growth analysis was observed at the vegetative (emergence, V6) and reproductive (silking (R1), dough stage (R4) and physiological maturity (R6)) stages and recorded when 50% and 75% of the plants reached the stages, respectively, as described in [35,36] (Table 2 and Table A1). The biomass from all the subplots was harvested at the recommended growth stages. The dry weight of each fraction of biomass harvested at each vegetative and reproductive stage was recorded before and after drying the samples at 70 °C for three days. Four plants were sampled from each subplot at V6, R1 and R4 for the biomass and leaf area measurements. The maize leaf area measurement was calculated by multiplying the manually measured length and maximum width by 0.75, reported as the maize calibration factor [37]. Accurate prediction of the leaf area measurement is critical for the precise forecasting of crop growth and yield by most crop simulation models [9,36]. The final maize aboveground biomass was harvested using standard procedures [36].

2.8. Modelling Approach

2.8.1. Description of the APSIM-Maize Module

APSIM version 7.9 [21] used in this study includes the APSIM-Maize, SoilWater, SoilN and SurfaceOM modules. It is a process-based model that runs at a daily time step and combines biophysical and management modules within a central engine to simulate the crop growth, development, soil water, phosphorus, carbon and nitrogen dynamics using the current and future climate scenarios [22,26,39]. It contains more than 80 modules that simulate plant, animal, soil, climate and management scenarios. It can simulate more than 20 crops and forests, such as alfalfa, eucalyptus, cowpea, pigeon pea, peanuts, cotton, lupin, maize, wheat, barley, sunflower, sugarcane, chickpea and tomato [7].
The APSIM-Maize module [40] consists of a phenology model, structure model and an arbitrator for nitrogen and biomass partitioning to various plant organs in the model. It simulates maize growth, leaf area and development on a daily time step. The growth of maize responds to precipitation, temperature, radiation, soil nitrogen (SoilN) and soil water supply (SoilWat) [41,42]. The phenology model simulates the growth stages, and each stage is divided into subphases. The duration of each phase is regulated by the temperature and photoperiod and may be affected by water and nitrogen (N) stresses.

2.8.2. Description of the DSSAT-CERES-Maize Model

The Decision Support System for Agrotechnology Transfer (DSSAT) v4.7 comprises 42 crop simulation models (CSMs), and this model is the CERES-Maize [36,43,44]. The DSSAT-CERES (Crop Environmental Resource Synthesis) maize model [45,46] is a process-oriented, management level model that simulates soil water and nitrogen dynamics, crop growth, yield and management scenarios at the field level on a daily time step from inputs of climate data, genotype, soil and management [3,45,47]. The soil water balance in DSSAT-CERES-Maize is based on Ritchie’s model, where the concept of a drained upper limit (DUL) and drained lower limit (LL) of the soil is used as the basis of the available soil water [48,49]. It is a widely used maize model and a reference for comparing new advances in maize simulations by farmers, researchers and planners in many countries [9].

2.8.3. Input Dataset into the APSIM-Maize and DSSAT-CERES-Maize Models

The minimum input data required by the APSIM-Maize and DSSAT-CERES-Maize models are: the minimum input data required to run APSIM-Maize and DSSAT CERES-Maize models includes daily records of rainfall, solar radiation, maximum and minimum air temperature, soil characterization data (physical and chemical properties and morphological properties for each layer), soil drained upper limit, soil lower limit, soil saturation and genetic cultivar coefficients characterizing the cultivar being grown (plant development and grain biomass and crop management information). The soil water balance is simulated to evaluate the potential yield reduction caused by soil water deficits. The data presented in Table 1 and Table 3 were used as inputs into the crop simulation models.

2.9. Parameterization, Calibration and Validation of the APSIM-Maize and DSSAT-CERES-Maize Models

2.9.1. Parameterization and Calibration of the APSIM-Maize and DSSAT-CERES-Maize Models

The soil analysis and aboveground biomass data (Table 1 and Table 3) collected from the 2016/2017 rainfed field experiment was used as inputs to calibrate the APSIM-Maize and DSSAT-CERES-Maize models. In APSIM, U and CONA parameters similar to the DSSAT-CERES-Maize model are used to determine the first and second evaporation stages. The U and CONA were set at 6 mm and 3.5 mm day−1, respectively and these parameters are crop-specific and determine the rate of root extension and lower limit of crop water extraction. The Air-dry (mm mm−1) soil moisture limit at which soil dries by evaporation was calculated as 0.5 × LL15 (015 cm soil layer), 0.9 × LL15 (1530 cm soil layer) and at deeper depths same as LL15 [50]. The crop lower limit (CLL in mm mm−1) was taken as equal to LL15. The soil water characteristics used as model inputs are shown in Table 1 above.
The growth and development module of APSIM-Maize uses a set of different coefficients as presented in Table 3 to define the phenology, crop growth and yield [51,52]. In the APSIM-Maize model, the B_130 cultivar was used as the starting point for calibrating the ZMS606, PHB 30G19 and PHB 30B50 maize cultivars. The Cultivar Specific Parameters (CSPs) for the three maize cultivars were obtained step-by-step based on the observed phenological stages (days after planting (DAP) to anthesis and maturity; leaf area index (LAI); biomass and grain yield (grain growth rate, grain yield, grain unit wt. and grain number)) using the manual trial and error method until the simulated values were within 9–20% of the observed values [53].
The CSPs in the DSSAT-CERES-Maize model were calibrated automatically using the GLUE (Generalized Likelihood Uncertainty Estimation) program embedded within DSSAT [43]. The CSPs describe the phenology and grain yield. Six parameters (P1, P2, P5, G2, G3 and PHINT), as shown in Table 2, are used by the DSSAT-CERES-Maize model [3,23]. The final leaf number in DSSAT-CERES-Maize is calibrated using the PHINT CSP. Using the soil and weather data, calibration procedure was done for the phenology stages, followed by grain development. The “P” CSPs (P1, P2 and P5) and “G” CSPs (G2 and G3) were computed using the phenological dates and grain yield, respectively. The derived CSPs for the cultivars used in the field experiment were added to the DSSAT’s genotype file to be used in the simulations. They control the growth and development of the crop and, by affecting the root growth and root water uptake, indirectly impact the soil water content [54]. The APSIM-Maize and DSSAT-CERES-Maize models were calibrated using phenology, physiological maturity, maximum leaf area index (mLAI), grain yield, aboveground biomass yield, grain unit wt., grain number m−2 and root soil water content to the depth of one meter.

2.9.2. Validation of the APSIM-Maize and DSSAT-CERES-Maize Models

The APSIM-Maize and DSSAT-CERES-Maize models were validated using the phenology, physiological maturity, mLAI, grain yield, aboveground biomass yield, grain unit wt. and grain number m−2. Cultivar specific parameters (CSPs) derived from the calibration of the two crop models were used to evaluate the robustness of the APSIM-Maize and DSSAT-CERES-Maize models in a subtropical environment of Zambia. The observed yield and yield components from the irrigated field experiment in 2016 were used to validate the two crop models. The validation of the APSIM-Maize and DSSAT-CERES-Maize models was essential for model application.

2.10. Data Analysis and Model Evaluation

Analysis of Variance

The treatment effects of the sowing date (SD) on the maize cultivar and N rate on the maize yield and yield components collected from the rainfed field experiment was computed using the analysis of variance for the split–split plot design. The analysis of variance was calculated using the ssp.plot function in the R programming software Agricolae package and means separated at p ≤ 5 using Tukey’s HSD test [55]. On the other hand, the treatment effect of the nitrogen fertilizer rate and cultivar on the maize yield and yield components observed at the irrigated experimental site was evaluated by a split–plot analysis of variance using the spp.plot function in the Agricolae package [55] in R Programming software.

2.11. Evaluation of the APSIM-Maize and DSSAT-CERES-Maize Models

The root mean square error (RMSE), normalized (RMSEn) root mean square error, d-stat (d) and 1:1 graphs (observed vs. simulated graphs) [45,56,57] were used to evaluate the crop models. In an adequately calibrated crop model, the EF should be above 0.70 [41]. The d-stat (d) ranges from 0 (the simulated is similar to the observed mean) to 1 (perfect model performance) [58]. The simulation is considered excellent with RMSEn <10%, good if 10–20%, acceptable or fair if 20–30% and poor >30% [59]. The d ≥ 0.70 is considered acceptable in crop modeling and evaluations [60].
RMSE = [ N 1 i = 1 n ( P i O i ) 2 ] 0.5 N 1
RMSEn = RMSE O ¯ 100 %
d stat = 1 P i O i 2 P i O i 2
EF = 1 ( P i O i ) 2 ( O i O ¯ ) 2
where: P i , O i ,   O ¯ and n are the simulated, observed, average observed values and the number of observations for the studied variables, respectively.

3. Results and Discussions

3.1. Analysis of Variance for the Sowing Dates, Cultivars, N Rate and Yield Components

3.1.1. Rainfed Condition

The treatment effect of the sowing date (SD), cultivar and nitrogen fertilizer application rates on the yield and yield components of maize are shown in Table 3 and Table 4. The treatment effect of SD and the cultivar were significant on the 100-grain wt., stover wt., cob wt., stem, wt., leaf blade wt., cob width and harvest index, as reported by [61,62]. The differences in the 100-grain wt. may have been caused by differences in the initial sizes of the spikelets, in growth rates during the exponential and linear phases of the grain-filling period [61]. The soil moisture stress after silking notably decreased the pooled values for the 100-seed weight (g) with an increasing N rate. The grain yield, aboveground biomass, grain number m−2, cob wt., cob width, harvest index and 100-grain wt. were reduced with a delay of 7 and 14 days in sowing, low soil fertility, lowering of the air temperature and reduction in cumulative rainfall during the duration of plant growth, as shown in Table 3 and Table 4. The cultivar treatment effect on the grain number m−2, 100-grain wt., stover wt., cob wt., husk wt., stem wt., leaf blade wt., harvest index and cob width was statistically significant. The plant LAI at any growth stage indicates the size of the assimilatory system, which contributes to dry matter accumulation and partitioning [63]. There was an interaction effect between the SD and cultivar on 100-grain wt. Furthermore, there was an interaction effect between the cultivar and nitrogen application rates on the stem wt. The nitrogen treatment effect on the yield and yield components was statistically nonsignificant. The application of N fertilizer did not significantly improve the maize yield [64]. However, deficiency or excess nitrogen application can reduce the maize yield [65].
The reduction in grain yield from SD1 to SD2 (1.91 t ha−1), SD1-SD3 (2.90 t ha−1) and SD2 to SD3 (0.99 t ha−1) were 21.04%, 31.83% and 13.83%, respectively. A reduction in the maize yield with a delay in SD was observed in Northern Nigeria [66]. Delayed SD reduces the duration of the grain filling process and yield [24]. The SD is a critical factor in capturing sufficient solar radiation with adequate soil moisture and fertility. The sowing date considerably affected the maize grain and biomass yield [2].

3.1.2. Irrigated Condition

The treatment effect of the cultivar and nitrogen fertilizer application rates on the yield and yield components of maize are shown in Table 5. The cultivar treatment effect under the irrigated condition significantly affected the 100-grain wt. and husk wt., as examined other authors [62]. The cultivar treatment effect on the 100-grain wt. and husk wt. was statistically significant. The treatment effect of the nitrogen fertilizer rate increased the grain wt., biomass wt., cob wt. and grain number m−2. The grain number m-2 increased with the higher N fertilizer application rates [67]. The highest husk yield was observed from ZMS 606 (383.94 kg ha−1), followed by PHB 30G19 (301.56 kg ha−1) and PHB 30B50 (267.27 kg ha−1). Pooled data for the husk yield showed that N1 (375.48 kg ha−1) had the highest, followed by N2 (289.70 kg ha−1) and N3 (287.59 kg ha−1), respectively. The grain number cob−1 and grain number m−2 were highly significantly affected by the N fertilizer treatment effect. The grain number m−2 increased with the higher N fertilizer rate. The maximum and a minimum numbers of grains per ear were recorded at N3 and N2, respectively. ZMS 606 recorded the highest mean grain number of 349.8, followed by PHB 30G19 (292.6) and PHB 30B50 (281.9). The highest grain number per ear was observed at the highest nitrogen application rate [68]. Similar results were obtained for maize yield response to the N fertilizer application [64]. However, nitrogen fertilizer application rate had no effect on the HI [69].
Stress as a result of water deficiency from the silking to maturity stages affected the ultimate sizes and yields of ears, and similar results were reported [28]. Adverse conditions such as water stress and nitrogen deficiency delay plant growth and slow silk development [70,71]. Nitrogen is a yield-limiting nutrient, and its application plays a significant role in improving the soil fertility [71].

3.2. Performance of APSIM-Maize and DSSAT-CERES-Maize Models in Simulating Growth and Yield for Three Maize Cultivars

3.2.1. Phenology

The comparison between the simulated days after planting (DAP) and observed days after planting (DAP) to anthesis, physiological maturity (DAP), maximum leaf area index (mLAI), grain unit wt. at maturity, grain number m−2, biomass and grain yield statistics are presented in Figure 3a–l and Table 6, Table 7, Table A2 and Table A3 for the APSIM-Maize and DSSAT-CERES-Maize models. The percent prediction (%PD) values for anthesis and maturity at SD1 (DSSAT-CERES-Maize model) were close to or equal to zero, and this indicted excellent agreement between the observed and simulated values. The %PD for anthesis and maturity was <20% of the observed values, and similar results were reported by [53]. The DSSAT-CERES-Maize model simulated maize cultivar anthesis, and physiological maturity was 62–68 and 122–134 DAP, respectively during the 2016/2017 season. The APSIM-Maize model simulated maize cultivar anthesis, and maturity was from 64 to 67 and 123 to 130 DAP, respectively. However, the observed maturity for the maize cultivars ranged from 120 to 130 DAP. Another study indicated that delayed sowing reduced the days after planting to anthesis and the number of leaves [24].
The results showed that there was good agreement between the observed and simulated anthesis (d = 0.47–0.49) and maturity (d = 0.71–0.80) for the cultivars using both models. The modeling of anthesis and maturity for the cultivars (ZMS 606, PHB 30G19 and PHB 30B50) was excellent (RMSEn ≤ 10). In the DSSAT-CERES-Maize model, the RMSE between the observed and simulated anthesis and maturity were 2.89 (d-stat = 0.47) and 3.13 (d-stat = 0.80) days, respectively. The anthesis and maturity RMSEn were 4.36 and 2.50 and considered excellent across the three cultivars. In APSIM-Maize, the RMSE for the simulation of DAP to anthesis and maturity were 1.91d (RMSEn = 2.89%, d-stat = 0.49) and 3.35d (RMSEn = 2.68%, d-stat = 0.71) days across the cultivars, respectively.
A multi-model ensemble of the APSIM, DSSAT-CERES and Aquacrop models have been used to simulate maize growth and yield under rainfed condition in Ethiopia [72]. The APSIM-Maize and DSSAT-CERES-Maize models simulated accurately days after planting to anthesis (RMSE = 1.73–4.09) and maturity (RMSE = 1.66–5.36). It was observed by [72] that the DSSAT-CERES-Maize model overestimated DAP to maturity for late-maturing cultivars. The results reported are comparable to this study (Table 6 and Table 7) [72]. APSIM-maize has been used to simulate the optimal sowing windows for maize [73]. The model simulated accurately anthesis (DAP) and maturity (DAP), grain yield and biomass for maize cultivars.

3.2.2. Biomass and Grain Yields

The grain yield results showed a wide range of %PD for simulated from observed values (−25.93 to +45.28% for APSIM-Maize and −21.74 to +49.06% for DSSAT-CERES-Maize). This indicated that grain yield was not well-simulated accurately for the cultivars. The evaluation statistics for the biomass yield, grain yield, grain unit wt. and grain number m−2 are shown in Table 8 and Table 9. The overall RMSE for grain (RMSE = 0.84 t ha−1, dstat = 0.92, R2 = 0.72) and biomass (RMSE = 2.87 t ha−1, d-stat = 0.89, R2 = 0.86) yield in the DSSAT-CERES-Maize model were 11.28% and 64.84% and considered as good and poor, respectively. The simulation of grain (RMSE = 1.38 t ha−1, d-stat = 0.57, R2 = 0.26) and biomass (RMSE = 3.19 t ha−1, d-stat = 0.90, R2 = 0.85) yield in APSIM-Maize was good and poor, with RMSEn of 18.42% and 71.64%, respectively. However, the simulated biomass in APSIM-Maize was acceptable, with a d-stat ≥ 0.70. In contrast, the grain and biomass yield evaluations (DSSAT-CERES-Maize model) for all three cultivars had a d-stat ≥ 0.70. The d-stat ≥ 0.70 has been reported as being acceptable in crop model evaluations [41,60]. The RMSE for a pooled biomass yield in the APSIM-Maize and DSSAT-CERES-Maize models were 3.19 t ha−1 and 2.87 t ha−1, respectively. The cultivar RMSE for the grain and biomass (Table 8 and Table 9) yield in DSSAT-CERES-Maize ranged from 0.68 to 1.11 t ha−1 and 2.65 to 3.03 t ha−1, respectively. The evaluation of the DSSAT-CERES-Maize model showed that the model was able to simulate the grain yield for the three cultivars accurately, as shown in Table 8 and Table 9. The modeling of the grain yield for ZMS 606 and PHB 30G19 and PHB 30B50 were considered excellent and good, respectively. The modeling of the biomass yield was considered good (PHB 30B50, RMSEn (10–20%)) and acceptable (ZMS 606 and PHB 30G19, RMSEn (20–30%)).

3.2.3. Grain Unit wt. and Grain Number per Square Meter

The simulation of the grain unit wt. and grain number m−2 in the DSSAT-CERES-Maize model was acceptable, with RMSE of 0.07 g and 379.81 grains m−2, respectively, as shown in Table 10 and Table 11. The underprediction of the grain number m−2 by the DSSAT-CERES-Maize model was also reported by [2,61]. In APSIM-Maize, the RMSE for the grain unit wt. and grain number m−2 were 0.05 g and 276.13 grains m−2, respectively. The simulation of the grain unit wt. and grain number m−2 in APSIM-Maize was acceptable, with a d-stat ≥ 0.70 and RMSEn of 10–20%. The sowing date analysis showed that delaying the sowing date from 12 December 2016 to 26 December 2016 and 9 January 2017 caused a decrease in the average yield for all cultivars. Other researchers reported that APSIM performs better in simulating the soil water, grain yield and nitrogen dynamics in different agricultural production systems in Australia [52]. In this study, DSSAT-CERES-Maize performed better in the simulation grain yield, grain unit wt. and grain number m−2. The DSSAT-CERES-Maize model explained 63%, 86% and 72% of the variation in the observed grain unit wt., biomass and grain yield, respectively. In APSIM-Maize, 85% of the biomass yield was explained by the model.

3.2.4. Leaf Area Index

The RMSE for the mLAI in the APSIM-Maize and DSSAT-CERES-Maize models were 0.69 and 2.04 m2 m−2, respectively. The simulations of the mLAI in APSIM-Maize and DSSAT-CERES-Maize models were good and poor, respectively. The maximum LAI in the models was underpredicted at all the treatment levels, as shown in Table 12. Similar results in the DSSAT-CERES-Maize model were reported [2]. The APSIM-Maize and DSSAT-CERES-Maize models simulate LAI, but in this situation, the functions performed poorly. The model failed to accurately simulate LAI for the three maize cultivars. LAI provides an index of the maize plant growth and is an important input into the APSIM-Maize and DSSAT-CERES-Maize models [36]. Accurate field or laboratory measurements of LAI are required to evaluate crop models properly. The simulated mLAI ranged from 3.02 to 3.35 m2 m−2 (APSIM-Maize) and 1.44 to 2.18 m2 m−2 (DSSAT-CERES-Maize), respectively. Typical LAI values of maize cultivars in dry land areas under rainfed production are 2.50–2.90 m2 m−2 [74] and agree with the results from the APSIM-Maize model.

3.2.5. Simulation of Root Soil Water Content in the Soil Layers

The statistics for pooled root soil water content under irrigated conditions are shown in Table 13 for the APSIM-Maize and DSSAT-CERES-Maize models. The predictive ability of the DSSAT-CERES-Maize model was excellent (RMSEn >10%, soil layer: 4), good (RMSEn 10–20%, soil layers: 2 to 3 and 5–7) and acceptable (RMSEn 20–30%, soil layers: 1 and 8). In the APSIM-Maize model, the predictive ability of the model was good (RMSEn 10–20%, soil layers: 1 and 6–8) and excellent (RMSEn >10%, soil layers 2–5). The modeling of the root soil water content may have been affected by unmeasured soil parameters, such as management scenarios (weeds, disease and pest control); LL; DUL and SAT. The APSIM soilwat module requires the determination of LL and DUL under field conditions [41]. The root soil water is important and should be modeled accurately, as it affects the crop growth and prediction of the grain and biomass yield [75]. Changes in the available root soil water content lead to variations in the biomass and grain yield for maize cultivars. Maize crops can only use root soil water that is present within the reach of their root system, and this could significantly influence the total available root soil water.
The analysis of SDs showed that the amount of root soil water available for plant growth and rainfall decreased during each maize growing period, and this affected the biomass and grain yield. Other researchers reported that root soil water profile layers of 5–15, 15–30 and 30–45-cm thickness are important for simulating the correct root soil water balance and plant water uptake [76]. Water stress in maize increased as the rainfall amount decreased during each SD duration, and this contributed to the grain yield reduction, similar to other findings [77].

3.3. APSIM-Maize and DSSAT-CERES-Maize Model Validation

3.3.1. Model Validation

Phenology (Anthesis and Maturity Days after Planting)

The observed and simulated phenological stages, grain unit wt., grain number m−2, grain and biomass yield. A number of validation statistics have been used to evaluate the models in examining their predictive ability, usefulness and robustness, as shown in Table 14 and Table 15. In the APSIM-Maize model, the mLAI, DAP to anthesis and maturity ranged from 3.26 to 3.37 m2 m−2, 87 to 88 days and 135 to 137 days under irrigated conditions, respectively. The DAP to anthesis and maturity were good (RMSEn = 10–20%), with RMSE being 16.69 and 17.36 days, respectively as reported by [59]. The mLAI, DAP to anthesis and maturity of the maize cultivars in the DSSAT-CERES-Maize model ranged from 1.61 to 2.11 m2 m−2, 80 to 86 days and 130 to 137 days, respectively. The RMSE for DAP to anthesis (RMSEn = 19.94%) and maturity (RMSEn = 12.77%) were 20.67 and 19.54 days. The vegetative and reproductive phases of the maize cultivars were delayed due to lowering of the temperature under irrigated conditions. Additionally, lowering temperatures during the growing season reduced the plant height through a decreased internode length, and similar results were reported [78].
The APSIM-Maize and DSSAT-CERES-Maize models simulate DAP to anthesis and maturity, but in this situation, the function performed poorly. The phenological stages (anthesis and maturity) took longer due to lower solar radiation and cooler air temperature during the field experiment. Plant stress during its growth period from silking to maturity affect the grain yield and yield components [28]. Adverse conditions such as lower temperature, soil water content and nitrogen deficiency delay plant growth and slows silk development and maturity [70,71].
Crop simulation models require accurate field observations for the timing of vegetative (V6) and reproductive stages (silking, dough and physiological maturity stages). This is important in crop simulation, so that partitioning of dry matter accumulation and growth duration is determined properly. The staging of phenological stages is vital in field experiments, since DAP to anthesis and physiological maturity are important determinants of productivity [43].

3.3.2. Biomass, Grain Yield and Leaf Area Index

The calibrated model was validated against the end-of-season mLAI, grain unit wt., grain number m−2, biomass and grain yield. The validation of the APSIM-Maize model was the grain unit wt. (RMSE = 0.02 g, RMSEn = 6.25% and d-stat = 0.75); grain number m−2 (RMSE = 456.21 grain m−2, RMSEn = 18.51% and d-stat = 0.59); biomass (RMSE = 6.87 t ha−1 and d-stat = 0.23) and grain yield (RMSE = 1.34 t ha−1, RMSEn = 20.18% and d-stat = 0.48). The simulation of the grain yield, grain unit wt. grain number m−2 and mLAI showed the robustness of the APSIM-Maize model considered acceptable in crop modeling [59].
Using the DSSAT-CERES-Maize model, the RMSE for grain unit wt. (RMSEn = 21.14%, R2 = 0.72 and d-stat = 0.55); grain number m−2 (RMSEn = 26.91% and d-stat = 0.52); grain (RMSEn = 45.08% and d-stat = 0.35) and biomass (RMSEn = 71.64% and d-stat = 0.24) yield were 0.06 g, 663.29 grains m−2, 2.99 t ha−1 and 6.65 t ha−1, respectively. The simulation of the grain unit wt. and grain number m−2 were evaluated to be good. However, the simulation of the grain (RMSEn = 45.08%) and biomass (RMSEn = 71.64%) yield were considered to be poor in crop modeling, as reported by [59]. The RMSE for mLAI in APSIM-Maize (RMSEn = 16.75%) and DSSAT-CERES-Maize (RMSEn = 49.56%) were 0.57 m2 m−2 and 1.67 m2 m−2 and considered acceptable and poor in the crop model evaluation, respectively. The accuracy of the simulated outcomes of the crop models using new cultivars depend on rigorous calibrations to minimize uncertainties.

4. Conclusions

The results showed that the maize grain yield and yield components are affected by varying sowing dates and cultivar and N application rates. Farmers could enhance maize yields by manipulating the sowing date, cultivar selection and application of nitrogen fertilizer in agriculture production systems. The treatment effect of the sowing date and cultivar are significant on the biomass yield, 100-grain weight, seed number m−2, cob length and cob width. The grain yield, aboveground biomass, seed number m−2 and 100-grain wt. were reduced with a delay in the sowing date of either 7 or 14 days. The reduction in grain yield from SD1-SD2 (1.91 t ha−1), SD1-SD3 (2.90 t ha−1) and SD2-SD3 (0.99 t ha−1) were 21.04, 31.83 and 13.83%, respectively under the rainfed experiment.
The APSIM-Maize and DSSAT-CERES-Maize model calibration parameters are quite demanding. The models require complete data on the phenology, biomass, grain yield, grain number m−2, grain unit wt., rate of grain filling in milligram day−1 and leaf area index. The two crop simulation models provided acceptable simulation of DAP to anthesis, physiological maturity, biomass and grain yield. The models appear to be similar in their predictive ability in providing acceptable simulation. The crop models could be used as decision support tools by agricultural research stations, policymakers, farmers and planners for forecasting crop yield. Crop growth and yield can be simulated successfully using crop growth models before establishing costly and labor-intensive field experiments. Farmers and researchers can use crop models in investigating the effects of different management scenarios (plant density, planting time, fertilizer applications, irrigation regime and scheduling) and crop yield forecasting using the calibrated and validated maize cultivars.
The APSIM-Maize and DSSAT-CERES-Maize model simulates DAP to anthesis and maturity, but in this study, the function performed poorly. The plant growth duration took longer than expected due to lower solar radiation and cooler air temperature under irrigated condition. The key areas for APSIM-Maize and DSSAT-CERES-Maize models improvement includes the simulation of phenological stage duration and LAI prediction during winter and summer in the subtropical environment.

Author Contributions

C.B.C. and E.P. developed the study, designed the field experiments and selected the study area. C.B.C. monitored the field experiments, collected and prepared, and analyzed the data. C.B.C. and E.P. prepared the articles. C.B.C. created and prepared the tables and figures. V.R.N.C. proofread the articles and also checked the analysis and conclusion. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data can be requested from the corresponding author.

Acknowledgments

The researchers wish to thank Zambia Agricultural Research Institute (ZARI) for providing research facilities for both irrigated and rainfed experiments.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Summary of the observed data collected from the rainfed experiment site. Reproduced with permission from [33].
Table A1. Summary of the observed data collected from the rainfed experiment site. Reproduced with permission from [33].
SD1SD2SD3
CultivarsZMS 60630G19PHB 30B50ZMS 606PHB 30G19PHB 30B50ZMS 606PHB 30G19PHB 30B50
N rate123123123123123123123123123
Land preparation29 November 2016
Basal dressing/planting12 December 201626 December 20169 January 2017
Top dressing30 January 201717 February 20173 March 2017
Herbicides14 December 2016
Herbicides23 December 2016 & 18 January 2017
Weeding17 January 2017
Pesticides29 December 2016
Phenological stages
Emergence19 December 201619 December 201689 December 20164 January 20174 January 20173 January 201717 January 201716 January 201717 January 2017
V66 January 20176 January 20176 January 201720 January 201720 January 201719 January 20176 February 20176 February 20175 February 2017
R115 February 201715 February 201713 February 20174 March 20172 March 20174 March 201719 March 201719 March 201717 March 2017
R414 March 201719 March 201712 March 201728 March 201728 March 201726 March 201712 April 201712 April 201710 April 2017
R614 April 201715 April 201713 April 201728 April 201727 April 201727 April 201718 May 201718 May 201719 May 2017
Biomass sampling
V66 January 20176 January 20176 January 201720 January 201720 January 201720 January 20176 February 20176 February 20176 February 2017
R115 February 201715 February 201713 February 20174 March 20174 March 20172 March 201721 March 201721 March 201721 March 2017
R416 March 201716 March 201716 March 201730 March 201730 March 201728 March 201713 April 201713 April 201713 April 2017
Final harvest3 May 201715 May 20171 January 2017
1 = N1; 2 = N2; 3 = N3; Pesticide = Monocrotophos, Fustac; Herbicide = Nicosulfuron; Termites: Terminator (Imidacloprid 30.5% SC) 350 g of Imidacloprid/liter.
Table A2. Comparison between the observed and simulated phenology, mLAI, biomass, grain and grain size in the CERES-Maize model under the rainfed condition. Reproduced with permission from [61].
Table A2. Comparison between the observed and simulated phenology, mLAI, biomass, grain and grain size in the CERES-Maize model under the rainfed condition. Reproduced with permission from [61].
AnthesisMaturitymLAIBiomassGrainUnit wtGrain no
Sowing Date 1
TrtObsSim%PDObsSim%PDObsSim%PDObsSim%PDObsSim%PDObsSim%PDObsSim%PD
V1N16564−1.541231230.003.701.54−58.3811.4013.7020.189.208.80−4.350.370.382.15240123226.91
V1N26564−1.541231230.003.291.55−52.8910.9013.7025.698.708.700.000.380.380.53227623282.28
V1N36564−1.541231230.003.151.56−50.489.8013.8040.827.908.8011.390.370.370.002128234110.01
V2N165650.001241240.003.791.96−48.2811.2014.8032.148.408.804.760.370.382.70225122972.04
V2N265650.001241240.003.681.97−46.4710.5014.7040.007.908.8011.390.410.38−7.321922229919.61
V2N365650.001241240.004.451.98−55.5113.1014.8012.989.708.80−9.280.380.380.0025192305−8.50
V3N163641.591221220.004.221.54−63.5112.5014.4015.209.5010.106.320.510.510.00186719966.91
V3N263641.591221220.004.561.55−66.0112.7014.5014.179.6010.206.250.500.512.00195220012.51
V3N363641.591221220.005.001.56−68.8014.2014.502.1110.8010.20−5.560.520.51−1.9220632012−2.47
Sowing Date 2
V1N16764−4.481231273.253.571.56−56.309.1012.0031.876.907.6010.140.320.3612.50211321953.88
V1N26764−4.481231273.253.451.56−54.7810.4011.9014.428.207.50−8.540.340.340.0024212194−9.38
V1N36764−4.481231273.253.261.60−50.928.6012.0039.536.507.6016.920.300.3413.33219022181.28
V2N16665−1.521221284.923.751.91−49.0710.2012.8025.497.307.604.110.340.365.88214621460.00
V2N26665−1.521221284.923.581.93−46.099.3012.7036.566.507.5015.380.300.3516.67212521541.36
V2N36665−1.521221284.924.081.94−52.4511.2012.9015.187.907.70−2.530.340.365.8823112158−6.62
V3N16862−8.821221230.823.641.47−59.6210.5012.1015.247.907.900.000.430.42−2.33186119082.53
V3N26862−8.821221230.823.201.47−54.067.3012.0064.385.307.9049.060.350.4117.141497190727.39
V3N36862−8.821221230.823.181.50−52.8311.0012.3011.828.008.101.250.370.4213.51190819291.10
Sowing Date 3
V1N16665−1.521291300.783.181.60−49.6910.1010.806.936.005.80−3.330.320.21−34.382300279321.43
V1N26665−1.521291300.783.341.59−52.408.9010.7020.226.505.80−10.770.310.21−32.262425275513.61
V1N36665−1.521291300.783.341.59−52.408.7010.6021.846.105.60−8.200.320.20−37.502068278934.86
V2N16968−1.451291343.883.612.18−39.619.5011.7023.166.105.30−13.110.310.19−38.712043276835.49
V2N26968−1.451291343.884.022.18−45.779.6011.8022.926.105.50−9.840.330.20−39.392243276523.27
V2N36968−1.451291343.883.402.17−36.1810.1011.9017.826.105.60−8.200.330.20−39.552331274517.76
V3N16762−7.46130127−2.313.551.45−59.158.4010.3022.625.605.50−1.790.350.24−31.431568229146.11
V3N26762−7.46130127−2.313.591.44−59.8910.0010.101.006.905.40−21.740.330.24−27.271628226639.19
V3N36762−7.46130127−2.313.721.44−61.299.3010.108.606.305.30−15.870.290.23−20.691521229851.08
DAP = Days after planting (anthesis, physiological maturity); mLAI (m2 m−2) = maximum leaf area index, biomass (t ha−1), unity weight (g); Negative deviations indicate under prediction while positive deviations indicate over prediction; %PD = percentage prediction deviation, Sim = Simulated; Obs = observed; DAP = Days after planting; mLAI = maximum leaf area index; V1 = ZMS 606; V2 = PHB 30G19; V3 = PHB 30B50; Trt = Treatment.
Table A3. Comparison between the observed and simulated phenology, mLAI, biomass, grain and grain size in the APSIM-Maize model under the rainfed condition. Reproduced with permission from [61].
Table A3. Comparison between the observed and simulated phenology, mLAI, biomass, grain and grain size in the APSIM-Maize model under the rainfed condition. Reproduced with permission from [61].
AnthesisMaturitymLAIBiomassGrainUnit wtGrain no
Sowing Date 1
TrtObsSim%PDObsSim%PDObsSim%PDObsSim%PDObsSim%PDObsSim%PDObsSim%PD
V1N165650.001231240.813.703.18−14.0511.4015.4035.099.208.90−3.260.370.394.8424012261−0.86
V1N265650.001231240.813.293.19−3.0410.9015.4041.288.708.902.300.380.393.1722762262−0.62
V1N365650.001231240.813.153.191.279.8015.4057.147.908.9012.660.370.395.41212822626.30
V2N165650.001241240.003.793.19−15.8311.2015.5038.398.408.905.950.370.382.70225123243.24
V2N265650.001241240.003.683.19−13.3210.5015.5047.627.908.9012.660.410.38−7.321922232420.92
V2N365650.001241240.004.453.19−28.3113.1015.5018.329.708.90−8.250.380.380.0025192324−7.74
V3N163653.171221230.824.223.14−25.5912.5015.3022.409.508.00−15.790.510.43−15.6918671851−0.86
V3N263653.171221230.824.563.14−31.1412.7015.3020.479.608.00−16.670.500.43−14.0019521852−5.12
V3N363653.171221230.825.003.14−37.2014.2015.307.7510.808.00−25.930.520.43−17.3120631852−10.23
Sowing Date 2
V1N167670.001231305.693.573.34−6.449.1015.4069.236.908.7026.090.320.4025.00211321752.93
V1N267670.001231305.693.453.34−3.1910.4015.8051.928.208.706.100.340.4017.6524212175−10.16
V1N367670.001231305.693.263.342.458.6015.8083.726.508.7033.850.300.4033.3321902175−0.68
V2N166671.521221284.923.753.34−10.9310.2015.8054.907.308.8020.550.340.3811.76214622906.71
V2N266671.521221284.923.583.34−6.709.3015.8069.896.508.8035.380.300.3826.67212522907.76
V2N366671.521221284.924.083.35−17.8911.2015.8041.077.908.8011.390.340.3811.7623112290−0.91
V3N16865−4.411221252.463.643.29−9.6210.5015.6048.577.907.70−2.530.430.430.0018611790−3.82
V3N26865−4.411221252.463.203.292.817.3015.60113.705.307.7045.280.350.4322.861497179019.57
V3N36865−4.411221252.463.183.293.4611.0015.6041.828.007.70−3.750.370.4316.2219081790−6.18
Sowing Date 3
V1N16665−1.521291300.783.183.05−4.0910.1013.5033.666.007.5025.000.320.29−9.382300258612.43
V1N26665−1.521291300.783.343.05−8.688.9013.5051.696.507.5015.380.310.29−6.45242525866.64
V1N36665−1.521291300.783.343.05−8.688.7013.5055.176.107.5022.950.320.29−9.382068258625.05
V2N16966−4.351291290.003.613.05−15.519.5013.5042.116.107.3019.670.310.28−9.682043258226.38
V2N26966−4.351291290.004.023.05−24.139.6013.5040.636.107.3019.670.330.28−15.152243258215.11
V2N36966−4.351291290.003.403.05−10.2910.1013.5033.666.107.3019.670.330.28−15.152331258210.77
V3N16764−4.48130128−1.543.553.02−14.938.4013.6061.905.607.2028.570.350.350.001568205330.93
V3N26764−4.48130128−1.543.593.02−15.8810.0013.6036.006.907.204.350.330.356.061628205326.11
V3N36764−4.48130128−1.543.723.02−18.829.3013.6046.246.307.2014.290.290.3520.691521205334.98
DAP = Days after planting (anthesis, physiological maturity); mLAI (m2 m−2) = maximum leaf area index, biomass (t ha−1), unity weight (g); Negative deviations indicate under prediction while positive deviations indicate over prediction; %PD = percentage prediction deviation, Sim = Simulated; Obs = observed; DAP = Days after planting; mLAI = maximum leaf area index; V1 = ZMS 606; V2 = PHB 30G19; V3 = PHB 30B50; Trt = Treatment.

References

  1. Lukeba, J.-C.L.; Vumilia, R.K.; Nkongolo, K.C.K.; Mwabila, M.L.; Tsumbu, M. Growth and leaf area index simulation in maize (Zea mays L.) under small-scale farm conditions in a sub-Saharan African region. Am. J. Plant Sci. 2013, 4, 575–583. [Google Scholar] [CrossRef] [Green Version]
  2. Chisanga, C.B.; Phiri, E.; Chizumba, S.; Sichingabula, H. Evaluating CERES-maize model using planting dates and nitrogen fertilizer in Zambia. J. Agric. Sci. 2015, 7, 79–97. [Google Scholar] [CrossRef] [Green Version]
  3. Jones, J.W.; Hoogenboom, G.; Porter, C.H.; Boote, K.J.; Batchelor, W.D.; Hunt, L.A.; Wilkens, P.W.; Singh, U.; Gijsman, A.J.; Ritchie, J.T. The DSSAT cropping system model. Eur. J. Agron. 2003, 18, 235–265. [Google Scholar] [CrossRef]
  4. Tsimba, R. Development of a Decision Support System to Determine the Best Maize (Zea mays., L ) Hybrid—Planting Date Option under Typical New Zealand Management Systems; Massey University: Palmerston North, New Zealand, 2011. [Google Scholar]
  5. Soltani, A.; Hoogenboom, G. Minimum data requirements for parameter estimation of stochastic weather generators. Clim. Res. 2003, 25, 109–119. [Google Scholar] [CrossRef] [Green Version]
  6. Keating, B.A.; Carberry, P.S.; Hammer, G.L.; Probert, M.E.; Robertson, M.E.; Holzworth, N.I.H.; Hargreaves, J.N.G.; Meinke, H.; Hochman, Z.; McLean, G.; et al. An overview of APSIM, a model designed for farming systems simulation. Eur. J. Agron. 2003, 18, 267–288. [Google Scholar] [CrossRef] [Green Version]
  7. McCown, R.L.; Hammer, G.L.; Hargreaves, J.N.G.; Holzworth, D.P.; Freebairn, D.M. APSIM: A novel software system for model development, model testing and simulation in agricultural systems research. Agric. Syst. 1996, 50, 255–271. [Google Scholar] [CrossRef]
  8. Hoogenboom, G.; Jones, J.W.; Wilkens, P.W.; Porter, C.H.; Batchelor, W.D.; Hunt, L.A.; Boote, K.J.; Singh, U.; Uryasev, O.; Bowen, W.T.; et al. Decision Support System for Agrotechnology Transfer; Version 4.0 [CD-ROM]; University of Hawaii: Honolulu, HI, USA, 2004. [Google Scholar]
  9. Lizaso, J.I.; Boote, K.J.; Jones, J.W.; Porter, C.H.; Echarte, L.; Westgate, M.E.; Sonohat, G. CSM-IXIM: A New Maize simulation model for DSSAT version 4.5. Agron. J. 2011, 103, 766–779. [Google Scholar] [CrossRef]
  10. Jones, D.S.; Creel, R.C.; Rios, B.A. Carbon isotope stratigraphy and correlation of depositional sequences in the Upper Ordovician Ely Springs Dolostone, eastern Great Basin, USA. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2016, 458, 85–101. [Google Scholar] [CrossRef]
  11. Timsina, J.; Godwin, D.; Humphreys, E.; Yadvinder-Singh Bijay-Singh Kukal, S.S.S.; Smith, D. Evaluation of options for increasing yield and water productivity of wheat in Punjab, India using the DSSAT-CSM-CERES-Wheat model. Agric. Water Manag. 2008, 95, 1099–1110. [Google Scholar] [CrossRef]
  12. Jones, J.W.; Tsuji, G.Y.; Hoogenboom, G.; Hunt, L.A.; Thornton, P.K.; Wilkens, P.W.; Imamura, D.T.; Bowen, W.T.; Singh, U. Underst. Options Agric. Prod. Syst. Approaches Sustain. Agric. Dev. In Decision Support System for Agrotechnology Transfer: DSSAT v3.; Tsuji, G.Y., Hoogenboom, G., Thornton, P.K., Eds.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1998; Volume 7, pp. 157–177. [Google Scholar]
  13. Sadras, V.; Rodriguez, D. Modelling the nitrogen-driven trade-off between nitrogen utilisation efficiency and water use efficiency of wheat in eastern Australia. Field Crop. Res. 2010, 118, 297–305. [Google Scholar] [CrossRef]
  14. Gaydon, D.S.; Wang, E.; Poulton, P.L.; Ahmad, B.; Ahmed, B.; Akhter, S.; Ali IAmarasingha, R.; Chaki, A.K.; Chen, C.; Choudhury, B.U.; et al. Evaluation of the APSIM model in cropping systems of Asia. Field Crop. Res. 2017, 204, 52–75. [Google Scholar] [CrossRef]
  15. Snow, V.O.; Houlbrooke, D.J.; Huth, N.I. Predicting soil water, tile drainage, and runoff in a mole-tile drained soil. N. Z. J. Agric. Res. 2007, 50, 13–24. [Google Scholar] [CrossRef]
  16. Wang, N.; Wang, E.; Wang, J.; Zhang, J.; Zheng, B.; Huang, Y.; Tan, M. Modelling maize phenology, biomass growth and yield under contrasting temperature conditions. Agric. For. Meteorol. 2018, 250–251, 319–329. [Google Scholar] [CrossRef]
  17. Monsi, M.; Saeki, T. On the factor light in plant communities and its importance for matter production. Ann. Bot. 2004, 95, 549–567. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Spitters, C.J.T.; Toussaint, H.A.J.M.; Goudriaan, J. Separating the diffuse and direct component of global radiation and its implications for modeling canopy photosynthesis Part I. Components of incoming radiation. Agric. For. Meteorol. 1986, 38, 217–229. [Google Scholar] [CrossRef]
  19. Ritchie, J.T. Model for predicting evaporation from row crops with incomplete cover. Water Resour. Res. 1972, 8, 1204–1213. [Google Scholar] [CrossRef] [Green Version]
  20. Achieng, J.; Kananpiu, F.; Chauhan, Y.; Rodriguez, D. Sustainable intensification of maize-bean production among smallholder farmers in western Kenya. In Proceedings of the 5th World Congress on Conservation Agriculture, Brisbane, Australia, 26–29 September 2011. [Google Scholar]
  21. Holzworth, D.P.; Huth, N.I.; DeVoil, P.G.; Zurcher, E.J.; Herrman, N.I.; McLean, G.; Chenu, K.; van Oosterom, E.J.; Show, V.; Murphy, C.; et al. APSIM—Evolution towards a new generation of agricultural systems simulation. Environ. Model Softw. 2014, 62, 327–350. [Google Scholar] [CrossRef]
  22. Gaydon, D. The APSIM Model—An overview. In The SAARC-Australia Project-Developing Capacity in Cropping Systems Modelling for South Asia; Gaydon, D.S., Saiyed, I., Roth, C.R., Eds.; SAARC Agriculture. Centre: Dhaka, Bangladesh, 2014; pp. 15–31. [Google Scholar]
  23. Jones, P.G.; Thornton, P.K. The potential impacts of climate change on maize production in Africa and Latin America in 2055. Glob. Environ. Chang. 2003, 13, 51–59. [Google Scholar] [CrossRef]
  24. Tsimba, R.; Edmeades, G.O.; Millner, J.P.; Kemp, P.D. The effect of planting date on maize grain yields and yield components. Field Crop. Res. 2013, 150, 145–155. [Google Scholar] [CrossRef]
  25. Bationo, A.; Tabo, R.; Kihara, J.; Hoogenboom, G.; Traore, P.C.S.; Boote, K.J.; Jones, J.W. Building capacity for modeling in Africa. In Improving Soil Fertility Recommendations in Africa Using the Decision Support System for Agrotechnology Transfer (DSSAT); Springer: Berlin/Heidelberg, Germany, 2012; pp. 1–7. [Google Scholar]
  26. Ahmed, M.; Aslam, M.A.; Hassan, F.U.; Asif, M.; Hayat, R. Use of APSIM to model nitrogen use efficiency of rain-fed wheat. Int. J. Agric. Biol. 2014, 16, 461–470. [Google Scholar]
  27. Soil Survey Staff. Keys to Soil Taxonomy, 12th ed.; USDA-Natural Resources Conservation Service: Washington, DC, USA, 2014.
  28. Du Plessis, J. Maize production. Resour. Cent. Dir. Agric. Inf. Serv. 2003, 1–38. [Google Scholar]
  29. Araya, A.; Hoogenboom, G.; Luedeling, E.; Hadgu, K.M.; Kisekka, I.; Martorano, L.G. Assessment of maize growth and yield using crop models under present and future climate in southwestern Ethiopia. Agric. For. Meteorol. 2015, 214, 252–265. [Google Scholar] [CrossRef]
  30. Sato, J.H.; de Figueiredo, C.C.; Marchão, R.L.; Madari, B.E.; Benedito, L.E.C.; Busato, J.G.; de Souza, D.M. Methods of soil organic carbon determination in Brazilian savannah soils. Sci. Agric. 2014, 71, 302–308. [Google Scholar] [CrossRef] [Green Version]
  31. Saxton, K.E.; Willey, P.H.; Rawls, W.J. Field and pond hydrologic analyses with the SPAW Model. In Proceedings of the 2006 ASABE Annual International Meeting, Portland, OR, USA, 9–12 July 2006; pp. 1–14. [Google Scholar]
  32. Saxton, K.E.; Rawls, W.J. Soil water characteristic estimates by texture and organic matter for hydrologic solutions. Soil. Sci. Soc. Am. J. 2006, 70, 1569–1578. [Google Scholar] [CrossRef] [Green Version]
  33. Chisanga, C.B.; Phiri, E.; Chinene, V.R.N.; Chabala, L.M. Projecting maize yield under local-scale climate change scenarios using crop models: Sensitivity to sowing dates, cultivar, and nitrogen fertilizer rates. Food Energy Secur. 2020, 9, e231. [Google Scholar] [CrossRef]
  34. Sentek. Calibration Manual: For Sentek Soil Moisture Sensors Version 2.0; Sentek Pty Ltd.: Stepney, Australia, 2011; p. 5069. [Google Scholar]
  35. Asseng, S.; Keating, B.A.; Fillery, I.R.P.; Gregory, P.J.; Bowden, J.W.; Turner, N.C.; Palta, J.A.; Abrecht, D.G. Performance of the APSIM-wheat model in Western Australia. Field Crop Res. 1998, 57, 163–179. [Google Scholar] [CrossRef]
  36. Hoogenboom, G.; Wilkens, P.W.; Tsuji, G.Y. DSSAT v3; University of Hawaii: Honolulu, HI, USA, 1999; Volume 4. [Google Scholar]
  37. Karuma, A.N.; Gachene, C.K.K.; Gicheru, P.T.; Mtakwa, P.W.; Amuri, N. Effects of tillage and cropping systems on maize and beans yield and selected yield components in a semi-arid area of Kenya. Trop. Subtrop. Agroecosyst. 2016, 19, 167–179. [Google Scholar]
  38. NSW Department of Primary Industries. Maize Growth and Development; State of New South Wales: Sydney, Australia, 2009.
  39. Keating, B.A.; Gaydon, D.; Huth, N.I.; Probert, M.E.; Verburg, K.; Smith, C.J.; Bond, W. Use of modelling to explore the water balance of dryland farming systems in the Murray-Darling Basin, Australia. Eur. J. Agron. 2002, 18, 159–169. [Google Scholar] [CrossRef]
  40. Brown, H.E.; Huth, N.I.; Holzworth, D.P.; Teixeira, E.I.; Zyskowski, R.F.; Hargreaves, J.N.G.; Moot, D.J. Plant modelling framework: Software for building and running crop models on the APSIM platform. Environ. Model Softw. 2014, 62, 385–398. [Google Scholar] [CrossRef] [Green Version]
  41. Archontoulis, S.V.; Miguez, F.E.; Moore, K.J. A methodology and an optimization tool to calibrate phenology of short-day species included in the APSIM PLANT model: Application to soybean. Environ. Model Softw. 2014, 62, 465–477. [Google Scholar] [CrossRef]
  42. Dietzel, R. A Comparison of Carbon Storage Potential in Corn- and Prairie-Based Agroecosystems; Iowa State University: Ames, IA, USA, 2014. [Google Scholar]
  43. Jones, J.W.; Hoogenboom, G.; Wilkens, P.W.; Porter, C.H.; Tsuji, G.Y. Decision Support System for Agrotechnology Transfer Version 4.0; DSSAT v4.5: Crop Model Documentation; University of Hawaii: Honolulu, HI, USA, 2010; Volume 4. [Google Scholar]
  44. Jones, J.W.; Hoogenboom, G.; Wilkens, P.W.; Porter, C.H.; Tsuji, G.Y. Decision Support System for Agrotechnology Transfer Version 4.5; DSSAT v4.5: ICASA Tools; University of Hawaii: Honolulu, HI, USA, 2010; Volume 3. [Google Scholar]
  45. Jones, C.A.; Kiniry, J.R. CERES-Maize: A Simulation Model of Maize Growth and Development; Texas A&M University Press: College Station, TX, USA, 1986. [Google Scholar]
  46. Kiniry, J.R.; Williams, J.R.; Vanderlip, R.L.; Atwood, J.D.; Reicosky, D.C.; Mulliken, J.; Cox, W.J.; Mascagni, H.J., Jr.; Hollinger, S.E.; Wiebold, W.J. Evaluation of two maize models for nine U.S. locations. Agron. J. 1997, 89, 421–426. [Google Scholar] [CrossRef]
  47. Keating, B.A.; Wafula, B.M.; McCown, R.L. Simulation of plant density effects on maize yield as influenced by water nitrogen limitations. In Proceedings of the International Plant Physiology Conference, New Delhi, India, 15–20 February 1988; pp. 547–559. [Google Scholar]
  48. Ritchie, J.T.; Porter, C.H.; Suleiman, A.A.; Judge, J.; Jones, J.W. Extension of an existing model for soil water evaporation and redistribution under high water content conditions. Soil Sci. Soc. Am. J. 2009, 73, 792. [Google Scholar] [CrossRef]
  49. Valeriano, A.; Ines, M.; Droogers, P.; Makin, I.W.; Das Gupta, A.; Loof, R.; Clemente, R.S.; Kyoshi, H. Crop Growth and Soil Water Balance Modeling to Explore Water Management Option; International Water Management Institute: Colombo, Sri Lanka, 1993. [Google Scholar]
  50. Dalgliesh, N.; Hochman, Z.; Huth, N.; Holzworth, D. Field Protocol to APSoil characterizations. Version 4; CSIRO: Canberra, Australia, 2016. [Google Scholar]
  51. Mesfin, T.; Moeller, C.; Parsons, D.; Meinke, H. Maize (Zea mays, L.) productivity as influenced by sowing date and nitrogen fertiliser rate at Melkassa, Ethiopia: Parameterisation and evaluation of APSIM-Maize. In Proceedings of the Building Productive, Diverse and Sustainable Landscapes, Proceedings of 17th ASA Conference, Hobart, Australia, 20–24 September 2015; pp. 20–24. Available online: https://www.semanticscholar.org/paper/Maize-(-Zea-mays-L.)-productivity-as-influenced-by-Mesfin-Moeller/8f7fb30537061873615af59f9ce90514424ec832 (accessed on 23 September 2021).
  52. Wang, E.; Chen, C.; Yu, Q. Modeling the response of wheat and maize productivity to climate variability and irrigation in the North China Plain. In Proceedings of the 18th World IMACS/MODSIM Congress, Cairns, Australia, 13–17 July 2009; pp. 13–17. [Google Scholar]
  53. Ahmed, M. APSIM and DSSAT models as decision support tools. In Proceedings of the 19th International Congress on Modelling and Simulation, Perth, Australia, 12–16 December 2011; pp. 1174–1180. [Google Scholar]
  54. Dokoohaki, H.; Gheysari, M.; Mousavi, S.-F.; Hoogenboom, G. Effects of different irrigation regimes on soil moisture availability evaluated by CSM-CERES-Maize model under semi-arid condition. Ecohydrol. Hydrobiol. 2017, 17, 207–216. [Google Scholar] [CrossRef] [Green Version]
  55. De Mendiburu, F. Agricolae: Statistical Procedures for Agricultural Research. R Package Version 1.2-4. 2016. Available online: https://cran.r-project.org/web/packages/agricolae/index.html (accessed on 27 July 2017).
  56. Harrison, S.R. Regression of a model on real-system output: An invalid test of model validity. Agric. Syst. 1990, 34, 183–190. [Google Scholar] [CrossRef]
  57. Willmott, C.J. Some comments on the evaluation of model performance. Bull. Am. Meteorol. Soc. 1982, 63, 1309–1313. [Google Scholar] [CrossRef] [Green Version]
  58. MacCarthy, D.S.; Kihara, J.; Adiku, S.G.K. Decision support tools for site-specific fertilizer recommendations and agricultural planning in selected countries in sub-Sahara Africa. Nutr. Cycl. Agroecosyst. 2017, 110, 265–289. [Google Scholar]
  59. Zhang, Y.; Feng, L.; Wang, E.; Wang, J.; Li, B. Evaluation of the APSIM-Wheat model in terms of different cultivars, management regimes and environmental conditions. Can. J. Plant. Sci. 2012, 92, 937–949. [Google Scholar] [CrossRef] [Green Version]
  60. Jing, Q.; Shang, J.; Huffman, T.; Qian, B.; Pattey, E.; Liu, J.; Dong, T.; Drury, C.F. Using the CSM–CERES–Maize model to assess the gap between actual and potential yields of grain maize. J. Agric. Sci. 2017, 155, 239–260. [Google Scholar] [CrossRef]
  61. Chisanga, C.B. Interactive Effects of N Fertilization Rate, Cultivars and Planting Date under Climate Change on Maize (Zea mays L.) Yield Using Crop Simulation and Statistical Downscaling of Climate Models; University of Zambia: Lusaka, Zambia, 2019. [Google Scholar]
  62. Chisanga, C.B.; Phiri, E.; Chinene, V.R.N. Evaluation of sowing date and fertilization with nitrogen in maize cultivars in rainy conditions in Zambia. Afr. J. Plant Sci. 2019, 13, 221–230. [Google Scholar]
  63. Ahmad, M.; Khaliq, A.; Ahmad, R.; Ranjha, A.M. Allometery and productivity of autumn planted maize hybrids under narrow row spacing. Int. J. Agric. Biol. 2010, 12, 661–667. [Google Scholar]
  64. Jindo, K.; Schut, A.G.T.; Langeveld, J.W.A. Sustainable intensification in Western Kenya: Who will benefit? Agric. Syst. 2020, 182, 102831. [Google Scholar] [CrossRef]
  65. Gul, S.; Khan, M.H.; Khanday, B.A.; Nabi, S. Effect of Sowing Methods and NPK Levels on Growth and Yield of Rainfed Maize (Zea mays L.); Hindawi Publishing Corporation: London, UK, 2015; pp. 1–7. [Google Scholar]
  66. Adnan, A.A.; Jibrin, J.M.; Kamara, A.Y.; Abdulrahman, B.L.; Shaibu, A.S.; Garba, I.I. CERES-Maize model for determining the optimum planting dates of early maturing maize varieties in Northern Nigeria. Front. Plant Sci. 2017, 8, 1118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Sharifi, R.S.; Namvar, A. Effects of time and rate of nitrogen application on phenology and some agronomical traits of maize (Zea mays, L.). Biologija 2016, 62, 35–45. [Google Scholar] [CrossRef]
  68. Zeidan, M.S.; Amany, A.; Bahr El-Kramany, M.F. Effect of N fertilizer and plant density on yield and quality of maize in sandy soil. Res. J. Agric. Biol. Sci. 2006, 2, 156–161. [Google Scholar]
  69. Arif, M.; Jan, M.T.; Khan, N.U.; Akbar, H.; Khan, S.A.; Khan, M.J.; Khan, A.; Munir, I.; Saeed, M.; Iqbal, A. Impact of plant populations and nitrogen levels on maize. Pak. J. Bot. 2010, 42, 3907–3913. [Google Scholar]
  70. Sangoi, L. Understanding plant density effects on maize growth and development: An important issue to maximize grain yield. Ciência Rural St. Maria 2001, 31, 159–168. [Google Scholar] [CrossRef] [Green Version]
  71. Singh, J.; Hadda, M.S. Phenology and thermal indices of maize (Zea mays L.) influenced by subsoil compaction and nitrogen fertilization under semi-arid irrigated conditions. J. Appl. Nat. Sci. 2014, 6, 349–355. [Google Scholar] [CrossRef] [Green Version]
  72. Feleke, H.G.; Savage, M.J.; Tesfaye, K. Calibration and validation of APSIM-Maize, DSSAT CERES-Maize and AquaCrop models for Ethiopian tropical environments. S. Afr. J. Plant Soil. 2021, 38, 36–51. [Google Scholar] [CrossRef]
  73. Beah, A.; Kamara, A.Y.; Jibrin, J.M.; Akinseye, F.M.; Tofa, A.I.; Ademulegun, T.D. Simulation of the optimum planting windows for early and intermediate-maturing maize varieties in the Nigerian savannas using the APSIM model. Front. Sustain. Food Syst. 2021, 5, 624886. [Google Scholar] [CrossRef]
  74. Kisaka, O.M.; Mucheru-Muna, M.; Ngetich, F.K.; Mugwe, J.N.; Mugendi, D.N.; Mairura, F.; Muriuki, J. Using APSIM-model as a decision-support-tool for long-term integrated-nitrogen-management and maize productivity under semi-arid conditions in Kenya. Exp. Agric. 2016, 52, 279–299. [Google Scholar] [CrossRef]
  75. Groenendyk, D.; Kaleita, A.; Thorp, K. Assimilating in situ soil moisture measurements into the DSSAT-CSM using a Kalman filter. In Proceedings of the 2011 ASABE Annual International Meeting, Louisville, KY, USA, 7–10 August 2011; pp. 4319–4349. [Google Scholar]
  76. Nangia, V.; Ahmad M ud, D.; Jiantao, D.; Changrong, Y.; Hoogenboom, G.; Xurong, M.; Wenqing, H.; Shuang, L.; Qin, L. Modeling the field-scale effects of conservation agriculture on land and water productivity of rainfed maize in the Yellow River Basin, China. Int. J. Agric. Biol. Eng. 2010, 3, 5–17. [Google Scholar]
  77. Ollenburger, M.H. Modeling Integrated Soil Fertility Management Options in Malawi; Michigan State University: East Lansing, MI, USA, 2012. [Google Scholar]
  78. Tsimba, R.; Edmeades, G.O.; Millner, J.P.; Kemp, P.D. Optimising maize planting date and hybrid selection using simulation modeling. Agron. N. Z. 2013, 43, 55–62. [Google Scholar]
Figure 1. Daily weather data collected during the rainfed field experiment at Mt Makulu during the 2016/2017 season.
Figure 1. Daily weather data collected during the rainfed field experiment at Mt Makulu during the 2016/2017 season.
Nitrogen 02 00027 g001
Figure 2. Daily weather data (June–November 2016) collected during the irrigated field experiment at Mount Makulu.
Figure 2. Daily weather data (June–November 2016) collected during the irrigated field experiment at Mount Makulu.
Nitrogen 02 00027 g002
Figure 3. Overall calibration statistics between the observed and simulated values for the APSIM-Maize and DSSAT-CERES-Maize models using 1:1 graphs under rainfed conditions. The graphs show simulated vs. observed values for the anthesis (DAP), maturity (DAP), grain yield (t ha−1), aboveground biomass (t ha−1), mLAI (m2 m−2) and grain number m−2.
Figure 3. Overall calibration statistics between the observed and simulated values for the APSIM-Maize and DSSAT-CERES-Maize models using 1:1 graphs under rainfed conditions. The graphs show simulated vs. observed values for the anthesis (DAP), maturity (DAP), grain yield (t ha−1), aboveground biomass (t ha−1), mLAI (m2 m−2) and grain number m−2.
Nitrogen 02 00027 g003
Table 1. Physical and chemical soil analyses under rainfed conditions. Reproduced with permission from [33].
Table 1. Physical and chemical soil analyses under rainfed conditions. Reproduced with permission from [33].
Depth (Cm)0–2020–4040–6060–8080–100Analysis Method
pH (water)7.307.207.507.707.601:5 soil water
Total N (%)0.0310.0420.0540.0610.036Modified Kjeldahl method
NO3N29.9048.7056.4070.1042.80
NH4N18.0029.2033.9042.1025.70
P extractable (mg kg−1)10.0011.0010.0018.0012.00Bray 1
K (mg kg1)1.050.991.120.590.89Ammonium acetate
Ca (cmol(+) kg−1)11.009.303.402.903.20Ammonium acetate
Mg (cmol(+) kg−1)3.502.702.301.001.30Ammonium acetate
OC (%)0.350.570.660.820.50Walkley & Black method
OM (%)0.6020.9801.1351.4100.860
CEC (cmol(+) kg−1)15.5713.026.854.525.42Ammonium acetate
Bulk density (g cm−3)1.431.411.411.461.36SPAW
Silt (%)12.8016.8012.8018.802.80Hydrometer method
Sand (%)39.6035.6037.6041.6037.60
Clay (%)47.6047.6049.6039.6059.60
Soil textureclayclayclayclayclaySPAW
LL0.2870.2870.2990.2440.350SPAW
DUL0.4070.4090.4190.3630.470
SAT0.4590.4670.4680.4470.487
SHC (mm h−1)0.3500.5000.2901.4800.010
Phosphorus; LL: lower limit (Wilting point); DUL: drained upper limit (Field Capacity); SAT: saturation; SHC: Saturated hydraulic conductivity; SPAW: Soil–Plant–Air–Water.
Table 2. Growth and development stages of maize. Reproduced with permission from [38].
Table 2. Growth and development stages of maize. Reproduced with permission from [38].
Vegetative StageReproductive Stage
Emergence (VE)silking (R1)
first leaf with collar (V1)blister (R2)
second leaf with full collar (V2)milk (R3)
third leaf with full collar (V3)dough (R4)
nth leaf with full collar (V(n))dent (R5)
tasseling (VT)maturity (R6)
Table 3. Treatment effect of the sowing date, cultivar and N on the grain wt., grain number m−2, grain wt., stover wt., aboveground biomass wt. and cob wt. under rainfed conditions.
Table 3. Treatment effect of the sowing date, cultivar and N on the grain wt., grain number m−2, grain wt., stover wt., aboveground biomass wt. and cob wt. under rainfed conditions.
Treatment/CultivarGrain wt. (g m−2)Grain No m−2Grain Unit wt. (g)Stover wt. (g m−2)Biomass wt. (g m−2)Cob wt. (g m−2)
SD1907.6 a2153 a42.35 a277.5 a1185.0 a158.0 a
SD2716.3 b2064 a34.27 b258.0 b974.3 b140.8 ab
SD3617.7 c2014 a31.86 b314.5 b932.2 b128.4 b
Significance***ns************
Tukey HSD 5%90.7817.901.2531.30115.2624.95
CV %18.5257.84.616.8417.055.1
ZMS 606732.6 a2258 a33.62 b234.5 b967.1 a109.10 b
P30B19733.9 a2210 a34.52 b318.0 a1052.0 a157.8 a
P30G50775.2 a1763 b40.34 a297.4 a1073.0 a160.3 a
Significancens*********ns***
Tukey HSD 5%61.78211.382.803.3524.2080.0111.43
CV %1515.5214.1215.514.13.3
Nitrogen (N) rate (N1)742.3 a2061 a36.70 a278.2 a1020.0 a137.37 a
Nitrogen (N) rate (N2)729.2 a2054 a35.91 a272.6 a1002.0 a137.74 a
Nitrogen (N) rate (N3)770.2 a2116 a35.87 a299.2 a1069.0 a152.07 a
Significancensnsnsnsnsns
Tukey HSD 5%90.7913.743.4331.88124.616.18
CV %18.544.21618.5220.45.6
Interaction (SD ∗ V) significancensns****nsns
Interaction (V ∗ N) significancensnsnsnsnsns
Interaction (V ∗ SD ∗ N) significancensnsnsnsnsns
Means sharing the same letter in the table do not differ statistically at p < 0.05; LSD = Least Mean Differences; ** = highly significant at the 5% level; *** = very highly significant at the 5% level; ns = Nonsignificant; At R6, weight = g m−2 square meter (g m−2 * 10 = t ha−1); wt. = weight.
Table 4. Treatment effect of the sowing date, cultivar and N on the husk wt., stem wt., leaf blade wt., harvest index, cob width and cob length under rainfed conditions.
Table 4. Treatment effect of the sowing date, cultivar and N on the husk wt., stem wt., leaf blade wt., harvest index, cob width and cob length under rainfed conditions.
Treatment/CultivarHusk wt. (g m−2)Stem wt. (g m−2)Leaf Blade wt. (g m−2)HICob Width (cm)Cob Length (cm)
SD133.30 b47.77 b20.78 b0.77 a5.49 a21.43 a
SD228.06 b47.59 b25.20 b0.73 b5.29 b21.3 a
SD355.06 a68.17 a35.03 a0.66 c5.04 c18.43 b
Significance****************
Tukey HSD 5%13.4039.786.450.020.171.94
CV %45.7012.6231.604.645.1011.41
ZMS 60637.73 ab51.23 b19.89 b0.75 a5.17 b19.48 a
P30B1944.72 a62.15 a31.59 a0.69 c5.51 a20.78 a
P30G5033.97 b50.15 b29.53 a0.72 b5.14 b20.91 a
Significance************ns*
Tukey HSD 5%2.1831.565.770.020.161.85
CV %32.0013.7236.003.83.307.8
Nitrogen (N) rate (N1)21.83 a54.88 a271.2 a0.72 a5.27 a20.33 a
Nitrogen (N) rate (N2)19.53 a51.56 a279.1 a0.73 a5.24 a20.33 a
Nitrogen (N) rate (N3)20.49 a57.08 a310.3 a0.72 a5.31 a20.51 a
Significancensnsnsnsnsns
Tukey HSD 5%2.3437.9354.340.020.100.94
CV %26.2015.3227.604.15.6016.5
Interaction (SD ∗ V) significancensnsnsnsns***
Interaction (V∗ N) significancens*nsnsnsns
Interaction (V ∗ SD ∗ N) significancensnsnsnsnsns
Means sharing the same letter in the table do not differ statistically at p < 0.05; LSD = Least Mean Differences; * = significant at the 5% level; ** = highly significant at the 5% level; *** = very highly significant at the 5% level; ns = Nonsignificant; At R6, weight = g m−2 square meter (g m−2 * 10 = t ha−1); wt. = weight.
Table 5. Treatment effect of the cultivar and N on the grain wt., grain number m−2, 100-grain wt., stover wt., aboveground biomass wt., grain No. m−2, cob wt., ear wt., husk wt., stem wt. and leaf blade wt. under irrigated conditions.
Table 5. Treatment effect of the cultivar and N on the grain wt., grain number m−2, 100-grain wt., stover wt., aboveground biomass wt., grain No. m−2, cob wt., ear wt., husk wt., stem wt. and leaf blade wt. under irrigated conditions.
Treatment/CultivarGrain wt.(g m−2)Grains No. m−2100-Grain wt.(g m−2)Stover wt. (g m−2)Biomass(g m−2)Cob wt.(g m−2)Husk wt.(g m−2)Stem wt.(g m−2)HILeaf Blade wt.(g m−2)
ZMS 606717.16 a2798.67 a27.53 c266.62 a983.98 a117.2138.3947.680.72 a84.79
PHB 30B19611.70 b2341.10 b29.33 b259.05 a870.75122.4930.1643.00 a0.80 a80.84
PHB 30B50662.00 b2254.94 b31.39 a272.39 a934.40126.4126.7349.850.71 a91.54
Significancensns**nsnsns*nsnsns
LSD 5%104.39341.361.7242.191127.7917.2810.8212.770.0423.53
CV %15.3123.004.4718.7021.824.314.5019.002.7023.2
N rate
N1565.69 b2132.55 b29.04 a250.98 a816.67105.0637.5544.000.69 a81.53
N2711.47 a2562.81 a30.26 a274.09 a985.55129.0328.9749.80.72 a89.75
N3713.71 a2699.35 a28.96 a273.20 a986.91132.03 a28.7646.73 a0.72 a85.91
Significance****nsns**nsnsnsns
LSD 5%104.39341.361.7242.19127.7917.7310.8212.970.0423.53
CV %15.3013.504.4715.4413.4014.1433.1727.005.1226.73
Interaction (V ∗ N)
Significancensnsnsnsnsnsnsnsnsns
Means sharing the same letter in the table do not differ statistically at p < 0.05; LSD = Least Mean Differences; * = significant at the 5% level; ** = highly significant at 5%; NS = Nonsignificant; sqm = square meter; wt. = weight.
Table 6. APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the days after planting (DAP) to anthesis for maize cultivars under rainfed conditions.
Table 6. APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the days after planting (DAP) to anthesis for maize cultivars under rainfed conditions.
APSIM-MaizeDSSAT-CERES-Maize
CultivarRMSEnd-StatRMSEMER2RMSEnd-StatRMSEMER2
ZMS 6060.870.890.58−0.330.752.900.421.91−1.670.00
PHB 30G192.740.531.83−0.670.061.220.930.82−0.670.94
PHB 30B504.100.352.71−1.330.116.890.284.55−3.330.96
Table 7. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the days after planting (DAP) to maturity for maize cultivars under rainfed conditions.
Table 7. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the days after planting (DAP) to maturity for maize cultivars under rainfed conditions.
APSIM-MaizeDSSAT-CERES-Maize
CultivarRMSEnd-StatRMSEMER2RMSEnd-StatRMSEMER2
ZMS 6063.300.634.123.000.252.320.732.892.330.82
PHB 30G192.770.653.462.000.183.610.714.513.670.60
PHB 30B501.730.862.160.670.841.460.911.83−0.670.96
Table 8. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the biomass yield for the maize cultivars under rainfed conditions.
Table 8. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the biomass yield for the maize cultivars under rainfed conditions.
APSIM-MaizeDSSAT-CERES-Maize
CultivarRMSEnd-statRMSEMER2RMSEnd-statRMSEMER2
ZMS 60653.310.255.215.090.1268.080.892.862.220.87
PHB 30G1943.300.304.564.410.1963.740.893.032.400.89
PHB 30B5043.160.434.604.180.1555.310.922.651.900.87
Table 9. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the grain yield for the maize cultivars under rainfed conditions.
Table 9. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the grain yield for the maize cultivars under rainfed conditions.
APSIM-MaizeDSSAT-CERES-Maize
CultivarRMSEnd-statRMSEMER2RMSEnd-statRMSEMER2
ZMS 60617.410.661.281.030.609.160.910.670.130.84
PHB 30G1917.340.661.271.000.589.270.930.68−0.040.76
PHB 30B5020.130.441.56−0.130.6214.310.911.110.080.69
Table 10. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the grain unit wt. for the maize cultivars under rainfed conditions.
Table 10. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the grain unit wt. for the maize cultivars under rainfed conditions.
APSIM-MaizeDSSAT-CERES-Maize
CultivarRMSEnd-statRMSEMER2RMSEnd-statRMSEMER2
ZMS 60615.050.610.050.020.1919.250.550.06−0.020.57
PHB 30G1912.390.680.040.000.2222.100.520.08−0.030.33
PHB 30B5015.020.720.060.000.5114.280.910.06−0.020.79
Table 11. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the grain number m2 for the maize cultivars under rainfed conditions.
Table 11. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the grain number m2 for the maize cultivars under rainfed conditions.
APSIM-MaizeDSSAT-CERES-Maize
CultivarRMSEnd-statRMSEMER2RMSEnd-statRMSEMER2
ZMS 60610.210.42230.5582.890.0014.660.35330.92176.890.00
PHB 30G1912.810.37283.20188.560.0016.410.30362.67194.000.00
PHB 30B5017.520.22308.84135.440.3224.840.30437.91304.780.31
Table 12. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for mLAI for the maize cultivars under rainfed conditions.
Table 12. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for mLAI for the maize cultivars under rainfed conditions.
APSIM-MaizeDSSAT-CERES-Maize
CultivarRMSEnd-statRMSEMER2RMSEnd-statRMSEMER2
ZMS 6067.270.480.24−0.170.1051.540.141.73−1.730.15
PHB 30G1918.150.400.69−0.620.0347.830.201.83−1.840.06
PHB 30B5024.300.420.94−0.700.0662.910.292.42−2.370.60
Table 13. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the root soil water content under rainfed conditions.
Table 13. The APSIM-Maize and DSSAT-CERES-Maize model calibration statistics for the root soil water content under rainfed conditions.
APSIM-MaizeDSSAT-CERES-Maize
Soil Layer (cm3 cm−3)RMSEnRMSERMSEnRMSE
soil layer 116.080.0428.950.09
soil layer 28.370.0319.530.07
soil layer 37.690.0310.250.04
soil layer 47.340.039.790.04
soil layer 59.560.0411.950.05
soil layer 611.740.0511.740.05
soil layer 711.690.0511.690.05
soil layer 814.120.0621.180.09
Table 14. Results of the validation statistics for the APSIM-Maize model under irrigated conditions.
Table 14. Results of the validation statistics for the APSIM-Maize model under irrigated conditions.
AnthesisMaturityGrain YieldBiomass YieldGrain SizeGrain No m−2mLAI
NRMSE16.0012.7720.1873.986.2518.5116.75
RMSE16.6919.541.346.870.02456.210.57
MAE20.6717.331.016.770.01354.440.50
CRM0.200.11−0.11−0.73−0.04−0.090.01
Pearson0.63NA−0.010.200.890.32−0.33
d-stat0.04NA0.480.230.750.590.11
Table 15. Results of the validation statistics for the CERES-Maize model under irrigated conditions.
Table 15. Results of the validation statistics for the CERES-Maize model under irrigated conditions.
AnthesisMaturityGrain YieldBiomass YieldGrain SizeGrain No m−2mLAI
NRMSE19.9411.3445.0871.6421.1426.9149.56
RMSE20.8017.362.996.650.06663.291.67
MAE20.6719.332.796.510.05568.891.59
CRM0.200.13−0.42−0.70−0.18−0.220.47
Pearson0.63NA0.210.270.850.320.38
d-stat0.04NA0.350.240.550.520.36
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chisanga, C.B.; Phiri, E.; Chinene, V.R.N. Evaluating APSIM-and-DSSAT-CERES-Maize Models under Rainfed Conditions Using Zambian Rainfed Maize Cultivars. Nitrogen 2021, 2, 392-414. https://0-doi-org.brum.beds.ac.uk/10.3390/nitrogen2040027

AMA Style

Chisanga CB, Phiri E, Chinene VRN. Evaluating APSIM-and-DSSAT-CERES-Maize Models under Rainfed Conditions Using Zambian Rainfed Maize Cultivars. Nitrogen. 2021; 2(4):392-414. https://0-doi-org.brum.beds.ac.uk/10.3390/nitrogen2040027

Chicago/Turabian Style

Chisanga, Charles B., Elijah Phiri, and Vernon R. N. Chinene. 2021. "Evaluating APSIM-and-DSSAT-CERES-Maize Models under Rainfed Conditions Using Zambian Rainfed Maize Cultivars" Nitrogen 2, no. 4: 392-414. https://0-doi-org.brum.beds.ac.uk/10.3390/nitrogen2040027

Article Metrics

Back to TopTop