Abstract

Soil formation occurs through numerous physical and chemical weathering processes acting to alter the parent rock on the Earth’s surface. Samples of surface soils were collected over a range of elevations (2000–3600 m) from profiles directly overlying basaltic to more felsic parent rocks, over a region in NW Ethiopia. The soils were investigated to determine their chemical composition and X-ray diffraction was used to identify and quantify individual mineral phases. The data set was analyzed using non-parametric statistics (Spearman’s Rank and Mann-Whitney U tests) to compare the soils forming over the two parent rocks. Principal component analysis (PCA) was used to identify the mineral alteration assemblage and formation during pedogenesis. The extent of alteration was quantified using several chemical weathering indices (Chemical Index of Alteration = CIA; Chemical Index of Weathering = CIW), including an index calculated by multivariate analyses of the soil chemical composition data (weathering “W” index). Further to this we devised and tested a new weathering index (Wmin) using multivariate analysis of the soil mineralogy, to estimate the extent of weathering and physico-chemical proprieties of the parent rock from which the soil formed.

The soils present a fair to advanced stage of alteration, with abundant iron (Fe) oxides (up to 40 wt%) and phyllosilicates (up to 57 wt%), including kaolinite-smectite (K-S) mixed-layer phases. The K-S was composed of either 30–50% kaolinite or 94–98% kaolinite layers. Discrete kaolinite was also present. The bimodal K-S mineralogical composition is likely due to two precursor phases: feldspar for the kaolinite-rich K-S and volcanic glass for the smectite-rich K-S. K-S with intermediate composition (50–94% kaolinite) was rare, due to its instability. Statistical analysis showed significant differences between the chemical compositions of the soils developed on the two different parent volcanic compositions. The soils overlying the more felsic parent rocks were less altered than those overlying the flood basalt. When comparing the weathering indices calculated in this study, we conclude that while the CIA and CIW may be more readily determined, the W and Wmin indices can elucidate information on the composition of the original rock from which they formed. The W index is more sensitive to certain variables when compared with the newly derived mineralogical Wmin index; however the Wmin index takes into account mineral phases within the sample, which provides a more detailed interpretation of weathering rates than chemistry alone. In addition the Wmin index correlated with meteorological variables, such as elevation (and consequently temperature and precipitation), known to influence the degree of pedogenesis. The Wmin index can be used to enhance our understanding of the processes that occur during weathering processes to supplement information gained from traditional chemical weathering indices.

Introduction

The formation of soils (pedogenesis) is a process by which weathering alters constituents within the parent deposit through the loss of more mobile (i.e., soluble) elements, concurrent enrichment of less mobile elements, combined with alteration and formation of new secondary minerals and accumulation of organic matter. The rate and extent of pedogenesis is highly dependent on a range of factors, including the lithology of the parent material, climatic conditions, time, presence of organisms within the strata, vegetation, and topography (e.g., Jenny 1941).

Extrusive igneous rocks are thought to be particularly sensitive to chemical erosion (Berner and Berner 1996), weathering 5–10 times faster than intrusive crystalline rocks (such as granite; Dessert et al. 2001). Furthermore, the weathering of mafic (basaltic) minerals is more rapid than that of felsic minerals (e.g., Sigfusson et al. 2008). Basaltic rocks are composed of primary minerals, mesostasis, and glass that, over time, develop into soils with comparably more stable mineral phases at surface conditions. Geochemical weathering of basalt begins with the loss of non-hydrolyzing cations (such as calcium, magnesium, and sodium) and the concomitant enrichment of silicon, aluminum, and iron ions in both crystalline and non-crystalline phases (Chorover et al. 2004). In terms of mineralogy, weathering typically removes the primary silicate phases (and glass), and replaces them with poorly crystalline aluminosilicates. Weathering can eventually produce iron oxide and oxyhydroxide phases that evolve toward more crystalline varieties (e.g., Rasmussen et al. 2010; Dahlgren et al. 2004; Righi et al. 1999).

Large igneous provinces are regions of extensive crustal emplacements of extrusive igneous rock, and one such example is the Ethiopia Volcanic plateau (EVP): a vast area (~0.81 × 106 km2; Dessert et al. 2003) igneous deposit in modern day Ethiopia. The EVP was formed when extensive eruptions deposited ~350 000 km3 of basalt (flood basalt; marking the appearance of the Ethiopian-Afar plume superswell) over a <2 Ma period (Mohr 1983; Mohr and Zanettin 1988). Several large shield volcanoes subsequently developed on the surface of the flood basalt (Mohr and Zanettin 1988), which in turn produced ~1/5 of the total volume of basaltic lava distributed in discrete outcrops across the flood basalt. Deposits from these (younger) shield volcanoes are generally more felsic, thinner, and less continuous than the extensive flood basalts that they cover. Mt. Choke is a shield volcano located in the center of the EVP, 215 km NW of Addis Ababa (Fig. 1), in the Gojam region. Mt. Choke (elevation 4052 m) stands ~1200 m above the surrounding flood basalts and the upper 300–400 m of the Mt. Choke sequence contains deposits that are more felsic than the surrounding flood basalts (Fig. 1). 40Ar/39Ar dating of a rhyolitic sample near the summit has provided an estimate of the age of the summit deposits of 22.4 (± 0.3) Ma (Kieffer et al. 2004). The Mt. Choke shield volcano is divided into six main deposits, differentiated by their igneous characteristics, from the summit to the base: Choke peak basalt, Choke rhyolite, Kutye basalt, Choke trachyte plug, Arat Mekeraker basalt and Rob Gebaya basalt (Fig. 1). The older, surrounding flood basalts are 25–29 Ma old (Hofmann et al. 1997), and are divided into four main deposits: Lumame, Debre Markos, Yejube, and Arero Gidabo basalts.

The soils near the summit of Mt. Choke are typically lighter and grayish brown in color, whereas soils at lower elevations of Mt. Choke and on the flood basalts are characteristically reddish brown. The soil units covering Mt. Choke have been classified as Haplic Alisols (deep soils with a predominantly clay or silt-clay texture), Eutric Leptosols (shallow soils with loam or clay-loam texture), and Eutric Vertisols (deep soils with clayey texture and angular/sub-angular blocky structure) (BCEOM 1999). The soil cover on Mt. Choke region is heavily influenced by morphological activity (e.g., erosion and land sliding), and calculations from remote sensing estimate that the soil thickness (i.e., depth to bedrock) varies from nonexistent to several meters (Teferi et al. 2010). The endemic natural forest cover in the Mt. Choke region has been predominantly replaced with grasslands and moorland, with patches of woody plant cover and ever-increasing plantations of Eucalyptus (Eucalyptus globulus). Isolated family-owned and worked fields are widespread, growing wheat, teff (Eragrostis tef), or other agricultural crops. Erosion rates are relatively high, and it is estimated that the farmed highland areas above the central and eastern parts of the basin are a major sediment source into the Blue Nile (Bewket and Teferi 2009). Soil degradation, both erosion and decreasing fertility, is a recognized problem in the Ethiopian Highlands (e.g., Van de Wauw et al. 2008). The majority (80–85%) of the Ethiopian population work within the agricultural sector (Alemayehu 2003), and the quality of the soil is continuously declining due to the intensification of cultivation driven by factors such as population increase and a need to increase food yields in times of drought (Devereux 2000).

The mineralogy and geochemistry of the surface soil can potentially offer an insight into the weathering process and various chemical alterations that have occurred to change the parent rock and regolith over time on the EVP. The parent material is found below the soil as unaltered, pristine rock. Increasing alteration is observed toward the upper soil layer, within the O and A horizons, typically subjected to the most intensive weathering processes. Geochemical-based weathering indices are commonly used to measure and compare the relative extent and intensity of soil pedogenesis based on the chemistry of the surface soils. In theory, the ionic potential of elements in the rock or soil determines whether the element is soluble and thus the extent to which the element will be mobile during weathering. Conventional chemical weathering indices (e.g., Chemical Index of Alteration = CIA; Nesbitt and Young 1982, Chemical Index of Weathering = CIW; Harnois 1988) are based on the chemical composition of the sample and aim to quantify the selective removal of mobile and soluble elements.

The use of standard weathering indices (CIA, CIW, etc.), however, may be restrictive and often not able to account for the changes occurring during real weathering scenarios (e.g., Duzgoren-Aydin and Aydin 2003; Price and Velbel 2003; Duzgoren-Aydin et al. 2002), as they do not consider the subtle variations (for example mineral or chemical heterogeneity) within the parent rock or the variable (often unpredictable) behavior of some components within the soil (for example iron). In situations where the parent rock is heterogeneous over a relatively small volume, conventional chemical weathering indices may not be sensitive enough to recognize subtle compositional differences within the soil that are associated with the parent rock. To derive a value to represent the degree of weathering within soils more accurately, Ohta and Arai (2007) used principal component analysis (PCA: a multivariate analysis tool, used primarily for exploratory data analysis), based on eight main oxides of both mobile and immobile elements (Al2O3, CaO, Fe2O3, K2O, MgO, Na2O, SiO2, TiO2). They derived a weathering index “W” for soils and sediments forming over igneous rocks. The W index has been shown to be a good indicator of the bulk geochemical alterations that occur during weathering, irrespective of igneous parent rock (e.g., Szilas and Garde 2013; Ohta and Arai 2007), which has the advantage of allowing comparative measurement of the degree of weathering between unrelated samples. In addition to igneous rocks, the W index has been applied to 3.75 Ga metamorphic mafic-ultramafic rocks with good results (Nuvvuagittuq Greenstone Belt, Hudson Bay, Canada; O’Neil et al. 2011). Calculation of the W index, however, has not proven to be entirely suitable for sedimentary parent rocks (e.g., Zhang et al. 2014), and further testing of this index is justified.

The main objective of this study was to investigate the extent of the weathering and alteration that has contributed to pedogenesis in one region of the EVP. Samples of surface soils were collected and the inorganic components were analyzed, using weathering indices, to reveal the processes of soil formation and extent of alteration over a range of elevations and differing parent rocks. As a second goal, the surface soils over the differing parent rocks were compared to identify whether the parent rock influenced the chemistry and mineralogy of the surface soils. To address these objectives, a new weathering index was developed, derived from the statistical theory given in Ohta and Arai (2007) although based on mineralogy rather than chemistry, which can be used to gain information on both the relative extent of weathering and the igneous nature of the parent rock. When investigating the formation of soil from igneous deposits the use of a weathering index based on soil mineralogy as well as the more conventional weathering indices derived from soil chemistry, can provide additional information that enhances our understanding of pedogenesis.

Experimental methods

Meteorological data

Meteorological data was sourced from http://www.WorldClim.org (accessed: February 2014; Hijmans et al. 2005), which generates data layers of interpolated average monthly climate information from weather stations on a 30 s resolution grid (from data gathered between 1950 and 2000). The mean annual precipitation (MAP) and mean annual temperature (MAT) were determined for every sample location using ArcGIS (ArcMap 10.1) (Table 1).

Sample collection

Surface soil samples (n = 135) were collected over the Mt. Choke region of the EVP, selected to best represent the elevation changes and major geologic formations of the region. Soils were collected from areas that are underlain by two main formations: the Mt. Choke shield deposits (Kutye, Arat Mekeraker, and Rob Gebaya) and the surrounding flood basalt deposits (Lumame, Debre Markos, and Yejube; Fig. 1). Soils collected from the Mt. Choke shield volcano were sampled at elevations 2874–3597 m, and those collected on the flood basalts at 2147–2943 m. At each site, approximately 1 kg of the surface soil (to a depth of ~5 cm) was collected using a plastic trowel and placed into sterile, airtight sample bags, then double bagged and labeled. The GPS waypoint and elevation was recorded for each sample site.

The soil was homogenized and sub-samples were separated at Addis Ababa University. From this point, samples were shipped to, processed, and analyzed at the Natural History Museum (NHM) in London (unless otherwise stated). Permission to export samples was obtained from Addis Ababa University and the Ethiopian Federal Ministry of Mines and Energy, under a DEFRA soil license granted to the NHM. Once in the U.K., further sub-samples of the soil were separated. Aliquots of the soil were dried at 60 °C for at least 48 h, lightly hand ground with a mortar and pestle, and sieved to <2 mm. Sub-samples for chemical analysis were ground in a Reich agate ball mill and sub-samples for X-ray diffraction (XRD) were further hand ground and passed through a 63 μm sieve.

In addition to the soil samples, whole-rock geochemical data from unaltered parent rock material from across the EVP, collected by Kieffer et al. (2004), were incorporated into the study to allow comparison of the weathered surface soils with parent rock samples. In this study we used data labeled “fresh” to “moderately weathered” (by description in Kieffer et al. 2004) that best represented the Mt. Choke shield and upper flood basalt deposits (the latter taken from the Lima Lion section) (for details see Appendix Table 11). Samples from Kieffer et al. (2004) were chosen to represent unaltered, pristine parent rock samples and will be referred to as “unweathered parent rock” hereafter. Whole rock geochemical data were also sourced from igneous “Certified Reference Materials” (CRMs) to compare the data that we collected with type-examples of igneous rocks. The data from Kieffer et al. (2004) and the CRMs used in the study are presented in Appendix Table 11, and plotted onto a TAS diagram in Appendix Figure 11 to ascertain their igneous origin.

Chemical analysis

The major elemental composition (oxides) of soil samples was determined (n = 135) using a lithium metaborate flux fusion and dissolution, followed by analysis using inductively coupled plasma–atomic emission spectroscopy (ICP-AES) (detailed in Le Blond et al. 2008). At least 1 blank was run for every 5 samples to check for background/cross-contamination. The accuracy and precision of the ICP-AES analysis was measured by regularly analyzing synthetic solutions containing known levels of the elements of interest from the following CRMs: BEN basalt, BHVO-1 basalt, BCR-2 basalt, AGV-1 andesite, GSP-1 granidiorite, JG-1 granidiorite, biotite, NIM-G granite, NIM-P pyroxenite, and STM-1 syenite. Limits of detection (LOD) and limits of quantification (LOQ; in parentheses) are as follows (in wt%): Al2O3 0.3(1.1); CaO 0.4(1.3); Fe2O3 0.3(1.0); K2O 0.02(0.06); MgO 1.1(0.4); MnO 0.004(0.01); Na2O 0.2(0.7); P2O5 0.01(0.04); SiO2 1.5(4.7); TiO2 0.3(1.0).

Mineralogy

X-ray diffraction (XRD) analysis enables identification and quantification of crystalline components in a bulk, powdered sample. Initially, for phase identification, the soil samples (n = 100) were close-packed into aluminum side-loading sample holders (to favor random orientation of the particles) and analyzed in an Enraf-Nonius X-ray diffractometer with an INEL curved position-sensitive detector (PSD). This apparatus collects data from 2–120 °2θ continuously throughout the measurement. The apparatus used CuKα radiation (λ = 1.5418 Å), at a voltage of 45 kV and current of 40 mA, with a primary slit of 0.1 mm by 5 mm. Samples were analyzed for 60 min.

The phases identified in the soils were quartz, feldspars, poorly crystalline Fe oxides, smectite, mica, chlorite, amorphous silica, and kaolin-group minerals. The quantification method is detailed in Cressey and Schofield (1996) and Cressey (1999), but briefly, standard mineral samples (from the collections at the NHM) representing those in the soils, were analyzed under identical preparation and instrument conditions. These standards were: a mixture of anorthite, albite, and orthoclase in equal proportions to represent the feldspars; goethite for the poorly crystalline Fe oxides; quartz; and volcanic glass for the amorphous silica phase. The proportion of each of the above mineral phases was determined by matching the contributions of the mono-mineralic standards to the experimental patterns and then correcting the values using the individual mineral X-ray absorption coefficients. The sum of the phyllosilicates (smectite, kaolinite, mica, and chlorite) was determined by difference. The proportion of the individual phyllosilicates was investigated in oriented mounts (see below). Repeated preparation and analysis of a sub-sample set (n = 15) produced quantification values within 5% for all mineral phases.

To investigate the octahedral configuration of the clays (trioctahedral vs. dioctahedral), the 060 peaks were studied. The samples used in the bulk identification XRD analysis were analyzed in the 59–64 °2θ angle range on a PANalytical X’Pert-PRO scanning XRD (240 mm radius) with a step size of 0.0167 °2θ, for a total counting time of 60 min. CuKα radiation was used, with a voltage of 45 kV and current of 40 mA. The incident beam passed through a 0.02 rad Soller slit, a 1° divergence slit, a 10 mm fixed mask and a 2° fixed anti-scatter slit. The instrument was fitted with a Ge monochromator and the diffracted beam was detected by an X’Celerator RTMS detector.

For the detailed investigation of the phyllosilicates, an aliquot of each of the samples was prepared as an oriented mount on a glass slide. This technique encourages preferred orientation of the clay particles parallel to the glass slide surface and enhances the intensity of the basal reflections (00l). A suspension of soil was prepared (4 mL deionized water was added to 80 ± 0.003 mg sample), sonicated for 1 min and 0.5 mL dropped onto a clean glass slide and left to dry. The PANalytical X’Pert-PRO diffractometer was used, set up as in the analysis of the 060 peaks, but the scan range was 2–60 °2α, the step size 0.0334 °2α and a total counting time of 60 min. Once the slides had been analyzed in air-dry conditions, they were glycolated (in a desiccator with ethylene glycol; in an oven at 60 °C for 24 h) to investigate the presence of swelling clays (smectite or vermiculite). The glycolated slides were analyzed using identical conditions as adopted for the air dried oriented mounts.

The XRD patterns from the glycolated samples were used for modeling and quantifying the phyllosilicate phases within the samples, using CLAYSim (version 1.1.6) from MDI. This program calculates the profiles of the 00l peaks of mixed-layer clays with two components, using chemical and experimental input variables. Up to 10 phases can be added to simulate mineral mixtures within the sample. The calculated patterns were matched to the experimental patterns using the whole pattern profile: peak position and shape, as well as background. The glycolated patterns were chosen for this investigation because their peak positions are more reliable than those from air dried patterns due to changes in relative humidity.

Scanning electron microscopy (SEM)

To further investigate the mechanisms of mineral formation within the soils, two samples (SO 5137 and DM-CM SO 7A) were selected that had high quantities of kaolinite and amorphous silica (alteration products). A LEO VP 1455 scanning electron microscope (SEM) with energy-dispersive X-ray (EDX) detector (Oxford Instruments X-Max detector) was used. The samples were set in a resin block, polished, and carbon coated (~12 nm thickness). Samples were analyzed in high vacuum (15 kV, working distance of 12 mm) with backscattered electrons. EDX analyses were acquired for 1 min and the spectra generated were used to identify the mineral phases.

Statistical data processing

The bulk mineralogy and major oxide composition of the soil samples were compared using the non-parametric statistical analysis test Spearman’s Rank correlation. A preliminary assessment of the data was carried out using multivariate regression analysis, to determine whether the individual mineral phases within the soils could be represented as functions of the bulk chemical oxides (e.g., mineral = a+b[oxide1]+c[oxide2]+d[oxide3].... etc.). Each individual mineral phase was assessed as a dependent variable using the chemical composition as predictor variables; collinear variables were excluded from the statistical analysis to avoid the problems of multi-collinearity in the regression analysis.

Two indices of chemical weathering, CIA (Al2O3 × 100/[Al2O3 + CaO* + Na2O + K2O]; Nesbitt and Young 1982) and CIW (Al2O3 × 100/[Al2O3 + CaO* + Na2O]; Harnois 1988); CaO* denotes the CaO exclusively from silicate minerals and can be applied in samples where the molar ratio of CaO/Na2O does not exceed 1; McLennan 1993), were calculated from the chemical composition of the soils. In addition, the W index from Ohta and Arai (2007) was calculated from the chemical compositional data of our soil samples and from the data provided by Kieffer et al. (2004). As a further step, we adapted the statistical framework from Ohta and Arai (2007) to calculate a W index value from the soil mineralogy, which we refer to as the “Wmin” index. The data were first centered log ratio (clr)-transformed (a mathematical function applied to each data point to ensure that the data meet the assumptions of the statistical inference procedure to be used; (Aitchison 1986) using the “Hotelling” package R statistical computing software. Zero values were altered to 0.05 to represent the maximum potential error in XRD quantitative measurements. The IBM SPSS statistics program (versions 20 and 21) was used for the PCA calculations. The data from the PCA were used to create a ternary plot (using the 3-axis rule of Tolosana-Delgado et al. 2005) that identifies and distinguishes the geochemical variations in the parent rock, in addition to the extent of weathering.

The degree to which the soils were influenced by the parent rock was also investigated by comparing the chemistry, mineralogy, and derived-weathering indices of the soils collected over the two parent rocks: the Mt. Choke shield basalts and surrounding flood basalt, using a non-parametric Mann-Whitney-U test.

Normative mineralogical data for the unweathered parent rock samples collected on Mt. Choke and the CRMs were not available, and therefore whole rock mineralogy was determined using Cross, Iddings, Pirsson, and Washington (CIPW) norm calculations (Kelsey 1965; Johannsen 1931; Cross et al. 1902). CIPW norm calculations work on the basis that a defined set of minerals (mode) within volcanic rocks can be used to quantitatively estimate the mineralogical composition based on total elemental analysis and defined mineral formulas (Appendix Table 11). The samples used for the CIPW norm calculations, were considered to be fresh or had only minor alteration.

Results

Climate

There was a negative correlation between MAP and MAT (r = −0.85, p = 0.00; Appendix Table 21). On average MAP in the sample sites was 1446.5 mm/yr, and MAT was 15.3 °C. Elevation correlated negatively with MAT (r = −0.82) and positively with MAP (r = 0.83).

Soil chemistry

The values for the major element chemical analyses and quantification of the mineral phases within the soil are shown in Table 1, and Spearman’s Rank correlations between the data are shown in Appendix Table 21. Strong correlations (where r ≥ 0.5 and p ≤ 0.01) involving wt% Fe2O3 were observed, negatively with SiO2 and K2O (r = −0.71 and −0.82, respectively) and positively with TiO2 (r = 0.77). In addition, Na2O3 correlates negatively with Al2O3 (r = −0.82) and positively with MgO (r = 0.53), while K2O correlates negatively with TiO2 (r = −0.66) and positively with SiO2 (r = 0.71). Finally, SiO2 correlates negatively with TiO2 (r = −0.62). The correlations observed in the surface soils approximately reflect the chemical characteristics of the volcanic parent rock, where the proportion of elements more commonly found in mafic rocks (Fe, Mg, Ca, Ti) were generally negatively correlated with those found in more felsic rocks (Si, Na, K).

SiO2 negatively correlated with MAP (r = −0.52) and MgO negatively correlated with MAT (r = −0.656) (Appendix Table 21). No other bulk oxide was found to have a strong correlation (where r ≥ 0.5) with either temperature or precipitation (MAT or MAP).

The data from the characterization of the soil samples were grouped according to the parent rock (i.e., Mt. Choke shield volcano and flood basalt). A non-parametric Mann-Whitney U test was used to determine whether the two sets of data were significantly different (where p <0.05). The outcomes of the statistical analysis are detailed in Table 2, which shows that the soils formed on the Mt. Choke shield volcano were significantly different in their chemistry (for all elements) compared with those formed on the flood basalt.

Soil mineralogy

The soils were generally iron-oxyhydroxide-rich (goethite; 16–40 wt%; Table 1), with varying proportions of amorphous silica (9–26 wt%; indicated by the characteristic XRD pattern showing a “hump” centered at ~21 °2θ; e.g., Fig. 2a), quartz (0–45 wt%), feldspars (either in concentrations 11–30 wt% or completely absent) and phyllosilicates (totaling 9–57 wt%). Other minor phases were occasionally observed in the XRD identification, for example pyroxene was observed and some minor crystalline iron oxide phases (e.g., hematite; Fig. 2b), but at levels <<5 wt%, and hence they were not included in the quantification. The phyllosilicates included mica, chlorite, kaolinite, and kaolinite-smectite (K-S), the latter with either high smectite (50 to 70%) or high kaolinite (95 to 97%) contents. The high smectite K-S was detected in the diffraction patterns of the oriented mounts, where the air-dry samples displayed a broad peak at 14.3 Å that was displaced after glycolation to 17.4 Å (Fig. 2c).

The 060 peak is broad and centered ~1.49–1.50 Å (Fig. 2b) confirming that the kaolinite and K-S were dioctahedral (octahedral cation abundance is Al >> Fe, Mg). The broad 060 peak may also mask the mica (dioctahedral) peak at 1.50 Å (similarly the quartz peak at 1.54 Å overlies and hides the influence of the chlorites (at ~1.53–1.55 Å).

The K-S phases were divided into their kaolinite and smectite components to produce an overall quantity of each layer type within the soil. For example, the amount of kaolinite is the sum of the percentage of the discrete kaolinite plus the proportion of kaolinite layers in the K-S (mixed-layer) phases. Figure 2d provides an example of match between a calculated and experimental XRD pattern of an oriented, glycolated mount preparation. The 001 peak of the calculated XRD pattern is displaced to a lower angle due to the combined effect of the low intensity of the peak and the steep background (Fig. 2d). The ClaySIM program, used in the quantification of the phyllosilicates, typically produces a steeper pattern background than the experimental pattern in the low-angle region (Plançon 2002). The simulations, however, rely on the position of all 00l peaks in the experimental patterns and the calculated pattern in Figure 2d can be considered to be a true representation of the experimental pattern.

Significant negative correlations were observed in the soils between the quantity of quartz and amorphous silica (r = −0.54; Appendix Table 21). A negative correlation between the amount of mica and kaolinite (r = −0.61) was also noted. Chlorite was found to positively correlate with both MAP and MAT (Appendix Table 21). No other mineral phase was found to have a strong correlation (where r ≥ 0.5) with either MAT or MAP.

The specific alteration products of some of the most abundant minerals in the original rock and the potential mechanisms of clay formation were further investigated using SEM. Specifically, two issues were explored: first, the alteration process by which kaolinite phases were formed; second, the form of the amorphous silica within the soil. For the SEM work, two soil samples (SO 5137 and DM-CM SO 7A) with high quantities of kaolinite and amorphous silica were fixed in resin and polished. Kaolinite was observed as an alteration product of both feldspars and perhaps quartz, forming distinct weathering rims. An example of a quartz grain is shown in Figure 3, where the outer rim of the quartz grain has been transformed into kaolinite. The similarities between the shapes of the grain and the coating also point to this theory of transformation rather than an inherited coating, however, more work is need to confirm this. Kaolinite-rich phases were ubiquitous within the soil samples analyzed by SEMEDS, but in most cases the K-S was found in the soil matrix and there was no preserved relation to the original mineral grains. In addition, no discrete regions of pure amorphous silica were observed in the resin soil blocks under the SEM (with EDS). It is therefore surmised that the amorphous silica must occur as fine groundmass (intergrown with clay minerals and Fe oxides) in the soil, which is in agreement with other studies (e.g., Singh and Gilkes 1993).

In terms of the soil mineralogy over the two different parent rocks, the quantity of iron oxides, quartz, amorphous silica, kaolinite, and chlorite was significantly different in soils over the flood basalt when compared with the Mt. Choke shield (Table 2). However, the content of feldspars, smectite, and mica measured in the soils were not significantly different over the two parent rocks.

Multivariate regression analysis was used to investigate associations between bulk chemical composition of the soils and their mineralogical phases (r ≥ 0.5, p ≤ 0.01 was considered significant). The most significant relationship was observed between the Fe2O3 concentration and quantities of FeOOH, mica, and chlorite from XRD quantification (R2 = 0.754, R2adjusted = 0.746):

 

wt%Fe2O3=3.235-(0.021×wt%mica)+(0.077×wt%chlorite)+(0.529×wt%Fe oxide)
(1)

Geochemical weathering indices

The CIA values for the soil samples varied in the range 68.4–96.7 (n = 135, average = 87.4), and the CIW values were 78.7–98.5 (n = 135, average = 93.1) (not shown). The average CIA and CIW values calculated for the unweathered parent rock were 40.8 and 41.8, respectively. The average CIA and CIW values calculated for the CRMs were 37.2 and 39.8, respectively.

The CIA and CIW values correlated positively (r = 0.75; Appendix Table 21), and the CIA correlated positively with kaolinite (r = 0.52) and iron oxide (r = 0.59). Both CIA and CIW correlated positively with Al2O3 and Fe2O3, and negatively with Ca-, Na-, and K-oxides.

The W index calculations (derived from Ohta and Arai 2007) were plotted onto a ternary diagram, to represent the three latent variables: mafic source, felsic source and relative extent of weathering (shown as the three vertices M-F-W on the ternary plot in Fig. 4). The majority of the soil samples concentrate in the region of the M-F-W plot that indicates high weathering, derived from basaltic rocks. The parent rocks and CRMs plot between the M-F vertices, the solid arc on which they plot indicating no alteration. The unweathered parent rock, however, show that there is variation in the chemistry of the samples collected from the flood basalt and Mt. Choke regions. The soils collected above the Kutye basalt deposit (near the summit of Mt. Choke shield) plot relatively closer to the “F” vertex, however, the majority of the soil samples are clustered together and it is challenging to resolve the detail in the mass of points at the “W” apex.

Deriving the Wmin index

PCA of the mineralogy data (Table 3) suggests that the first principal component (PC1) is positively correlated with the amorphous silica, kaolinite, amorphous iron oxides, mica, and smectite, and negatively with feldspars, crystalline iron oxide, pyroxene, and olivine. The second principal component (PC2) shows a strong positive correlation with chlorite, a weak negative relationship with olivine, and also a weak but positive correlation with quartz. PC1 and PC2 represent 70.1 and 13.7% of the variance (respectively) and collectively the three components explain 84.8% of the variance within the samples. The resultant biplot using PC1 and PC2 (e.g., Aitchison and Greenacre 2002) is shown in Figure 5. The values for the individual soils are plotted as discrete points on the biplot, and the “mineral” phase loadings are represented by the arrows radiating from the central point of the biplot (representing the PCA loading scores in Table 3). The mineral phases (variables) are represented in this bi-plot by a vector, and the direction and length of the vector indicate how each variable contributes to the two principal components in the plot. Values for the mineralogy of the CRMs and unweathered parent rock (Appendix Table 11) are from CIPW norm calculations, and are plotted on the biplot, in addition to the soil samples, for comparison. The unweathered parent rocks and CRM samples have negative PC1 scores (Fig. 5), whereas the soil samples have comparatively more positive PC1 scores. There is also separation of the samples according to the PC2 scores, from those samples that are smectite-rich (with negative PC2 scores), to those that are quartz and/or chlorite-rich (with positive PC2 scores).

The orthogonal components (calculated from PCA), representing the main controls on the mineralogy of the soils (PC1 and PC2, Fig. 5), were reformatted to plot on a ternary diagram using Equation 2 (see Ohta and Arai 2007, for further details):

 

V1=exp(-16PC 1+-12PC 2),V2=exp(-16PC 1+-12PC 2),V3=exp(26PC 1)
(2)

Then, V1, V2, and V3 were normalized and used to create the “mafic”, “felsic,” and “Wmin” vertices in the ternary diagram (Fig. 6a). The unweathered parent rock and CRMs, plot along a line delimiting the furthest distance from the “Wmin” apex on the ternary diagram, next to the mafic-felsic joining line (Fig. 6a). In a similar way to the W index, the parent rock samples plot along the range of mineralogy delimiting felsic and mafic samples accordingly. The soil samples collected in this study plot as an arc progressing away from the felsic-mafic axis toward the “Wmin” apex (Fig. 6a). The extent of weathering within the soils is shown by the variable distance between the “Wmin” apex and the datapoints. It is worth noting that the scales on the ternary diagrams used to represent the W (Fig. 4) and Wmin (Fig. 6a) indices do not represent a complete scale from 0 to 100%, but are specific to the compositional and weathering ranges included in the analysis.

The Wmin index correlates positively with MAT (r = 0.57), and negatively with both MAP and elevation (r = −0.5 and r = −0.58, respectively; Appendix Table 21).

Discussion

Soil formation

Pedogenesis occurs over time to alter the chemistry and mineralogy of the parent rock to create soils that contain mineral phases considered stable at the Earth’s surface. The surface soils characterized in this investigation were composed of minerals that typically result from the alteration of basaltic rock over time, such as iron oxides, phyllosilicates, and amorphous silica. During weathering, the most mobile elements (Mg, Ca, Na, and K) become depleted via various processes including congruent and incongruent dissolutions, oxidation, and hydration reactions. The mechanisms occurring during weathering that alter the primary minerals to secondary minerals may be either sub-solidus replacement or dissolution-precipitation, and can form either crystalline structures or non-crystalline phases (e.g., Viti et al. 2007). The breakdown of the primary minerals in our soil samples (pyroxene, olivine, feldspars) and volcanic glass has resulted in the loss of non-hydrolyzing cations (at least Ca2+, Na+) and the accumulation of Si4+, Al3+, and Fe3+ cations preserved in crystalline (e.g., kaolinite) and amorphous/poorly crystalline phases (e.g., FeOOH). Hence, with increasing weathering the total amount of phyllosilicates and proportion of Fe oxides increases within the soil (Delvaux et al. 1989).

Kaolinite and mixed-layer K-S clays were present in most of our soil samples. These minerals have been found to be a common constituent in soils with a basaltic parent rock across tropical, subtropical, and semi-arid regions (Vingiani et al. 2004; Righi et al. 1999; Delvaux et al. 1990). Primary minerals in igneous rocks eventually weather to clay minerals (Taylor and McLennan 1985; Taylor et al. 1983) and in regions with high annual precipitation, as in our study area, kaolinite is typically the final product of weathering. The soils analyzed in our study contained both a proportion of discrete kaolinite and K-S that was either kaolinite-rich (with 95 to 97% kaolinite layers) or smectite-rich (with 50 to 70% smectite layers). No K-S mixed-layer phases were found with intermediate kaolinite contents (i.e., between 50 and 94%). As smectite transforms into kaolinite in highly weathered soil environments, one would expect to find the entire range of K-S compositions present in our soils (e.g., Herbillon et al. 1981). Our modeling procedure does not totally exclude the existence of K-S of “intermediate” composition; however, if this phase is present, it is in much lower concentrations (relative to pure kaolinite, or kaolinite- and smectite-rich K-S forms). We conclude that this gap in K-S composition indicates a degree of alteration heterogeneity within the soils that needs to be explored further.

The multivariate analysis, comparing the soil mineralogy with the bulk chemistry, indicated that the majority of the iron in the soils is held within the amorphous/poorly crystalline iron oxide phases. There was also a positive correlation with chlorite indicating that the chlorite contains iron. The negative correlation of the mica in the regression equation implies that the mica must be in a form that is Al rich and Fe poor. When kaolinite and smectite were included in the regression analysis (Eq. 2), their contributions to the R2 value were insignificant. We recall that the kaolinite and smectite included the proportion of layers in K-S as well as end-member kaolinite; i.e., smectite and kaolinite indicates layers in K-S to a large extent. The lack of correlation between kaolinite and iron was surprising, because one would have expected a significant negative correlation between them. Our interpretations are that both smectite and kaolinite layers contain iron and that the quantity of iron in these clays must be approximately proportional to the iron in the glass and primary minerals from which they formed. Hence, no correlation emerged between these minerals and iron content. The incorporation of iron into smectite, at variable proportions, is well known; however, it is more difficult to justify the presence of iron in kaolinite. Nonetheless, the presence of iron was confirmed within kaolinite, as SEM-EDS analysis consistently showed iron within the spectra of kaolinitic grains (e.g., Fig. 3). The incorporation of Fe3+ ions into the octahedral sites of kaolinite has previously been documented (Weaver et al. 1967; Malden and Meads 1967), and Iriarte et al. (2005) synthesized kaolinite with up to ~1/3 of octahedral sites occupied by iron (Si2Al1.4Fe0.6O5[OH]4). In addition, in our soil samples the majority of the kaolinite was found in a mixed-layered (K-S) form, where kaolinite layers have been shown to incorporate significant amounts of iron (Dudek et al. 2006; Ryan and Huertas 2009).

We propose that the kaolinite- and smectite-rich K-S phases are the products of alteration from different mineral phases in the parent rock. More specifically, feldspar would alter into kaolinite and kaolinite-rich K-S, while volcanic glass and perhaps additionally olivine and pyroxene, would alter to smectite (which would subsequently alter to K-S of progressively kaolinitic composition). These processes were compatible with our SEM observations. In this study, kaolinite and kaolinite-rich K-S was observed surrounding grains of feldspar and less frequently, but notably, around quartz (Fig. 3), in what appeared to be alteration rims. This observation is of interest because, in systems similar to this, smectite is perceived to be an intermediate stage of weathering toward K-S that ends in the formation of kaolinite. However, our observations suggest that K-S can also form without a smectite precursor. The alteration of quartz to kaolinite (or kaolinite-rich K-S) may be made possible by the direct contact with other mineral grains containing Al, also in the process of weathering or dissolution. The smectite-rich K-S was predominantly found in the soil matrix and without spatial or morphological reference to original primary minerals (i.e., relic features) and hence may indicate that smectite-rich K-S formed mainly from the alteration of volcanic glass in the parent rock matrix.

Weathering indices

Weathering indices are conventionally used to quantify the extent of weathering within a sample soil (or other unconsolidated sediment or rock) and help to explain processes occurring during pedogenesis. CIA values calculated for pristine basalts collected on the Deccan Traps in India were ~35 (Babechuk et al. 2014), whereas the CIA values for residual clays lie between 85 and 100 (Taylor and McLennan 1985) and pure kaolinite approaches 100 (Nesbitt and Young 1982). The CIW values calculated for pristine basalts collected in Iceland were between 37 and 38, and andosol soils (developed over volcanic rock) ranged from 37–45 for weakly weathered soils to 50–77 for more strongly weathered soils (Óskarsson et al. 2012). Similarly, CIW values for strongly weathered volcanic soils (laterites) developing in warmer climates in Brazil and China have been calculated to be between 70 and 100 (e.g., Ma et al. 2007; Truckenbrodt et al. 1991). The CIA and CIW values calculated in this study reveal that the extent of weathering within the soil samples is relatively high (CIA average = 90.5, CIW average = 96.1). The unweathered parent rock samples and the CRMs had significantly lower values for CIA and CIW, of <42, which indicated that they effectively represent the unaltered, parent rock.

The ternary diagram for the W index shows the chemical composition of the parent rock samples, relative to the soil samples (Fig. 4). The CRMs plotted as expected, with the rhyolite RGM-1 (Appendix Table 11), plotting near the “F” (felsic) axis, the basaltic samples plotting near the “M” (mafic axis) and the andesite AGV-1 midway between the “F” and “M” axis. The unweathered parent rocks, however, have varied compositions: Almost all of the rock samples collected on the flood basalt plot near the “M” axis, with the exception of sample 139, which plots near to the “F” axis. However, this sample was classified as “trachyte” and was recognized as being more felsic than the other samples from the flood basalt deposits. The majority of the parent rock samples collected on the Mt. Choke shield were classified as “basalts” and plot accordingly near to the “M” axis on the ternary diagram, with the exception of sample 240. This sample was identified as a “trachyandesite” and plots accordingly on the ternary plot, indicating that there is a distinct felsic character to the deposits on the upper elevations of Mt. Choke. The chemical characteristics of the soil reflect that there may be an influence by the felsic parent rock because the soils collected on the Kutye basalt at the summit of Mt. Choke (~3700 m) plot relatively closer to the “F” axis. In addition, these soils from the summit of Mt. Choke appear to be less weathered (i.e., further from the “W” axis) when compared with the other soil samples.

Information from the Wmin index

The biplot calculated for the Wmin index (Fig. 5) suggests that PC1 corresponds to the degree of weathering, which was defined as the amount of feldspars, crystalline iron oxide, pyroxene, and olivine, vs. the amount of amorphous silica, amorphous iron oxide and phyllosilicates within the sample. This is exemplified by the placement of the parent rock samples relative to the soil samples along the PC1 axis. Mica is not expected to be a weathering product but perhaps its resistance to weathering makes it more concentrated in the soils and thus positively correlated with weathering. One other, alternative, explanation for the positive correlation with weathering is that here we did not differentiate between mica and illite (the latter being a product of weathering in soils).

Chlorite had a strong positive loading on the PC2 axis and quartz a weaker positive loading. Olivine, and to a lesser extent smectite, had a moderate negative loading on the PC2 axis. Chlorite is a group of phyllosilicates that can be formed by hydrothermal alteration of igneous rocks and the alteration of ferromagnesian minerals (Meunier 2005). Chlorite was present in all the soil samples collected on the Kutye and Arat Makeraker basalts deposit, had a negative correlation with elevation, but was only sporadically present in the soils collected at lower elevations. In addition, the amount of chlorite was also negatively associated with MAT and positively correlated with MAP.

When considering the placement of olivine and quartz, the PC2 axis is likely to represent the chemistry of the parent rock composition, with a mafic to felsic transition from negative to positive PC2 values. If the unweathered parent rocks are taken as an example, negatively loaded PC2 samples have increased amounts of olivine and are from a more mafic rock, whereas positively loaded PC2 samples have more quartz and are more felsic. It must be recognized, however, that PC2 does not represent a completely regular scale from mafic to felsic, as the more mafic parent rocks extend from PC2 = −2.5 to 0.5, and the felsic samples plot in the range = >0.5. The very high loading of chlorite in PC2 may be related to chemistry (Al-rich) or to the higher efficiency of the more felsic magma in producing chlorite because of the higher water content. The soils showed a wide degree of separation, particularly between those collected on the Kutye and Arat Makeraker basalts and the Yejube basalt deposits. The soils collected on the Rob Gebaya, Lumame, and Debre Markos basalts, however spanned the range of the PC2 axis, which is representative of the wide variation in parent rock composition. The degree to which the soils have been weathered (calculated as “Wmin”) was correlated with the increase in mafic mineralogy of the soil (Fig. 6b), and Wmin correlated positively with the amount of smectite (r = 0.55).

Although the PCA has differentiated the samples according to both the extent of weathering and the nature of the unweathered parent rock, we must consider that mineralogy of the parent rocks was derived from CIPW norm calculations and not actual measurements (for example quantification by XRD). It must be acknowledged that the CIPW norm calculations have a distinct error when used in this context, as the volcanic glass likely present in the unweathered parent rock matrix is not accounted for. This would influence the amorphous silica content, and potentially reduce the loading in the positive region on the PC1 axis. Mineral phases quantified by the CIPW norm calculations, for example pyroxene and crystalline Fe oxide phases, were observed in the SEM and XRD analysis (respectively) of the soil samples. Hence, the mineral phases predicted by CIPW are likely to be present in the original parent rock and thus the CIPW norm calculations have provided a good estimate of the rock mineralogy. It is also worth noting that the iron oxide derived from the XRD quantification was predominantly in a poorly crystalline form. However, we observed trace (i.e., <<5%) amounts of crystalline iron oxide, hematite, in the XRD patterns of the soil samples (Fig. 2b), that is negligible in the quantification of iron phases that was carried out.

Comparison of weathering indices

In principle, a correlation is expected between the chemical weathering indices and the proportion of secondary minerals formed by weathering. In this study, there was a positive correlation of CIA values with both amorphous Fe oxide (r = 0.59) and kaolinite (r = 0.52) (Appendix Table 21), which are both weathering products. However, no significant correlation existed between either the CIW or W index and mineral content.

In addition to the extent of weathering, the W and Wmin index ternary diagrams also discriminated between samples according to the nature of the parent rock from which they originated. There was a good correlation between the extent of weathering and the nature of the parent rock for the soil samples analyzed in this study (r = 0.62, Fig. 6b). Soil samples that were more felsic, in terms of their mineralogy, were comparatively less weathered. This compositional effect on weathering extent can also be observed in the ternary plots (Figs. 4 and 6a). For example, the Kutye basalts were relatively more felsic than the other basalts and the corresponding soils were shown to be comparatively less weathered (Figs. 4 and 6a). This is an expected result because minerals in progressively more felsic rocks, such as quartz and mica, are more resistant to weathering than those typically found in more mafic rocks (pyroxene and amphibole); likewise K-feldspar is more resistant than Na and Ca feldspars (e.g., Retallack 2001; Albarede and Michard 1986; Holland 1984; Edmond et al. 1979). Small chemical differences in magma can result in significantly different mineral components within the parent rock, which impacts the overall resistance to weathering. Hence, a weathering index based on mineralogy, has the advantage over a chemical index that it will more readily reveal the cause of differential weathering linked to mineralogy.

Comparison of the Wmin and W indices shows the following differences. (1) In the W index ternary diagram (Fig. 4), the unweathered parent rocks plotted along a concave arc spanning the mafic-felsic apices and extend into the plot, whereas in the Wmin index diagram (Fig. 5) these same data points plotted next to the straight line joining the mafic-felsic apices. (2) The pathways of weathering denoted on the W index diagram (dashed lines, Fig. 4) are not necessarily intuitive and need to be identified for each parent rock type, whereas the weathering pathways in the Wmin index ternary diagram can be determined by the relative distance from the W apex. From these two observations the Wmin index configuration is more intuitive as the parent rocks are at the furthest distance away from the W axis, and soil samples simply plot relative to each other without conforming to a pathway that must be specified. (3) The data points in the W index plot are compressed toward the W axis, whereas the plotting range is extended on the Wmin plot. Thus, differences between samples may be lost or masked when using the W plot, whereas the Wmin index ternary diagram may reveal comparatively greater resolution within sample data points.

The W and Wmin weathering indices provide a method of establishing both the degree of weathering, and an indication as to the nature of the parent rock from which the soil sample was formed. The CIA and CIW chemical indices of weathering, when used alone, can provide an estimate of the extent of weathering but details are lost within the complexity of the weathering processes. This is in general agreement with the conclusions from others (e.g., Price and Velbel 2003; Duzgoren-Aydin et al. 2002), who have questioned the validity of using a selective array of elements to explain the complex processes occurring within weathering profiles. Although it is often less time consuming to determine the bulk chemical composition and calculate the W index for a suite of samples, we would additionally advocate the use of the Wmin index.

Implications

The extent of weathering within the soils, and samples within a soil profile, is traditionally quantified by chemical indices based on the relative proportions of common immobile to mobile cations (CIA, CIW). An advancement on these chemical indices is the multivariate statistical analysis of the entire cation suit that has been used to derive the W index (detailed in Ohta and Arai 2007), which not only calculates the relative extent of weathering but also uses the data to identify the composition of the parent rock. Relating the surface soil to the parent rock from which it was formed is a valuable addition in understanding the evolution of the soils. Following the multivariate mathematical approach, we derive a new method of using mineralogy to determine the degree of weathering within soils forming over igneous parent rocks and the Wmin index, which also attempts to identify the nature of the parent rock. In this study, we noted that the soils appear to maintain a weak signature of the composition of the original magmatic rock from which they were formed. We identified the benefits of calculating the Wmin index and used this index to explain the mineralogical processes that have occurred within the soil profile to create the surface soil over time. Briefly, the Wmin index data points plotted in a more intuitive manner and with greater resolution than the W index. Nonetheless, the W index is possibly more accurate, due to the greater limitations of mineral phase quantification and is potentially less time-consuming. However, when the W index and Wmin index are used in conjunction, information on both the geochemical and mineralogical processes that occur during the weathering of an igneous parent rock can be acquired.

To better understand and explain the processes occurring during soil formation over igneous deposits, we advocate the use of multivariate statistical analysis of soil data to compare the relative weathering extent between soil samples. Although there are other factors that will inevitably influence the characteristics of soils, such as erosion and aeolian deposit, the W and Wmin indices provide a single value for a surface soil that identifies the comparative extent of weathering, and if quantitative data on soil geochemistry or mineralogy is available, then calculation of the W and Wmin indices can be carried out retrospectively.

The analysis techniques outlined in this paper can be used on various soil samples, and can therefore provide a protocol that can be employed over wide study areas and on a global scale. The aim is to utilize this methodology for surface soil characterization to rapidly compare properties where there is pre-existing data on the chemistry and mineralogy of the soil. The W and Wmin indices can be used across other flood basalt provinces as well, to test the sensitivity of the indices in differentiating soils from parent rocks with subtle differences in their mineralogy. An additional suggestion is to utilize this methodology to compare the soils forming over very different igneous deposits (ranging from ultramafic to felsic) to better understand the impact that weathering has over different time-scales and the regeneration of soil profiles over time. Single values, obtained from the indices, can be mapped spatially and coupled with other categories of information (for example water flow direction and accumulation, climatic, vegetation, from in situ measurements or remote sensing) to conduct comprehensive studies of soils within different environments and weathering regimes. This has implications for agriculture and mineral prospecting, and can potentially be applied over large land areas.

Acknowledgments

This work was funded by a grant from the Wellcome Trust University Award (No. 091956) held at BSMS. The authors thank E. Humphreys-Williams for help with the ICP-AES data analysis and J. Najorka for XRD advice, both at the NHM, London. Thanks also to M. Riishuus for insightful comments on the manuscript in an earlier form.

Open access: Article available to all readers online.
1
Deposit item AM-15-115168, Appendix Figures and Tables. Deposit items are stored on the MSA web site and available via the American Mineralogist Table of Contents. Find the article in the table of contents at GSW (ammin.geoscienceworld.org) or MSA (www.minsocam.org), and then click on the deposit link.