Geospatial Analysis of Soil Properties and Their Effects on Maize and Cassava Production in Ohaji/Egbema Imo State, Nigeria ()
1. Introduction
It is envisaged that Africa’s soil resource will continue to deteriorate, probably as a result of changes in climate, land use and human activities in general [1] . The main problems for soils in the African continent are contamination and irreversible losses due to increasing erosion. Soil erosion, in particular, is regarded as one of the major and most widespread forms of land degradation, and as such, poses severe limitations to sustainable agricultural land use [1] .
However despite Nigerian government campaigns and slogans, farm production has not kept pace with food demand. Most food crops produced in the country come from the efforts of small-scale farmers who depend largely on traditional farming systems for their agricultural inputs [2] . The recurrent food crisis in Nigeria is partly due to high rate of population growth over the food production level and erratic amounts of food crops produced from year to year. This however, can be attributed to high susceptibility of the country to serious environmental hazards from low rainfall, extreme temperature, acid rain, gas flaring, oil spillage, poor soil quality, deforestation, continuous cropping and unhindered desert encroachment [3] . Arable crops such as cassava, maize, yam and cocoyam are the chief sources of dietary energy for the majority of the people living in the lowland tropics, and many of the sub-humid tropics of West and Central Africa [4] .
Therefore, their production and utilization must be given prime attention in food policy. Even though farmers have not yet attained the desired technical efficiency in their production as a result of weak access to external inputs such as fertilizers and herbicides [5] .
[6] pointed out that optimum soil nutrients are sine qua non for sustainable agriculture for food and nutrition security. [7] , in their study concluded that rapid decline in crop yield was associated with continuous cropping in Netherlands. Also [8] noted that the low-available water holding capacity of the soil results in poor crop growth in the tropics. In another study, [9] believes that rising temperature affects soil moisture, which in turn could affect soil fertility. Supporting the above remark, [10] observed that the major agro-climatic constraints on agricultural production are related to insufficient, excessive or irregular moisture supply, which in turn will affect the length of growing period (LGP) of crops [11] . While assessing the characteristics of soils in the Guinea Savanna zone of Nigeria it is reported that bulk density increased down the profile, in accordance with increase in clay content.
According to [12] , soil physical properties play a key role in soil sustainability and crop production; that they determine how easily plant roots can grow to access soil nutrients, and how easily water can flow through the soil to deliver these nutrients.
[13] reported that organic matter was generally highest in the top most horizon of Imo River Basin while it was consequently highest in the topsoil because of the continuous addition of organic matter to the soil surface by plant residues. [14] , attributed this decreasing trend of organic matter with depth primarily to the effect of nutrient bio-cycling.
According to [15] ECEC is used to estimate the potential fertility of the soil, the possible response to fertilizer application and as a rough guide to the types of clay minerals present. Generally ECEC was reported to be low in Imo river basin [16] [17] [18] . Percentage Base Saturation (% Bsat) is another good indicator of soil fertility [15] . According to [14] , the low nitrogen in most of the pedons is a common phenomenon in the soils of Southeastern Nigeria and is as a result of the high nitrogen losses sustained in these soils through leaching of nitrates, as well as the rapid mineralization of organic matter under the isohyperthermic soil temperature regime.
Ultimately, physical factors such as climate, topography, land cover and specific soil characteristics have important effects on the processes of soil erosion and soil formation, which in turn affect crop production. Therefore, this study was investigating physic-chemical soil properties known as environmental factors among others, directly and indirectly influencing production of maize and cassava in Ohaji/Egbema in Imo State of southeastern Nigeria.
2. Materials and Methods
2.1. Location
The study area, Ohaji/Egbema in Imo State of Southeastern Nigeria is located on latitude 5˚15'0"N and 5˚36'0"N and longitude 6˚35'0"E and 7˚0'0"E. It is an oil-rich Local Government Area of Imo State with Its headquarters at Mmahu-Egbema. The local council has an area of 890 km2 and lies in the south-western part of Imo State. Besides Oil and Gas exploration, Ohaji/Egbema is in the rich agricultural zone of Imo State. Hence, agriculture is the mainstay of the people with faming, fishing, palm oil processing, hunting, and animal husbandry vigorously practiced. The area is naturally humid as the entire southeastern Nigeria, invariably, soils from the area humid which is good for crop production.
Moreover, Ohaji/Egbema is typically rainforest vegetation (Figure 1). Arable crops grow very well in the area which includes cassava, yam, maize, cocoa yam, etc. Therefore, soil moisture (helping chemical and biological activities of the soil) was expected to be high, while soil temperature (affecting physical and chemical processes in the soil) was expected to be low in the area. Therefore, for the purpose of this study, it was important to test characteristics of the soils from the area to know their potentials for production of the arable crops.
Figure 1. Rainforest vegetation at Mmahu, in Ohaji/Egbema in Imo State. Source: Field work, 2019.
2.2. Primary Data Collection
A free survey was employed for the study after the ten locations in the area was purposely selected for representation of the soil sampling. The free survey was guided by size and extent of farming activities in the area. Two types of sampling were conducted. Auger sampling was conducted at 0 - 20 and 20 - 40 cm depths for arable crop production purposes while soil profile sampling was done for in-depth characterization and scientific classification of soils. Soil profile was dug, described and sampled using standard procedures as recommended by [19] . Soil samples were collected based on horizon differentiation (see Figure 2 and Figure 3). Sampling started from the deepest horizon upwards in each soil profile. Soil profiles were geo-referenced using hand held Global Positioning System (GPS) Receiver. Core soil samples were collected for bulk density determinations. Soil samples were bagged and taken to the laboratory for analysis.
2.3. Secondary Data Collection
The secondary data used was annual crop yield data on maize and cassava (in tons/hectare), for Imo State from 1988 to 2017 (30 years). The data was collected from Agriculture Development Programme (ADP) office in Owerri, Imo State. The crop yields data (see Table 1), for the entire state was used to represent crop yields for the study area, which was correlated with physicochemical soil properties to determine their relationship. The values of Maize and cassava in metric tons/hectare were obtained from the archrivals of the Agriculture Development Program (ADP), Owerri office, Imo State.
Figure 2. Soil samples collection at Obokafia Ohaji/Egbema, Imo State.
Figure 3. Field officer collecting soil samples at Mmahu, Ohaji/Egbema, Imo State.
3. Results and Discussions
3.1. Results
The primary and secondary data collected for this research were subjected to series of analytical methods including laboratory soil analyses, geospatial data analysis, and statistics models analysis.
3.1.1. Laboratory Studies
Soil samples were air-dried and sieved using 2-mm sieve.
Particle size distribution was determined by hydrometer method [20] while
Table 1. Crop yield for Imo State in Metric tons (1988-2017).
Source: ADP Owerri, Imo State.
bulk density was measure by core procedure [21] . Results from particle size analysis (sand, silt and clay) values were used to obtained textural class using textural triangle. In the process of the analysis, Robinson’s pipette method (Pipette) and Bouyoucos hydrometer method were used.
Bulk density values were used to calculate total porosity of soils given a relationship between bulk density and particle density [22] .
where BD = determined bulk density;
PD = particle density assured to be 2.65 mg∙m−3 (2.65 g/cm−3).
Gravitaional moisture content (
) was measured using the procedure as outlined in [23] .
where
= gravimetric moisture content;
Ws = weight of wet soil sample;
Ds = weight of dry soil sample.
Soil pH water and pH KCL were determined electronically in 1:2.5 soil solutions:water-ratio [24] .
Soil organic Carbon was measured by wet digestion using the procedure outlined in [25] .
Soil organic matter was calculated by multiplying organic carbon value by 1.724.
Exchangeable basic cations were extracted using ammonium acetate at pH7. Therefore, exchangeable calcium and magnesium were determined using ethylene diamine-tetraascetic acid (EDTA) titration, and exchangeable potassium and sodium were determined by flame photometry [25] .
Exchangeable acidity (exchangeable hydrogen and aluminum) were measured by apparent titration [26] .
Cation exchange capacity (ACEC) was measured at pH of 7.0 (neutral) [25] .
Total nitrogen (TN) was determined by micro-kjedahl apparatus [27] . Available phosphorous was estimated by Bray 2 method according to procedure of [28] .
Base saturation was computed as a sum of exchangeable basic cations (Ca, Mg, K, Na) divided by Cation Exchange Capacity, multiplied by 100% [25] .
3.1.2. Date Analytical Techniques
In probability theory and statistics, coefficient of variation is a standardize measure of dispersion of a probability distribution or frequency distribution. It is a statistical measure of the distribution of data points in a data series around the mean, and it represents the ratio of the standard deviation to the mean. While variation is a measure of how far from the mean data set varies, coefficient of variation compares the standard deviation to the mean of each sample.
In this study, soil data from auger samples were subjected to analysis of variance (ANOVA). Therefore, variation in properties among soil horizons were obtained using coefficient of variation. Also, the correlation coefficient (r) and the coefficient of determination for simple linear regression (r2) were calculated to determine the degree of association or relationship among some physico-chemical properties of the soils from selected profiles around the study area. Similarly, soil parameters (independent variables) were correlated with crop yield (dependent variables) toestablish relationship between soil properties and production of maize and cassava in the area.
3.1.3. Geospatial Data Analysis
From the coordinates captured in the fields by the global positioning system device at the soil samples location points, Arc-GIS software was used in the GIS laboratory to produce GIS vector “soil sample collection location map” of the study area. Also digital elevation model DEM of the study area was produced using the coordinates (Figure 4 and Figure 5).
Figure 4. Area study map (Ohajie/Egbema) showing soil sampling locations.
Figure 5. Study area map showing soil sampling points and elevations.
3.1.4. Physical Properties of Ohaji/Egbema Soils in Imo State
The physical properties of soil of Ohaji/Egbema Local Government Area of Imo State are presented in Table 2. Statistical formula for standard error was used for the calculation regarding the soil physical properties (variables) shown in Table 2.The formula is given as;
.
Results therefore, showed that the mean value of CS, FS, TS, Si and Cl from 0 - 190 cm depth were 534, 132, 559.8, 26.4 and 307.6 g/kg respectively. There were also variations in the textural classes of the soil from A to Bg3 horizons. The bulk density increased from A to Bg3 horizon due to higher organic matter content (Table 2). There were low variability in soil bulk density, coarse sand, total porosity and AWC and medium variability in clay, fraction, field capacity, and permanent wilting point. High variability occurred in fine sand, total sand, silt and moisture content.
3.1.5. Chemical Properties of Ohaji/Egbema Soils in Imo State
The chemical properties of soils in Ohaji/Egbema Local Government Area Imo State are shown in Table 3. The soils are strongly acidic with low variability. The mean values of organic carbon, organic matter, and total nitrogen were 11.74, 20.1 and 1.20 g∙kg respectively. These values were low using [19] rating and the highest values were recorded at the epipedon (0 - 18 cm depth). This could be as a result of litter fall and higher biodiversity within this layer of soil.
Exchangeable Ca, Mg, K and Na were low. Effective cation exchange capacity and base saturation were low according to [29] . Available phosphorus was low with mean value of 14.68 mg/kg. There was low variability in soil pH. Exchangeable K, and Na, and base saturation while medium variation occurred in exchangeable Mg, total exchangeable bases, effective cation exchange capacity, ACEC, and available P. Organic carbon, organic matter, total nitrogen, exchangeable Ca and H had high variability.
Table 2. Physical properties of Ohaji/Egbema soils in Imo State.
Table 3. Chemical properties of Ohaji/Egbema soils in Imo State.
3.1.6. Relationship between Soil Properties and Crop Yield in the Study Area
Results of the relationship between soil properties and crop yield are presented in Table 4. Results showed that there was significant positive relationship between moisture content and cassava yield (r = 0.6743) and significant negative correlation between moisture content and maize yield (r = −0.6490). Exchangeable Al significantly correlated positively with cassava yield and negatively with maize yield. Bulk density, CS, exchangeable K and TS significantly correlated negatively with cassava yield and positively with maize yield. There was significant positive relationship between Cl, effective cation exchange capacity, exchangeable Mg, Si, total porosity and total N with cassava yield and negative relationship with maize yields. The results showed that soil properties influence the yield of cassava and maize in the area and Imo State in general. Therefore improving the fertility status of soil and improving soil quality is a panacea for better cassava and maize yield. Similarly some soil properties that favour cassava yield may not favour maize yield and hence understanding these factors will help in policy managing agricultural lands for crop production.
The “r” means the correlation coefficient value for the relationships between soil properties and crop yields. Hence, soil data from auger sample were statistically correlated with crop yields using coefficient of variation at both the 0.01 and 0.05 probability levels.
3.2. Discussion
Generally, this research supports previous studies, which revealed that physical and chemical properties of soils influence crop yields in different regions of the world. Soil physical properties determine how easily plant roots can grow to access soil nutrients, and how easily water can flow through the soil to deliver these nutrients. Supporting this, studies in Southern Cameroon by [30] [31] ; in Nigeria by [32] and [33] ; in the US by [34] ; show that soil physical properties affected maize grain yields [35] , indicated that inadequate moistures in the soil
Table 4. Relationship between soil factors and yields of maize and cassava in the study area.
* and ** = sig at 0.05 and 0.01 probability levels respectively. Source: Field work (2019).
reduce maize production.
[36] reported, very low to moderate silt/clay ratios (0.04 to 0.029) that generally decreased with depth while working on soils at Isienyi Ibeku, which is within Ikwuano soil zone in Abia State. This study however reported high to very high silt/clay ratios (1:2.99) in Ohaji/Egbema in general, and the ratios increased with depth, which was in agreement with the research by [37] , which reported that the bulk density of sub surface horizons were usually significantly greater than those of surface horizons. The value of bulk density in the study area at 0 - 20 cm depth was 1.33 g/cm3 while the value at 20 - 40 cm depth was 1.39 g/cm3.
This supports that BD increases with increase in depth down the profile. This research further observed that there is irregular silt to sand ratio in the entire Southeast Nigeria.
[38] reported high variations in silt fraction and total sand and fine sand in southeast Nigeria. This study showed high and irregular variations these soil properties which supports the previous study by [38] .
This research showed that pH in water {pH(H2O)} was higher than pH in KCl {pH(KCl)}. For example, mean of pH in water was 5.12 and pH(KCl), 4.24, and generally, the soil of the area was highly acidic with low pH values. Hence [14] , reported that such reaction is characteristic of soils in southeastern Nigeria and it is the result of the acidic nature of the parent rocks, coupled with the influence of the leached profile under high annual rainfall condition.
[39] and [40] reported that total nitrogen is usually highest in the topsoil, decreasing down the profile. In support, this research reported high values of total nitrogen at the topsoil in Ohaji/Egbema, decreasing down the profile, from 2 g/kg at A horizon to 0.3 g/kg at Bg3. These results generally showed that the nitrogen values are decreasing down the profile.
[13] reported, that organic matter was generally highest in the top most horizon of Imo River Basin and highest in the topsoil because of the continuous addition of organic matter to the soil surface by plant residues. This study, however, is in consonance with the reports of the previous studies by indicating that organic matter contents in the study area are generally higher in the topsoil. In the area, OM decreased from total value of 124.46 kg in the topsoil to 67.56 kg in the sub soil.
ECEC was reported to be low in Imo river basin [16] [17] [18] . The present study supports previous studies by reporting also a low ECEC in the area with total value of 11.63 cmol/kg in the top soil and 10.63 cmol/kg in the subsoil. However, the values of ECEC decrease with increase in depth, which is a common condition with most soil properties in the southeast Nigeria as reported in other studies.
4. Conclusions
The study established that there are variations in physical and chemical properties of soils in Ohaji/Egbema area of Imo State, Nigeria and the variability in the soil properties which affect crops growing in the area. It further indicated that the soil properties negatively influenced maize and cassava yields, particularly in the area, and it is susceptible to future decrease. From the study also, it is observed that Ohaji/Egbema soil contains less plant nutrient elements that will support crop production. Therefore, the chemical properties of soils in the area are low according to [19] and [29] soil fertility rating. Hence, there is therefore need to boost the productivity of soils in the area and entire Imo State for maximum crop production, in order to achieve sustainable food security.
However, apart from soil characteristics, there are other factors that can cause variations in cassava and maize yield in the area. These factors are, not limited to;
・ Farm machines or technologies in use;
・ Climate change and variability;
・ Level of fertilizer and insecticide application;
・ Adaptive practices to rainfall variation and climate change;
・ Educational level and farming experience of the farming population.
Recommendations
The data for this research allow suggesting recommendations to achieve increased productivity and resources use efficiency for these major food crops in the area.
Therefore, the government will need to design policies to counteract the effects of climate change in relation to agriculture. This will allow farmers to adjust their management strategies like early planting and to minimize the risk of crop failure.
Also, since this study ascertained that some soil properties negatively influence production of maize and cassava in the area, therefore, improving the fertility status of soil and improving soil quality will be a panacea for better maize and cassava yield. Similarly, some soil properties that favour maize yield may not favour cassava yield and thus, understanding these factors will help in policy managing agricultural lands for crop production. Therefore, farmers in this region can reduce the potential damage through adaptation such as the practice of organic agriculture, which is one of the most important strategies for adaptation to climate change and variability.
Furthermore, using different crop varieties with varying planting dates is significant at farm level. Since different crops are affected differently by climatic changes, farmers in the area can grow varieties of crops on the same farmland or on a different lands and this will reduce the risk of complete crop failure.
They are as well recommended to practice crop adaptation using possible existing genetic diversity and biotechnology.
Finally, farmers in the region should resort to some more management practices such as contour ploughing, ridges and cross-bars between ridges, which will definitely reduce erosion and runoff thereby preventing leaching of soil nutrients.
Acknowledgements
Nigerian Tertiary Education Trust Fund (TETFUND) is solely acknowledged for sponsoring this research through Institutional Base Research (IBR) programme. Also, Imo State University Owerri (our Institution), is highly acknowledged for attracting the research funding. Franklyn Okoro and two other persons are acknowledged for taking the research team round the study area during field studies.