Open access peer-reviewed chapter

Multivariate Assessment of California Bearing Ratio with Contrasted Geotechnical Properties of Soils in Ilorin-Lokoja Highway

Written By

Attah Fakeye, Olusegun Ige and Olufemi Ogunsanwo

Submitted: July 8th, 2020 Reviewed: July 30th, 2020 Published: April 28th, 2021

DOI: 10.5772/intechopen.93523

Chapter metrics overview

277 Chapter Downloads

View Full Metrics


California Bearing Ratio (CBR) is an important parameter used in designing pavement layers in road construction but testing this parameter requires time, labor, and huge cost. The study therefore applies multivariate approach to evaluate CBR based on contrasted geotechnical parameters along Ilorin-Lokoja highway. The results obtained showed that the migmatite-gneiss-derived soils are slightly more fines (< 0.075 mm; 7.4–59.6%), more plastic (PI; 1.6–39%), and have low strength (MDD = 1.8 mg/m3; CBR = 29.0%) than the metasediments (11–57.7%, 2.0–30%, 1.6 mg/m3, 23.6%) and older granite soils (8.2–32.7%, 2.6–13.4%, 1.7 mg/m3, 27.8%), respectively. The principal component analysis (PCA) revealed three major components (eigenvalues >1) which accounted for 83.8% of the total variance at the rate of 33.4, 14.7, and 11.4%. Major contributing variables for the components were fines (R = 0.87), plasticity index (R = 0.7), and coarse sand (R = 0.67%). Spatial distribution of these groups established interplay of sediment-gradation and moisture-connection evident in hierarchical cluster analysis that revealed patterns of homogeneity and soil relationships. Regression analysis established five models from predictor variables such as fines, activity, free swell, liquid and plastic limits, weighted plasticity index, optimum moisture content, and maximum dry density with the coefficient of determination (R2 = 0.33) and root mean square error (RMSE) of 7.80.


  • multivariate
  • principal component analysis
  • regression
  • hierarchical analysis
  • geotechnical properties

1. Introduction

Identification and quantitative characterization of soils are of dire importance in geotechnical assessment despite the difficulties experienced using conventional approach. Index properties are important parameters in the analysis of geotechnical engineering problems, particularly to estimate strength of the soil material. Conversely, laboratory test takes 2–4 days to measure compaction and California Bearing Ratio (CBR) values for pavement design. As a result, they are expensive and time-consuming. Also due to lack of specialized personnel, these tests are oftentimes avoided in many soil investigation programs. Thus, the need to incorporate statistical approach in predicting soil properties becomes inevitable.

Several authors have applied this approach in relating and predicting soil properties. One to one relationship was presented among soil properties [1] such as liquid limit (LL), plastic limit (PL), plasticity index (PI), optimum moisture content (OMC), and maximum dry density (MDD). Furthermore, Carter and Bentley explained that soil type, density, moisture content play an important role in soil relationship [2] and correlated soil expansion index and plasticity index, fine fraction and weighted plasticity index (i.e., product of PI and percentage passing 0.425 mm). Apart from index properties, some researchers like Owoseni et al. [3] and Yildrin and Gunaydin [4] observed that California Bearing Ratio depends on other factors such as type of soils, permeability of soil, maximum dry density and optimum moisture content. To correct overlapping problem and uncertainty in prediction, Yitagesu et al. applied multiple regressions to improve the ability of predicting soil properties, and better model the extent of their relationship [5].

This paper attempts to identify geotechnical characteristics of soils developed on different rocks and establish relationships among various properties in order to estimate soil strength capability in three lithological units. Multivariate approach using principal component analysis (PCA) and hierarchical classification methods are used to identify patterns, detect and classify new parameters into groups; and further propose regression models to determine CBR values in view of huge cost and labor.


2. Methodology

Soil samples (130 samples) were collected along the Ilorin-Lokoja highway (>300 km length) which spans across latitude 7°25′N-8°40′N and longitude 4°30′E-6°45′E. Simultaneously, Garmin GPS was used to record coordinates of sample locations. The topography ranges from a relatively flat to hilly, undulating terrain with elevation ranging between 100 and 700 m above sea level. The highway is overlying the Precambrian Basement rock of South Western Nigeria (Figure 1) and cut across three geologic units: the migmatite-gneiss complex (denoted by PCB), the metasediments/volcanic series (PCM) and the older granite series (PCG) [6].

Figure 1.

Geology of Nigeria showing the study highway overlain by sampling points (red dots).

Majority of the rock is the migmatite—gneiss essentially made up of migmatite and banded gneiss. Others are flaggy quartzite with biotite gneiss, undifferentiated schist, porphyritic granite (porphyroblastic), and medium-coarse grained biotite and hornblende granite. Temperature ranges from 25 to 35°C. Climate is dry to wet, with a mean annual rainfall of 1200 mm. Due to heavy rainfall, considerable moisture change occur in the soils which dries up at prolonged dry season. This induces soils susceptibility to volume changes.

2.1 Laboratory analysis

Geotechnical tests were carried out on air-dried (35–40°C) soil samples at the Soil Geotechnical Laboratory of Nigerian Building and Road Research Institute following the British Standard [7] Part 2: Clause 9.2, 4.5, 5.3, 5.4, and Part 4, Clause 3.3 and 3.4 methods. The soil engineering parameters obtained include natural moisture content, Atterberg limits, particle size distribution, free swell, compaction test and California Bearing Ratio.

2.2 Statistical analysis

In SPSS statistical software, 20 soil parameters were explored and their relationships examined. Data transformation was applied to ensure equal influence on the model thus, fulfilling the linear model assumptions. The strength and relationship trends on the dataset were examined from Pearson correlation matrix and quantitative measures of linear associations determined. Principal component analysis (PCA) approach was incorporated to reduce the data with many variables, identify clusters, and transform the soil variables into new uncorrelated variables that preserve most of the information [8]. Components with eigenvalues >1 were retained and subjected to varimax rotation to maximize correlation between the factors and measured variables. Thereafter, Agglomerative Hierarchical Cluster (AHC) analysis was computed to identify analogous behavior among different soil characteristics and soil individuals using Ward’s method and squared Euclidean distance as a measure of similarity between soils [9].


3. Results and discussions

The statistical summary of the laboratory test is shown in Table 1. The soils exhibited wide variations of data clustering around the mean value (1.08–88.6%) and high coefficients of variation (1.7–147%). The median of some parameters was lower than the mean value, indicating a low effect of abnormality on sampling values.

3.1 Particle size characteristics

In the migmatite-gneiss derived soils (PCB), gravel and coarse sand varied with coefficient of variation (CV) from 18.3 to 100% (<23.5%), medium to fine sand was between 8.0 and 86% (>32%) while the percentage of silt and clay were 3.2–50.8% (52.4%) and 0.9–34.6% (58.8%) respectively. However, the percentage of fines (<0.075 mm) ranged between 7.4 and 59.6% (48.0%). This proportion of fines is similar to those reported by Ige et al. [10]. In the metasediment derived soils (PCM), content of gravel and coarse sand were higher from 43 to 100% (<20.6%), while medium to fine sand was between 15.2 and 84.2% with 36% CV. Similarly, the percentage of silt and clay ranged between 6.1–44.3% (41.2%) and 1.8–28.7% (46.3%). The proportion of fines (11.7–57.7%; 36.6% CV) is relatively as high as the PCB origin. Similarly, the older granite rock (PCG) exhibited a wide range of gradation with gravel and coarse sand ranging between 19.4 and 100% (<33.3%). Medium to fine sand content was lower (14–92%) (CV = 28–34.9%) while percentage of silt and clay varied between 2.4–24.5% and 1.2–18.8% (CV = 48.6–60.8%), respectively. The amount of fines (8.2–32.7%) and CV (34.6%) are very low in this area.

On one hand, this granularity is similar to the work of Nwaiwu et al. [11] where the lateritic soils are enriched with gravel and sands ranging between 28.2–40% and 42.2–48% resp. However, the high percentage passing through No. 200 (0.075 mm) BS sieve suggests the soil is predominantly of fine materials and classified according to Unified Soils Classification System (USCS) system as clayey sand (SC), silty sand (SM) and silty, clayey sands (SC-SM). Other soil classes obtained include poorly graded sand with silt or clay (SP-SM, SP-SC), poorly graded gravel with clay or silty clayey gravel (GP-GC, GC-GM), silty gravel (GM), sandy lean or fat clay, (CL, CH), and sandy silt or elastic silt (ML, MH) that occurred in low percentage. Similarly, according to American Association of State Highways and Transportation Officials (AASHTO) system, the most dominating classes are A-2 and A-7 soils, hence rated as excellent to good and fair to poor materials for road use.

3.2 Consistency limits

A wide range of plasticity (Figure 2) characterized the inorganic silty clayey soils in the area. The liquid limit varied between 13.4 and69% with a lower range experienced in PCG derived soils (<32.4%), while plastic limit and plasticity index at the PCB unit ranged between 2.2–50% and 1.62–39%, with mean values of 28.4, 19.0, and 9.4%, respectively (Table 1). The Casagrande plasticity chart revealed majority of the soils from the migmatite-gneiss origin placed above the A-line, indicating that they are composed of inorganic clay material and exhibited low to medium plasticity, implying low to medium swelling and compressibility. The moderate plasticity suggests low to medium dry strength, which could easily crumble under load thus leading to pavement failure and possible erosion under climatic threat. The distribution of the soil samples on the chart portrayed the variability in soil plasticity characteristics.

Figure 2.

Casagrande chart of plasticity-liquid limit relationship.

Moreover, free swell (Fsw) varied from 1.8 to 28.4% in PCB, 2.8–19.6% in PCM and 3.2–4.67% in PCG with mean values ranging from 6.4, 7.1, and 4.04%, respectively; while soil activity with mean values oscillated between 0.09–8.5 (1.1), 0.17–7.85 (0.9), and 0.2–5.08 (1.2) within the 3 units. The weighted plasticity index (wPI) value ranged between 0.25–31.9% (5.1%), 0.92–25.3% (5.8%), and 0.87–7.3% (3.5%) with mean from the 3 units. In PCB soils, activity tends to be higher than normal (8.5), high weighted plasticity index (wPI), plasticity ratio (PIr), and swelling potential (SP) indicating that the soils are active. The result of natural moisture content (NMC) (2.6–22.7%) is fairly high, considering the time of sample collection. This indicates the soil potential for water retention, which is a property of fine-grained soils. The high water content also suggests the presence of high water table earlier reported by Adams et al. [12]. These observations correspond with Bayamack et al. [13]. The derived plasticity parameters (wPI, PIr, SP, and LLr) represent the effective contribution of the plasticity of fines to the performance of the entire soil materials, depending on the amount of fines.

3.3 Compaction and California Bearing Ratio

The maximum dry density (MDD) of the soils from PCB area (Table 1) increases with mean to 2.6 mg/m3 (1.77 mg/m3) at 25% (13%) optimum moisture contents (OMC). These values are higher than those obtained in metasediment (PCM) and older granite (PCG) units with 2.1 mg/m3 (1.6 mg/m3) MDD and 22.5% (15%) OMC. The low density-moisture relationship implies low strength instigated by loose soils that are susceptible to erosion. The interaction of the subgrade with water greatly reduces strength and therefore promotes continuous failure of the overlying pavement. Few examples of soil compaction curves (Figure 3) illustrate distinct peak of maximum dry density at optimum moisture content.

Figure 3.

Compaction curves of selected soil samples.

The CBR values at 95% OMC after 48 hours of immersion varied between 10 and 56.4% for PCB, 11 and 45% for PCM and 12.1 and 37.2% for PCG soils (Table 1). The mean values within the three lithological units varied between 28.8, 23.6, and 27.8%, respectively. For unsoaked condition, the CBR varied in a higher rate from 12.5 to 75.0% within the 3 units. The result showed a reduction in strength due to soaking suggesting a probable drastic reduction in strength by more than half during wet condition and the penetration resistance becomes reduced due to excessive moisture. These values are similar to those found along Ado Ekiti-Akure road (27–100%) by Adams and Adetoro [14]. The low mean CBR value (<30%) suggests that the soils may not withstand ground vibrations when vehicular load is applied and reinforces its susceptibility to erosion. Soil improvement measures are therefore, envisaged for the stability of soils for adequate strength.

3.4 Simple linear regression

High statistically significant correlation (R > 0.70) is recorded among 13 soil attributes pairs (Table 2) such as gravel, coarse sand (CS), medium sand (MS), fine sand (FS), silt, clay, fines, swelling potential (SP), free swell (Fsw), liquid limit (LL), plasticity index (PI), dry density (DD) and maximum dry density (MDD) which raises the issue of multi-collinearity. However, other parameters exhibit low correlations (R < 0.50) including sand, activity (Ac), plastic limit (PL), moisture content (MC) and optimum moisture content (OMC).

Migmatite-Gneiss (PCB) = 87Metasediment/Metavolcanic (PCM) = 33Older Granite (PCG) = 11

Table 1.

Statistical summary of soil properties.

Parameter abbreviations are described in the Abbreviations section

PropertiesR valuesPropertiesR values
FS and fines0.97Fsw and CBRu−0.62
Gravel and CS0.96Fsw and PI0.61
BD and MDD0.95LL and PI0.61
DD and MDD0.95NMC and clay0.60
Silt and clay0.94MS and silt0.60
Silt and fines0.92LL and MC0.59
BD and DD0.90PL and MC0.56
PI and wPI0.89PI and MC0.56
Clay and fines0.88FS and sand−0.56
Fsw and LL0.88NMC and LL0.56
FS and silt0.80MS and clay0.56
LL and PL0.80DD and OMC−0.55
FS and clay0.77MDD and OMC−0.54
CS and MS0.75Gravel and sand0.52
MC and OMC0.72fines and OMC0.52
Gravel and MS0.70FS and LL0.52
LL and CBRu−0.70FS and PL0.51
MS and FS0.68FS and OMC0.51
MS and fines0.68NMC and Fsw0.51
Fsw and wPI0.68Clay and Ac−0.51
Fsw and PL0.66fines and LL0.51
LL and wPI0.64NMC and fines0.50
PL and CBRu−0.62Sand and fines−0.50

Table 2.

Pearson significant correlation of soil properties.

This could be attributed to the presence of high fine fractions and potential influence from environmental factors. The result corroborates with the observations obtained on gneiss derived laterite in Central Cameroun [15] and reaffirms the views of the earlier scholars that geotechnical properties of laterites depends on the parent materials, climate, vegetation, topography and duration of the laterization phenomenon [16].

3.5 Multivariate analysis of soil properties

3.5.1 Principal component analysis

Among the multivariate analysis techniques, principal component analysis is the most frequently used because it is the starting point in data mining which aims at minimizing the dimensionality of the data. Seven principal components (PCs) were extracted with eigenvalues >1 which accounted for 83.8% of the total variance of data (Table 3).

Initial eigenvalues
PCsEigen values% of VarianceCum % of variance

Table 3.

Eigenvalues and proportions of variance explained by PCA.

However, the first five PCs accounted for >70% of variability in measured soil properties. While PC1 explained 33.4% of the total variance with fines as the major contributing variable (R = 0.87), PC2 accounted for additional 14.7% of the total variance with plasticity index (PI) as the second major contributing variable (R = 0.70). In PC3, 11.4% was accounted for, with coarse sand (CS) contributing more (R = 0.67). Other components accounted for <15% and as such were removed as they explained less variance than individual variable in the dataset [8].

Based on the communality estimates, the five factors explained more than 90% of variance in MDD, PI, LL, DD, BD, FS, CS, SP, fines and gravel; > 80% in wPI, PL, Fsw, Wr, LLr, NMC, sand, silt and clay; > 70% in MS, MC, CBRu, and Dr.; above 60% in PIr, OMC, and Ac; and 53% in CBRs (Table 4). According to Johnson and Wichern [17], a high communality suggests that a high proportion of the variability is explained by the factor with a higher preference over a low communality estimate. By implication, the factors fairly explained the variance in soaked CBR and as such required a regression model to predict the property. The values obtained are similar to those obtained by Shukla et al. [18].

Soil variablesPC1PC2PC3PC4PC5Communalities

Table 4.

Proportion of variance and communality estimates of soil variables.

The coefficient of linear correlation between the variables and their factors (Table 4) give a meaning to the principal components. The parameters are well represented and explained by the factorial axes on the correlation circle (Figure 4). This graph shows three groups of variables, suggesting the existence of correlation between them.

Figure 4.

Score plot and correlation circle obtained with PCA.

PC1 positively correlates (> 0.84) with SP, PI, LL, Fsw, wPI and MC, NMC, PL, CBRu (Table 4) and is termed plasticity parameters. PC2 demonstrated very high positive correlation with soil densities (MDD, DD, BD) (> 0.93) and negatively correlated with moisture contents (MC, OMC) (<−0.6) and is termed moisture-density or compaction parameters since the variables are important functions of soil moisture density. It also showed moderate positive loading from CBRs (0.37) resulting from significant correlation between MDD and OMC. Similarly, PC3 defined as fine gradation parameters showed highest positive correlation (0.72) with clay and NMC; FS, silt and fines (0.53, 0.61, 0.67); and negatively correlated (>0.75) with activity and liquid limit ratio (LLr). These variables are a function of fine soil texture. PC4 and PC5 are positively correlated (>0.70) with coarse materials (gravel, coarse and medium sand), referred to as coarse soil texture.

3.5.2 Agglomerative hierarchical analysis

Agglomerative hierarchical clustering (AHC) of the principal components showed a representation of the soil variables in homogenous classes where soils of same class exhibited similar values. The three major classes distinctively categorized by the AHC are depicted by dendrogram (Figure 5) which also displayed the cohesiveness of clusters formed. Class 1 soils were characterized by fines (silt, clay), plasticity and moisture contents (MC, OMC) parameters. Approximately, 47% of these soils were of the migmatite-gneiss derived origin (PCB) with a p-value <0.001. Similarly, Class 3 showed significant clustering with densities (BD, DD, and MDD) parameters of mainly metasediment origin (PCM). Class 2 soils performed poorly owing to few parameters and p-value >0.05, which might be from the older granite derived origin (PCG).

Figure 5.

Dendrogram of the studied soil variables.

3.5.3 Multiple regression analysis

Following stepwise regression method, five models were generated in (Eqs. 15). The result indicated that between 25 and 33% of the variation in soil properties was explained by the combination of these predictors. In Eq. 5, 70% training dataset accounted for 33% variance with coefficient of determination (R2 = 0.33) and root mean square error of performance (RMSE = 7.8). Given the p-value <0.001 computed by analysis of variance (ANOVA), the significance level (5%) and the low bias (0.05), the prediction by the explanatory variables is significant.


4. Conclusion

All the variables exhibited a large variation of data clustering around the mean value and high coefficients of variation. The soils within the area are predominantly very fine sands with a high percentage passing through No. 200 (0.075 mm). The natural moisture content and soil activity are moderately high which may be due to soil potential for retaining water. The result of Atterberg limits shows most of the soils classified as inorganic clayey soils under A-2 and A-7 groups. Their variability is locational and between lithology which reflects the influence of parent materials and pedogenic activities. Application of correlation analysis has allowed for the determination of the relationship between index properties, compaction and CBR and for deriving multivariate relationships for the assessment of CBR based on these parameters. Strong correlation existed among index properties but showed weak relationship with CBR under soaked condition. Principal component analysis categorized the measured soil parameters into five major groups for which first three components explained more than half of the total variance. Hierarchical classification validated the clustering of different individuals/variables based on the parent material. Five empirical models were obtained between soaked CBR and index parameters. The coefficient of determination (R2) and root mean square error (RMSE) revealed that the models obtained were able to predict the target variable to a good degree of accuracy. The pavement challenges witnessed on the highway are attributable to the poor subgrade, the influence of geology and lack of drainage. It is therefore recommended to stabilize the soil with cement or lime in order to improve the strength and provide drainage along the road corridor.



My appreciation goes to my Supervisors, Prof. Olusegun Ige and Prof. Olufemi Ogunsanwo for their support and indebt supervision of the research as well as Nigerian Building and Road Research Institute for assisting in the laboratory test.


Conflict of interest

The authors declare no conflict of interest.



Grain size distribution parameters according to their diameter
Ggravel size (4.0–13.2 mm)
CScoarse sand size (0.6–2 mm)
MSmedium sand size (0.2–0.425 mm)
FSfine sand size (0.075–0.2 mm)
Ssilt size (0.002 mm)
Cclay size (0.001 mm)
Atterberg limit parameters
LLliquid limit
PLplastic limit
PIplasticity index
wPIweighted plasticity index [(PI × % passing 0.425 mm)/100]
NMCnatural moisture content
Fswfree swell
Compaction parameters
MDDmaximum dry density
DDdry density
BDbulk density
OMCoptimum moisture content
MCmoisture content
CBRuunsoaked California Bearing Ratio
CBRssoaked California Bearing Ratio
Other derived parameters
LLrliquidity ratio (LL/MC)
PIrplasticity ratio (LL/PL)
Wrmoisture ratio (MC/OMC)
Drdensity ratio (DD/MDD)
Acactivity (PI/clay fraction)
SPswelling potential (PI/PIr)
SDstandard deviation
CVcoefficient of variation [(SD/mean) × 100]
Rcorrelation coefficient
R2coefficient of determination
RMSEroot mean square error


  1. 1. Patel RS, Desal MD. CBR predicted by index properties for alluvial soils of South Gujarat. In: Proceedings of the Indian Geotechnical Conference. Mumbai: Geotrendz. IGS Mumbai Chapter & IIT Bombay; 16-18 December 2010. pp. 79-82
  2. 2. Carter M, Bentley SP. Shrinkage and swelling characteristics. In: Soil Properties and Their Correlations. Chichester, UK: John Wiley and Sons; 2016. DOI: 10.1002/9781119130888.ch8
  3. 3. Owoseni JO, Adeyemi GO, Asiwaju-Bello YA, Anifowose AYB. Engineering geological assessment of some lateritic soils in Ibadan, southwestern Nigeria using bivariate and regression analyses. African Journal of Science and Technology (AJST), Science and Engineering Series. 2012;12(1):59-71
  4. 4. Yildirin B, Gunaydin O. Estimation of California bearing ratio by using soft computing systems. Expert Systems with Applications. 2011;38(5):6381-6391
  5. 5. Yitagesu FA, Van Der Meer F, Van Der Werff H, Seged H. Evaluation of soil expansion index from routinely determined geotechnical parameters. Soil Science Society of America Journal. 2011;75:1640-1651
  6. 6. Oluyide PO, Nwajide CS, Oni AO. The Geology of Ilorin area with explanations on the 1:250,000 series, Sheet 50 (Ilorin). Geological Survey of Nigeria Bulletin. 1998;42:1-84
  7. 7. British Standards Institution 1377: Methods of Test for Soils for Civil Engineering Purposes. General Requirements and Sample Preparation. London, United Kingdom: British Standards Institution; 1990. p. 6. ISBN: 0-580-18030
  8. 8. Jollife IT. Principal Component Analysis. 2nd ed. Springer series in Statistics. New York, USA: Springer-Verlag; 2002
  9. 9. Murtagh F, Legendre P. Ward’s hierarchical clustering method: Clustering criterion and agglomerative algorithm. Journal of Classification. 2014;31(3):274-295
  10. 10. Ige OO, Fakeye AM, Ogunsanwo O. Geotechnical assessment of some along Ilorin-Lokoja highway—Implication on suitability for road construction. Bulletin of the Science Association of Nigeria. 2018;29:149-162
  11. 11. Nwaiwu CM, Afolayan JO, Osinubi KJ. Predicting hydraulic conductivity of compacted lateritic soils, a reliability approach. Electronic Journal of Geotechnical Engineering. 2006;24(1):45-56
  12. 12. Adams JO, Aderinola OS, Akinwamide JT. Geotechnical study of pavement indices influencing failures along ado-Ajabandele-Ikere road, South Western, Nigeria. International journal of Engineering Innovation and Research. 2015;4(4):567-572
  13. 13. Bayamack JF, Onana VL, Ndzie AT, Ngo’oze A, Ohandja HN, Eko RM. Assessment of the determination of California bearing ratio of laterites with contrasted geotechnical properties from simple physical parameters. Transportation Geotechnics. 2019;19:84-95
  14. 14. Adams JO, Adetoro AE. Analysis of road pavement failure caused by soil properties along Ado-Ekiti–Akure Road, Nigeria. International Journal of Novel Research in Engineering and Sciences. 2014;1(1):1-7
  15. 15. Ndzie AT, Onana VL, Ngo’oze A, Nyassa OH, Ekodeck GE. Influence of hydromorphic conditions in the variability of geotechnical parameters of gneiss derived lateritic gravels in a savannah tropical humid area (Centre Cameroon), for road construction purposes. Transportation Geotechnics. 2017;12:70-84
  16. 16. Adeyemi GO, Olarewaju VO, Akintunde OB, Meshida EA. Mineralogical and geotechnical characteristics of some subgrade soils in a section of the Ibadan/Ile-Ife expressway, Southwestern Nigeria. Journal of Applied Sciences. 2003;6(2):3536-3547
  17. 17. Johnson RA, Wichern DW. Applied Multivariate Statistical Analysis. Eaglewood Cliffs, New Jersey: Prentice-Hall; 1992
  18. 18. Shukla MK, Lal R, Ebinger M. Soil quality indicators for the north Appalachian experimental watersheds in Coshocton, Ohio. Soil Science. 2005;169:195-205

Written By

Attah Fakeye, Olusegun Ige and Olufemi Ogunsanwo

Submitted: July 8th, 2020 Reviewed: July 30th, 2020 Published: April 28th, 2021