Abstract

Soil quality (SQ) degradation continues to challenge sustainable development throughout the world. One reason is that degradation indicators such as soil quality index (SQI) are neither well documented nor used to evaluate current land use and soil management systems (LUSMS). The objective was to assess and identify an effective SQ indicator dataset from among 25 soil measurements, appropriate scoring functions for each indicator and an efficient SQ indexing method to evaluate soil degradation across the LUSMS in the Mai-Negus catchment of northern Ethiopia. Eight LUSMS selected for soil sampling and analysis included (i) natural forest (LS1), (ii) plantation of protected area, (iii) grazed land, (iv) teff (Eragrostis tef)-faba bean (Vicia faba) rotation, (v) teff-wheat (Triticum vulgare)/barley (Hordeum vulgare) rotation, (vi) teff monocropping, (vii) maize (Zea mays) monocropping, and (viii) uncultivated marginal land (LS8). Four principal components explained almost 88% of the variability among the LUSMS. LS1 had the highest mean SQI (0.931) using the scoring functions and principal component analysis (PCA) dataset selection, while the lowest SQI (0.458) was measured for LS8. Mean SQI values for LS1 and LS8 using expert opinion dataset selection method were 0.874 and 0.406, respectively. Finally, a sensitivity analysis (S) used to compare PCA and expert opinion dataset selection procedures for various scoring functions ranged from 1.70 for unscreened-SQI to 2.63 for PCA-SQI. Therefore, this study concludes that a PCA-based SQI would be the best way to distinguish among LUSMS since it appears more sensitive to disturbances and management practices and could thus help prevent further SQ degradation.

1. Introduction

Globally, declining in soil quality (SQ) has posed a tremendous challenge to increasing agricultural productivity, economic growth, and healthy environment [1, 2]. The underlying causes for SQ degradation are largely related to inappropriate land use and soil management, erratic and erosive rainfall, steep terrain, deforestation, and overgrazing [24]. Most of the causes are resulted from a desperate attempt by farmers to increase production for the growing population which aggravate SQ degradation more in the developing countries, which mainly depend on natural resources (agriculture) [1, 4]. Misuse of natural resources that leads to degradation can also be stimulated by socioeconomic and political issues, for example, land tenure, capital, and infrastructure [5]. SQ degradation by soil erosion such as soil nutrient depletion and changes in soil physical indicators is largely recognized as a principal cause aggravated by the effect of inappropriate land use and soil management in the developing countries like Ethiopia [6, 7].

Interest in the evaluation of soil degradation particularly the quality of soil resources has been increased as soil is critically important component of the Earth’s biosphere, functioning not only in the production of food and fibers but also in the maintenance of environmental quality [8]. In normal conditions, the soil can maintain equilibrium by pedogenetic processes [911]. However, this equilibrium is easily disturbed by anthropogenic activities (e.g., agricultural practices, deforestation, and overgrazing), and such effects are mainly noticed in the developing countries with poor technical and financial resources to manage natural resources [10, 11]. In order to make sound decisions regarding sustainable land use systems, knowledge of SQ related to different land use scenarios is essential [12]. It is therefore most important to assess SQ degradation of different land use and soil management systems using soil quality index (SQI) since many of the factors that influence sustainable productivity are related to SQ. Information on SQI can support to further prioritization and then device management strategies that improve soil resources sustainably [11]. To do so, applying the concept of SQI is desirable as individual soil properties in isolation may not be sufficient to quantify changes in SQ related to land use and soil management systems [11, 13]. In line to this, many studies reported that indexing SQ indicators based on a combination of soil properties could better reflect the status of SQ degradation as compared to individual parameters [6, 1315].

Despite the importance of SQI in describing SQ degradation or aggradations, there is no universally accepted dataset selection, scoring, and SQ indexing method for field conditions. Previous studies reported that different methods of minimum dataset selection (MDS), scoring, and SQ indexing have been applied but SQI results varied even for the same conditions [11, 13, 16]. The most widely reported MDS methods of SQ indicators are expert opinion and statistical tools (e.g., regression, principal component analysis (PCA)) [11, 13]. An expert can generate a list of appropriate SQ indicators on the basis of ecosystem processes and functions and other decision rules such as management goals for a site associated with soil functions as well as other site-specific factors, like region or crop sensitivity as selection criteria [13, 16].

The transformation of the datasets into scores (scoring function) can be done using linear and nonlinear scoring techniques [11, 13]. Studies elsewhere compared the two scoring methods to represent soil system function but the value of nonlinear scoring method was reported higher than the linear method [11, 16, 17]. There are different types of linear and nonlinear scoring functions, even though none of the previous studies have evaluated them all simultaneously [11, 16]. Different SQ indexing methods have been also used by different researchers [13, 1618]. The same authors have reported that there are differences in SQI values among the various SQ indexing methods (e.g., additive, weighted, and max-min objective functions). Despite the fact that there is diversity in data selection, scoring, and SQ indexing methods, previous studies have limitation in evaluating the methods using the same data simultaneously in a similar field conditions.

Regardless of the above limitation, having SQI of long-term land use and soil management systems is necessary in order to locate areas to be carefully managed for sustainable development. The use of site-specific SQI can help planners and decision makers to evaluate which land use and management system is most sustainable and vice-versa in a given situation [18, 19]. These authors also noted that SQI can reflect the extent of SQ degradation and thereby give support to suggest appropriate remedial measures such as optimum fertilizer rates and planning of other suitable land management practices considering potentials and constraints of different fields at large scale such as a catchment.

In general, SQI is a useful assessment tool that may help move soil conservation and resource management beyond assessments of soil erosion and changes in productivity [13]. SQI can thus provide the necessary information for planners and decision makers to make informed decisions against SQ degradation using the introduction of appropriate interventions. Despite such importance of SQI in combating SQ degradation, only few studies have been reported in relation to various land use and soil management systems. This indicated that research on SQI has been mostly neglected for unknown reasons, with the most probable reason which could be technical and financial limitations.

Many approaches assessing SQ degradation using the concept of SQI have been already developed and applied elsewhere [6, 11, 13, 1519]. In this study, such concepts are adopted and evaluated to narrow the knowledge/information gap of SQI across different land use and soil management systems in the northern Ethiopia. The objective of this study was to assess and identify an effective SQ indicator dataset among 25 soil measurements, appropriate scoring functions for each indicator, and an efficient SQ indexing method to evaluate soil degradation across the LUSMS in the Mai-Negus catchment of northern Ethiopia.

2. Materials and Methods

2.1. Study Area

This evaluation was conducted in the Mai-Negus catchment in Tigray regional state, northern Ethiopia (Figure 1). The catchment covers an area of 1240 ha, with a landscape consisting of generally rugged terrain at altitudes ranging from 2060 to 2650 m above sea level. Land use is dominantly arable with teff (Eragrostis tef) being the primary crop on >80% of the land area. The remainder of the catchment is pasture with scattered patches of mixed tree, bush, and shrub cover. The major rock types are lava pyroclastic and metavolcanic. According to FAO-UNESCO Soil Classification System, soils are dominantly Leptosols at very steep positions, Cambisols on middle to steep slopes, and Vertisols on flat areas [20]. Annual rainfall averages 700 mm but is very erratic in amount and distribution throughout the catchment. Mean annual temperature was 22°C.

2.2. Selection of Land Use and Soil Management Systems

Eight land use and soil management systems (LUSMS) were selected on the basis of three steps. First, information on historical and existing LUSMS in the catchment was collected and described. Soil sampling units were then identified across each LUSMS. Finally, composite soil samples were collected, processed, and analyzed for several SQ indicators using standard laboratory procedures. The first step is described below, with details for the second and third steps given in Sections 2.3 and 2.4.

Field reconnaissance surveys and informal group discussions were conducted in June 2009 by the author, two development agents, and six farmers who are knowledgeable about the catchment and local farming systems. The six farmers were selected purposively because the large group size made it impractical for all the household heads to participate and doing so would have been problematic for discussion and consensus building. The dominant cropping history and soil management practices for each LUSMS were identified and described by the team. In addition, terrain characteristic and soil factors were documented for each rain-fed agricultural LUSMS. All eight LUSMS were selected as much as possible to be from similar soil type (Cambisols) and a range in slope gradient in the catchment. Topographical characteristics of each sampling units are presented in Table 1.

Based on land use information acquired in the study, eight LUSMS that represent the best and worst management practices being used throughout the study catchment were identified and are described (Table 2). The LUSMS selected for SQ evaluation were (i) natural forest (LS1), (ii) plantation on protected areas (LS2), (iii) grazed land (LS3), (iv) teff (Eragrostis tef)-faba bean (Vicia faba) rotations (LS4), (v) teff-wheat (Triticum vulgare)/barley (Hordeum vulgare) rotations (LS5), (vi) teff monocropping (LS6), (vii) maize (Zea mays) monocropping (LS7), and (viii) uncultivated marginal land (LS8). The various LUSMS were in place for various amounts of time ranging from 5 to 6 years for teff monocropping and 20 to 30 years for maize monocropping. Average age for the other systems was about 10 years except for the plantation, grazed land, and uncultivated marginal land with which was in place for more than 15 years. To assess the impact of LUSM on SQ indicators, it is either necessary to have a baseline against which human induced differences can be measured [21] or to measure the same systems repeatedly in time [13]. For this study, the researcher chose to use the natural forest (LS1) as a reference, assuming the soil in those areas is less disturbed than in cultivated or grazed areas.

2.3. Soil Sampling, Processing, and Analysis

After identifying the eight LUSMS locations, three soil sampling units and their corresponding areas were selected by the researcher considering representativeness and uniformity of the fields. Soil samples were collected from 24 areas () in June 2009 and analyzed for 24 potential SQ indicators. Each soil sampling unit ranged from 50 to 80 m2. Five to eight soil samples from 0 to 20 cm depth (plow layer) were collected randomly from each unit and mixed to form a composite sample. The number of composite samples was determined by the size and homogeneity (hydrologic conditions) of each sampling unit. Fewer samples were collected from homogenous, small fields compared to large, heterogeneous fields. The sampling focused on the plow layer because this is where most SQ changes are expected to occur due to long-term land use and soil management practices. Each composite soil sample was mixed thoroughly in a bucket before taking a 500 g subsample that was air dried and sieved to pass a 2 mm mesh before analysis.

The soil samples were analyzed for selected physical, chemical, and biological SQ indicators. Soil texture was determined using the Bouyoucos hydrometer method [22] and soil bulk density (BD) by the core method [23]. Percent pore space (total porosity) was computed from BD and average particle density (PD) of 2.65 g cm−3 as [24]. Soil aggregate stability (SAS) was measured using the wet sieve method [25] and maximum water holding capacity (MWHC) determined by equilibrating the soil with water through capillary action in a KR box [26]. A-horizon depth was directly measured in pits opened to a 60 cm depth within each LUSMS.

Soil pH was determined using a 1 : 2.5 soil to water ratio with a combined glass electrode [27]. Soil organic (OC) was determined by the Walkley-Black method [28], available phosphorus (Pav) by the Olsen method [29], total nitrogen (TN), and total phosphorus (TP) by the Kjeldahl digestion method [30]. Cation exchange capacity (CEC) was determined by ammonium acetate buffered at pH 7 [31]. Exchangeable bases (calcium, Ca; magnesium, Mg; potassium, K) were analyzed after extraction using 1 M ammonium acetate at pH 7.0. Iron and zinc were determined using 0.005 M diethylene triamine pentaacetic acid (DTPA) extraction as described in Baruah and Barthakur [26].

Exchangeable sodium percentage (ESP) was calculated by dividing exchangeable Na+ by CEC. Base saturation percentage (BSP) was calculated by dividing the sum of base forming cations by CEC, multiplied by 100% [32]. A 25th SQ indicator, earthworm population, was monitored monthly as a biological indicator throughout the wet season (mid-June–mid-September 2009). Three randomly collected soil samples ( cm) from a 1 m2 area were passed through a 10 mm sieve to separate and then count the average number of earthworms in each LUSMS sampling unit.

2.4. Soil Quality Index Computations

After measuring the 25 potential SQ indicators using field and laboratory analysis techniques, different methods for computing SQI values were evaluated. Although the type of data used for each SQI may differ, the process of SQ indexing follows the same three basic steps regardless of the method used (Figure 2). These steps are indicator selection, interpretation/scoring, and integration into index value [11, 33]. The detailed descriptions of each step are given below.

2.4.1. Indicator Selection (Step  1)

Potential SQ indicators were selected based on their sensitivity to management practices, ability to describe major soil processes, ease and cost of sampling and laboratory analysis, and significance of increasing productivity (agronomic) and protecting environmental soil functions. Two methods for selecting a minimum dataset (MDS), expert opinion (EO), and principal component analysis (PCA) were compared. For the EO approach MDS, variables were chosen from the 25 potential SQ indicators based on researcher knowledge and literature recommendations [3436].

2.4.2. Indicator Transformation/Scoring Function (Step  2)

After selecting the MDS using EO or PCA (Step  1), each value was transformed using two different scoring techniques [11, 33]. Both linear and nonlinear scoring functions were compared as described below.

(I) Linear Scoring Functions. Three linear scoring functions (LSF) were identified and evaluated. For the first LSF approach, indicators were ranked in ascending for “more is better” or descending order for “less is better” in terms of soil functions [37]. Each “more is better” indicator was divided by the highest value in the group such that it received a score of 1. For “less is better” indicator with the lowest observed value, it was divided by itself so that it received a score of 1. Threshold values were identified and used as outlined by Liebig et al. [37]. This approach was therefore described as the “Liebig Linear Scoring Function.”

For the second LSF approach, SQ indicator values were transformed to a common range between 0.1 and 1.0 using homothetic transformation equations (1) and (2) [38]. In the context of this study, this is termed as “homothetic transformation method of LSF” : where and are values of the variables after transformation. is value of the variable to transform and and are the maximum and minimum threshold values of the variable (Table 3). Equation (1) is used for “more is better” scoring function, (2) for “less is better,” and a combination of both equations for “optimum is better” scoring function. This LSF approach is therefore described as the “homothetic transformation method” .

The third LSF approach was adopted from Masto et al. [11], Glover et al. [39], and Masto et al. [40]. It is described as where is the linear score, is the soil property value, and and are lower and upper threshold values (Table 3). For values below and above the threshold, the score is zero. Equation (3) was used for “more is better,” whereas (4) for “less is better” and a combination of both for “optimum is better.” This approach is therefore described as “Glover LSF method” .

(II) Nonlinear Scoring Functions. Two nonlinear scoring function (NLSF) approaches were also evaluated. The first NLSF approach transformed the indicators using curves constructed with CurveExpert version 1.3 shareware (http://www.flu.org.cn/en/download-79.html) as described by Andrews et al. [13]. The shape of each curve, that is, bell-shaped (midpoint is optimum), sigmoid with an upper asymptote (more is better) or sigmoid with a lower asymptote (less is better), was determined according to agronomic and environmental soil functions using data from undisturbed fields (natural ecosystem), literature values, and knowledge of experts. To develop indicator curve using CurveExpert, it was assumed that levels of activity found in undisturbed soil (natural forest) would have a score at or near 1.0. Polynomial fit model was applied after examining the data by curve finder and then drawing the curve based on the observed values of SQ indicators from undisturbed ecosystem. After interpolated by polynomial Lagrangian interpolation method in the CurveExpert version 1.3, the corresponding transformed value was analyzed for each untransformed indicator value in each LUSMS. The -axis for the functions represented a site-specific expected range of values of the soil properties. The -axis, ranging from 0 to 1, was the transformed score. This approach is termed as the “CurveExpert method of NLSF” .

The second NLSF used (5) to normalize SQ indicators suggested by Masto et al. [11]; Glover et al. [39]; Masto et al. [40] and is therefore termed as the “GloverNLSF method” : where is the indicator value, is the baseline or value of the soil property where the score equals 0.5 or about the midpoint between the upper and lower threshold value, and is the slope. Baseline values are generally regarded as the minimum target value. If soil indicator values are located within the control limits, the system is considered to be in an acceptable state. Conversely, if the values lie outside the threshold limits, the system is considered to be in a state of degradation [11, 40, 44].

Critical values or thresholds were established based on a range of values measured in natural ecosystems, best managed systems, and values adopted from literatures and personal experiences of the researcher to better fit to the local conditions within the catchment (Table 3). This table also presents the indicator scoring curve, baseline, and threshold values used to transform selected SQ indicators. For detailed descriptions of these standard scoring functions the reader is referred to Masto et al. [11]; Glover et al. [39]; Masto et al. [40]; Karlen et al. [45].

Measured values are transformed into unitless scores ranging from 0 to 1 so that scores can be combined and averaged into a single value such that a score of 1 represents the highest potential function for that system; that is, the indicator is nonlimiting with regard to the pertinent soil functions and processes [46]. An advantage of indexing is that important information can often be captured by scoring that might otherwise go undetected when examining only the observed values [13].

2.4.3. Soil Quality Indexing (Step  3)

Two SQ indexing methods are commonly found in literature, for example, Masto et al. [11]; Glover et al. [39]; Masto et al. [40]. The first is unscreened transformation (the additive index) and the second uses principal component analysis (PCA). Both methods were applied to data from each LUSMS. The SQI values were then compared with those from natural forest land systems to assess the degree of soil degradation or improvement. Details of the two SQ indexing methods are summarized below.

(I) Unscreened Transformation Based Soil Quality Indexing (Unscreened-SQI). SQ indicators are integrated into an index (SQI) by summing the scores from individual indicators and dividing by the total number of indicators (i.e., an additive model) as described in Masto et al. [40]: where SQI is the soil quality index, is the linear or nonlinear scored value of individual indicators, and is the number of indicators included in the dataset.

(II) Principal Component Analysis Based Soil Quality Indexing (PCA-SQI). A standard PCA was conducted using all SQ indicators that showed significantly differences among the LUSMS. Under each principal component (PC), only the variables having high factor loadings and eigenvalues >1 that explained at least 5% of the data variations were retained for indexing. Among well-correlated variables within each PC, the variable with the highest correlation coefficient (absolute value) and loading factor was chosen. If the highly weighted variables (≥±0.7 eigenvector) were not well correlated (), each was considered important and retained in the PC for SQ indexing. As each PC explains a certain amount of variation within the total dataset, this provides a “weight” for the variables chosen under a given PC. The final PCA based SQI (7) described by Masto et al. [40] was used for this study: where PCA-SQI is principal component analysis (PCA) based soil quality index, is the PCA weighing factor equal to the ratio of variance of each factor to total cumulative variance coefficients in the equation, and is scored value of each SQ indicator.

2.5. Evaluation of Soil Quality Indexing Methods

The SQ indexing methods were evaluated using sensitivity analysis described by Masto et al. [11] as where and are the maximum and minimum SQI observed under each scoring procedure using each dataset selection methods. The SQ indexing method with higher value of sensitivity is more preferable as this is sensitive to perturbations and management practices [21].

2.6. Data Analysis

Data were analyzed using statistical software package of SPSS 18.0 [47]. One-way analysis of variance (ANOVA) was performed to determine the effects of LUSMS on SQI. Mean SQI was tested for its level of significance at probability level . Data were also analyzed using correlation and factor analysis. A PCA was used to examine the relationship among the 24 SQ indicators by statistically grouping them into four PC factors through the varimax rotation procedure. Varimax rotation with Kaiser Normalization was used because this results in a factor pattern that highly loads into one factor [48]. Communalities estimate the portion of variance of each soil attribute that explains for the factors. A high communality for a soil attribute indicates a high proportion of the variance explained in the factor. Less importance should be ascribed to soil attributes with low communalities when interpreting the factors [48].

3. Results and Discussion

3.1. Grouping Soil Quality Indicators

A moderate to strong correlation () among many SQ indicators within the different LUSMS was observed, indicating a multicollinearity effect (data not shown). Factor analysis can help reduce the number of indicators analyzed needed for indexing by identifying components that best account for the variability and thus minimizing data redundancy (multicollinearity effect). The 25 SQ indicators analyzed to evaluate the eight LUSMS were grouped using a PCA. This resulted in four principal component (PC) groups that best explained variability in the data (Table 4). Communalities of the SQ indicators (Table 4) show that individual indicators accounted for 64 to 97% of the variance. Indicators with high communality get preference over those with low communality [49]. When combined, the first four PCs factors with eigenvalues >1 explained about 88% of the soil variability among the eight LUSMS. The first two PCs accounted for about 52% of the variance, indicating that these would be potential components to assess SQ effects within the LUSMS.

Eigenvectors for the first four PCs (Table 4) show that OC, TN, EW, porosity, and Zn were retained for PC1. The correlation coefficient between TN and the other variables was less than 0.6, so TN was retained in PC1. However, the high correlation coefficient () of OC with the other high loading variables suggested that OC was most important so it too was retained in PC1 for SQ indexing. Several literature references show that EW populations are highly influenced by organic matter availability [50, 51]. Furthermore, according to Jenkinson [52], soil microbial biomass comprises 1 to 4% of the total OC and 2 to 6% of the total organic nitrogen. Based on this information, we chose to exclude EW data from PC1 since its contribution is explained by soil organic matter. For such reasons, PC1 is referred to as the “soil organic matter factor.”

Soil CEC, TP, exchangeable K and Ca were the highly loaded factors attributed to PC2. The CEC values were strongly correlated () with K, and Ca so CEC was selected as the preferred variable for retention in PC2. The correlation coefficient of TP with the other variables was <0.6 and since the cutoff is , TP was also retained in PC2. Phosphorus is frequently a limiting factor for crop production in northern Ethiopia soils, so its inclusion in a SQ index is logical for assessing SQ degradation within the various LUSMS. Based on these two factors, PC2 is referred to as the “soil macro-nutrient factor.”

The only highly loaded variables in PC3 were silt, DBD, and pH. The correlation of silt with the highly loaded variables was less than 0.6, so it was retained in PC3 for index development. A similar analysis showed a higher partial correlation for DBD (), indicating that it was also a preferred indicator for retention in PC3. Both silt and DBD also showed higher communality effect when compared to pH. These indicators also influenced the SQ in an opposite direction (Table 4); their inclusion in any SQ indexing is crucial to assess variability associated with the various LUSMS. Again, based on the critical indicators, PC3 is referred to as the “soil physical property factor.” The fourth PC is referred to as the “soil micro-nutrient factor,” because Fe is the highly loaded variable (Table 4).

Based on the four PCs, a composite PCA-SQI consisting of soil OC, TN, CEC, TP, silt, DBD, and Fe was chosen to assess SQ variability among the LUSMS. Weighting factors were developed based on the percent variation explained by the first four PCs (Table 4), resulting in a final normalized PCA based SQI equation: where PCA-SQI is a PCA based soil quality index and is the score (linear or nonlinear) for each variable (Table 4), with coefficients based on the variance accounted for by each PC.

3.2. Integration of Soil Quality Indicators into a Soil Quality Index
3.2.1. Unscreened Transformation Based Soil Quality Index (Unscreened-SQI)

An unscreened-SQI using the PCA and EO selected minimum datasets showed the highest SQI values for the (homothetic transformation method) and (CurveExpert method of NLSF) models used to compare the eight LUSMS (Tables 5 and 6). For all comparisons, the nonlinear unscreened-SQI showed higher values than the linear scored functions. In addition, unscreened-SQI values for both LSF and NLSF methods were greater for EO minimum datasets than for PCA selected datasets.

Among LUSMS, LS1 had a significantly higher () unscreened-SQI value whereas a lower value in LS8 using either the PCA or EO selected datasets. The general SQI pattern among the LUSMS was LS1 > LS2 > LS7 > LS4 > LS3 > LS5 > LS6 > LS8, (Tables 5 and 6). Mean unscreened-SQI differences between the LS1 reference soils and average unscreened-SQI values for PCA and EO datasets showed the following soil degradation levels: LS2 (−8%), LS3 (−20%), LS7 (−24%), LS4 (−26%), LS6 (−38%), LS5 (−45%), and LS8 (−58%). Generally, this study demonstrated that LS2 followed by LS3 and LS7 is more advantageous in maintaining SQ than the other LUSMS in the catchment.

3.2.2. PCA Based Soil Quality Indices (PCA-SQI)

For both PCA and EO databases, NLSF methods resulted in higher SQI values than LSF (Tables 5 and 6). In contrast to unscreened-SQI values, PCA-SQI values derived from PCA selected datasets were higher than those derived using EO datasets. Once again, the and indexing methods resulted in the highest PCA-SQI values for both data selection methods (Tables 5 and 6). Among the LUSMS, the PCA-SQI value for LS1 was significantly higher () than the other land uses although it has been just slightly better than LS2 and LS3. LS8 had the lowest PCA-SQI value followed closely by LS5 and LS6 (Tables 5 and 6). Soil degradation based on PCA-SQI differences between each LUSMS and the LS1 reference averaged −6, −21, −23, −24, −34, −36, and −59% for LS2, LS3, LS7, LS4, LS5, LS6, and LS8, respectively, when computed using average values for PCA and EO derived datasets.

3.3. Dataset Selection, Scoring, and Indexing Methods Comparison

Using NLSF, PCA-SQI values were greater than unscreened-SQI values but with LSF methods unscreened-SQI values were higher. Comparing MDS selection methods, unscreened-SQI values derived using an EO dataset were higher for both LSF and NLSF methods (Tables 5 and 6). The and scoring methods resulted in higher SQI values than the other methods. Overall, for almost all LUSMS, NLSF values were generally greater than LSF values. Mean SQI values using PCA selected datasets ranged from 0.458 (LS8) to 0.932 (LS1) (Table 5). The overall average using EO selected datasets ranged from 0.406 (LS8) to 0.874 (LS1) (Table 6). This implies that mean SQI values calculated from both indexing methods using EO datasets resulted in lower index values than with PCA selected datasets. Therefore, this study concludes that adopting the PCA dataset selection method is the best for evaluating the soil degradation status of LUSMS in northern Ethiopia and under similar environmental conditions. The use of a PCA dataset selection method also minimizes any disciplinary bias that may be associated with the EO selection method.

3.4. Evaluation of Indexing Methods

The unscreened-SQI values varied from 0.410 to 0.933 and 0.397 to 0.954 for the PCA and EO dataset selection methods, respectively. Likewise, PCA-SQI values ranged from 0.347 to 0.934 and 0.335 to 0.801 for the two dataset selection methods, respectively (Table 7). A sensitivity analysis showed that both indexing and dataset selection methods influenced the SQI values for the various LUSMS throughout the catchment (Table 7). It also confirmed that indexing, scoring, and MDS selection can be used to evaluate soil degradation among LUSMS in the Mai-Negus catchment of northern Ethiopia. However, the sensitivity analyses do not agree with findings reported by Masto et al. [11] and Andrews et al. [53] who favored NLSF over LSF. One reason for this difference is that the previous reports were based on single scoring functions from LSF and NLSF methods and their sensitivity was not evaluated. Also, many of the LSF and NLSF approaches evaluated in this project were not addressed at the same time and place in the previous studies.

3.5. Synthesis of Dataset Selection, Scoring, and SQ Indexing Methods

Dataset selection methods (PCA versus EO), scoring functions, and indexing approaches all influenced the SQI values used to compare LUSMS in northern Ethiopia (Tables 5, 6 and 7). This study results were consistent with previous studies that also showed these factors contributed to variability in SQI values [11, 16, 34, 35]. Arguments favoring EO method for selection of a MDS stress its focus on sustainable management goals [11, 36, 45], but many researchers have also relied on statistical techniques such as PCA (e.g., [11, 16, 53]). This study confirmed PCA selection could reduce expert biases as compared to an EO method. Andrews et al. [16] reported that EO and PCA dataset selection methods produced almost similar SQI values, but PCA-SQI values using PCA and NLSF methods in this study were higher. This was not true for unscreened-SQI which was obtained by summing scores for individual indicators.

Many studies have shown that NLSF resulted in higher SQI values than LSF [11, 16], but a sensitivity analysis showed that the method explained variation better than the NLSF method for both indexing methods. This was followed by another LSF method denoted as . Several strategies for integrating scored indicators values into an overall SQI have been proposed, but no approach has received universally acceptance [11].

SQI values were segregated into physical, chemical, and biological components (Figure 3) to show the component effects on SQ and thus soil degradation within each LUSMS in the catchment. The individual factors also permit comparisons among the scoring methods that help illustrate differences among them. Finally, although such comparisons are interesting, they are not adequate for assessing SQ effects of LUSMS because numerous interactions also contribute to the overall SQI values [13, 15]. This is also a major reason that use of a single SQ indicator regardless of how its scored or interpreted is not suitable for comparing overall soil degradation dynamics among LUSMS.

Several studies have shown that agroforestry can be a viable option for restoration of degraded land as it improves soil indicators such as organic matter, nitrogen cycling, soil structure, and biological activity [5456]. Others have added that extensive rooting systems and protective canopies of the Leucaena leeucocephala and Sesbania sesban trees can protect the soil from erosion and create favorable conditions for plant and microbial growth [56, 57]. Litter-fall and fine roots decomposition both contribute to these processes. This study shows higher SQI value for LS1 and LS2 because of better vegetation cover (natural and plantation), which is consistent with previous reports [5457]. Overall, this study confirmed that SQI values can be developed and used as tools for early identification of land degradation particularly the soil component. They can also be used to support decision making processes with regard to sustainable soil management planning in the context of different LUSMS at catchment scale.

4. Conclusion

This study demonstrated that natural forest land systems (LS1) have relatively good soil quality (SQ), whereas uncultivated marginal land systems (LS8) have a seriously degraded soil within the Mai-Negus catchment of northern Ethiopia. Areas being managed in a teff-barley/wheat rotation land system (LS5) and teff monocropping (LS6) also had not reasonably good SQ. SQI values for LS8 and LS6 provided an early warning regarding the severity of soil degradation since more than 50% of the original soil is degraded when compared to LS1 (reference soil). The PCA-SQI method resulted in a higher SQI value than the unscreened-SQI for NLSF when datasets were selected using a PCA method. However, the reverse was true for the datasets selected using an EO method. The NLSF methods resulted in higher SQI values than LSF methods for all LUSMS regardless of the SQ indexing method. Among NLSF, the CurveExpert method showed the highest SQI for the dataset selected by the PCA method. However, a sensitivity analysis indicated that PCA-SQI with the LSF (denoted as or the “Liebig method”) and a PCA derived dataset had the highest values. The lowest sensitivity values were associated with an unscreened-SQI using NLSF method denoted as (Glover method NLSF) using PCA derived datasets. This study also indicated that the PCA-SQI method with a PCA selected dataset and a LSF designated as was the most sensitive for assessing differences among the eight LUSMS in northern Ethiopia. Based on this study result, the researcher concludes that use of relevant datasets, scoring functions, and indexing methods are all important factors influencing SQI values within the eight LUSMS. Also, for improving degraded soils and maintaining SQ, appropriate interventions within the various LUSMS should be identified, prioritized, and monitored using SQI values as decision aids.

Conflict of Interests

The author does not have a direct financial relation with the commercial identity mentioned in this paper that might lead to a conflict of interests.

Acknowledgments

The author gratefully acknowledges the financial support by DAAD/GIZ (Germany) through the Center for Development Research (ZEF), University of Bonn (Germany), and field work supported by Aksum University (Ethiopia). The author also highly appreciates the cooperation of the participant farmers and assistance offered by the local administration and extension agents during the field work. The author is also grateful to the anonymous reviewers for their comments which helped in improving this paper.