# Soil Quality Assessment Strategies for Evaluating Soil Degradation in Northern Ethiopia

Authors: Gebreyesus Brhane Tesfahunegn
Journal: Applied and Environmental Soil Science
 (2014)

---

## 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.


---

## Body

## 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 [2–4]. 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 [9–11]. 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, 13–15].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, 16–18]. 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, 15–19]. 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 (Figure1). 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.Figure 1
Locations of the study area: Ethiopia (a), Tigray Region (b), and Mai-Negus catchment (study site) (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 Sections2.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 Table1.Table 1
Topographic characteristics of each sampling unit in the eight LUSMS selected for soil quality assessment within the Mai-Negus catchment in northern Ethiopia.
LUSMS
Sampling Unit 1
Sampling Unit 2
Sampling Unit 3
Elevation(m)
Slope(%)
UTMa
Elevation(m)
Slope(%)
UTM
Elevation(m)
Slope(%)
UTM
Latitude
Longitude
Latitude
Longitude
Latitude
Longitude
LS1
2175
7
463657
1561482
2170
5.5
463856
1561592
2173
6.2
464039
1561523
LS2
2124
4.5
463383
1560411
2152
5
463615
1560503
2138
5.5
463248
1560503
LS3
2136
5
463425
1561193
2165
5.5
463489
1561466
2151
5.0
463535
1561614
LS4
2168
6
463592
1561250
2145
6.5
462538
1561454
2157
6.0
463168
1561832
LS5
2141
6.5
463606
1561338
2163
7
463317
1561763
2152
6.8
463707
1561912
LS6
2165
6
463734
1561459
2136
5.5
463363
1561500
2151
5.8
464124
1561328
LS7
2162
6.5
462810
1561314
2155
7
462641
1561614
2159
6.8
464027
1561202
LS8
2149
7
462795
1561152
2138
5
462745
1561718
2144
6.0
463913
1562153
aUniversal Transverse Mercator 37 North (UTM-37N) in meters is the projection system.
LUSMS: land use and soil management systems; LS1: natural forest (reference); LS2: plantation of protected area; LS3: grazed land; LS4: teff (Eragrostis tef)-faba bean (Vicia faba) rotation; LS5: teff-wheat (Triticum vulgare)/Barley (Hordeum vulgare) rotation; LS6: teff monocropping; LS7: Maize (Zea mays) mono-cropping; LS8: uncultivated-marginal land system.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 (Table2). 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.Table 2
Description of the eight land use and soil management systems (LUSMS) evaluated within the Mai-Negus catchment in northern Ethiopia.
Serialnumber
Land use and soil management systems
Description
1
Natural forest (LS1)
Less disturbed land, used as a reference in the system, which has native trees, vegetation, and grass cover.
2
Plantation of protected area (LS2)
Sesbania (Sesbania sesban) and Leucaena (Leucaena leucocephala) trees plantation was established 16 years ago and grass was used by cuting and carrying during the dry season for livestock, protected throughout the year from livestock interferences; no fertilizer application and less intensive soil conservation measures exist.
3
Grazed  land (LS3)
Open grazed land appeared during the dry season 16 years ago with no inclusion of any improved management practices, for example, soil and water conservation and enrichment of plant/grass species. It overstocked in the dry months (November–June) but it is a swampy area for the rest of months.
4
Teff (Eragrostis tef)-faba bean (Vicia faba) rotation (LS4)
Fields were harvested of teff (Eragrostic tef (Zucc) Trot) crop before soil samples were collected and rotated with faba bean (Vicia faba  L.) for more than 5-6 years. Urea and diammonium phosphate (DAP) fertilizers were applied each for teff at 50 kg ha−1 y−1 but sometimes reduced by half depending on resource availability and quality of the soil. Teff needed 4–6 times tillage and at least one time weeding. For faba bean, 2-3 times tillage was sufficient with no or one time weeding and addition of manure around homestead is common practice but is urea or DAP rarely used. Soil and water conservation was used at field borders.
5
Teff-wheat (Triticum vulgare)/Barley (Hordeum vulgare) rotation (LS5)
The fields were planted wheat (Tritium vulgare L.) before teff (Eragrostic  tef (Zucc) Trot) and soil samples were collected after teff was harvested. Wheat (Triticum vulgare L.)/barley (Hordeum vulgare  L.) with teff (Eragrostic tef) was rotated for more than 6 years. 50 kg ha−1 y−1 of each urea and DAP fertilizers was used for teff field but the amount varies with crop color at vegetative stage and soil quality condition for fertilizer rate application on wheat/barley crop fields. The fertilizer rate used for wheat/barley is lower than teff.
6
aTeff mono-cropping land system (LS6)
For more than 5-6 years, teff  (Eragrostic tef (Zucc) Trot) was sown continuously with 50 kg ha−1 y−1 urea and 100 kg ha−1 y−1 DAP fertilizers; 5–7 tillage frequency; at least one time hand weeding. Manure and intensive SWC for fields around homestead was applied.
7
Maize (Zea mays L.)-mono cropping land system (LS7)
Maize (Zea mays L.) was planted for more than 20–30 years continuously. The addition of manure was estimated at 6–12 t ha−1 y−1, which varies with manure and labor availability for transportation; 2-3 tillage frequency with at least one times hand hoeing and weeding. This is always practiced around homestead fields. Conservation measures are also well executed.
8
Uncultivated marginal land system (LS8)
It was terraced with wide spacing since the last 20 years but most of it is broken and used as an open grazing land throughout the year; on some spot areas it had very few to few naturally growing but over grazed grass species like Bermuda grass (Cynodon dactylon L.) management practices such as fertilizer and species enrichment were not introduced. Farmers considered it as the most degraded soil, for example, abandoned land.
aTeff is the dominant crop in the study catchment and other parts of northern Ethiopia. It is an annual cereal crop (belonging to the grass family) which has sparse crop canopies and provides little cover to the soil against erosion. It has very fine seeds that require repeated plowing of fields and preparation of fine seedbeds, which increases the vulnerability of the soil to erosion.


### 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 (8LUSMS×3samplingunits) 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 {1-BD/PD}×100 [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 (25×25×20 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 (Figure2). These steps are indicator selection, interpretation/scoring, and integration into index value [11, 33]. The detailed descriptions of each step are given below.Figure 2
A conceptual model for computing soil quality indices. (From Andrews [33].)


#### 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 [34–36].


#### 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 (IL) was therefore described as the “Liebig Linear Scoring Function.”For thesecond 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”(IIL):
(1)Y=0.1+((X-b)(a-b))*0.9,(2)Z=1-((X-b)(a-b))*0.9,
where Y and Z are values of the variables after transformation. X is value of the variable to transform and a and b 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”(IIL).Table 3
Soil quality indicators, scoring curve, and threshold and baseline limits used for evaluating eight LUSMS in northern Ethiopia.
Indicator
Scoring curve
Threshold
Baselinec
Optimumd
Slope at baseline
Source of limits
Lowera
Upperb
Sand (%)
Optimum
0
60
L, 30
36
0.440
Natural ecosystem, best managed soils
U, 50
Silt (%)
More is better
0
38
19
—
0.249
Best managed soil, natural ecosystem
Clay (%)
More is better
0
31
15
—
0.268
Natural ecosystem
SAS (%)
More is better
10
60
30
—
0.296
Harris et al. [17]; natural system
BD (g cm−3)
Less is better
1
2.2
1.5
1.2
−0.323
Harris et al. [17];  natural ecosystem
A-horizon depth (cm)
More is better
0
20
10
—
0.318
Natural ecosystem
MWHC (%)
More is better
20
58
30
—
0.198
Gregory et al. [41]; natural system
OCe (%)
More is better
1
6.5
3.5
—
1.046
Kay and Angers [42]; natural ecosystem
CEC (cmolc kg−1)
More is better
6
46
20
—
0.245
Natural ecosystem, best managed soil
TN (%)
More is better
0.05
0.54
0.30
—
25.408
Natural ecosystem,
TP (mg kg−1)
More is better
200
1316
650
—
0.005
Natural ecosystem
Pav (mg kg−1)
More is better
5
29
15
—
0.433
Mausbach and Seybold [43]; natural ecosystem
Zn (mg kg−1)
Optimum
2
20
L, 10
14
0.855
Mausbach and Seybold [43]; natural system
U, 18
Fe (mg kg−1)
Optimum
10
50
L, 20
26
0.577
Harris et al. [17]; natural   system
U, 40
Earth worm countper m2
More is better
0
11
5
—
0.470
Natural ecosystem
aSoils at or below the threshold values are prone to structural destabilization, erosion, and low productivity; so the scoring value is 0.
bSoils at or beyond this values  no further increase in productivity or decrease in erosion rate are achieved the upper threshold; values at and above this level thus receive a score of 1.0.
cValues receive a score of 0.5 and are generally regarded as the minimum target values.
dThe value is given a score of 1.0 if the desired relationship is bell-shaped.
eAccording to Kay and Angers [42], irrespective of soil type if SOC contents are below 1%, it may not be possible to obtain potential yields.
L: lower; U: upper; SAS: soil aggregate stability; BD: bulk density; MWC: maximum water holding capacity; OC: organic carbon; CEC: cation exchangeable capacity; TN: total nitrogen; TP: total phosphorous; Pav: available phosphorous; Zn: available Zinc; Fe: available iron; —: implies not applicable.Thethird LSF approach was adopted from Masto et al. [11], Glover et al. [39], and Masto et al. [40]. It is described as
(3)(Y)=(x-s)(t-s),(4)(Y)=1-(x-s)(t-s),
where Y is the linear score, x is the soil property value, and s and t 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”(IIIL).(II) Nonlinear Scoring Functions. Two nonlinear scoring function (NLSF) approaches were also evaluated. The firstNLSF 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 Lagrangianinterpolation method in the CurveExpert version 1.3, the corresponding transformed value was analyzed for each untransformed indicator value in each LUSMS. The X-axis for the functions represented a site-specific expected range of values of the soil properties. The Y-axis, ranging from 0 to 1, was the transformed score. This approach is termed as the “CurveExpert method of NLSF”(IN).Thesecond 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”(IIN):
(5)NLSF(Y)=1[1+e-b(x-A)],
where x is the indicator value, A 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 b 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 (Table3). 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]:
(6)unscreenedSQI=(∑i=1nSin),
where SQI is the soil quality index, S is the linear or nonlinear scored value of individual indicators, and n 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 (r<0.60), 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:
(7)PCA-SQI=∑i=1nWiSi,
where PCA-SQI is principal component analysis (PCA) based soil quality index, Wi is the PCA weighing factor equal to the ratio of variance of each factor to total cumulative variance coefficients in the equation, and Si 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
(8)Sensitivity(S)=SQI(max)SQI(min),
where SQI(max) and SQI(min) 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 (P)≤0.05. 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].


## 2.1. Study Area

This evaluation was conducted in the Mai-Negus catchment in Tigray regional state, northern Ethiopia (Figure1). 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.Figure 1
Locations of the study area: Ethiopia (a), Tigray Region (b), and Mai-Negus catchment (study site) (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 Sections2.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 Table1.Table 1
Topographic characteristics of each sampling unit in the eight LUSMS selected for soil quality assessment within the Mai-Negus catchment in northern Ethiopia.
LUSMS
Sampling Unit 1
Sampling Unit 2
Sampling Unit 3
Elevation(m)
Slope(%)
UTMa
Elevation(m)
Slope(%)
UTM
Elevation(m)
Slope(%)
UTM
Latitude
Longitude
Latitude
Longitude
Latitude
Longitude
LS1
2175
7
463657
1561482
2170
5.5
463856
1561592
2173
6.2
464039
1561523
LS2
2124
4.5
463383
1560411
2152
5
463615
1560503
2138
5.5
463248
1560503
LS3
2136
5
463425
1561193
2165
5.5
463489
1561466
2151
5.0
463535
1561614
LS4
2168
6
463592
1561250
2145
6.5
462538
1561454
2157
6.0
463168
1561832
LS5
2141
6.5
463606
1561338
2163
7
463317
1561763
2152
6.8
463707
1561912
LS6
2165
6
463734
1561459
2136
5.5
463363
1561500
2151
5.8
464124
1561328
LS7
2162
6.5
462810
1561314
2155
7
462641
1561614
2159
6.8
464027
1561202
LS8
2149
7
462795
1561152
2138
5
462745
1561718
2144
6.0
463913
1562153
aUniversal Transverse Mercator 37 North (UTM-37N) in meters is the projection system.
LUSMS: land use and soil management systems; LS1: natural forest (reference); LS2: plantation of protected area; LS3: grazed land; LS4: teff (Eragrostis tef)-faba bean (Vicia faba) rotation; LS5: teff-wheat (Triticum vulgare)/Barley (Hordeum vulgare) rotation; LS6: teff monocropping; LS7: Maize (Zea mays) mono-cropping; LS8: uncultivated-marginal land system.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 (Table2). 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.Table 2
Description of the eight land use and soil management systems (LUSMS) evaluated within the Mai-Negus catchment in northern Ethiopia.
Serialnumber
Land use and soil management systems
Description
1
Natural forest (LS1)
Less disturbed land, used as a reference in the system, which has native trees, vegetation, and grass cover.
2
Plantation of protected area (LS2)
Sesbania (Sesbania sesban) and Leucaena (Leucaena leucocephala) trees plantation was established 16 years ago and grass was used by cuting and carrying during the dry season for livestock, protected throughout the year from livestock interferences; no fertilizer application and less intensive soil conservation measures exist.
3
Grazed  land (LS3)
Open grazed land appeared during the dry season 16 years ago with no inclusion of any improved management practices, for example, soil and water conservation and enrichment of plant/grass species. It overstocked in the dry months (November–June) but it is a swampy area for the rest of months.
4
Teff (Eragrostis tef)-faba bean (Vicia faba) rotation (LS4)
Fields were harvested of teff (Eragrostic tef (Zucc) Trot) crop before soil samples were collected and rotated with faba bean (Vicia faba  L.) for more than 5-6 years. Urea and diammonium phosphate (DAP) fertilizers were applied each for teff at 50 kg ha−1 y−1 but sometimes reduced by half depending on resource availability and quality of the soil. Teff needed 4–6 times tillage and at least one time weeding. For faba bean, 2-3 times tillage was sufficient with no or one time weeding and addition of manure around homestead is common practice but is urea or DAP rarely used. Soil and water conservation was used at field borders.
5
Teff-wheat (Triticum vulgare)/Barley (Hordeum vulgare) rotation (LS5)
The fields were planted wheat (Tritium vulgare L.) before teff (Eragrostic  tef (Zucc) Trot) and soil samples were collected after teff was harvested. Wheat (Triticum vulgare L.)/barley (Hordeum vulgare  L.) with teff (Eragrostic tef) was rotated for more than 6 years. 50 kg ha−1 y−1 of each urea and DAP fertilizers was used for teff field but the amount varies with crop color at vegetative stage and soil quality condition for fertilizer rate application on wheat/barley crop fields. The fertilizer rate used for wheat/barley is lower than teff.
6
aTeff mono-cropping land system (LS6)
For more than 5-6 years, teff  (Eragrostic tef (Zucc) Trot) was sown continuously with 50 kg ha−1 y−1 urea and 100 kg ha−1 y−1 DAP fertilizers; 5–7 tillage frequency; at least one time hand weeding. Manure and intensive SWC for fields around homestead was applied.
7
Maize (Zea mays L.)-mono cropping land system (LS7)
Maize (Zea mays L.) was planted for more than 20–30 years continuously. The addition of manure was estimated at 6–12 t ha−1 y−1, which varies with manure and labor availability for transportation; 2-3 tillage frequency with at least one times hand hoeing and weeding. This is always practiced around homestead fields. Conservation measures are also well executed.
8
Uncultivated marginal land system (LS8)
It was terraced with wide spacing since the last 20 years but most of it is broken and used as an open grazing land throughout the year; on some spot areas it had very few to few naturally growing but over grazed grass species like Bermuda grass (Cynodon dactylon L.) management practices such as fertilizer and species enrichment were not introduced. Farmers considered it as the most degraded soil, for example, abandoned land.
aTeff is the dominant crop in the study catchment and other parts of northern Ethiopia. It is an annual cereal crop (belonging to the grass family) which has sparse crop canopies and provides little cover to the soil against erosion. It has very fine seeds that require repeated plowing of fields and preparation of fine seedbeds, which increases the vulnerability of the soil to erosion.


## 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 (8LUSMS×3samplingunits) 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 {1-BD/PD}×100 [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 (25×25×20 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 (Figure2). These steps are indicator selection, interpretation/scoring, and integration into index value [11, 33]. The detailed descriptions of each step are given below.Figure 2
A conceptual model for computing soil quality indices. (From Andrews [33].)


### 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 [34–36].


### 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 (IL) was therefore described as the “Liebig Linear Scoring Function.”For thesecond 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”(IIL):
(1)Y=0.1+((X-b)(a-b))*0.9,(2)Z=1-((X-b)(a-b))*0.9,
where Y and Z are values of the variables after transformation. X is value of the variable to transform and a and b 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”(IIL).Table 3
Soil quality indicators, scoring curve, and threshold and baseline limits used for evaluating eight LUSMS in northern Ethiopia.
Indicator
Scoring curve
Threshold
Baselinec
Optimumd
Slope at baseline
Source of limits
Lowera
Upperb
Sand (%)
Optimum
0
60
L, 30
36
0.440
Natural ecosystem, best managed soils
U, 50
Silt (%)
More is better
0
38
19
—
0.249
Best managed soil, natural ecosystem
Clay (%)
More is better
0
31
15
—
0.268
Natural ecosystem
SAS (%)
More is better
10
60
30
—
0.296
Harris et al. [17]; natural system
BD (g cm−3)
Less is better
1
2.2
1.5
1.2
−0.323
Harris et al. [17];  natural ecosystem
A-horizon depth (cm)
More is better
0
20
10
—
0.318
Natural ecosystem
MWHC (%)
More is better
20
58
30
—
0.198
Gregory et al. [41]; natural system
OCe (%)
More is better
1
6.5
3.5
—
1.046
Kay and Angers [42]; natural ecosystem
CEC (cmolc kg−1)
More is better
6
46
20
—
0.245
Natural ecosystem, best managed soil
TN (%)
More is better
0.05
0.54
0.30
—
25.408
Natural ecosystem,
TP (mg kg−1)
More is better
200
1316
650
—
0.005
Natural ecosystem
Pav (mg kg−1)
More is better
5
29
15
—
0.433
Mausbach and Seybold [43]; natural ecosystem
Zn (mg kg−1)
Optimum
2
20
L, 10
14
0.855
Mausbach and Seybold [43]; natural system
U, 18
Fe (mg kg−1)
Optimum
10
50
L, 20
26
0.577
Harris et al. [17]; natural   system
U, 40
Earth worm countper m2
More is better
0
11
5
—
0.470
Natural ecosystem
aSoils at or below the threshold values are prone to structural destabilization, erosion, and low productivity; so the scoring value is 0.
bSoils at or beyond this values  no further increase in productivity or decrease in erosion rate are achieved the upper threshold; values at and above this level thus receive a score of 1.0.
cValues receive a score of 0.5 and are generally regarded as the minimum target values.
dThe value is given a score of 1.0 if the desired relationship is bell-shaped.
eAccording to Kay and Angers [42], irrespective of soil type if SOC contents are below 1%, it may not be possible to obtain potential yields.
L: lower; U: upper; SAS: soil aggregate stability; BD: bulk density; MWC: maximum water holding capacity; OC: organic carbon; CEC: cation exchangeable capacity; TN: total nitrogen; TP: total phosphorous; Pav: available phosphorous; Zn: available Zinc; Fe: available iron; —: implies not applicable.Thethird LSF approach was adopted from Masto et al. [11], Glover et al. [39], and Masto et al. [40]. It is described as
(3)(Y)=(x-s)(t-s),(4)(Y)=1-(x-s)(t-s),
where Y is the linear score, x is the soil property value, and s and t 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”(IIIL).(II) Nonlinear Scoring Functions. Two nonlinear scoring function (NLSF) approaches were also evaluated. The firstNLSF 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 Lagrangianinterpolation method in the CurveExpert version 1.3, the corresponding transformed value was analyzed for each untransformed indicator value in each LUSMS. The X-axis for the functions represented a site-specific expected range of values of the soil properties. The Y-axis, ranging from 0 to 1, was the transformed score. This approach is termed as the “CurveExpert method of NLSF”(IN).Thesecond 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”(IIN):
(5)NLSF(Y)=1[1+e-b(x-A)],
where x is the indicator value, A 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 b 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 (Table3). 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]:
(6)unscreenedSQI=(∑i=1nSin),
where SQI is the soil quality index, S is the linear or nonlinear scored value of individual indicators, and n 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 (r<0.60), 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:
(7)PCA-SQI=∑i=1nWiSi,
where PCA-SQI is principal component analysis (PCA) based soil quality index, Wi is the PCA weighing factor equal to the ratio of variance of each factor to total cumulative variance coefficients in the equation, and Si is scored value of each SQ indicator.


## 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 [34–36].


## 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 (IL) was therefore described as the “Liebig Linear Scoring Function.”For thesecond 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”(IIL):
(1)Y=0.1+((X-b)(a-b))*0.9,(2)Z=1-((X-b)(a-b))*0.9,
where Y and Z are values of the variables after transformation. X is value of the variable to transform and a and b 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”(IIL).Table 3
Soil quality indicators, scoring curve, and threshold and baseline limits used for evaluating eight LUSMS in northern Ethiopia.
Indicator
Scoring curve
Threshold
Baselinec
Optimumd
Slope at baseline
Source of limits
Lowera
Upperb
Sand (%)
Optimum
0
60
L, 30
36
0.440
Natural ecosystem, best managed soils
U, 50
Silt (%)
More is better
0
38
19
—
0.249
Best managed soil, natural ecosystem
Clay (%)
More is better
0
31
15
—
0.268
Natural ecosystem
SAS (%)
More is better
10
60
30
—
0.296
Harris et al. [17]; natural system
BD (g cm−3)
Less is better
1
2.2
1.5
1.2
−0.323
Harris et al. [17];  natural ecosystem
A-horizon depth (cm)
More is better
0
20
10
—
0.318
Natural ecosystem
MWHC (%)
More is better
20
58
30
—
0.198
Gregory et al. [41]; natural system
OCe (%)
More is better
1
6.5
3.5
—
1.046
Kay and Angers [42]; natural ecosystem
CEC (cmolc kg−1)
More is better
6
46
20
—
0.245
Natural ecosystem, best managed soil
TN (%)
More is better
0.05
0.54
0.30
—
25.408
Natural ecosystem,
TP (mg kg−1)
More is better
200
1316
650
—
0.005
Natural ecosystem
Pav (mg kg−1)
More is better
5
29
15
—
0.433
Mausbach and Seybold [43]; natural ecosystem
Zn (mg kg−1)
Optimum
2
20
L, 10
14
0.855
Mausbach and Seybold [43]; natural system
U, 18
Fe (mg kg−1)
Optimum
10
50
L, 20
26
0.577
Harris et al. [17]; natural   system
U, 40
Earth worm countper m2
More is better
0
11
5
—
0.470
Natural ecosystem
aSoils at or below the threshold values are prone to structural destabilization, erosion, and low productivity; so the scoring value is 0.
bSoils at or beyond this values  no further increase in productivity or decrease in erosion rate are achieved the upper threshold; values at and above this level thus receive a score of 1.0.
cValues receive a score of 0.5 and are generally regarded as the minimum target values.
dThe value is given a score of 1.0 if the desired relationship is bell-shaped.
eAccording to Kay and Angers [42], irrespective of soil type if SOC contents are below 1%, it may not be possible to obtain potential yields.
L: lower; U: upper; SAS: soil aggregate stability; BD: bulk density; MWC: maximum water holding capacity; OC: organic carbon; CEC: cation exchangeable capacity; TN: total nitrogen; TP: total phosphorous; Pav: available phosphorous; Zn: available Zinc; Fe: available iron; —: implies not applicable.Thethird LSF approach was adopted from Masto et al. [11], Glover et al. [39], and Masto et al. [40]. It is described as
(3)(Y)=(x-s)(t-s),(4)(Y)=1-(x-s)(t-s),
where Y is the linear score, x is the soil property value, and s and t 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”(IIIL).(II) Nonlinear Scoring Functions. Two nonlinear scoring function (NLSF) approaches were also evaluated. The firstNLSF 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 Lagrangianinterpolation method in the CurveExpert version 1.3, the corresponding transformed value was analyzed for each untransformed indicator value in each LUSMS. The X-axis for the functions represented a site-specific expected range of values of the soil properties. The Y-axis, ranging from 0 to 1, was the transformed score. This approach is termed as the “CurveExpert method of NLSF”(IN).Thesecond 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”(IIN):
(5)NLSF(Y)=1[1+e-b(x-A)],
where x is the indicator value, A 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 b 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 (Table3). 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]:
(6)unscreenedSQI=(∑i=1nSin),
where SQI is the soil quality index, S is the linear or nonlinear scored value of individual indicators, and n 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 (r<0.60), 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:
(7)PCA-SQI=∑i=1nWiSi,
where PCA-SQI is principal component analysis (PCA) based soil quality index, Wi is the PCA weighing factor equal to the ratio of variance of each factor to total cumulative variance coefficients in the equation, and Si 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
(8)Sensitivity(S)=SQI(max)SQI(min),
where SQI(max) and SQI(min) 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 (P)≤0.05. 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 (r>0.7) 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.Table 4
Principal component analysis results using 25 potential soil quality indicators to evaluate eight LUSMS in the Mai-Negus catchment, northern Ethiopia.
Eigenvector
Principal component, PC
Communalities
1
2
3
4
Organic carbon, OC
0.86
0.29
0.12
0.28
0.92
Total nitrogen, TN
0.86
0.21
0.04
0.36
0.92
Earthworm per m2, EW
0.79
0.27
0.28
0.23
0.82
Porosity
0.74
0.38
0.40
0.30
0.94
Zinc, Zn
0.70
0.46
0.31
0.16
0.82
Water-holding capacity, WHC
0.69
0.47
0.47
0.20
0.96
Available phosphorus, Pav
0.69
0.34
0.31
0.11
0.71
A-horizon depth, AHD
0.64
0.37
0.48
0.41
0.94
Soil aggregate stability, SAS
0.60
0.45
0.52
0.33
0.94
Exchangeable potassium, K
0.27
0.83
0.28
0.17
0.86
Cation exchangeable capacity, CEC
0.42
0.76
0.51
0.29
0.95
Exchangeable calcium, Ca
0.41
0.72
0.40
0.28
0.93
Total phosphorous, TP
0.45
0.70
0.38
0.11
0.85
Exchangeable magnesium, Mg
0.53
0.68
0.38
0.24
0.94
Sum base forming cations, SBF
0.45
0.63
0.40
0.26
0.97
Silt
0.19
0.39
0.84
0.04
0.90
Dry bulk density, DBD
−0.40
−0.38
−0.78
−0.30
0.94
pH
0.22
0.26
0.74
0.26
0.82
Sand
−0.34
−0.36
−0.67
−0.30
0.91
EC
0.38
0.44
0.55
−0.03
0.64
Exchangeable sodium, Na
−0.18
−0.18
0.01
−0.57
0.82
Iron, Fe
−0.30
−0.08
−0.13
0.86
0.85
Exchangeable sodium percentage, ESP
−0.17
−0.25
−0.49
−0.59
0.97
Clay
0.51
0.13
0.21
0.63
0.72
Eigenvalue
7.26
5.16
4.79
3.82
—
Variance (%)
30.26
21.51
19.96
15.90
—
Cumulative variance (%)
30.26
51.77
71.73
87.63
—
Boldface eigenvector values correspond to the PCs highly weighted variables examined for the index.
Bold-italic factors correspond to the indicators retained in the SQ index. The weight of the variables included in the index was decided using the variance of each factor.
PCA-SQI = 0.303OC + 0.0303TN + 0.215CEC + 0.215TP + 0.200silt + 0.200DBD + 0.159Fe.
Normalized PCA-SQI =(0.303OC + 0.0303TN + 0.215CEC + 0.215TP + 0.200silt + 0.200DBD + 0.159Fe)/(1.595) = 0.190OC + 0.190TN + 0.135CEC + 0.135TP + 0.125silt + 0.125DBD + 0.100Fe.Eigenvectors for the first four PCs (Table4) 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 (r>0.87) 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 (r=0.85) 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 r=0.7, 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 (r>0.80), 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 (Table4), resulting in a final normalized PCA based SQI equation:
(9)PCA-SQI=0.190OC+0.190TN+0.135CEC+0.135TP+0.125silt+0.125DBD+0.100Fe,
where PCA-SQI is a PCA based soil quality index and S 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 theIIL(homothetic transformation method) and IN (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.Table 5
Soil quality indices developed using a PCA selected dataset to compare eight LUSMS in the Mai-Negus catchment of northern Ethiopia.
LUSMS
Unscreened SQI
PCA-SQI
LSFa
NLSFa
LSF
NLSF
Meanb
IL
IIL
IIIL
IN
IIN
IL
IIL
IIIL
IN
IIN
LS1
0.913a
0.941a
0.899a
0.956a
0.953a
0.913a
0.925a
0.881a
0.973a
0.964a
0.932a
LS2
0.823a
0.843a
0.793b
0.920a
0.916a
0.826b
0.845a
0.795b
0.928a
0.921a
0.861a
LS3
0.621bc
0.659b
0.626c
0.747bcd
0.741bcd
0.619c
0.658b
0.593c
0.773bc
0.754bc
0.679b
LS4
0.648b
0.673b
0.627c
0.780bc
0.777bc
0.618c
0.608bc
0.587c
0.774bc
0.772bc
0.686b
LS5
0.539cd
0.582bc
0.522d
0.689cd
0.686cd
0.484d
0.535cd
0.470d
0.696cd
0.692cd
0.589c
LS6
0.506d
0.532cd
0.510d
0.655de
0.652de
0.463d
0.514d
0.453d
0.656d
0.653d
0.559c
LS7
0.639b
0.676b
0.660c
0.787b
0.784b
0.597c
0.638b
0.591c
0.791b
0.781b
0.694b
LS8
0.410e
0.475d
0.405e
0.562e
0.559e
0.347e
0.418e
0.341e
0.513e
0.546e
0.458d
mean
0.637
0.672
0.630
0.762
0.759
0.609
0.642
0.589
0.763
0.760
0.682
LSD (P=0.05)
0.095
0.099
0.091
0.097
0.096
0.086
0.082
0.076
0.089
0.089
0.083
Means followed by different letters in the same column are significantly different atP=0.05; LSD: least significance difference.
aDetails of LSF and NLSF descriptions can be found in Section 2.4 (Step 2).
bThis shows overall mean of the different scoring and indexing methods across the columns for the same LUSMS as tested statistically.
LSF: linear scoring function; NLSF: nonlinear scoring function; IL: Liebig method LSF; IIL: Homothetic transformation method of LSF; IIIL: Glover method LSF; IN: CurveExpert method NLSF; IIN: Glover method NLSF; LUSMS: land use and soil management systems; SQI: soil quality index; PCA: principal component analysis; LS1: natural forest (reference); LS2: plantation of protected area; LS3: pasture land system; LS4: teff-faba bean rotation; LS5: teff-barley/wheat rotation; LS6: teff-monocropping; LS7: Maize monocropping; LS8: uncultivated-marginal land soil system.Table 6
Soil quality indices developed using an EO selected dataset to compare eight LUSMS in the Mai-Negus catchment of northern Ethiopia.
LUSMS
Unscreened SQI
PCA-SQI
LSFa
NLSFa
LSF
NLSF
Meanb
IL
IIL
IIIL
IN
IIN
IL
IIL
IIIL
IN
IIN
LS1
0.937a
0.943a
0.925a
0.966a
0.964a
0.796a
0.799a
0.779a
0.818a
0.816a
0.874a
LS2
0.828b
0.831b
0.786b
0.920a
0.918a
0.724b
0.723b
0.690b
0.783a
0.781a
0.799b
LS3
0.707c
0.720c
0.673c
0.830b
0.828b
0.612c
0.622c
0.584c
0.708b
0.666b
0.695c
LS4
0.652c
0.662c
0.610c
0.775b
0.773b
0.556d
0.565d
0.522d
0.651c
0.650b
0.642d
LS5
0.529d
0.552d
0.491d
0.658c
0.656c
0.448e
0.468e
0.417e
0.548d
0.546c
0.531e
LS6
0.501d
0.526d
0.461d
0.621c
0.619c
0.428e
0.450e
0.397e
0.497e
0.517c
0.502e
LS7
0.657c
0.670c
0.622c
0.773b
0.771c
0.545d
0.558d
0.514d
0.653c
0.651b
0.641d
LS8
0.397e
0.434e
0.361e
0.511d
0.509d
0.335f
0.369f
0.307f
0.420f
0.418d
0.406f
mean
0.651
0.667
0.616
0.756
0.754
0.556
0.569
0.526
0.635
0.631
0.636
LSD (P=0.05)
0.066
0.059
0.067
0.060
0.060
0.044
0.038
0.042
0.042
0.051
0.048
Means followed by different letters in the same rows are significantly different atP=0.05; LSD: least significance difference.
aDetails of LSF and NLSF description can be found in Section 2.4 (Step 2).
bThis shows the overall mean of the different scoring and indexing methods across the columns for the same LUSMS tested statistically.
Explanations of abbreviations are similar to footnotes under Table5.Among LUSMS, LS1 had a significantly higher (P≤0.05) 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 (Tables5 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 IN and IIL 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 (P=0.05) 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 (Tables5 and 6). The IIL and IN 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 (Table7). 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.Table 7
Dataset selection, scoring function, and indexing method comparisons for eight LUSMS in the Mai-Negus catchment of northern Ethiopia.
SQ indexing method
Linear scoring function
Nonlinear scoring function
IL
IIL
IIIL
IN
IIN
Range
S
Range
S
Range
S
Range
S
Range
S
Unscreened-SQIa
0.410–0.913
2.227
0.475–0.941
1.981
0.405–0.899
2.220
0.562–0.956
1.701
0.559–0.953
1.705
PCA-SQIa
0.347–0.913
2.631
0.418–0.923
2.208
0.341–0.881
2.584
0.513–0.973
1.898
0.546–0.964
1.766
Unscreened-SQIb
0.397–0.937
2.360
0.434–0.943
2.173
0.361–0.925
2.562
0.511–0.966
1.890
0.508–0.964
1.897
PCA-SQIb
0.335–0.796
2.376
0.369–0.799
2.165
0.307–0.779
2.537
0.420–0.818
1.952
0.418–0.816
1.948
Mean
—
2.399
—
2.132
—
2.476
—
1.859
—
1.830
adenotes minimum dataset chosen using principal component analyses (PCA).
bdenotes minimum dataset chosen by expert opinion (EO).
IL: Liebig method LSF; IIL: Homothetic transformation method of LSF; IIIL: Glover method LSF; IN: CurveExpert method NLSF; IIN: Glover method NLSF; S: Sensitivity analysis; LSF: linear scoring function; NLSF: nonlinear scoring function.


### 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 (Tables5, 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 IL method explained variation better than the NLSF method for both indexing methods. This was followed by another LSF method denoted as IIIL. 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 (Figure3) 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.LUSMS effects on mean physical, chemical, and biological soil quality indicator scores for different indexing methods within the Mai-Negus catchment of northern Ethiopia.
(a)
(b)
(c)
(d)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 [54–56]. 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 [54–57]. 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.


## 3.1. Grouping Soil Quality Indicators

A moderate to strong correlation (r>0.7) 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.Table 4
Principal component analysis results using 25 potential soil quality indicators to evaluate eight LUSMS in the Mai-Negus catchment, northern Ethiopia.
Eigenvector
Principal component, PC
Communalities
1
2
3
4
Organic carbon, OC
0.86
0.29
0.12
0.28
0.92
Total nitrogen, TN
0.86
0.21
0.04
0.36
0.92
Earthworm per m2, EW
0.79
0.27
0.28
0.23
0.82
Porosity
0.74
0.38
0.40
0.30
0.94
Zinc, Zn
0.70
0.46
0.31
0.16
0.82
Water-holding capacity, WHC
0.69
0.47
0.47
0.20
0.96
Available phosphorus, Pav
0.69
0.34
0.31
0.11
0.71
A-horizon depth, AHD
0.64
0.37
0.48
0.41
0.94
Soil aggregate stability, SAS
0.60
0.45
0.52
0.33
0.94
Exchangeable potassium, K
0.27
0.83
0.28
0.17
0.86
Cation exchangeable capacity, CEC
0.42
0.76
0.51
0.29
0.95
Exchangeable calcium, Ca
0.41
0.72
0.40
0.28
0.93
Total phosphorous, TP
0.45
0.70
0.38
0.11
0.85
Exchangeable magnesium, Mg
0.53
0.68
0.38
0.24
0.94
Sum base forming cations, SBF
0.45
0.63
0.40
0.26
0.97
Silt
0.19
0.39
0.84
0.04
0.90
Dry bulk density, DBD
−0.40
−0.38
−0.78
−0.30
0.94
pH
0.22
0.26
0.74
0.26
0.82
Sand
−0.34
−0.36
−0.67
−0.30
0.91
EC
0.38
0.44
0.55
−0.03
0.64
Exchangeable sodium, Na
−0.18
−0.18
0.01
−0.57
0.82
Iron, Fe
−0.30
−0.08
−0.13
0.86
0.85
Exchangeable sodium percentage, ESP
−0.17
−0.25
−0.49
−0.59
0.97
Clay
0.51
0.13
0.21
0.63
0.72
Eigenvalue
7.26
5.16
4.79
3.82
—
Variance (%)
30.26
21.51
19.96
15.90
—
Cumulative variance (%)
30.26
51.77
71.73
87.63
—
Boldface eigenvector values correspond to the PCs highly weighted variables examined for the index.
Bold-italic factors correspond to the indicators retained in the SQ index. The weight of the variables included in the index was decided using the variance of each factor.
PCA-SQI = 0.303OC + 0.0303TN + 0.215CEC + 0.215TP + 0.200silt + 0.200DBD + 0.159Fe.
Normalized PCA-SQI =(0.303OC + 0.0303TN + 0.215CEC + 0.215TP + 0.200silt + 0.200DBD + 0.159Fe)/(1.595) = 0.190OC + 0.190TN + 0.135CEC + 0.135TP + 0.125silt + 0.125DBD + 0.100Fe.Eigenvectors for the first four PCs (Table4) 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 (r>0.87) 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 (r=0.85) 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 r=0.7, 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 (r>0.80), 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 (Table4), resulting in a final normalized PCA based SQI equation:
(9)PCA-SQI=0.190OC+0.190TN+0.135CEC+0.135TP+0.125silt+0.125DBD+0.100Fe,
where PCA-SQI is a PCA based soil quality index and S 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 theIIL(homothetic transformation method) and IN (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.Table 5
Soil quality indices developed using a PCA selected dataset to compare eight LUSMS in the Mai-Negus catchment of northern Ethiopia.
LUSMS
Unscreened SQI
PCA-SQI
LSFa
NLSFa
LSF
NLSF
Meanb
IL
IIL
IIIL
IN
IIN
IL
IIL
IIIL
IN
IIN
LS1
0.913a
0.941a
0.899a
0.956a
0.953a
0.913a
0.925a
0.881a
0.973a
0.964a
0.932a
LS2
0.823a
0.843a
0.793b
0.920a
0.916a
0.826b
0.845a
0.795b
0.928a
0.921a
0.861a
LS3
0.621bc
0.659b
0.626c
0.747bcd
0.741bcd
0.619c
0.658b
0.593c
0.773bc
0.754bc
0.679b
LS4
0.648b
0.673b
0.627c
0.780bc
0.777bc
0.618c
0.608bc
0.587c
0.774bc
0.772bc
0.686b
LS5
0.539cd
0.582bc
0.522d
0.689cd
0.686cd
0.484d
0.535cd
0.470d
0.696cd
0.692cd
0.589c
LS6
0.506d
0.532cd
0.510d
0.655de
0.652de
0.463d
0.514d
0.453d
0.656d
0.653d
0.559c
LS7
0.639b
0.676b
0.660c
0.787b
0.784b
0.597c
0.638b
0.591c
0.791b
0.781b
0.694b
LS8
0.410e
0.475d
0.405e
0.562e
0.559e
0.347e
0.418e
0.341e
0.513e
0.546e
0.458d
mean
0.637
0.672
0.630
0.762
0.759
0.609
0.642
0.589
0.763
0.760
0.682
LSD (P=0.05)
0.095
0.099
0.091
0.097
0.096
0.086
0.082
0.076
0.089
0.089
0.083
Means followed by different letters in the same column are significantly different atP=0.05; LSD: least significance difference.
aDetails of LSF and NLSF descriptions can be found in Section 2.4 (Step 2).
bThis shows overall mean of the different scoring and indexing methods across the columns for the same LUSMS as tested statistically.
LSF: linear scoring function; NLSF: nonlinear scoring function; IL: Liebig method LSF; IIL: Homothetic transformation method of LSF; IIIL: Glover method LSF; IN: CurveExpert method NLSF; IIN: Glover method NLSF; LUSMS: land use and soil management systems; SQI: soil quality index; PCA: principal component analysis; LS1: natural forest (reference); LS2: plantation of protected area; LS3: pasture land system; LS4: teff-faba bean rotation; LS5: teff-barley/wheat rotation; LS6: teff-monocropping; LS7: Maize monocropping; LS8: uncultivated-marginal land soil system.Table 6
Soil quality indices developed using an EO selected dataset to compare eight LUSMS in the Mai-Negus catchment of northern Ethiopia.
LUSMS
Unscreened SQI
PCA-SQI
LSFa
NLSFa
LSF
NLSF
Meanb
IL
IIL
IIIL
IN
IIN
IL
IIL
IIIL
IN
IIN
LS1
0.937a
0.943a
0.925a
0.966a
0.964a
0.796a
0.799a
0.779a
0.818a
0.816a
0.874a
LS2
0.828b
0.831b
0.786b
0.920a
0.918a
0.724b
0.723b
0.690b
0.783a
0.781a
0.799b
LS3
0.707c
0.720c
0.673c
0.830b
0.828b
0.612c
0.622c
0.584c
0.708b
0.666b
0.695c
LS4
0.652c
0.662c
0.610c
0.775b
0.773b
0.556d
0.565d
0.522d
0.651c
0.650b
0.642d
LS5
0.529d
0.552d
0.491d
0.658c
0.656c
0.448e
0.468e
0.417e
0.548d
0.546c
0.531e
LS6
0.501d
0.526d
0.461d
0.621c
0.619c
0.428e
0.450e
0.397e
0.497e
0.517c
0.502e
LS7
0.657c
0.670c
0.622c
0.773b
0.771c
0.545d
0.558d
0.514d
0.653c
0.651b
0.641d
LS8
0.397e
0.434e
0.361e
0.511d
0.509d
0.335f
0.369f
0.307f
0.420f
0.418d
0.406f
mean
0.651
0.667
0.616
0.756
0.754
0.556
0.569
0.526
0.635
0.631
0.636
LSD (P=0.05)
0.066
0.059
0.067
0.060
0.060
0.044
0.038
0.042
0.042
0.051
0.048
Means followed by different letters in the same rows are significantly different atP=0.05; LSD: least significance difference.
aDetails of LSF and NLSF description can be found in Section 2.4 (Step 2).
bThis shows the overall mean of the different scoring and indexing methods across the columns for the same LUSMS tested statistically.
Explanations of abbreviations are similar to footnotes under Table5.Among LUSMS, LS1 had a significantly higher (P≤0.05) 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 (Tables5 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 IN and IIL 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 (P=0.05) 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.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 theIIL(homothetic transformation method) and IN (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.Table 5
Soil quality indices developed using a PCA selected dataset to compare eight LUSMS in the Mai-Negus catchment of northern Ethiopia.
LUSMS
Unscreened SQI
PCA-SQI
LSFa
NLSFa
LSF
NLSF
Meanb
IL
IIL
IIIL
IN
IIN
IL
IIL
IIIL
IN
IIN
LS1
0.913a
0.941a
0.899a
0.956a
0.953a
0.913a
0.925a
0.881a
0.973a
0.964a
0.932a
LS2
0.823a
0.843a
0.793b
0.920a
0.916a
0.826b
0.845a
0.795b
0.928a
0.921a
0.861a
LS3
0.621bc
0.659b
0.626c
0.747bcd
0.741bcd
0.619c
0.658b
0.593c
0.773bc
0.754bc
0.679b
LS4
0.648b
0.673b
0.627c
0.780bc
0.777bc
0.618c
0.608bc
0.587c
0.774bc
0.772bc
0.686b
LS5
0.539cd
0.582bc
0.522d
0.689cd
0.686cd
0.484d
0.535cd
0.470d
0.696cd
0.692cd
0.589c
LS6
0.506d
0.532cd
0.510d
0.655de
0.652de
0.463d
0.514d
0.453d
0.656d
0.653d
0.559c
LS7
0.639b
0.676b
0.660c
0.787b
0.784b
0.597c
0.638b
0.591c
0.791b
0.781b
0.694b
LS8
0.410e
0.475d
0.405e
0.562e
0.559e
0.347e
0.418e
0.341e
0.513e
0.546e
0.458d
mean
0.637
0.672
0.630
0.762
0.759
0.609
0.642
0.589
0.763
0.760
0.682
LSD (P=0.05)
0.095
0.099
0.091
0.097
0.096
0.086
0.082
0.076
0.089
0.089
0.083
Means followed by different letters in the same column are significantly different atP=0.05; LSD: least significance difference.
aDetails of LSF and NLSF descriptions can be found in Section 2.4 (Step 2).
bThis shows overall mean of the different scoring and indexing methods across the columns for the same LUSMS as tested statistically.
LSF: linear scoring function; NLSF: nonlinear scoring function; IL: Liebig method LSF; IIL: Homothetic transformation method of LSF; IIIL: Glover method LSF; IN: CurveExpert method NLSF; IIN: Glover method NLSF; LUSMS: land use and soil management systems; SQI: soil quality index; PCA: principal component analysis; LS1: natural forest (reference); LS2: plantation of protected area; LS3: pasture land system; LS4: teff-faba bean rotation; LS5: teff-barley/wheat rotation; LS6: teff-monocropping; LS7: Maize monocropping; LS8: uncultivated-marginal land soil system.Table 6
Soil quality indices developed using an EO selected dataset to compare eight LUSMS in the Mai-Negus catchment of northern Ethiopia.
LUSMS
Unscreened SQI
PCA-SQI
LSFa
NLSFa
LSF
NLSF
Meanb
IL
IIL
IIIL
IN
IIN
IL
IIL
IIIL
IN
IIN
LS1
0.937a
0.943a
0.925a
0.966a
0.964a
0.796a
0.799a
0.779a
0.818a
0.816a
0.874a
LS2
0.828b
0.831b
0.786b
0.920a
0.918a
0.724b
0.723b
0.690b
0.783a
0.781a
0.799b
LS3
0.707c
0.720c
0.673c
0.830b
0.828b
0.612c
0.622c
0.584c
0.708b
0.666b
0.695c
LS4
0.652c
0.662c
0.610c
0.775b
0.773b
0.556d
0.565d
0.522d
0.651c
0.650b
0.642d
LS5
0.529d
0.552d
0.491d
0.658c
0.656c
0.448e
0.468e
0.417e
0.548d
0.546c
0.531e
LS6
0.501d
0.526d
0.461d
0.621c
0.619c
0.428e
0.450e
0.397e
0.497e
0.517c
0.502e
LS7
0.657c
0.670c
0.622c
0.773b
0.771c
0.545d
0.558d
0.514d
0.653c
0.651b
0.641d
LS8
0.397e
0.434e
0.361e
0.511d
0.509d
0.335f
0.369f
0.307f
0.420f
0.418d
0.406f
mean
0.651
0.667
0.616
0.756
0.754
0.556
0.569
0.526
0.635
0.631
0.636
LSD (P=0.05)
0.066
0.059
0.067
0.060
0.060
0.044
0.038
0.042
0.042
0.051
0.048
Means followed by different letters in the same rows are significantly different atP=0.05; LSD: least significance difference.
aDetails of LSF and NLSF description can be found in Section 2.4 (Step 2).
bThis shows the overall mean of the different scoring and indexing methods across the columns for the same LUSMS tested statistically.
Explanations of abbreviations are similar to footnotes under Table5.Among LUSMS, LS1 had a significantly higher (P≤0.05) 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 (Tables5 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 IN and IIL 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 (P=0.05) 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 (Tables5 and 6). The IIL and IN 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 (Table7). 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.Table 7
Dataset selection, scoring function, and indexing method comparisons for eight LUSMS in the Mai-Negus catchment of northern Ethiopia.
SQ indexing method
Linear scoring function
Nonlinear scoring function
IL
IIL
IIIL
IN
IIN
Range
S
Range
S
Range
S
Range
S
Range
S
Unscreened-SQIa
0.410–0.913
2.227
0.475–0.941
1.981
0.405–0.899
2.220
0.562–0.956
1.701
0.559–0.953
1.705
PCA-SQIa
0.347–0.913
2.631
0.418–0.923
2.208
0.341–0.881
2.584
0.513–0.973
1.898
0.546–0.964
1.766
Unscreened-SQIb
0.397–0.937
2.360
0.434–0.943
2.173
0.361–0.925
2.562
0.511–0.966
1.890
0.508–0.964
1.897
PCA-SQIb
0.335–0.796
2.376
0.369–0.799
2.165
0.307–0.779
2.537
0.420–0.818
1.952
0.418–0.816
1.948
Mean
—
2.399
—
2.132
—
2.476
—
1.859
—
1.830
adenotes minimum dataset chosen using principal component analyses (PCA).
bdenotes minimum dataset chosen by expert opinion (EO).
IL: Liebig method LSF; IIL: Homothetic transformation method of LSF; IIIL: Glover method LSF; IN: CurveExpert method NLSF; IIN: Glover method NLSF; S: Sensitivity analysis; LSF: linear scoring function; NLSF: nonlinear scoring function.


## 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 (Tables5, 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 IL method explained variation better than the NLSF method for both indexing methods. This was followed by another LSF method denoted as IIIL. 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 (Figure3) 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.LUSMS effects on mean physical, chemical, and biological soil quality indicator scores for different indexing methods within the Mai-Negus catchment of northern Ethiopia.
(a)
(b)
(c)
(d)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 [54–56]. 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 [54–57]. 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, theCurveExpert method(IN) 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 IL 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 IIN (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 IL 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.



---

*Source: 646502-2014-02-04.xml*