The purpose of this study was to simulate the potential yields of maize under optimal growth conditions across Ghana using simulation modeling based on plant physiological processes as affected by weather and crop characteristics, and to determine the potential yield gap with observed yields reported by farmers, as well as the factors that influence that gap. A total of 676 observed data points were used to validate a calibrated a Light Interception and Utilization (LINTUL-1) crop model. The validated model was used to simulate the potential yields at 2374 randomly sampled points from existing cropland locations to create a spatial overview of the yield levels through mapping techniques.
This analysis was driven by the following questions and objectives:
The Study Area is Ghana. The map below shows maize field observations obtained from peer-review papers in green dots, 10 random sample points per district based on land cover in red dots (the landcover layer MODIS MCD12Q1 helped to sample from crop lands), district boundaries, and agroecological zones in Ghana.
The LINTUL model for maize was initially calibrated based on existing research findings and expert knowledge. The parameters are as shown in the image.
The calibrated model was used to simulate the potential maize yields at all the random sampled point locations (2374 points), and mapped using Ordinary Kriging geostatistical method as shown in the image.
The Ordinary Kriging results was further plotted to visualize the distribution by agroecological zones as shown in the figure below. Guinea savanna has the highest potential whilst Moist evergreen has the lowest.
Data on observed maize yield and production for year 2017 was downloaded from the Ghana open data website. This data was geocoded because it had no coordinates assigned. As there were fewer districts than the shapefile used to perform geocoding, the results had some empty districts shown in the figure below. Yield gaps were calculated based on the observed data.
The formula for calculating potential yield gap is: Potential Yield Gap = Avg. Pot. Yield Per District - Obs. Yield Per District. The yield gap was calculated for each district and mapped as shown in the maps below. The first map from the left is the Observed Yields, the map in the middle is the Potential Yield, and the third map on the far right is the Potential Yield Gap. It can be seen that there are higher gaps (6-8 t/ha) in the northern part of Ghana than the south.
To visualize the yield gaps for the districts that have empty data, Ordinary Kriging geostatistical method was used to predict and visualize the yield gaps covering the entire country. The statistical results show that there is a spatial autocorrelation range of 117,172m, a nugget (variation at very small distances) of 0.22, and a sill (maximum variation) of 0.37. The wide range implies that the potential yield gap values are similar in a very large geographical area which is evident in the smoothness of the output map shown
Based on the potential yield per district, the potential production per district was calculated. The formula for calculating potential production per district is: Potential Production Per District = Harvested Area per District (Ha) * Avg. Pot. Yield Per District (T/Ha)
To understand the driving factors for the potential maize yield gaps, Random Forest machine learning model was used to explain the variability (68%) of the potential yield gaps and to determine the important variables which were used in a regression model to estimate their effects on the yield gaps. The figure below shows the formula used and the variable importance plot of the Random Forest model. The selected variables used in regression modelling contributes to 79.1% of the explained variability.
The figure below shows the regression analysis results. The red arrows indicate variables that have a positive influence on the yield gaps whilst the green arrows show variables that have a negative influence on the yield gaps. Based on the results, it could be estimated that a 1 unit increase in Organic Nitrogen may reduce the potential yield gaps by about 1.4 t/ha.