Soil Science Society of America Journal 64:1423-1433 (2000)
© 2000 Soil Science Society of America
DIVISION S-5-PEDOLOGY
Quantifying Soil Morphology in Tropical Environments
Methods and Application in Soil Classification
Anne Gobin,
Paul Campling,
Jozef Deckers and
Jan Feyen
Institute for Land and Water Management, Katholieke Universiteit Leuven, Vital Decosterstraat 102, 3000 Leuven, Belgium
anne.gobin{at}agr.kuleuven.ac.be
 |
ABSTRACT
|
|---|
We tested the hypothesis that readily observed and easily measured morphological variables can be used to characterize the soils sampled and described in southeastern Nigeria for purposes of land use and management. Field tests were developed for estimating soil texture and amount of ironstone nodules. Two new soil color indices provided an immediate means of diagnosing the soil drainage regime in case of the color index (CI) and soil forming processes in tropical soils in case of the redness index (RI). The indices correlated negatively with organic C content
and positively with dithionite-extracted Fe2O3 (0.44) and Al2O3 (0.51). Inexpensive field tests for color, texture, and ironstone can be quantified using color indices and laboratory measurements. The local soil classification was quantified by means of color indices (RI, CI) and percentages of ironstone, sand, silt, and clay measured in the A horizon. A classification based on soil texture, ironstone, and color was used to define classes for the B horizon. The two first principal components (PC) extracted from soil morphological variables measured on the upper three horizons of 72 pedons explained 64.7% of the total variance. Nonhierarchical clustering performed on the two PCs produced seven clusters that compare well with the great groups of U.S. soil taxonomy. Principal component analysis on 20 soil chemical and morphological variables confirmed that soil texture, ironstone, and soil color account for most of the variation of the soils and provide an efficient means of characterizing tropical soils derived from sedimentary parent material.
Abbreviations: BS, base saturation CEC, cation-exchange capacity CI, color index CV, coefficient of variation EA, exchangeable acidity PC, principal component RI, redness index RR, redness rating SD, standard deviation
 |
INTRODUCTION
|
|---|
THE INCREASING USE of geographical information systems and earth observation techniques in land resources analysis has highlighted the need for quantitative data on the spatial distribution pattern of soil characteristics. However, many soil surveys have concentrated on the vertical sequence of horizons within pedons, paying less attention to the spatial distribution (FAO, 1998). Parametric soil surveys concentrate on measuring single soil characteristics, and provide, in concert with process modeling, a suitable paradigm for spatial prediction across low-surveyed regions (McKenzie and Austin, 1993; Moore et al., 1993). The use of readily observed and quantifiable morphological characteristics to distinguish between different soils as a first step to determining their spatial extent is our major concern here.
Soil morphological descriptions are commonly recorded and mostly used to aid classification purposes (Soil Survey Staff, 1998; FAO, 1998). Standard soil morphological descriptions have been quantified and combined in a soil profile development index for evaluating soil development (Bilzi and Ciolkosz, 1977; Meixner and Singer, 1981; Harden, 1982). In addition, McKenzie and Jacquier (1997) describe the use of soil morphological descriptions together with inexpensive field tests to derive soil hydraulic properties in low-surveyed regions. Last, soil color indices present a quantified approach for assigning drainage class to a soil (Megonigal et al., 1993; Thompson and Bell, 1996) and assessing the Fe oxide mineralogy (Torrent et al., 1980; Mokma, 1993).
Farmers often describe soils in combinations of single morphological characteristics (e.g. red sand or stone) and often relate their decision-making on land use and management to these soil descriptions (Gobin et al., 1998, 1999). Quantifying these morphological characteristics opens new perspectives for incorporating farmers' knowledge into land resources information systems, and enables statistical modeling to be used. Farmers' knowledge could benefit scientific understanding of soils, contribute to international agricultural development, and facilitate exchange between farmers and researchers (Sandor and Furbee, 1996; Alexander, 1996; Habarurema and Steiner, 1997; Norton et al., 1998).
The objectives of this study were to quantify field observations of soils made by scientists and local farmers, relate field observations to laboratory analysis, and identify soil variables for distinguishing soils in a tropical area. Our method for quantifying field observations uses laboratory techniques and color indices measured on 72 pedons. We developed two soil color indices that incorporate a weighting factor for matrix and mottling colors, and we identified soil variables that distinguish the soils using exploratory methods, analysis of variance, and multivariate analysis.
 |
Materials and methods
|
|---|
Study area
The 589-km2 study area comprises the River Ebonyi headwater catchment and a part of the Udi-Nsukka Plateau located in Enugu State, Nigeria (Fig. 1)
. The study area traverses from west to east a part of the Udi-Nsukka Plateau and Escarpment and a section of the Cross River Plains, characterized by a complex sedimentary history (Benkhelil, 1988). In general, argillaceous rocks underlie the Plains, whereas the Plateau and Escarpment are formed by sandstone. Flat-topped hills or ridges with surface ironstone are common on the Plateau and east of the Escarpment and represent remnants of the Coal Measure formations (Jungerius, 1964). The fluvial landscape farther on the Plains is characterized by a meandering river channel bordered by narrow river banks, with seasonally submerged backswamps from combined river flooding and upland run-on (Gobin et al., 1998). The upper interfluve is undulating to rolling, whereas the lower interfluve is flat to undulating.
The region has a humid tropical climate with a distinct dry season between November and March and an annual rainfall averaging
1500 mm yr-1 (Köppen climate classification is Aw). The area is situated in the transition zone between lowland, Guinea-Congolian, wetter-type rainforest and Guinea savanna, resulting in a mosaic vegetation pattern (Hopkins, 1979; White, 1992). Luxuriant evergreen forest fringes the river and perennial streams on the regularly flooded soils of the valley bottom, whereas along seasonal streamlines corridors of semideciduous trees and bushes are found. Moist semideciduous forest occurs on the Plateau and shale dominant interfluve, whereas drought-tolerant tree species and tall grasses mainly occupy the Escarpment and denudated gravelly interfluve areas.
Location of Pedons
Generalized information on geology was derived from the 1:250000 geological maps of Nigeria (Shell-BP, 1957). Vegetation, drainage, and relief were obtained from stereoscopic analysis of 1:40000 aerial photographs from 1962 and topographic maps at a scale of 1:50000. The study area was stratified according to combined geology, drainage, relief, and vegetation characteristics. Based on this stratification, 40 pedons were located along the steepest environmental gradients following the catena concept (Sommer and Schlichting, 1997) to sample the full range of soils in the study area. Soil pits were dug (Fig. 1) at representative locations of each landscape along the toposequences. These locations were identified using a clinometer, a compass, and a global positioning system (Trimble Pathfinder Basic and Software; Trimble, 1992). The catena concept was also employed in former soil surveys within the study area: 32 pedons were sited in the field and on the topographic map (Akamigbo and Opara, 1977; Federal Department of Agricultural Land Resources, 1985; Asadu, 1986; Nwadialor, 1989) (Fig. 1). A relational database (Prague and Irwin, 1997) was set up to process all information obtained on the 72 pedons.
Field Methods
Soil profile observations were fully described according to the FAO guidelines (1990). Soil texture of the fine earth was approximated using a manual field test (International Land Development Consultants, 1981). The extent to which a tablespoon (
15 mL) of moist fine earth can be shaped is indicative of its texture (Table 1)
. On the Interfluve, the test was improved with the knowledge that loam is powdery when dry and leaves dirt on the skin when moist, and clay displays shining faces when augering. The occurrence of ironstone nodules, usually between 2 and 20 mm in size, was expressed in percentage surface covered: few (05%), common (515%), many (1540%), abundant (4080%), and dominant (>80%). Soil colors were determined per horizon using the revised standard soil color charts (Takehara, 1992), which are based on the Munsell soil color chart (Munsell, 1954). The surface percentage of mottling and of each soil matrix color was estimated in the field. The mottling size was recorded in millimeters as fine (<5 mm), medium (515 mm), or coarse (>15 mm), and the mottling abundance was described as few (<2%), common (220%), and many (2050%).
The term local soil classification refers to the classification that local farmers use to name the soils in the Igbo language. Techniques borrowed from participatory rural appraisals (Chambers, 1992) were conducted along the toposequence with the aid of two local village guides and trained interpreters and facilitators from outside the village to elicit the local soil classification scheme. Open-ended and semistructured interviews (Mettrick, 1993) were administered to villagers cultivating or owning fields within 50 m from a profile pit. The soil adjacent to each field was augered and described to verify its similarity with the reference soil pit and to facilitate the interview. Villagers were asked to compare their descriptions, knowledge, and name of a particular soil with soils located farther along the toposequence. Accounts of within-field variation were not incorporated into the classification scheme. Schematic diagrams were constructed for each of the toposequences, and the soil classification was verified during a group discussion involving farmers and village elders.
Laboratory Methods
Samples were taken per horizon, air-dried at an ambient temperature of
25°C, crushed, and passed through a 2-mm sieve. The percentage of the fraction >2 mm was weighed, and the ironstones were weighed separately. Soil chemical analysis was carried out on the sieved fraction (Table 2)
. Soil classifications according to U.S. soil taxonomy (Soil Survey Staff, 1998) were inserted into the database. Local soil descriptions were compared with soil variables measured by standard laboratory techniques and expressed as numerical values.
Soil Color Indices
Soil color based on the soil Munsell color chart consists of hue, value, and chroma. The color notations were converted into color indices to arrive at a single numerical value. Previously published work suggests that color indices are useful for quantifying differences in soil morphology. Hurst (1977) defined the Redness Rating (RR) as an assigned numerical value for hue, multiplied by chroma and divided by value. The RR proved a good discriminator between hematite and goethite and correlated highly with the fine earth hematite content in red Mediterranean soils (Torrent et al., 1980). Chroma, augmented with or without an arbitrary value for hue, and in both cases corrected for all mottle and matrix fractions, corresponded well with the hydric conditions of seasonally saturated soils (Thompson and Bell, 1996). The decrease in Munsell value along a toposequence was found to be highly correlated with increased organic C content (Fernandez et al., 1988). A decreasing color number, which was based on value and chroma added to a numerical value for hue, was highly correlated with increasing organic C and to a lesser extent with Fe and Al content in Spodosols (Mokma, 1993). Harden (1982) and Harden et al. (1991) also developed Melanization, rubification, color lightening and color paling for desert soils.
Two new indices, a redness index (RI) and a color index (CI), combined and modified components from previous studies to suit soil colors typical for active tropical weathering and to incorporate a weighting for matrix and mottling colors. Hues of all soil horizons in the study area ranged from 10YR (Yellowish Red) to 10R (Red) and were transformed to a numerical value (HT) so that redder hues would correspond with higher values:
 | (1) |
Where "hue" is the Munsell hue. Numerical values, based on the mean of the possible abundance range, were assigned to each mottling fraction depending on the combination of the recorded abundance and mottling size. A numerical value of 0.01 was assigned to few-fine and few-medium mottling; 0.11 corresponded with the combinations few-coarse, common-medium, common-fine, and many-fine; 0.35 was reserved for common-coarse, many-medium, and many-coarse. The sum of all mottling fractions was not allowed to exceed the field-estimated total fraction of mottle colors. The fraction of each matrix color (Pi) is:
 | (2) |
Where pi is the percentage of each matrix color estimated in the field description, m is the number of mottle colors, Sj is the fraction of each mottle color. A redness index (RI) and a color index (CI) were designed according to:
 | (3) |
 | (4) |
Where n is the number of matrix colors, V is the Munsell value, C is the Munsell chroma. The second terms of Eq. [2] and [3] are referred to as the mottling redness index (MoRI) and mottling color index (MoCI), respectively.
Statistical Analysis
The statistical analysis aimed at identifying soil variables that can be used to distinguish the different soils of the study area.
Exploratory and Correlation Analysis
Exploratory data analysis included descriptive statistics and a check for normality (SAS Institute, 1990) on all soil characteristics per horizon. Most soil variables were transformed to normal distribution using a natural logarithm or square root. Soil texture was transformed according to Eq. [5]; the amount of ironstone was transformed according to Eq. [6] (Webster and Olliver, 1990).
 | (5) |
 | (6) |
Where x is sand, silt, or clay fraction expressed in percentage of fine earth, and y is percentage of ironstone. Correlation between the soil variables was examined.
Analysis of Variance
Mean, standard deviation (SD) and coefficient of variation (CV) were used to evaluate the within- and between-classes homogeneity in the different soil characteristics. Classes of the A horizon from all 72 pedons were based on the local classification, and classes of the B horizon were based on soil texture, ironstone, and color. A Fisher test was used to assess whether the class means are different from each other at the 0.05 significance level (SAS Institute, 1990). Multiple comparison was performed on each normalized soil variable, and subsequent statistical grouping was accomplished through Duncan-Waller multiple range tests at the 0.05 significance level (SAS Institute, 1990).
Multivariate Analysis
Multivariate analysis was used to identify interrelationships among soil chemical and morphological variables and to study differentiation patterns among the different soils using two matrices. A 72 by 27 matrix, X, consisted of nine normalized morphological variables per horizon for the upper three horizons (p) of 72 pedons (n), namely horizon thickness, sand, silt, clay, percentage of ironstone, and four soil color indices (RI, CI, MoRI, and MoCI). A second 72 by 60 matrix, Y, was created from 72 pedon observations on 20 normalized soil chemical, textural, and morphological variables for the upper three horizons. The calculations explained hereafter for X were also performed on Y. All normalized soil variables were standardized to zero mean and unit variance:
 | (7) |
Where x are the normalized observation values of the n by p matrix,
is the mean vector of each normalized variable, and sj the corresponding variance. The corresponding standardized correlation matrix (RX) was defined:
 | (8) |
RX was examined according to three criteria (Jobson, 1992; SPSS, 1994):
Test for Zero Correlation
The correlation matrix equals the identity matrix when u has a
2 distribution with 0.5p(p - 1) degrees of freedom, whereby u equals:
 | (9) |
Where
j are the eigenvalues of the correlation matrix. The test statistic is highly significant for matrices that are not orthogonal and thus appropriate for factoring.
Kaiser-Mayer-Olkin (KMO) Measure of Sampling Adequacy.
 | (10) |
Where rij is the correlation coefficient and aij is the partial correlation coefficient between variables i and j. Values for KMO below 0.5 are unacceptable.
Individual Measure of Sampling Adequacy (MSAi) for Each Variable i
A value close to one is excellent, whereas values smaller than 0.5 require elimination of the particular variable.
 | (11) |
Principal components (SAS Institute, 1990) were extracted from the correlation matrix, using the equation Z = XV, where Z is the matrix of unstandardized PCs, X the data matrix, V the matrix of eigenvectors with V'V = VV' = I, and I is the identity matrix (Jobson, 1992). The solution is an eigenvalue problem according to:
 | (12) |
Where
is a diagonal matrix of eigenvalues of RX (Eq. [8]), and V is orthogonal to RX. As the PCs are standardized, the values of the standardized components for n observations were determined using:
 | (13) |
Where Z* = Z
-1/2 and V* =
1/2V.
The p PCs are standardized linear combinations of the p original soil variables, with coefficients equal to the eigenvectors of the correlation matrix. The total variance (s
) is given by:
 | (14) |
Where
i is the eigenvalue or variance of the ith PC (PCi). Larger eigenvalues correspond with a larger explained portion of the total variance. All PCs were given a rank (eigenvalue number) according to their eigenvalue. The mean eigenvalue was used as a criterion to retain a reduced set of q components:
 | (15) |
The correlation between the original variables and PCs (loadings) are used for interpreting the PCs. Scatterplots of the observations in the plane of the PCs with the largest eigenvalues (scores) (Eq. [14]) provide a good estimate of the pattern of variation. The same procedure (Eq. [7][15]) was employed to extract PCs from matrix Y.
Nonhierarchical cluster analysis following the k means algorithm was performed on the first two unstandardized PCs of matrix Z derived from matrix X (Eq. [12] and [13]). The Euclidean distance between group centroids was used to measure the proximity between groups:
 | (16) |
Where d is the Euclidean distance between class r and s, and
j contain the coordinates of the centroids. The dispersion within classes, expressed as the sums of squares of deviations from the group means, was minimized through subsequent iterations in order to arrive at an optimal number of classes. Records are reassigned until they are located in the group with the nearest centroid. The resulting clusters were compared with classification based on soil texture, ironstone, and soil color and with classification according to U.S. soil taxonomy.
 |
Results
|
|---|
Field Observations
In general, all pedons examined are closely related to parent material, which in turn is strongly associated with landform and landscape. On the Plateau, soils formed from sandstones and colluvial material derived from Upper Coal Measures occur on nearly flat slopes and dry valleys. They consist of reddish brown sandy loam over crumbly red sandy clay loam (2.5YR 4/810R 4/8). At Escarpment, soils are formed on sandstone and contain high amounts of sand throughout the profile; the subsoil contains weak crumb-structured sand with a very friable consistence. Soils rich in ironstone nodules dominate the Residual Hills. The Interfluve is characterized by soils derived from shale and having an argillic horizon. On the ridges and residual hills of the Interfluve, ironstone material occurs at the surface as lag gravel or at shallow depth, mostly in concretionary form. Fairly well-drained soils of the Upper Interfluve display argillic subsurface horizons without pronounced mottling. Soils displaying red Fe mottles (2.5YR 4/6) within the 2-m profile depth occur on the slopes and River Terraces of the Upper Interfluve and on the somewhat better-drained areas of the Lower Interfluve. Seasonally waterlogged soils of the Lower Interfluve contain plinthite and a yellow mottling (10YR 6/8) in a gray soil matrix. Sandy to very sandy soils have developed on the riverbanks, whereas the soils of the backswamp have a finer texture and display redoximorphic features.
Farmers of the Nsukka Agricultural Zone can identify major soil types according to morphological characteristics of the soil to the depth they till or position in the landscape (Gobin et al., 1998, 1999). Although the local soil descriptions are not the same for every village, the underlying concepts or characteristics are similar. The characteristics used are easy to recognize and include the occurrence of ironstone, texture, and color, to which secondary soil properties such as workability, drainage, and water-holding capacity are attributed. The local soil names are based on the soil morphological characteristics, and a translation of the local names was used to reconstruct the local classification (Fig. 2)
. The soil names are directly linked to decision-making in land use and management. For example, soils of the residual hills and ridges that are named stone are considered marginal for agricultural uses.

View larger version (29K):
[in this window]
[in a new window]
|
Fig. 2 Local soil classification scheme. Local farmers refer to the part of the Cross River Plains within the study area as Ebonyi or Isi Uzo. Ebonyi is the main river of the study area; Isi Uzo is the district named after the old main road crossing the study area. Dark soil is a synonym for black soil; clay is often referred to as sticky
|
|
Quantifying Field Observations
The huge differences in soils enabled the field estimation method to be tested for a wide range of soil textures. Compared with the laboratory textural analysis, the accuracy was 100% in detecting USDA soil textural classes. Field estimates of ironstone content corresponded with the following classes by weight percentage: few (<10%), common (1030%), many (3050%), abundant (5070%), and dominant (>70%). RI and CI correlated positively with dithionite-extracted Fe2O3 (R = 0.44) and Al2O3 (0.51) and negatively with organic C content (-0.39) (Fig. 3)
. However, soil color is also dependent on the mineralogical composition of parent materials. Values for RI and CI were calculated for dominant mottle and matrix colors (Table 3)
. The highest values for RI were found in red soils dominated by sesquioxides (oxic), whereas low values for RI were found on gray and yellow soils (argillic). High values for CI corresponded with well-drained soils; intermediate values occurred on fairly well drained clayey soils (argillic); and the lowest values were associated with seasonally saturated soils (redoximorphic). Both indices were strongly correlated, but the difference in physical meaning, that is, relation to drainage (CI) and mineralogy (RI), justified keeping them both in as independent measurements.

View larger version (38K):
[in this window]
[in a new window]
|
Fig. 3 Scatterplots and correlation coefficients between CI, RI, and soil chemical properties (** indicates significance at the 0.01 level)
|
|
Summary statistics for each soil variable showed that laboratory determination of texture, ironstone content, and the two color indices were useful in distinguishing the local soil classes (Fig. 2, Table 4)
. Despite its frequent use in soil science (Webster and Olliver, 1990), the CV is difficult to interpret, particularly in the case of proportional data. The classes "sticky sand/red" (1 [Class number; Table 4]), clay (5) and "heavy clay" (6) had the same SD for ironstone content, which made their CV, ranging from 92.3 to 323.7, entirely dependent on the mean (Table 4). However, the Duncan-Waller grouping provided an easy way of recognizing differences between soil classes for a particular soil variable; classes with the same letter for a group were statistically not significant from each other for that variable. "Stone" and "stony" related to the presence of ironstone nodules (>2-mm fraction). The ironstone of the "stone" (3) and "stony clay" (4) class content differed significantly from each other and from the other classes. The fine earth sand fraction (50 µm2 mm) significantly separated the "sand" classes (1, 2, and 7) from the other classes. The local term "stickiness" of a soil, often translated as clay, compared best with a combination of the fine earth silt and clay fraction (<50 µm) (Table 4). The composite variable "clay and silt" was significantly different for the "heavy clay" class (6) and the "clay" class (5); the former also contained significantly less sand. The term "heavy" reflected the workability rather than the textural composition per se. The two color indices also related to the local observations (Table 4). Farmers discerned a yellow to red color range and a degree of whiteness and blackness when describing a soil. The yellow to red range related to soil Fe dynamics. The "sand" classes (1, 2, and 7) displayed significant differences for the two color indices. The color white was mainly used in conjunction with clay, and referred to a thin light-colored powdery surface layer that occurred on well-drained clays. The black color was linked to the presence of organic matter. Organic C and N values were highest for soils rich in ironstone and lowest for the Escarpment soils; a somewhat similar tendency was observed for the cation-exchange capacity (CEC) (Table 4).
Similar to the local classification, a classification based on soil texture, color, and ironstone was made for the B horizon (Table 5)
. The amount of ironstone significantly differed for Class 3, soils having ironstone throughout the profile, and Class 4, soils having ironstone in subsoil, from other classes. Classes 3 and 4 both had higher values for the color indices, and their morphological similarity implied that they could only be separated on the basis of chemical variables (pH, base saturation [BS]). However, land use and management are similar on the two soil classes. Sand content and the RI showed significant differences between Classes 1, 2, and 7. The classes 5B, 6A, and 6B displayed significant differences in clay content; Class 6A had a higher ironstone and lower sand content. Class 6B could be further isolated using the CI and RI. Both RI and CI increased with profile depth; the differentiation between the soil classes was pronounced for CI and even more so for RI values of the B horizon (Tables 4 and 5). Differences between the top- and subsoil are often related to organic C and eluviation processes. Fluctuating water tables influence redox processes in low-lying areas of the landscape. Compared with the topsoil classes, the subsoil classes also displayed a higher level of differentiation in the soil chemical variables. Values for the CEC were highest for Classes 6A and 6B, followed in descending order by Classes 5A, 4, 3, 5B, 1, 2, and 7.
Comparison with U.S. Soil Taxonomy
There is a distinct relation between landform and soil according to the local classification and U.S. soil taxonomy (Table 6)
. Quantifying the local descriptions in terms of soil texture, ironstone, and color indices enabled application of the classification to the subsoil and examining the separation of soil classes, as defined in Table 4 and 5, using multivariate analysis. Cluster analysis was performed on the two first PCs extracted from RX. RX had a KMO of 0.813, a value of the test statistic for zero correlation of 3128 (P < 0.001), and values for MSAi >0.7. The PC method of factor extraction yielded five independent components according to criterion Eq. [14], explaining 85.5% of the total variance. Principal Component 1 seemed to be a measure of mineralogy since the most significant loadings are sand, silt, RI, and subsoil ironstone (Table 7)
. Principal Component 2 had high positive loadings for CI, surface ironstone, and subsoil RI; this component seemed to be a measure of drainage. Prinicipal Component 3 related to mottling color and intensity, PC4 mainly expressed the stoniness, and PC5 was a measure for thickness of the upper horizons. The projection of component scores on the diagram of PC1 vs. PC2 showed a wide scatter with an uneven distribution that could be separated into seven distinct clusters (Fig. 4, I)
. Along the PC1 axis, soils having a high value for the RI or sand content could be distinguished from soils with high silt contents or high subsoil ironstone content or mottling. Along the PC2 axis, negative values were for soils with a high sand content or poor drainage characteristics (low values for CI), whereas positive values indicated high values for CI or surface ironstone. As expected, the clusters presented similarities with the soil morphological classification (Fig. 4, II). Some of the soil morphological classes separated well into great groups or subgroups of U.S. soil taxonomy, particularly where the distinctions are associated with soil geomorphic relationships (Fig. 4, III). Cluster 1 represented the red, slightly sticky, sand soils on the Plateau, which are Rhodic Kandiustox. Cluster 2 contained red, sand soils or Ustoxic Quartzipsamments of the Escarpment. Cluster 3 grouped soils high in subsoil ironstone content, whereas Cluster 4 consisted of soils containing plinthite and some subsurface ironstone. Soils having a high silt content and mottling belonged to Cluster 5. These seasonally waterlogged soils of the lowland were classified as "heavy clay" in the local classification and Plinthaquults, Plinthohumults or Plinthustults in U.S. soil taxonomy. Cluster 6 grouped the riverside sand and corresponded with the Oxyaquic Quartzipsamments. Cluster 7 represented soils that are dominated by clay and are fairly well drained; they are clay soils of the upland Interfluve and were classified as Kandiustults.
View this table:
[in this window]
[in a new window]
|
Table 7 Loading values (between parentheses) of soil variables for principal component (PC) analysis performed on submatrix X containing only soil morphological variables
|
|

View larger version (45K):
[in this window]
[in a new window]
|
Fig. 4 Scatterplots for pedons on the first two principal components according to analysis on X. (I) nonhierarchical multivariate classification, (II) soil morphological classification according to Table 5, (III) great groups of U.S. soil taxonomy
|
|
Results from the PC analysis performed on Y, including the soil chemical variables, showed that the first diagram, PC1 vs. PC2, was more successful in separating soil classes or great groups (Fig. 5)
. RY had a KMO value of 0.735, a test statistic for zero correlation of 5184 (P < 0.001), and large values for MSAi. Ten PCs explaining 85% of the total variance were extracted (Table 8)
. Principal Component 1 contributed to 31.1% of the variance and was strongly dominated by soil texture, CEC, subsurface ironstone, and RI (Table 8). Prinicipal Component 2, representing 15.36% of the variance expressed the soil color (CI, RI), organic C content, and pH. The soil chemical variables pH, BS, and exchangeable acidity (EA), dominated PC3, whereas mottling colors, organic C, and BS correlated strongly with PC4. The loading values of the remaining components were dominated by soil chemical variables, apart from PC7 and PC8 (Table 8). The soil chemical variable CEC contributed to distinguishing the soil classes (Fig. 5, I.A and II.A). The PC3 vs. PC4 diagram (Fig. 5, I.B and II.B) was not successful in separating soil classes. The Duncan-Waller grouping (Tables 4 and 5) provided information that the soil chemical variables, having a high loading on PC3 and PC4, were not significantly different between the soil classes.

View larger version (42K):
[in this window]
[in a new window]
|
Fig. 5 Scatterplots for pedons projected in the plane of the first four principal components according to analysis on Y. (I) great groups of U.S. soil taxonomy, (II) soil morphological classification according to Table 5
|
|
View this table:
[in this window]
[in a new window]
|
Table 8 Loading values (in parentheses) of soil variables for principal component (PC) analysis performed on the full data set (Y)
|
|
 |
Conclusions
|
|---|
The local soil classification is based on soil texture, ironstone, and soil color of the topsoil to tillage depth, and relates directly to land use and management. Field tests developed for estimating soil texture and amount of ironstone nodules were in accordance with conventional laboratory measurements. Two new color indices, the redness index (RI) and color index (CI), were developed on the basis of Munsell color notation including weighting factors for soil matrix and mottling colors. Analysis of variance techniques and multivariate analysis of soil chemical and morphological variables showed that soil color, soil texture, and ironstone content of both top- and subsoil explained most of the variation of the soils in this tropical case study area. However, the limitation of the soil morphological classification is that it does not have the sufficient chemical information to precisely distinguish between the great groups for exploring potential land use and management strategies, nor does it provide evidence of the generic origin of soils.
Soil texture, occurrence of ironstone, the RI, and the CI vary with parent material, which in turn is strongly associated with position in the landscape. The implications of our results suggest that digital terrain analysis and subsequent soil-landscape modeling could be useful for predicting the spatial distribution of these easily discernible soil characteristics that relate directly to local land use and management.Bray Kurtz 1945; FAO 1990; SPSS Inc 1994
 |
ACKNOWLEDGMENTS
|
|---|
Funding for this research was provided by the Belgian Agency for Development Co-operation (BADC) through the Inter-University Project on "Water Resources Development for domestic use and small scale irrigation in the rural areas of southeastern Nigeria". The authors thank the Belgian Soil Service for analyzing the soil samples. Special appreciation is extended to the project staff and farmers who contributed to this particular study. We are indebted to Prof. Akamigbo and Dr. Igwe from the Department of Soil Science of the University of Nigeria Nsukka, who kindly provided documents by the DSS (1977) and the FDALR (1985), and facilitated access to the Department's library. The authors are grateful to the anonymous referees for their constructive comments on an earlier version of this paper.
Received for publication March 1, 1999.
 |
REFERENCES
|
|---|
- Akamigbo, F., and N. Opara. 1977. Natl. Soil Correlation Conf., 5th. Nsukka, Nigeria. 12 Sept. 1977. Dep. of Soil Sci., Univ. of Nigeria, Nsukku.
- Asadu, C.L.A. 1986. Pedological studies of six soil toposequences on three sedimentary formations in Anambra State, Nigeria. M.S. thesis. Dep. of Soil Science, Univ. of Nigeria, Nsukka.
- Alexander M.J. The effectiveness of small-scale irrigated agriculture in the reclamation of mine land soils on the Jos Plateau of Nigeria. Land Degrad. Dev. 1996;7:77-84.
- Benkhelil, J. 1988. Structure et évolution géodynamique du bassin intracontinental de la Bénoué (Nigéria). Bull. Centres Recherche Exploration - Production Elf - Aquitaine 12:30128.
- Bilzi A.F., Ciolkosz E.J. A field morphology rating scale for evaluating pedological development. Soil Sci. 1977;124:45-48.
- Bray R.H., Kurtz L.T. Determination of total, organic and available forms of phosphorus in soils. Soil Sci. 1945;59:39-45.
- Chambers R. Methods for analysis by farmers: The professional challenge. UK: IIED, 1992.
- FAO. Guidelines for soil profile description. AGLS, Rome: FAO, 1990.
- FAO. 1998. The world reference base for soil resources. World Soil Resources Rep. 84. ISSS/AISS/IBG/ISRIC/FAO, Rome.
- Federal Department of Agricultural Land Resources. 1985. The reconnaissance Soil Survey of Anambra State Nigeria (1:250,000). FDALR, Anambra State, Nigeria.
- Fernandez R.N., Schulze D.G., Coffin D.L., Van Scoyoc G.E. Color, organic matter, and pesticide adsorption relationships in a soil-landscape. Soil Sci. Soc. Am. J. 1988;52:1023-1026.[Abstract/Free Full Text]
- Gobin A., Campling P., Deckers J., Feyen J. Integrated toposequence analysis at the confluence zone of the River Ebonyi headwater catchment (southeastern Nigeria). Catena 1998;32:173-192.
- Gobin A., Campling P., Deckers J., Poesen J., Feyen J. Soil erosion assessment at the Udi - Nsukka Cuesta (southeastern Nigeria). Land Degrad. Dev. 1999;10:141-160.
- Habarurema E., Steiner K.G. Soil suitability classification by farmers in southern Rwanda. Geoderma 1997;75:75-87.
- Harden J.W. A quantitative index of soil development from field descriptions: Examples from a chronosequence in central California. Geoderma 1982;28:1-28.[ISI]
- Harden J.W., Taylor E.M., Hill C., Mark R.K., McFadden L.D., Reheis M.C., Wells S.G. Rates of soil development from four soil chronosequences in the southern Great Basin. Quaternary Res. 1991;35:383-399.
- Hopkins B. Forest and savanna, 2nd ed London: Heinemann Education Books, 1979.
- Hurst V.J. Visual estimation of iron in saprolite. Geol. Soc. Am. Bull. 1977;88:174-176.[Abstract]
- International Land Development Consultants. 1981. Agricultural compendium for rural development in the tropics and subtropics. Elsevier, Amsterdam, the Netherlands.
- Jobson J. Applied multivariate data analysis. Volume II: Categorical and multivariate methods. New York: Springer Verlag, 1992.
- Jungerius P.D. Some aspects of the geomorphological significance of soil texture in Eastern Nigeria. Z. Geomorphol. 1964;9:332-345.
- McKenzie N.J., Austin M.P. A quantitative Australian approach to medium and small-scale surveys based on soil stratigraphy and environmental correlation. Geoderma 1993;57:329-355.
- McKenzie N.J., Jacquier D. Improving the field estimation of saturated hydraulic conductivity in soil survey. Aust. J. Soil Res. 1997;35:803-825.
- Megonigal J.P., Patrick W.H., Faulkner S.P. Wetland identification in seasonally flooded forest soils: soil morphology and redox dynamics. Soil Sci. Soc. Am. J. 1993;57:140-149.[ISI]
- Meixner R.E., Singer M.J. Use of a field morphology rating system to evaluate soil formation and discontinuities. Soil Sci. 1981;131:114-123.
- Mettrick H. Development oriented research in agriculture. An ICRA notebook. The Hague, the Netherlands: CIP, 1993.
- Mokma D.L. Color and amorphous materials in Spodosols from Michigan. Soil Sci. Soc. Am. J. 1993;57:125-128.
- Moore I.D., Gessler P.E., Nielsen G.A., Peterson G.A. Soil attribute prediction using terrain analysis. Soil Sci. Soc. Am. J. 1993;57:443-452.[Abstract/Free Full Text]
- Munsell. 1954. Munsell soil color charts. Munsell Color Company, Baltimore, MD.
- Norton J.B., Pawluk R.R., Sandor J.A. Observation and experience linking science and indigenous knowledge at Zuni, New Mexico. J. Arid Environ. 1998;39:331-340.
- Nwadialor B.E. Soil-landscape relationships in the Udi-Nsukka plateau. Catena 1989;16:111-120.
- Prague C.N., Irwin M.R. The Access 97 Bible. Foster City, CA: IDG Books Worlwide, 1997.
- Sandor J.A., Furbee L. Indigenous knowledge and classification of soils in the Andes of Southern Peru. Soil Sci. Soc. Am. J. 1996;60:1502-1512.[Abstract/Free Full Text]
- SAS Institute. 1990. SAS/STAT user's guide. Version 6. 4th ed. Vol. 1. SAS Inst., Cary, NC.
- Shell-BP. 1957. Geological maps, 1:250,000 Sheets: Enugu (72), Geol. Survey Nigeria.
- Sommer M., Schlichting E. Archetypes of catenas in respect to matterA concept for structuring and grouping catenas. Geoderma 1997;76:1-33.
- SPSS Inc. SPSS Professional Statistics 6.1. Chicago: SPSS Inc, 1994.
- Soil Survey Staff. 1998. Keys to soil taxonomy. 8th ed U.S. Gov. Print. Office, Washington, DC.
- Takehara O.H. Revised standard soil color charts, 08.11. Giesbeek, the Netherlands: Eijkelkamp Agrisearch Equipment, 1992.
- Thompson J.A., Bell J. Color index for identifying hydric conditions for seasonally saturated Mollisols in Minnesota. Soil Sci. Soc. Am. J. 1996;60:1979-1988.[Abstract/Free Full Text]
- Torrent J., Schwertmann U., Schulze D.G. Iron oxide mineralogy of some soils of two river terrace sequences in Spain. Geoderma 1980;23:191-208.
- Trimble. 1992. Operating the GPS Pathfinder Basic Receivers version 5.29. Pathfinder software reference version 2.3. Trimble Navigation Ltd., Sunnyvale, CA.
- Walkley A., Black I.A. Determination of organic carbon in soil. Soil Sci. 1934;37:29-38.
- Webster R., Olliver M. Statistical methods in soil and land resource survey. New York: Oxford Univ. Press, 1990.
- White, F. 1992. Vegetation map of Nigeria. In W. Lanson (ed.) Plant ecology in West Africa: Systems and processes. Bayero University, Kano, Nigeria.
This article has been cited by other articles:

|
 |

|
 |
 
P. J. Drohan, E. J. Ciolkosz, and G. W. Petersen
Soil Survey Mapping Unit Accuracy in Forested Field Plots in Northern Pennsylvania
Soil Sci. Soc. Am. J.,
January 1, 2003;
67(1):
208 - 214.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
P. Campling, A. Gobin, and J. Feyen
Logistic Modeling to Spatially Predict the Probability of Soil Drainage Classes
Soil Sci. Soc. Am. J.,
July 1, 2002;
66(4):
1390 - 1401.
[Abstract]
[Full Text]
[PDF]
|
 |
|