Next Article in Journal
Validation of Instruments for the Improvement of Interprofessional Education through Educational Management: An Internet of Things (IoT)-Based Machine Learning Approach
Previous Article in Journal
Methodology for the Identification of Vehicle Congestion Based on Dynamic Clustering
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Studies on Soil Quality Index Estimation of a Hilly-Zone Sub-Watershed in Karnataka

by
M. Bhargava Narasimha Yadav
1,*,
P. L. Patil
2 and
M. Hebbara
1
1
Department of Soil Science and Agricultural Chemistry, University of Agricultural Sciences, Dharwad 580 005, India
2
Yettinagudda Campus, University of Agricultural Sciences, Dharwad 580 005, India
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(24), 16576; https://doi.org/10.3390/su152416576
Submission received: 5 September 2023 / Revised: 4 November 2023 / Accepted: 6 November 2023 / Published: 6 December 2023

Abstract

:
The assessment of soil quality aims to evaluate the utility and health of soils. In agricultural studies, soil productivity can be likened to soil quality. Evaluating the Soil Quality Index (SQI) solely based on surface properties offers an incomplete picture because productivity is influenced by both surface and subsurface characteristics, with the latter associated with pedogenic processes. Additionally, relying on weighted averages of soil properties from a soil profile for the SQI may offer an overall summary, but it can occasionally obscure variations that manifest across different soil horizons. Therefore, the present study was conducted to assess the SQI in the Ganjigatti sub-watershed using data from 27 soil profiles and three different methods: (1) assessment of horizon-wise SQI by subjecting the soil properties of every horizon to principal component analysis (PCA), followed by the calculation of the weighted averages of the SQI for each soil profile (SQI-1); (2) calculation of the weighted averages of the soil properties for each soil profile, subjected to PCA, and followed by an SQI assessment (SQI-2); and (3) SQI assessment considering the properties of the Ap horizon for each soil profile (SQI-3). Additionally, to validate SQI methodologies, correlation studies were conducted against major crop yields in the sub-watershed. The results showed that cation exchange capacity (CEC) has the most significant weight and contribution to the SQI determined using MDS, followed by porosity, exchangeable sodium percentage (ESP), organic carbon (OC), CN ratio, and total N. SQI-1 was most strongly correlated with crop yield; the correlation coefficient ranged from 0.69 to 0.74. Among all the three methodologies, SQI-1 and -2 were better methods for assessment of SQI compared to SQI-3. In the SQI-1 method, the soil quality of pedons ranged from 0.26 (pedon-26) to 0.74 (pedon-11). The majority of the area in the sub-watershed (72.40%) fell within the medium category of SQI (0.35–0.55), followed by the high category of SQI (>0.55), which comprised 12.92%, and the low SQI (<0.35), which comprised 6.45% of the sub-watershed.

1. Introduction

Soil is an important non-renewable natural resource on which humans, plants, and animals rely for survival. Research is being conducted globally to identify the best strategies for protecting soils and using them to increase agricultural productivity while preserving environmental quality through improved management methods. On a global scale, available land resources are decreasing at an alarming rate, largely due to a high population pressure [1]. Approximately 10 hectares of land are lost every minute due to various degradation processes such as erosion, nutrient depletion, salinity, acidity, alkalinity, and compaction [2]. The United Nations Environmental Programme (UNEP) has estimated that during the second half of the 20th century, around 2 billion hectares of farmland had suffered degradation [3]. India has a total geographical area of about 328.8 Mha, of which 180 Mha is agricultural land with different types of soils. It supports 17.5% of the world’s population with only 2.4% of the world’s geographical area and 9% arable land. The demand for food, fuel, and energy has increased many folds, and the growing population needs to be fed with shrinking and deteriorating land and water resources [4,5]. Nearly 120 million hectares of arable land in India is reported to be degraded [6].
The detreating of soil quality is one of the main causes of agricultural productivity stagnation and a serious threat to food supply and environment security. An action plan is needed to minimise natural resource degradation and improve soil quality. It should adhere to sustainability principles to ensure that the soil can be passed on to the next generation in improved conditions compared to what was inherited from the previous generation. Therefore, management approaches must align with ecological integrity, economic viability, and social and political approval. At this time, watershed management is an accepted strategy for the development of rainfed agriculture. This approach is multidisciplinary, broad, and intensive [7,8]. A watershed is a geographical and hydrological region that drains to a common point, and it is regarded as a suitable physical entity for assessing, planning, and managing natural resources. Since the 1980s, watershed-focused development has been the primary approach in India’s rainfed regions, aiming to preserve natural resources, boost agricultural output, and enhance rural livelihoods. In this context, assessing the soil quality of soils in watersheds is crucial to managing soil resources for maximum effectiveness in the here-and-now without compromising their viability for the future [9].
The assessment of soil quality is a sensitive and dynamic method for documenting the status of the soil, as well as the soil’s response to management and its resilience to stress, whether that stress is imposed by natural forces or by human interventions. Soil quality is defined as “the capacity of a specific kind of soil to function, within natural or managed ecosystem boundaries, to sustain plant and animal productivity, maintain or enhance water and air quality, and support human health and habitation” [10]. Soil quality, in the context of agricultural output, is defined as the extent to which it can maintain productivity [11]. In essence, soil quality pertains to the functionality of the soil [12,13], whereas soil health portrays the soil as a dynamic, non-renewable, and finite living resource [14,15]. In order to evaluate the status of soil degradation and the shifting patterns that resulted from various land uses and smallholder management interventions, it is required to conduct a fundamental assessment of soil quality [16]. This is the main cause of Africa’s inability to produce as much food as needed to meet demand, and the continent’s per-capita food production is dropping [17,18], mostly as a result of a decline in the quality of its soil.
Soil quality is such a complicated and multifaceted functional concept. It is not possible to directly measure it in the field or in the laboratory; rather, it can only be inferred from the characteristics of the soil. In order to estimate soil quality, a variety of soil parameters or indicators have been identified. There is currently no precise standard or set of criteria that is commonly acknowledged for use in evaluating soil quality [19,20]. Considering the diverse and intricate nature of soil characteristics, numerous assessment approaches have been created, including soil quality models [21], methods for establishing soil quality indices (SQI) [22,23,24], soil health cards and testing kits [25], fuzzy association rules [26], and the soil multifunctionality index [27]. One of these approaches, the soil quality index (SQI), is typically employed for the assessment of soil quality [28] due to its flexibility and ease of use in quantifying changes in various types of soil.
Soil quality indicators capture variability resulting from land management practices, encompassing a range of chemical, physical, and biological soil characteristics. To evaluate the effects of various management practices, it is crucial to establish a baseline or reference values for these soil quality indicators [29]. In one study [24], multiple physical indicators such as potential available water capacity (AWC), soil penetration resistance, bulk density (BD), mean weight diameter (MWD), aggregate size distributions, the fraction of water-stable aggregates (WSA), and geometric mean diameter (GMD) were employed, alongside other chemical indicators, to assess soil quality. Another study [30] evaluated soil quality by considering various biological and physico-chemical soil quality indicators to examine the sustainability of different management and land-use systems. The study identified soil pH, porosity, cation exchange capacity (CEC), available phosphorus (P), BD, total organic carbon (TOC), earthworm population, and plant available water holding capacity (PAWC) as the most responsive indicators. However, total nitrogen (N), exchangeable potassium (K), total phosphorus (P), and K showed moderate sensitivity, while the percentage base saturation was a weaker indicator. In a separate investigation [31], various soil organic carbon (SOC) fractions and the activities of various soil enzymes like dehydrogenase, phosphatase, aryl sulfatase, and fluorescein diacetate hydrolases (FDAse) were employed as biological indicators to assess soil quality in diverse cropping systems in the northwestern Himalayas. Additionally, ref. [32] utilised macroporosity, microporosity, soil hydraulic conductivity (sHC), moisture saturation (MS), effective saturation, aggregate size distribution, aggregate stability index (ASI), exchangeable calcium (Ca) and magnesium (Mg), exchangeable acidity, potential acidity, aluminium saturation, basal respiration, carbon stock (C stock), and nitrogen stock (N stock) as potential soil quality indicators. Further, ref. [33] conducted a study to assess soil quality in the Lakkampura mini-watershed (Karnataka) based on the interaction of soil carbon stocks with other soil parameters. Similarly, ref. [34] assessed the soil quality of the sub-humid southern plains of Rajasthan by using fertility characters. In the Garhwal Himalayas of Uttarakhand, ref. [35] demonstrated the impact of landslides on soil quality. They utilised Principal Component Analysis (PCA) to determine that in resource-constrained situations, soil organic carbon (SOC), available phosphorus (P), and clay content should be selected as the key indicators for tracking fluctuations in soil quality.
Soil quality can be conceptualised in two aspects, viz., inherent and dynamic soil quality. The inherent soil quality shows little change over time whereas dynamic soil quality changes with respect to soil management. The changes in soil properties may occur within hours to a period of decades with respect to the response level of soil properties. However, the limits to which dynamic soil properties can change are dictated by inherent properties [20]. Previously, many evaluations of soil quality focused on the characteristics of the surface soil [24,31,33,34], and there has been limited research utilising information from soil profiles [36,37,38]. Based on the information provided above, a hypothesis was formulated suggesting that analysing the soil quality index (SQI) across different horizons within a profile yields more comprehensive and accurate insights compared to relying on a weighted average of properties for the SQI. While weighted averages offer a broad overview, they can sometimes obscure variations present in distinct horizons. The assessment of SQI horizon-wise enables the examination of soil quality at each individual horizon, facilitating a better grasp of soil quality nuances. The current study was executed with the objectives: (a) identifying soil quality indicators; (b) evaluating the SQI using horizon-wise soil properties, weighted averages of soil properties, and properties of the Ap horizon in each soil profile; and (c) establishing correlations between the SQI and crop yield.

2. Materials and Methods

2.1. Description of the Study Area

The study was conducted in the Ganjigatti sub-watershed, located 37 km away from Dharwad town, Karnataka, during 2021–2022. It is located between 15°10′10.114″ and 15°17′1.147″ N latitudes and 75°0′57.672″ and 75°4′50.525″ E longitudes. This area lies in a hilly zone (zone-9) in Karnataka and a hot semi-arid ecological sub region. It is characterized by ustic soil moisture and isohyperthermic soil temperature regimes (Figure 1) [39]. The study area covers 4323.84 ha and receives an annual average rainfall of 917 mm (average annual rainfall of the zone ranges from 539 to 1256 mm). Annual average kharif (June to September), rabi (October to January), and summer (February to May) rainfall are 616 mm, 139 mm and 162 mm, respectively. The study area’s soils originated from chlorite schist with shale being the predominant parent material. These soils exhibit a coarse texture and are relatively shallow on uplands, and in the lower elements of landscape, they become finer in texture and deeper. The predominant soil types in the area are black and red soils. Black soils are well-suited for cultivating crops such as cotton, wheat, finger millets, sorghum, and oilseeds. Among these, maize and soybean are the primary crops grown during the kharif season, while sorghum and greengram are the main crops cultivated during the rabi season. The length of the growing period (LGP) in this region typically spans from 120 to 150 days.

2.2. Soil Sampling and Analysis

A detailed soil survey of the study area was performed using a high-resolution satellite image (Worldview-2, 50 cm SR) and cadastral maps. The image and scanned topographic sheet were geocoded and subdivided using ArcGIS 10.8.2, resulting in subsets at a scale of 1:7920. During the comprehensive soil survey, landforms and physiographic divisions were identified based on geological characteristics, drainage patterns, surface features, slope attributes, and land usage. Once the landforms were delineated on the satellite image, we conducted thorough traversing of each landform to select representative areas for transect studies. These transects were positioned roughly perpendicular to the contours and aimed to include the variations within each landform. Within each selected transect, profiles were established at closely spaced intervals to account for changes in land features, such as alterations in slope, erosion, presence of gravels, stones, or variations in soil depth. A total of twenty-seven (27) soil profiles were examined in this manner. For each pedon, detailed morphological features were described as per a Soil Survey Manual [39] such as colour, texture, structure, consistency, etc.
Soil samplings were air-dried, ground, sieved (<2 mm), and analysed for soil physical, physico-chemical and chemical properties. Particle size analysis of soil samples was carried out using the international pipette method [40]. Bulk density of soil samples was determined using the clod method [41]. The particle density was assessed using the specific gravity bottle method as outlined in reference [42]. The total porosity was subsequently computed based on the bulk density and particle density. The maximum water-holding capacity and volume expansion were determined using Keen Raczkowski’s method [43]. Moisture retention was determined at −33 kPa and −1500 kPa using pressure plate apparatus and pressure membrane, respectively and available soil moisture content was determined from the difference in moisture held at −33 kPa and −1500 kPa [44].
The pH and electrical conductivity (EC) were measured employing 1:2.5 soil:water ratio [45]. Organic carbon (OC) was determined using the method of [46]. CaCO3 equivalent (%) was determined using the method described by [40]. Soil inorganic carbon was assessed using the formula in [40]. Exchangeable cations were assessed through a soil extraction process involving 1 N neutral ammonium acetate. The concentrations of calcium (Ca) and magnesium (Mg) in the extract were determined through Versenate titration, while sodium (Na) and potassium (K) concentrations were determined using flame photometry. The quantities of exchangeable cations, specifically Ca2+, Mg2+, K+, and Na+, were expressed in units of cmol (p+) per kilogram of soil, as detailed in [47]. Additionally, the cation exchange capacity (CEC) of the soils was determined utilising a method involving 1 N sodium acetate at pH 8.2, as outlined in [48]. Base saturation was calculated by summing the bases extracted by dividing by CEC (by 1 N sodium acetate) and multiplying by 100. Exchangeable acidity (Al+3 + H+) of soil samples was determined using the KCl extractant method, as described by [40]. Exchangeable sodium percentage (ESP) was estimated as the ratio of exchangeable sodium to CEC. Available nitrogen was determined using the modified alkaline permanganate method as described by [49]. The total nitrogen (%) was determined using the Macrokjeldahl digestion cum distillation method as described by [45]. C:N ratio was estimated as the ratio of total carbon (%) to total nitrogen (%). Descriptive statistics for 34 parameters were performed with SPSS ver. 26.

2.3. Assessment of Soil Quality Index

A statistical model based on principal component analysis (PCA) was employed to estimate Soil Quality Index (SQI), as described in [10,22,28,50,51,52,53]. The SQI calculation involved four key steps: (a) defining the management objective, (b) selecting specific indicators as Minimum Data Set (MDS), (c) assigning scores to the chosen indicators, and (d) computing the SQI.

2.3.1. Indicator Selection

The PCA model was applied to establish a Minimum Data Set (MDS) with the aim of reducing the burden of indicators in the model and preventing data duplication [22]. The PCA method is preferred due to its objectivity, as it incorporates various statistical tools such as multiple correlations, factor analysis, and cluster analysis. These tools help in mitigating any potential bias and data redundancy by employing mathematical formulas to select an MDS, as elaborated in references [11,28]. The SQI of the sub-watershed was assessed using three methodologies: (1) horizon-wise SQI assessment followed by weighted averages of the SQI of each soil profile; (2) weighted averages of the soil properties of each soil profile subjected to PCA and followed by SQI assessment; and (3) considering the soil properties of the Ap horizon for SQI assessment. The primary purpose of PCA is to reduce the complexity of a comprehensive dataset containing numerous interconnected variables while preserving as much of the inherent variability as possible. This is accomplished by transforming the data into a fresh set of variables known as principal components (PCs), which are both uncorrelated and arranged in a manner that ensures the initial ones capture the majority of the variation observed across all the original variables [54,55]. In essence, the PCA method was selected as a means of streamlining the data and identifying the most suitable indicator(s) for representing and estimating SQI [8,55].
Principal components (PCs) in a dataset are combinations of variables that capture the maximum variation within the dataset. They describe vectors that closely fit the n observations in a p-dimensional space while being mutually orthogonal. In this study, the approach outlined by [22] was followed. In previous studies, it has been common to assume that PCs with high eigenvalues and variables with high factor loadings best represent system attributes. Consequently, only PCs with eigenvalues greater than 1 were chosen. For each PC, each variable had a corresponding eigenvector weight value or factor loading. Only the variables with the highest weights were retained for inclusion in the Minimum Data Set (MDS). Variables were considered “highly weighted” if they had the highest weight under a specific PC and an absolute factor loading value within 10% of the highest values under the same PC [56]. However, when multiple variables were retained under a specific PC, a multivariate correlation matrix was used to calculate correlation coefficients between these parameters, as detailed in references [22,28]. If the parameters were significantly correlated (with a correlation coefficient r > 0.60 and p < 0.05), then the one with the highest loading factor was included in the MDS, while all others were excluded to prevent redundancy. Parameters that were not correlated under a particular PC were considered important and retained in the MDS, as explained by [28,57].

2.3.2. Scoring of Indicators

The indicators selected for the Minimum Data Set (MDS) were transformed into dimensionless values ranging from 0 to 1 using a linear scoring method [58]. The ranking of indicators was determined by whether a higher value was deemed favourable or unfavourable concerning soil function. For indicators where a higher value indicated better performance, each indicator value was divided by the highest value, ensuring that the highest value received a score of 1. Conversely, for indicators where lower values were preferable, the lowest value was divided by each data value so that the lowest value received a score of 1. In the case of indicators like pH, an “optimum” threshold value was considered. These indicators were initially scored as “higher is better” up to the threshold value (for example, pH 7.5), and beyond the threshold, they were scored as “lower is better”. This scoring method is consistent with the approach outlined in references [22,28].

2.3.3. Soil Quality Index Calculation

Following that, the chosen observations were converted into numerical scores within the range of 0 to 1. These scores were then combined into indices for each soil sample through a weighted additive method. To determine the specific weighting value assigned to each Principal Component (PC), the variance explained by each PC was divided by the maximum total variation among all PCs generated by the PCA.
SQI = Σ Principal Component Weight ∗ Individual soil parameter score
Principal Component Analysis (PCA) and descriptive statistics were carried out using SPSS version 26 and Microsoft Excel version 16 to compute the Soil Quality (SQ) indices. Crop yield information spanning an eight-year period (2014–2021) was sourced from the Department of Agriculture, Government of Karnataka, Dharwad. The average yield data for soybean, maize, sorghum, and greengram were collected on soil phases basis and then correlated with the Soil Quality Index (SQI).

3. Results

3.1. Physical, Physico-Chemical, and Chemical Properties

Descriptive statistics for the 34 soil properties used for soil quality assessment are given in Table 1, Table 2 and Table 3. Soil depth among the pedons varied from 20 to 200+. The soil texture of the Ganjigatti sub-watershed is mostly limited to clay, sandy clay, and sandy clay loam, with few exceptions. Clay illuviation is an important soil-forming process. The movement of clay has an important pedological significance as clay bears a good relationship with other properties. Coefficient of variation (CV) is used to interpret soil variability [59].

3.2. Assessment of Soil Quality Index of Ganjigatti Sub-Watershed

3.2.1. Selection of Soil Quality Indicators

To determine the minimum data set (MDS) and key factors affecting soil quality index, PCA was carried out on the total data set. The PCA method is a common and accepted approach for determining MDS and has been widely used in previous studies [22,60,61,62]. The SQI of this sub-watershed was assessed using three methodologies: (1) assessment of horizon-wise SQI by subjecting the soil properties of every horizon to PCA followed by calculating the weighted averages of the SQI of each soil profile (SQI-1); (2) weighted averages of the soil properties of each soil profile subjected to PCA and followed by SQI assessment (SQI-2); and (3) considering the Ap horizon properties of each soil profile for the SQI assessment (SQI-3).
In the SQI-1 method, the soil properties of all 27 profiles horizon-wise were subjected to PCA to reduce the data dimension. The PCA data for the sub-watershed showed that all the eigenvalues were >1, which explained 85.658% of the cumulative variance in the data. The MDS were chosen based on the highly weighted factor loading of variables. In total, eight PCs were extracted for the sub-watershed with eigenvalues >1 (Table 4). A representative scree plot showing the variation of eigenvalues with soil components is shown in Figure 2. The parameters in each PC were considered based on higher values of the factor loading. The soil parameters obtained from PCA under PC1 were the sum of exchangeable basic cations, exchangeable Ca, exchangeable Mg, pH, CEC, field capacity, volume expansion, permanent wilting point, CaCO3, and inorganic carbon. However, a multivariate correlation matrix was utilised to calculate the correlation coefficients between the parameters when more than one variable was retained under a given PC [22,28]. To avoid redundancy, only the parameter with the highest loading factor was kept in the MDS if there was a significant correlation between them (r > 0.60, p 0.05). The non-correlated parameters under a particular PC were considered important and retained in the MDS [22,57]. Among these highly weighted variables, CEC is a parameter that governs nutrients’ holding capacity and is an indicator of soil fertility. It is a very important soil property, influencing the soil’s structural stability, nutrient availability, and pH. Other parameters are highly correlated to each other, so CEC was retained for MDS in PC1 (Table 5). In PC2, porosity was retained under MDS between bulk density and porosity because of its higher factor loading and better correlation (Table 6). From PC3, exchangeable Na, ESP, and CN ratio were selected, and exchangeable Na and CN ratio were considered for MDS based on correlation (Table 7). Bulk density, CEC/clay ratio, exchangeable acidity, silt, and total carbon were selected as indicators from PC4, PC5, PC6, PC7, and PC8, respectively. Some of the selected parameters were not independent of each other. Exchangeable acidity is correlated with porosity, CN ratio, ESP, BD, and silt (Table 8). CEC/clay ratio is a derived parameter that is dependent on the CEC of the soil [63].
In the SQI-2 method, the weighted average of the soil properties of all the 27 profiles subjected to PCA was used to reduce the data dimension. The PCA data for the sub-watershed show the eight PCs with eigenvalues > 1, which explain 90.128% of the cumulative variance in the data (Table 9). A representative scree plot showing the variation of eigenvalues with soil components is shown in Figure 3. The soil parameters selected from PC1 were the sum of exchangeable basic cations, exchangeable Ca, exchangeable Mg, pH, CEC, CaCO3, and inorganic carbon. However, multivariate correlation between these parameters indicated a high correlation (Table 10), and only CEC was retained in the MDS because it is a very important soil property influencing soil structural stability, nutrient availability, and soil pH. In PC2, total N was retained under MDS between available N and total N because of its higher factor loading after correlation results. In PC4, bulk density was retained under MDS between bulk density and porosity because porosity is a derived parameter that is dependent on the BD and PD of the soil. CN ratio, exchangeable acidity, exchangeable Na, particle density, and total carbon were selected as indicators from PC3, PC5, PC6, PC7, and PC8, respectively, since they were the only highly weighted parameters. The correlation data between the selected indicators of MDS are presented in Table 11.
In the SQI-3 method for soil properties, all 27 soil profiles of the Ap horizon were subjected to PCA to reduce the data dimension. The PCA data for the sub-watershed show the eight PCs with eigenvalues > 1, which explain 87.265% of the cumulative variance in the data. The MDS were chosen based on the highly weighted factor loading of variables (Table 12). A representative scree plot showing the variation of eigenvalues with soil components is shown in Figure 4. The soil parameters selected from PC1 were the sum of exchangeable basic cations, exchangeable Ca, exchangeable Mg, and CEC. However, multivariate correlation between these parameters indicated a high correlation, and only CEC was retained in the MDS (Table 13). From PC2, organic carbon, available N, and total N were selected, and organic carbon was considered because of its higher factor loading after correlation results (Table 14). From PC4, bulk density, porosity, and CN ratio were selected, and bulk density and CN ratio were considered based on correlation (Table 15). ESP, porosity, particle density, total carbon, and EC were selected as indicators from PC3, PC5, PC6, PC7, and PC8, respectively, since they were the only highly weighted parameters. The correlation data between the selected indicators of MDS are presented in Table 16.

3.2.2. Normalisation of MDS Values

After the selection of parameters for the MDS, all selected observations were transformed using linear scoring functions (less is better, more is better, and optimum) as described in the Materials and Methods. In this study, CEC, CEC/clay ratio, organic carbon, total carbon, total nitrogen, silt, and porosity were considered as more being good from the point of view of soil quality when they were in increasing order; hence, the “more is better” approach was followed. ESP, exchangeable Na, exchangeable acidity, EC, particle density, and bulk density were considered as being less value is better for soil quality; hence, the “less is better” approach was adopted. In the case of the CN ratio, the “optimum is better” approach was followed. After converting the chosen observations into numerical scores within the 0 to 1 range, a weighted additive method was employed to combine them into indices for each soil sample [24,28]. Each Principal Component (PC) explained a specific portion of the variation in the dataset (Table 4, Table 9 and Table 12). To assign a particular weight value to each PC, the variation explained by each PC was divided by the total maximum variation among all PCs selected for the Minimum Data Set (MDS) [28]. For instance, in the case of SQI-1 % variance, which had a value of 48.237, it was divided by the cumulative total variance (85.658) to yield a weight value of 0.563 for PC-1 (Table 4).

3.2.3. SQI Calculation

Thereafter, to obtain the weighted additive SQI, the weighted MDS indicator scores for each observation were summed up using the following equation:
SQI = Σ Principal Component Weight ∗ Individual soil parameter score
The soil quality index of 27 soil profiles in the Ganjigatti sub-watershed was assessed using three methodologies and is presented in Table 17. In the SQI-1 method, the soil quality of pedons ranged from 0.26 (pedon-26) to 0.74 (pedon-11). Pedons 12, 22, and 26 were classified under low SQI (<0.35); pedons 4, 6, 9, 11, 19, 23, and 27 were classified under high SQI (>0.55); and the remaining pedons were under medium SQI (0.35–0.55). All three low SQI pedons were upland pedons, which are characterised by a shallow depth (<35 cm), a low CEC, a low clay content, a low CEC/clay ratio, an acidic pH, and a high exchangeable acidity. Interestingly, all three of these pedons were classified as Entisols. In the high SQI pedons, pedons 4 and 9 were moderately deep (75 to 100 cm), and the remaining pedons were deep to very deep (100 to 200 cm). These pedons were characterised by a high CEC, high clay content, high porosity, low BD, and low exchangeable acidity. Pedons 4 and 9 were located midland, and the remaining pedons (6, 11, 19, 23, and 27) were located in low-land topography positions. The SQI-1 values of the Ganjigatti sub-watershed depicted in maps using GIS software (Figure 5) showed that the majority of the area is under the medium category of SQI (0.35–0.55), comprising about 72.40% of the sub-watershed area. The high category of SQI (>0.55) comprised 12.92%, and low SQI (<0.35) comprised 6.45% of the sub-watershed.
In the SQI-2 method, the soil quality of pedons ranged from 0.32 (pedon-26) to 0.85 (pedon-11). The spatial distribution map of SQI using the SQI-2 method (Figure 6) revealed that the low, medium, and high categories of SQI covered 2.18, 54.11, and 35.49% of the total geographical area of the Ganjigatti sub-watershed, respectively. In the SQI-3 method, the soil quality ranged from 0.39 (pedon-26) to 0.89 (pedon-4), and the spatial distribution map of SQI-3 (Figure 7) indicated that 81.04% of the area was classified under the high category of SQI.

3.2.4. Pearson’s Correlation Coefficients of Soil Quality Index (SQI) versus Crop Yield

The primary aim of this research was to assess the validity of Soil Quality Index (SQI) methodologies in relation to crop yield, which is an aspect that has been relatively understudied with a few exceptions noted in references [22,24,36,38,50]. The challenge lies in establishing a correlation between crop yields and SQI, whether measured at different soil depths or as a combination of soil layers, as there is currently no standardised technique available, and this area has received limited attention until now. However, in order to evaluate the index performance in the current study, SQIs were linked to the end-point variable (i.e., yield). The estimated SQI values were correlated with the recorded yields of soybean, maize, sorghum, and greengram in the Ganjigatti sub-watershed. Results of the present study revealed that (Table 18) yields of all four crops were significantly correlated with all three SQI methodologies.

4. Discussion

SQI is a composite of a few selected soil indicator characteristics, and it necessitates the selection of the most relevant qualities with a dominant influence on soil functions. It is noteworthy that the MDSs derived from SQI-1, SQI-2, and SQI-3 were very similar, with the exception of the presence or absence of two or three different parameters and the order (or) weightage being different. The SQI computation and ranking of soils based on SQI values reflected these results. It has been argued that a more comprehensive data collection or a larger subset of indicators may more accurately indicate soil quality; nevertheless, this can lead to data duplication when there is a high correlation between the indicators used [64,65]. ESP and CN ratio were included under MDS from PC3 of the SQI-1 method; B.D and CN ratio were included under MDS from PC4 of the SQI-3 method since they were not correlated. Within the MDS of three distinct methods, Cation Exchange Capacity (CEC) had the highest influence on the Soil Quality Index calculation, followed by porosity, Exchangeable Sodium Percentage (ESP), organic carbon, CN ratio, and total nitrogen. These factors have been frequently recognised as significant and sensitive variables in the development of the SQI, as indicated in previous studies [24,30,36,38,62,66,67,68].
Each of these factors plays a pivotal role in influencing the amalgamation of soil physicochemical and biological attributes, as well as soil fertility, productivity, and the various components that contribute to crop yield [69]. CEC is a very important parameter for the Ganjigatti sub-watershed because it influences the nutrient-supplying capacity of soils as it depends on the quantity and type of clay, soil pH, and organic matter [38,70]. On the other hand, CEC provides indications about the clay mineralogy of the soil, which is also responsible for the quality of soil [37]. Porosity significantly influences soil quality by affecting water-holding capacity, soil aeration, nutrient availability, root growth, and microbial activity. Mainly, pore spaces facilitate the availability and movement of air or water within the soil environment. Soil porosity and its capacity to retain water are additional soil characteristics that contribute to creating a favourable environment for microorganisms and moisture retention in the soil. Consequently, this promotes the enhanced decomposition of soil organic carbon by microbial communities. Soil porosity offers space for microbial proliferation and enhances soil aeration, nutrient accessibility, as well as the soil’s capacity for both water drainage and retention [37,71]. Maintaining a well-structured and porous soil allows for better water and nutrient management, supports root development, and promotes the activity of beneficial soil organisms, ultimately leading to improved plant growth and productivity.
Organic carbon is interlinked to all measures of soil quality through a complex web of many factors. This factor has a vital function in the cycling and storage of essential soil nutrients, contributes to the formation of soil structure, and serves as the primary nutrient source for heterotrophic microorganisms within the soil [69]. The C:N ratio is an important parameter that affects soil quality and nutrient availability. It is a measure of the relative amounts of C and N in soil. A higher ratio can lead to nitrogen immobilisation and slower decomposition, while a lower ratio promotes faster decomposition and nutrient release. Achieving a balanced C:N ratio that meets the specific needs of plants can contribute to a healthier soil and improved crop productivity [72,73]. ESP is a measure of the sodium content relative to other cations in the soil. It is an important parameter for assessing soil quality, particularly in terms of soil structure, permeability, and fertility. High ESP in soil can degrade soil structure, reduce water infiltration and drainage, affect nutrient availability, and alter soil pH. Managing high ESP levels through appropriate soil management practices is crucial for maintaining soil quality and supporting healthy plant growth [36,74,75]. Total nitrogen (TN) is a critical component of soil quality and plays a fundamental role in supporting plant growth and overall ecosystem functioning. It is important to note that while nitrogen is crucial for plant growth and soil quality, its excessive application or imbalanced ratios with other nutrients can lead to environmental issues such as nitrogen leaching, eutrophication of water bodies, and air pollution. In summary, total nitrogen is a key factor in soil quality, affecting nutrient availability, plant growth and yield, soil fertility, organic matter decomposition, and microbial activity.
Results on the correlation between crop yields and SQI methods revealed that the yields of all four crops were significantly correlated with all three SQI methodologies. Several studies [38,62,68,76,77] have also found a significant correlation between the yield components of various agricultural products and SQIs with different coefficients of variations. A relatively higher correlation coefficient was observed with SQI-1 and -2 than with SQI-3 (Table 17), which indicates the pedological significance and the effect of soil subsurface phenomena on the physiological conditions of plant systems, which take nutrients and water from the subsurface. Since crop production is influenced by both surface and subsurface attributes, which are inherently linked to pedogenic processes, the soil quality index (SQI) assessment using solely surface soil parameters does not provide enough information [78]. However, inherent properties limit the extent to which dynamic soil properties can change [79]. The pedogenic processes have an impact on the inherent properties of the soil, and the changes are particularly noticeable in tropical climates because of physical and chemical weathering that is accelerated by high temperatures and precipitation. Many of the early assessments of soil quality [22,23] used surface (dynamic) soil parameters, while studies utilising soil profile data (dynamic and inherent) are scarce [36,37]. Although surface soil characteristics are simple to measure and assess, they only provide partial information because pedogenic processes in the soil control section drive soil functions. The soil properties that have the higher impact on soil functions can be determined by evaluating the soil quality utilising both surface and subsurface properties. Soil profile parameters that are determined through soil genesis and reflected by taxonomy are the only ones that should be used to evaluate soil quality [80,81]. Therefore, pedogenesis must be taken into consideration when assessing the soil quality in SAT soils, and both the surface and subsurface soil properties should be given the appropriate weight.
In order to improve and restore the soil’s health, ref. [36] used the SQI to evaluate the soil quality of the Indo-Gangetic Plains. There was a moderate association between SQI and yields in the rice-wheat system. They drew a conclusion using the trial-and-error method by taking into account the percentage contribution of SQI for each layer in the calculation. Finally, they came to the conclusion that a composite SQI value was produced by combining 70% of the surface SQI value and 30% of the subsurface SQI value. In the present study, a good and significant correlation between crop yields and the SQI-3 method was probably due to the fact that all four crops (soybean, maize, sorghum, and green gram) considered in this study are shallow-rooted, and their root density is in the top layer of the soil, which can drive the nutrients from the surface horizons. In view of this, we conclude that when subsurface soil variables are considered alongside dynamic surface properties while evaluating the SQI using the weighted index approach, a good correlation between the SQI and defined soil function is established [36,38,80,81,82].
Furthermore, assessing the horizon-wise soil quality index (SQI-1) of a profile can provide more detailed and accurate information compared to using a weighted average (SQI-2). Weighted averages may provide an overall summary measure but can sometimes mask variations that occur at different horizons. Assessing the SQI horizon-wise allows for a more granular understanding of the profile’s soil quality. Here are a few reasons why horizon-wise assessment is advantageous. (1) A soil profile consists of different layers or horizons, each with its unique characteristics and properties. By assessing the SQI horizon-wise, you can account for variations in soil quality across different horizons. This helps capture the heterogeneity within the profile and provides a more comprehensive understanding of soil quality. (2) Soil quality problems, such as compaction, erosion, or nutrient depletion, may occur at specific horizons rather than being evenly distributed across the entire profile. Evaluating the SQI horizon-wise enables the identification of specific horizons that may require attention or management interventions. (3) Horizon-wise assessment allows for more targeted and site-specific management decisions. For example, if the SQI is low in a particular horizon, you can focus on implementing appropriate soil management or nutrient management strategies specific to that horizon. (4) Soil horizons interact with each other, influencing the overall soil quality of the profile. By examining the SQI horizon-wise, you can gain insights into how different horizons interact and affect the overall soil quality. This knowledge is crucial for understanding the complex dynamics within the soil profile. (5) For the assessment of the SQI through PCA, it is better to consider the soil properties of every horizon (SQI-1) instead of the weighted average of the soil properties of a profile (SQI-2). Since PCA uses a number of statistical tools, it could avoid any bias and data redundancy by choosing an MDS using mathematical formulas from huge data sets [11,28]. In conclusion, while weighted averages (SQI-2) can be useful for providing an overall summary measure, horizon-wise assessment of the Soil Quality Index (SQI-1) provides a more detailed and accurate representation of soil quality by accounting for variations, specific issues, targeted management decisions, and interactions between different horizons. The current data suggest that the SQI values were relatively higher in the SQI-2 and SQI-3 methods than in SQI-1. However, SQI-1 appears to be the best method among the three, particularly under the long-term scenario, especially due to its objective approach and relatively higher correlation with crop yield.
According to the SQI-1 method, soil quality in the study area varied from low to high. The large variation in soil quality is due to soil heterogeneity and soil degradation caused by erosion. The SQI of the study area ranged from 0.26 to 0.74, and the spatial distribution of the SQI of the study area (Figure 5) shows that about 6.45% of the sub-watershed has a low category of SQI (<0.35). The low SQI of the Ganjigatti sub-watershed might be due to unfavourable microclimatic conditions, soil erosion, inappropriate land use and management practices, and the high removal of available nutrients by annual crops without replenishment. These soils are located in upland positions with moderate to steep slopes, and, hence, the leaching of basic cations and erosion of the soil to lowlands could be a possible reason for the low soil quality. The cultivation of deep-rooted crops like cotton and pigeon pea is not suitable in these soils because of the limitation imposed by shallow depths (<35 cm). The SQI of pedons 1, 2, 3, 5, 7, 8, 10, 13, 14, 15, 16, 17, 18, 20, 21, 24, and 25 was in the medium category (0.35–0.55), which is the bulk of the area covered in the sub-watershed and comprising about 72.40%. The medium SQI might be due to soil erosion, land use, and management practices fallowed in that area. On the other hand, a few mid- (4 and 9) and low-land pedons (6, 11, 19, 23 and 27) showed a high SQI (>0.55) because of the deposition of soil and nutrients from uplands through the soil erosion process. These soils occupy a sizable area of 12.92% of the sub-watershed.
Soil quality can be used to evaluate cropping systems and recommend alternatives in a particular region [83]. Climate change affects short-term and long-term soil processes, which must be considered when establishing management measures to maintain soil resources and sustain agricultural output. Most crop models, such as InfoCrop and CERES-Wheat, are driven by biophysical parameters, rainfall variability, water balance, and economic implications [84,85] and pay little attention to soil quality [86], especially in India. According to [87], alterations in soil quality over time serve as indicators of whether the soil condition is sustainable or not. Maintaining soil quality at the desired level presents a complex challenge due to the interplay of climatic conditions, soil properties, plant factors, and human activities. This challenge is particularly pronounced in lowland rice cropping systems, largely owing to practices involving soil puddling during preparation [88]. The SQI includes many soil properties as indicators of soil quality, so integrating them into simulation models to predict the effects of climate change on soil functions and crop yield will improve the knowledge and accuracy of models, allowing for better management practices.

5. Conclusions

The Soil Quality Index is a valuable tool for assessing and managing soil health, which is fundamental for sustainable agriculture and environmental preservation. The study highlights the need for a holistic approach that considers both surface and subsurface soil properties, factors that impact soil quality.
1. The current study, within the context of the MDS-determined SQI, has demonstrated that CEC holds the most substantial influence and impact, with porosity, ESP, OC, CN ratio, and total N following in order of importance.
2. By assessing the SQI horizon-wise, we can analyse the quality of the soil at each individual horizon separately. Furthermore, the horizon-wise assessment of SQI provided reliable frameworks for soil quality evaluation in Ganjigatti sub-watershed in the hilly zone of Karnataka. The soil quality within the study area showed significant spatial variability, ranging from low to high SQI. In the sub-watershed, the largest portion of the area, accounting for 72.40%, is categorised as having a medium SQI (ranging from 0.35 to 0.55), while 12.92% is classified as having a high SQI (greater than 0.55), and 6.45% is in the low SQI category (less than 0.35). This approach allows for a more granular understanding of the soil quality.
3. There is a significant correlation between crop yields and all three SQI methodologies, with SQI-1 and SQI-2 showing higher coefficients. This underscores the importance of considering both surface and subsurface soil properties when assessing the SQI.
Factors like soil erosion, microclimatic conditions, land use, and management practices play a crucial role in determining soil quality. SQI can guide management practices in agriculture by identifying areas with a low soil quality and suggesting alternative cropping systems.

Author Contributions

P.L.P. and M.H., designed the study, data validation, data interpretation, and final manuscript editing; M.B.N.Y., analysis, data validation and interpretation, and final manuscript preparation. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

No applicable.

Informed Consent Statement

No applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The study is part of the REWARD (Rejuvenating Watershed for Agricultural Resilience through innovative Development) project funded by World Bank. The authors duly acknowledge the financial support. The authors would also like to thank the anonymous reviewers for their valuable time, suggestions, and critical comments, which helped to improve the quality of this research paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Foley, J.A.; Ramankutty, N.; Brauman, K.A.; Cassidy, E.S.; Gerber, J.S.; Johnston, M.; Mueller, N.D.; O’Connell, C.; Ray, D.K.; West, P.C.; et al. Solutions for a cultivated planet. Nature 2011, 478, 337–342. [Google Scholar] [CrossRef] [PubMed]
  2. UNCCD. Land and Soil the Context of a Green Economy for Sustainable Development, Food Security and Poverty Eradication. The Submission of the UNCCD Secretariat to the Preparatory Process for Rio+ 20. 2011. Available online: http://www.unscd2012.org/content/documents/unccd.pdf (accessed on 5 September 2023).
  3. Oldemann, L.R.; Hakkeling, R.T.A.; Sombrock, W.G. World Map of the Status of Human Induced Land Degradation; An Explantory Note; International Soil Reference and Information Centre: Wageningen, The Netherlands, 1990. [Google Scholar]
  4. Sarkar, D.; Chaturvedi, A.; Mahapatra, S.K. New vistas in perspective of land use planning. J. Indian Soc. Soil Sci. 2009, 57, 587–600. [Google Scholar]
  5. Lal, R. Soil and sanskriti. J. Indian Soc. Soil Sci. 2013, 61, 267–274. [Google Scholar]
  6. Maji, A.K.; Obi Reddy, G.P.; Sarkar, D. Degraded and Wastelands of India—Status and Spatial Distribution; A Joint Publication of ICAR and NAAS: New Delhi, India, 2010. [Google Scholar]
  7. Narayana, V.V.D.; Sastry, G.; Patnaik, U.S. Watershed Management; Indian Council of Agricultural Research: New Delhi, India, 1990. [Google Scholar]
  8. Patil, P.L.; Yadav, M.B.N. Geospatial technologies for watershed development. In Proceedings of the 31st National Conference on Innovative Resource Management Approaches for Coastal and Inland Ecosystems to Sustain Productivity and Climate Resilience, Navsari Agricultural University, Navsari, Gujarath, India, 13–15 October 2022; Soil Conservation Society of Indian: New Delhi, India, 2022; pp. 16–35. [Google Scholar]
  9. Dhawan, A.S. Soil quality in relation to agricultural sustainability: A road map a head. J. Indian Soc. Soil Sci. 2014, 62, 70–77. [Google Scholar]
  10. Karlen, D.L.; Mausbach, M.J.; Doan, J.W.; Cline, R.G.; Harris, R.F.; Schuman, G.E. Soil quality: A concept, definition and framework for evaluation. Soil Sci. Soc. Am. J. 1997, 61, 4–10. [Google Scholar] [CrossRef]
  11. Doran, J.W.; Parkin, T.B. Quantitative indicators of soil quality: A minimum data set. Soil Sci. Soc. Am. J. Spec. Publ. 1996, 49, 25–38. [Google Scholar]
  12. Harris, R.F.; Bezdicek, D.F. Descriptive aspects of soil quality/health. In Defining Soil Quality for Sustainable Environment; Doran, J.W., Coleman, D.C., Bezdicek, D.F., Stewart, B.A., Eds.; Soil Science Society of America: Madison, WI, USA, 1994; pp. 23–35. [Google Scholar]
  13. Karlen, D.L.; Andrews, S.S.; Wienhold, B.J.; Doran, J.W. Soil quality: Human kind’s foundation for survival. J. Soil Water Conserv. 2003, 58, 171–179. [Google Scholar]
  14. Letey, J.; Sojka, R.E.; Upchurch, D.R.; Cassel, D.K.; Olson, K.; Payne, B.; Petrie, S.; Price, G.; Reginato, R.J.; Scott, H.D.; et al. Deficiencies in the soil quality concept and its application. J. Soil Water Conserv. 2003, 58, 180–187. [Google Scholar]
  15. Doran, J.W.; Zeiss, M.R. Soil health and sustainability: Managing the biotic component of soil quality. Appl. Soil Ecol. 2000, 15, 3–11. [Google Scholar] [CrossRef]
  16. Lal, R.; Stewart, B.A. Soil management: Experimental basis for sustainability and environmental quality. In Advanced Soil Science Series; CRC Press: Boca Raton, FL, USA, 1995; pp. 1–67. [Google Scholar]
  17. Lal, R.; Blum, W.E.H.; Valentin, C.; Stewart, B.A. (Eds.) Methods for Assessment of Soil Degradation; CRC Press: Boca Raton, FL, USA, 1997; p. 554. [Google Scholar]
  18. Lal, R. Basic concepts and global issues: Soil quality and agricultural sustainability. In Soil Quality and Agricultural Sustainability; Ann Arbor Science: Chelsea, MI, USA, 1998; pp. 3–12. [Google Scholar]
  19. Mukhopadhyay, S.; Masto, R.; Yadav, A.; George, J.; Ram, L.; Shukla, S. Soil quality index for evaluation of reclaimed coal mine spoil. Sci. Total Environ. 2016, 542, 540–550. [Google Scholar] [CrossRef]
  20. Biswas, S.; Hazra, G.C.; Purakayastha, T.J.; Saha, N.; Mitran, T.; Roy, S.S.; Basak, N.; Mandal, B. Establishment of critical limits of indicators and indices of soil quality in rice–rice cropping systems under different soil orders. Geoderma 2017, 292, 34–48. [Google Scholar] [CrossRef]
  21. Larson, W.E.; Pierce, F.J.; Doran, J.; Coleman, D.; Bezdicek, D.; Stewart, B. The dynamics of soil quality as a measure of sustainable management. In Defining Soil Quality for a Sustainable Environment; Doran, J.W., Coleman, D.C., Bezdicek, D.F., Stewart, B.A., Eds.; Soil Science Society of America, Inc.: Madison, WI, USA, 1994; pp. 37–51. [Google Scholar]
  22. Andrews, S.S.; Karlen, D.L.; Mitchell, J.P. A comparison of soil quality indexing methods for vegetable production systems in Northern California. Agric. Ecosyst. Environ. 2002, 90, 25–45. [Google Scholar] [CrossRef]
  23. Armenise, E.; Redmile-Gordon, M.A.; Stellacci, A.M.; Ciccarese, A.; Rubino, P. Developing a soil quality index to compare soil fitness for agricultural use under different managements in the Mediterranean environment. Soil Tillage Res. 2013, 130, 91–98. [Google Scholar] [CrossRef]
  24. Mukherjee, A.; Lal, R. Comparison of soil quality index using three methods. PLoS ONE 2014, 9, e105981. [Google Scholar] [CrossRef] [PubMed]
  25. Ditzler, C.A.; Tugel, A.J. Soil quality field tools: Experiences of USDA-NRCS soil quality institute. Agron. J. 2002, 94, 33–38. [Google Scholar] [CrossRef]
  26. Xue, Y.J.; Liu, S.G.; Hu, Y.M.; Yang, J.F. Soil quality assessment using weighted fuzzy association rules. Pedosphere 2010, 20, 334–341. [Google Scholar] [CrossRef]
  27. Singh, A.K.; Rai, A.; Banyal, R.; Chauhan, P.S.; Singh, N. Plant community regulates soil multifunctionality in a tropical dry forest. Ecol. Indic. 2018, 95, 953–963. [Google Scholar] [CrossRef]
  28. Andrews, S.S.; Mitchell, J.P.; Mancinelli, R.; Karlen, D.L.; Hartz, T.K.; Horwath, W.R.; Pettygrove, G.S.; Scow, K.M.; Munk, D.S. On-farm assessment of soil quality in California’s Central valley. Agron. J. 2002, 94, 12–23. [Google Scholar]
  29. Bünemann, E.K.; Bongiorno, G.; Bai, Z.; Creamer, R.E.; De Deyn, G.; de Goede, R.; Fleskens, L.; Geissen, V.; Kuyper, T.W.; M€ader, P.; et al. Soil quality—A critical review. Soil Biol. Biochem. 2018, 120, 105–125. [Google Scholar] [CrossRef]
  30. Ezeaku, P.I. Evaluation of agro-ecological approach to soil quality assessment for sustainable land use and management systems. Acad. J. 2015, 10, 501–512. [Google Scholar]
  31. Sofi, J.A.; Bhat, A.G.; Kirmai, N.A.; Wani, J.A.; Lone, A.H.; Ganie, M.A.; Dar, G.I. Soil quality index as affected by different cropping systems in North Western Himalayas. Environ. Monit. Assess. 2016, 188, 161. [Google Scholar] [CrossRef] [PubMed]
  32. Stefanoski, D.C.; de Figueiredo, C.C.; Santos, G.G.; Marchão, R.L. Selecting soil quality indicators for different soil management systems in the Brazilian Cerrado. Pesqui. Agropecu. Bras. 2016, 51, 1643–1651. [Google Scholar] [CrossRef]
  33. Sathish, A.; Madhu, L. Assessment of soil quality based on relationship of soil carbon stocks with different soil parameters in Lakkampura mini-watershed in Davanagere under semi-arid Alfisols in India. J. Indian Soc. Soil Sci. 2021, 69, 354–368. [Google Scholar] [CrossRef]
  34. Yadav, K.K.; Mali, N.L.; Kumar, S.; Surya, J.N.; Moharana, P.C.; Nogiya, M.; Meena, R.L. Assessment of soil quality and spatial variability of soil properties using geo-spatial techniques in sub-humid Southern Plain of Rajasthan, India. J. Indian Soc. Soil Sci. 2022, 70, 69–85. [Google Scholar] [CrossRef]
  35. Goyal, D.; Joshi, V.; Gupta, N.; Cabral-Pinto, M.M.S. Soil quality assessment in a landslide chronosequence of Indian Himalayan region. Land 2022, 11, 1819. [Google Scholar] [CrossRef]
  36. Ray, S.K.; Bhattacharyya, T.; Reddy, K.R.; Pal, D.K.; Chandran, P.; Tiwary, P.; Mandal, D.K.; Mandal, C.; Prasad, J.; Sarkar, D.; et al. Soil and land quality indicators of the Indo-Gangetic Plains of India. Curr. Sci. 2014, 107, 1470–1486. [Google Scholar]
  37. Moncada, M.P.; Gabriels, D.; Cornelis, W.M. Data-driven analysis of soil quality indicators using limited data. Geoderma 2014, 235–236, 271–278. [Google Scholar] [CrossRef]
  38. Vasu, D.; Singh, S.K.; Ray, S.K.; Duraisami, V.P.; Tiwary, P.; Chandran, P.; Nimkar, A.M.; Anantwar, S.G. Soil quality index (SQI) as a tool to evaluate crop productivity in semi-arid Deccan plateau, India. Geoderma 2016, 282, 70–79. [Google Scholar] [CrossRef]
  39. Soil Survey Staff. Keys to Soil Taxonomy by Soil Survey Staff; USDA Natural Resources Conservation Service: Ann Arbor, MI, USA, 2022; p. 410. [Google Scholar]
  40. Piper, C.S. Soil and Plant Analysis: A Laboratory Manual of Methods for the Examination of Soils and the Determination of the Inorganic Constituents of Plants; Reprint for Asia; Scientific Publishers: Jodhpur, India, 2010. [Google Scholar]
  41. Singh, R.A. Soil Physical Analysis; Kalyani Publishers: New Delhi, India, 1980; pp. 64–87. [Google Scholar]
  42. Black, G.R.; Hartge, K.H. Bulk density and Particle density. In Methods of Soil Analysis, Part-1; Monograph, A.K., Ed.; No.9, Agronomy Series; American Society of Agronomy Inc.: Madison, WI, USA, 1986; pp. 363–382. [Google Scholar]
  43. Sankaram, A. Keen Raczkowaski Box Measurements. In A Laboratory Manual for Agricultural Chemistry; Asia Publishing House: Bombay, India, 1960; p. 142. [Google Scholar]
  44. Black, C.A. Methods of Soil Analysis Part—I. Physical and Mineralogical Properties; Agronomy Monograph No. 9.; American Society of Agronomy, Inc.: Madison, WI, USA, 1965; p. 1188. [Google Scholar]
  45. Jackson, M.L. Soil Chemical Analysis, 1st ed.; Prentice Hall Ltd.: New Delhi, India, 1973. [Google Scholar]
  46. Walkley, A.; Black, I.A. An examination of the Degtjareff method for determining soil organic matter, and a proposed modification of the chromic acid titration method. Soil Sci. 1934, 37, 29–38. [Google Scholar] [CrossRef]
  47. Thomas, G.W. Exchangeable cations. In Methods of Soil Analysis Part 2. Chemical and Microbiological Properties; Agronomy Monograph No. 9.; American Society of Agronomy Inc.: Madison, WI, USA, 1982; pp. 159–164. [Google Scholar]
  48. Bower, C.A.; Reitemeier, R.F.; Fireman, M. Exchangeable cations analysis of saline and alkali soils. Soil Sci. 1952, 73, 251–261. [Google Scholar] [CrossRef]
  49. Sharawat, K.L.; Burford, J.R. Modification of alkaline permanganate method for assessing the availability of soil nitrogen in upland soils. Soil Sci. 1982, 133, 53–57. [Google Scholar] [CrossRef]
  50. Mandal, U.K.; Ramachandran, K.; Sharma, K.L.; Satyam, B.; Venkanna, K.; Udaya Bhanu, M.; Mandal, M.; Masane, R.N.; Narsimlu, B.; Rao, K.V.; et al. Assessing soil quality in a semiarid tropical watershed using a geographic information system. Soil Sci. Soc. Am. J. 2011, 75, 1144–1160. [Google Scholar] [CrossRef]
  51. Navas, M.; Benito, M.; Rodriguez, I.; Masaguer, A. Effect of five forage legume covers on soil quality at the Eastern plains of Venezuela. Appl. Soil Ecol. 2011, 49, 242–249. [Google Scholar] [CrossRef]
  52. Li, Q.; Xu, M.X.; Liu, G.B.; Zhao, Y.G.; Tuo, D.F. Cumulative effects of a 17-year chemical fertilization on the soil quality of cropping system in the Loess Hilly Region, China. J. Plant Nutr. Soil Sci. 2013, 176, 249–259. [Google Scholar] [CrossRef]
  53. Masto, R.E.; Chhonkar, P.K.; Singh, D.; Patra, A.K. Alternative soil quality indices for evaluating the effect of intensive cropping, fertilisation and manuring for 31 years in the semi-arid soils of India. Environ. Monit. Assess. 2008, 136, 419–435. [Google Scholar] [CrossRef] [PubMed]
  54. Dunteman, G.H. Principal Components Analysis; Sage: Newcastle upon Tyne, UK, 1989. [Google Scholar]
  55. Jolliffe, I. Principal Component Analysis; Wiley Online Library: Hoboken, NJ, USA, 2005. [Google Scholar]
  56. Wander, M.M.; Bollero, G.A. Soil quality assessment of tillage impacts in Illinois. Soil Sci. Soc. Am. J. 1999, 63, 961–971. [Google Scholar] [CrossRef]
  57. Andrews, S.S.; Carroll, C.R. Designing a soil quality assessment tool for sustainable agroecosystem management. Ecol. Appl. 2001, 11, 1573–1585. [Google Scholar] [CrossRef]
  58. Liebig, M.A.; Varvel, G.; Doran, J.W. A simple performance based index for assessing multiple ecosystem functions. In Agronomy and Horticulture; Faculty Publications, University of Nebraska: Lincoln, NE, USA, 2001; p. 336. [Google Scholar]
  59. Wilding, L.P. Spatial Variability: Its documentation, accommodation, and implication to soil surveys. In Soil Spatial Variability; Nielsen, D.R., Bouma, J., Eds.; PUDOC: Wageningen, The Netherlands, 1985. [Google Scholar]
  60. Govaerts, B.; Sayre, K.D.; Deckers, J. A minimum data set for soil quality assessment of wheat and maize cropping in the highlands of Mexico. Soil Tillage Res. 2006, 87, 163–174. [Google Scholar] [CrossRef]
  61. Mamehpour, N.; Rezapour, S.; Ghaemian, N. Quantitative assessment of soil quality indices for urban croplands in a calcareous. Geoderma 2021, 382, 114781. [Google Scholar] [CrossRef]
  62. Kafei, F.; Rezapour, S.; Dalalian, M.R.; Sabbaghtazeh, E.; Rafieyan, O. Soil quality index as affected by long-time continuous cultivation in a Mediterranean sub-humid region. Rend. Lincei 2023, 34, 563–575. [Google Scholar] [CrossRef]
  63. Hazelton, P.A.; Murphy, B.W. Interpreting Soil Test Results: What Do All the Numbers Mean? CSIRO Publishing: Melbourne, Australia, 2007. [Google Scholar]
  64. Jokela, W.; Posner, J.; Hedtcke, J.; Balser, T.; Read, H. Midwest cropping system effects on soil properties and on a soil quality index. Agron. J. 2011, 103, 1552–1562. [Google Scholar] [CrossRef]
  65. Qi, Y.; Darilek, J.L.; Huang, B.; Zhao, Y.; Sun, W.; Gu, Z. Evaluating soil quality indices in an agricultural region of Jiangsu Province, China. Geoderma 2009, 149, 325–334. [Google Scholar] [CrossRef]
  66. Karaca, S.; Dengiz, O.; Turan, İ.D.; Özkan, B.; Dedeoğlu, M.; Gülser, F.; Ay, A. An assessment of pasture soils quality based on multi indicator weighting approaches in semi-arid ecosystem. Ecol. Indic. 2021, 121, 107001. [Google Scholar] [CrossRef]
  67. Zhang, Y.; Chen, Z.; Wang, X.; Huang, B. Comprehensive assessments of soil fertility and environmental quality in plastic greenhouse production systems. Geoderma 2021, 385, 114899. [Google Scholar]
  68. Tian, K.; Zhang, B.; Zhang, H.; Huang, B.; Darilek, J.L.; Zhao, Y.; Yang, J. Evaluation of soil quality in major grain-producing region of the North China Plain: Integrating minimum data set and established critical limits. Ecol. Indic. 2020, 117, 106613. [Google Scholar] [CrossRef]
  69. Brady, N.C.; Weil, R.R. The Nature and Properties of Soils; Prentice Hall: Upper Saddle River, NJ, USA, 2016. [Google Scholar]
  70. Mesfin, D.; Assefa, E.; Simane, B. Variability of soil quality indicators along with the different landscape positions of Choke Mountain agroecosystem, upper Blue Nile Basin, Ethiopia. Heliyon 2022, 8, e09850. [Google Scholar] [CrossRef]
  71. Ramesh, T.; Bolan, N.S.; Kirkham, M.B.; Wijesekara, H.; Kanchikerimath, M.; Rao, C.S.; Sandeep, S.; Rinklebe, J.; Ok, Y.S.; Choudhury, B.U.; et al. Soil organic carbon dynamics: Impact of land use changes and management practices: A review. In Advances in Agronomy; Sparks, D.L., Ed.; Academic Press: Cambridge, MA, USA, 2019; pp. 1–107. [Google Scholar]
  72. Watson, C.; Atkinson, D.; Gosling, P.; Jackson, L.; Rayns, F. Managing soil fertility in organic farming systems. Soil Use Manag. 2002, 18, 239–247. [Google Scholar] [CrossRef]
  73. Flavel, T.C.; Murphy, D.V. Carbon and nitrogen mineralization rates after application of organic amendments to soil. J. Environ. Qual. 2006, 35, 183–193. [Google Scholar] [CrossRef]
  74. Rengasamy, P. Soil salinization. In Oxford Research Encyclopedia of Environmental Science; Oxford University Press: New York, NY, USA, 2016. [Google Scholar]
  75. Ayub, M.A.; Ahmad, H.R.; Ali, M.; Rizwan, M.; Ali, S.; Zia Ur Rehman, M.; Waris, A.A. Salinity and its tolerance strategies in plants. In Plant Life under Changing Environment: Responses and Management; Academic Press: Cambridge, MA, USA, 2020; pp. 47–76. [Google Scholar] [CrossRef]
  76. Mendes, I.C.; Sousa, D.M.G.; Dantas, O.D.; Lopes, A.A.C.; Junior, F.B.R.; Oliveira, M.I.; Chaer, G.M. Soil quality and grain yield: A win–win combination in clayey tropical Oxisols. Geoderma 2021, 388, 114880. [Google Scholar] [CrossRef]
  77. Roy, D.; Datta, A.; Jat, H.S.; Choudhary, M.; Sharma, P.C.; Singh, P.K.; Jat, M.L. Impact of long term conservation agriculture on soil quality under cereal based systems of North West India. Geoderma 2022, 405, 115391. [Google Scholar] [CrossRef]
  78. Carter, M.R. Concepts of soil quality. In Soil Quality Is in the Hands of the Land Manager; MacEwan, R.J., Carter, M.R., Eds.; Symposium Proceedings; University of Ballarat: Ballarat, Australia, 1996. [Google Scholar]
  79. Norfleet, M.L.; Ditzler, C.A.; Puckett, W.E.; Grossman, R.B.; Shaw, J.N. Soil quality and its relationship to pedology. Soil Sci. 2003, 168, 149–155. [Google Scholar] [CrossRef]
  80. Merril, S.D.; Tanaka, D.L.; Liebig, M.A.; Krupinsky, J.M.; Hansan, J.D.; Anderson, R.L. Sequence effects among crops on alluvial-derived soil compared to those on glacial till-derived soils in the northern Great Plains, USA. Agric. Syst. 2012, 107, 1–12. [Google Scholar] [CrossRef]
  81. Merril, S.D.; Liebig, M.A.; Tanaka, D.L.; Krupinsky, J.M.; Hanson, J.D. Comparison of soil quality and productivity at two sites differing in profile structure and topsoil properties. Agric. Ecosyst. Environ. 2013, 179, 53–61. [Google Scholar] [CrossRef]
  82. Hewitt, A.E. Soil Properties Relevant to Plant Growth: A Guide to Recognising Soil Properties Relevant to Plant Growth and Protection; Manaaki Wenua Press: Lincoln, New Zealand, 2004. [Google Scholar]
  83. Dadhwal, K.S.; Mandal, D.; Shrimali, S.S. Soil quality index for different land use systems in northwestern hilly region of India. J. Indian Soc. Soil Sci. 2011, 59, 169–176. [Google Scholar]
  84. Asseng, S.; Zhu, Y.; Wang, E.; Zhang, W. Crop modelling for climate impact and adaptation. In Crop Physiology Applications for Genetic Improvements and Agronomy, 2nd ed.; Sadras, V.O., Calderini, D., Eds.; Elsevier, B.V.: Amsterdam, The Netherlands, 2015. [Google Scholar]
  85. Dasgupta, P.; Bhattacharjee, D.; Kumari, A. Socio economic analysis of climate change impacts on food grain production in Indian states. Environ. Dev. 2013, 8, 5–21. [Google Scholar] [CrossRef]
  86. Bhattacharyya, T.; Pal, D.K.; Easter, M. Modelled soil organic carbon stocks and changes in the Indo-Gangetic Plains, India from 1980 to 2030. Agric. Ecosyst. Environ. 2007, 122, 84–94. [Google Scholar] [CrossRef]
  87. Doran, J.W. Soil health and global sustainability: Translating science into practice. Agric. Ecosyst. Environ. 2002, 88, 119–127. [Google Scholar] [CrossRef]
  88. Chaudhury, J.; Mandal, U.K.; Sharma, K.L.; Ghosh, H.; Biswapati, M. Assessing soil quality under long-term rice based cropping system. Commun. Soil Sci. Plant Anal. 2005, 36, 1141–1161. [Google Scholar] [CrossRef]
Figure 1. Location map of the Ganjigatti sub-watershed, Karnataka, India.
Figure 1. Location map of the Ganjigatti sub-watershed, Karnataka, India.
Sustainability 15 16576 g001
Figure 2. Representative scree plot showing the variation of eigenvalues with soil components for SQI-1.
Figure 2. Representative scree plot showing the variation of eigenvalues with soil components for SQI-1.
Sustainability 15 16576 g002
Figure 3. Representative scree plot showing the variation of eigenvalues with soil components for SQI-2.
Figure 3. Representative scree plot showing the variation of eigenvalues with soil components for SQI-2.
Sustainability 15 16576 g003
Figure 4. Representative scree plot showing the variation of eigenvalues with soil components for SQI-3.
Figure 4. Representative scree plot showing the variation of eigenvalues with soil components for SQI-3.
Sustainability 15 16576 g004
Figure 5. Spatial distribution of SQI (by SQI-1) of the Ganjigatti sub-watershed.
Figure 5. Spatial distribution of SQI (by SQI-1) of the Ganjigatti sub-watershed.
Sustainability 15 16576 g005
Figure 6. Spatial distribution of SQI (by SQI-2) of the Ganjigatti sub-watershed.
Figure 6. Spatial distribution of SQI (by SQI-2) of the Ganjigatti sub-watershed.
Sustainability 15 16576 g006
Figure 7. Spatial distribution of SQI (by SQI-3) of the Ganjigatti sub-watershed.
Figure 7. Spatial distribution of SQI (by SQI-3) of the Ganjigatti sub-watershed.
Sustainability 15 16576 g007
Table 1. Descriptive statistics of horizon-wise soil properties used for soil quality assessment.
Table 1. Descriptive statistics of horizon-wise soil properties used for soil quality assessment.
MinimumMaximumMeanSDCVSkewness
Coarse sand (%)1.6844.9015.478.2553.341.000
Fine sand (%)3.2535.9021.627.2833.65−0.180
Sand (%)4.9373.2537.0814.1838.230.099
silt (%)4.2032.3617.046.2236.480.359
clay (%)18.3576.4545.8811.8625.85−0.184
B.D (Mg m−3)1.201.651.380.095.930.362
P.D (Mg m−3)2.322.572.440.062.330.169
Porosity (%)30.6749.7943.123.427.93−0.868
MWHC (%)30.4778.8050.589.4418.650.083
VE (%)0.3839.469.577.4377.591.468
FC (%)11.4645.8827.997.4726.67−0.056
PWP (%)6.4329.1016.495.3532.400.340
AWC (%)3.6419.2111.503.0926.800.172
pH5.159.677.081.1215.790.023
EC (dS m−1)0.080.920.250.1558.282.122
OC (%)0.231.340.620.2539.141.029
Ex. Ca [cmol (p+) kg−1]2.1640.3215.519.3660.330.681
Ex. Mg [cmol (p+) kg−1]1.0819.328.435.1861.440.663
Ex. Na [cmol (p+) kg−1]0.138.060.701.15163.044.985
Ex. K [cmol (p+) kg−1]0.091.590.320.2060.823.674
Sum of Ex. Basic cations3.8964.6624.9614.8859.600.724
Ex. Al+3 [cmol (p+) kg−1]0.000.390.060.11172.291.604
Ex. H+ [cmol (p+) kg−1]0.101.780.540.2647.542.344
Exchangeable Acidity0.102.020.600.3355.031.832
CEC [cmol (p+) kg−1]5.2078.7131.5415.2348.290.732
BS (%)36.9398.0475.4515.3420.34−0.585
ESP (%)0.4412.772.091.8789.113.216
CEC/clay ratio0.281.120.670.2434.490.085
CaCO3 (%)0.4020.506.736.0189.291.018
IC (%)0.052.460.810.7389.291.018
TC (%)0.462.751.430.6343.670.774
Available N (kg/ha)70.00217.00137.3736.5426.60.328
Total N (%)0.040.120.080.0326.620.363
C:N ratio3.7512.847.961.9123.940.176
Table 2. Descriptive statistics of weighted average soil properties used for soil quality assessment.
Table 2. Descriptive statistics of weighted average soil properties used for soil quality assessment.
MinimumMaximumMeanSDCVSkewness
Coarse sand (%)4.6543.5816.728.4450.451.302
Fine sand (%)9.4334.3922.866.5928.81−0.228
Sand (%)14.2270.2239.5813.5834.300.150
silt (%)5.9928.0715.415.4035.010.545
clay (%)21.2360.2945.0211.0024.43−0.492
B.D (Mg m−3)1.251.581.390.085.220.591
P.D (Mg m−3)2.342.542.440.062.360.092
Porosity (%)33.5746.9142.842.926.82−1.310
MWHC (%)30.8569.7848.959.2618.92−0.019
VE (%)0.4526.248.956.4371.860.911
FC (%)12.6840.1027.047.1526.44−0.182
PWP (%)7.6423.6515.814.6829.590.171
AWC (%)4.1016.5111.232.9225.97−0.290
pH5.369.106.901.0214.760.376
EC (dS m−1)0.060.490.230.1249.240.796
OC (%)0.331.020.600.1524.790.890
Ex. Ca [cmol (p+) kg−1]2.6633.8213.588.7364.290.966
Ex. Mg [cmol (p+) kg−1]1.2117.867.454.8264.730.981
Ex. Na [cmol (p+) kg−1]0.164.560.600.82136.744.719
Ex. K [cmol (p+) kg−1]0.130.670.300.1342.971.489
Sum of Ex. Basic cations4.2552.1221.9213.8663.241.000
Ex. Al+3 [cmol (p+) kg−1]0.000.330.080.11135.121.117
Ex. H+ [cmol (p+) kg−1]0.211.570.560.2645.102.486
Exchangeable Acidity0.211.770.640.3351.441.756
CEC [cmol (p+) kg−1]6.1859.3328.8914.3549.660.882
BS (%)40.2292.9972.0915.2521.16−0.367
ESP (%)0.537.362.001.3668.112.698
CEC/clay ratio0.291.000.630.2233.840.398
CaCO3 (%)0.4517.105.555.3997.211.452
IC (%)0.052.050.670.6597.211.452
TC (%)0.622.481.280.5845.301.215
Available N (kg/ha)99.75197.53138.3131.7822.980.543
Total N (%)0.060.110.080.0223.110.554
C:N ratio4.9810.947.761.5920.410.063
Table 3. Descriptive statistics of Ap horizon soil properties used for soil quality assessment.
Table 3. Descriptive statistics of Ap horizon soil properties used for soil quality assessment.
MinimumMaximumMeanSDCVSkewness
Coarse sand (%)6.7040.2720.017.6638.270.504
Fine sand (%)11.4532.8823.516.4027.20−0.332
Sand (%)18.1562.6543.5212.3528.37−0.492
silt (%)5.5029.8816.366.2838.400.656
clay (%)21.8756.6540.149.9524.770.066
B.D (Mg m−3)1.211.521.350.085.750.474
P.D (Mg m−3)2.322.542.420.052.060.896
Porosity (%)37.3949.7944.303.167.12−0.351
MWHC (%)31.2959.8445.907.4116.14−0.155
VE (%)0.4710.954.983.6573.230.558
FC (%)13.9031.9323.445.4923.40−0.128
PWP (%)8.1420.4713.283.5726.860.465
AWC (%)4.6214.1210.172.4323.85−0.56
pH5.158.386.470.9314.240.374
EC (dS m−1)0.080.640.230.1253.031.834
OC (%)0.471.340.810.2529.990.616
Ex. Ca [cmol (p+) kg−1]3.2428.4811.386.9360.901.013
Ex. Mg [cmol (p+) kg−1]1.0817.886.344.1264.901.351
Ex. Na [cmol (p+) kg−1]0.130.930.380.2258.080.982
Ex. K [cmol (p+) kg−1]0.140.650.340.1542.940.705
Sum of Ex. Basic cations4.6641.8418.4310.8658.911.032
Ex. Al+3 [cmol (p+) kg−1]0.000.320.090.12124.730.917
Ex. H+ [cmol (p+) kg−1]0.281.100.620.1930.530.671
Exchangeable Acidity0.281.240.710.2737.290.364
CEC [cmol (p+) kg−1]7.3058.6925.3311.8046.581.18
BS (%)36.9395.7370.0815.9122.70−0.308
ESP (%)0.495.441.681.1467.911.892
CEC/clay ratio0.301.120.630.2335.220.378
CaCO3 (%)0.5012.003.713.0682.621.571
IC (%)0.061.440.450.3782.621.571
TC (%)1.420.622.040.3729.100.338
Available N (kg/ha)119.00217.00165.8029.3117.68−0.042
Total N (%)0.070.120.100.0217.99−0.057
C:N ratio4.9612.268.61.9122.21−0.028
Table 4. Principal components of soil quality parameters, eigenvalues, and component matrix variables for soil profiles (SQI-1).
Table 4. Principal components of soil quality parameters, eigenvalues, and component matrix variables for soil profiles (SQI-1).
Principal ComponentsPC1PC2PC3PC4PC5PC6PC7PC8
Eigenvalue16.4012.7482.3471.9921.6101.5591.4001.068
% variance48.2378.0816.9035.8604.7354.5844.1183.140
% cumulative variance48.23756.31863.22169.08173.81678.40082.51885.658
Weightage factor0.56310.09430.08060.06840.05530.05350.04810.0367
Factor loadings (Rotated component matrix)
Coarse sand−0.830−0.296−0.1090.2580.1320.105−0.0510.014
Fine sand−0.669−0.374−0.3120.1780.168−0.2030.1410.195
Sand−0.826−0.364−0.2240.2410.163−0.0430.0430.108
Silt0.5340.1550.0400.218−0.2380.403−0.5040.136
Clay0.7080.3540.246−0.403−0.071−0.1600.213−0.201
BD0.228−0.6220.247−0.4680.270−0.194−0.2900.011
PD0.3490.0780.3250.142−0.0630.038−0.315−0.104
Porosity−0.0570.660−0.1560.504−0.3340.2410.202−0.065
MWHC0.8940.1140.164−0.070−0.059−0.0530.1380.009
VE0.832−0.2010.2440.012−0.071−0.0320.1920.159
FC0.8980.123−0.064−0.1440.045−0.1240.160−0.004
PWP0.8730.088−0.058−0.053−0.0460.0040.089−0.030
AWC0.6620.145−0.054−0.2580.188−0.3070.2330.042
pH0.884−0.173−0.1480.182−0.058−0.105−0.012−0.074
EC0.529−0.3480.3480.3060.1570.0860.1190.212
OC−0.5540.3920.4330.2170.310−0.323−0.2360.044
Ex. Ca0.9240.099−0.1020.0050.2820.093−0.0430.035
Ex. Mg0.9040.149−0.041−0.0210.2340.003−0.024−0.039
Ex. Na0.563−0.2610.5440.3880.0250.0490.328−0.068
Ex. K0.3630.380−0.0570.2000.2200.041−0.147−0.189
Sum of Ex.cations0.9440.099−0.0370.0280.2640.063−0.0120.001
Ex. Al+3−0.6770.0330.369−0.1300.1340.3440.0510.162
Ex. H+−0.6080.1010.242−0.3460.2190.4780.099−0.103
Exchangeable Acidity−0.6890.0890.306−0.3100.2120.4810.093−0.028
CEC0.8920.141−0.031−0.0040.3280.1040.067−0.11
BS0.7490.051−0.1040.077−0.151−0.244−0.2180.262
ESP0.298−0.3530.5830.419−0.171−0.0040.36−0.052
CEC/clay ratio0.705−0.001−0.2850.2690.4190.236−0.0540.053
CaCO30.880−0.152−0.0760.0540.1570.224−0.1030.197
Inorganic carbon0.880−0.152−0.0760.0540.1570.224−0.1030.197
Total carbon0.2280.2480.352−0.284−0.363−0.004−0.1130.632
Available N−0.6550.4920.0770.1380.337−0.1900.1450.270
Total N−0.6600.4870.0830.1440.329−0.2000.1440.266
C:N ratio−0.1760.0540.5690.2420.112−0.293−0.475−0.275
Bold-face factor loadings were considered highly weighted, and the underlined were retained in MDS.
Table 5. Correlation coefficient (Pearson’s) for highly loaded parameters in PC1 (SQI-1).
Table 5. Correlation coefficient (Pearson’s) for highly loaded parameters in PC1 (SQI-1).
MWHCFCPWPpHEx.CaEx.MgSEBCCECCaCO3IC
MWHC1
FC0.84 **1
PWP0.80 **0.94 **1
pH0.75 **0.74 **0.73 **1
Ex. Ca0.78 **0.82 **0.79 **0.77 **1
Ex. Mg0.82 **0.80 **0.76 **0.77 **0.91 **1
SEBC0.82 **0.83 **0.80 **0.80 **0.98 **0.96 **1
CEC0.79 **0.80 **0.75 **0.72 **0.95 **0.93 **0.97 **1
CaCO30.73 **0.74 **0.73 **0.79 **0.87 **0.80 **0.87 **0.78 **1
IC0.73 **0.74 **0.73 **0.79 **0.87 **0.80 **0.87 **0.78 **1.00 **1
** Correlation is significant at the 0.01 level (2-tailed).
Table 6. Correlation coefficient (Pearson’s) for highly loaded parameters in PC 2 (SQI-1).
Table 6. Correlation coefficient (Pearson’s) for highly loaded parameters in PC 2 (SQI-1).
BDPorosity
BD1
Porosity−0.90 **1
** Correlation is significant at the 0.01 level (2-tailed).
Table 7. Correlation coefficient (Pearson’s) for highly loaded parameters in PC 3 (SQI-1).
Table 7. Correlation coefficient (Pearson’s) for highly loaded parameters in PC 3 (SQI-1).
Ex. NaESPCN
Ex. Na1
ESP0.87 **1
CN0.110.151
** Correlation is significant at the 0.01 level (2-tailed).
Table 8. Correlation coefficient (Pearson’s) for selected MDS parameters from PCA results (SQI-1).
Table 8. Correlation coefficient (Pearson’s) for selected MDS parameters from PCA results (SQI-1).
CECPorosityESPCNBDCEC/ClayEx. AciditySiltTC
CEC1
Porosity−0.051
ESP0.15−0.26 *1
CN−0.170.23 *0.031
BD0.11−0.12−0.040.22 *1
CEC/clay0.82 **0.150.00−0.120.081
Ex. Acidity−0.00−0.79 **0.27 *−0.33 **−0.46 **−0.181
Silt0.36 **0.31 **0.030.130.050.48 **−0.28 **1
TC0.060.020.03−0.030.06−0.09−0.060.28 *1
** Correlation is significant at the 0.01 level (2-tailed). * Correlation is significant at the 0.05 level (2-tailed).
Table 9. Principal components of soil quality parameters, eigenvalues, and component matrix variables for soil profiles (SQI-2).
Table 9. Principal components of soil quality parameters, eigenvalues, and component matrix variables for soil profiles (SQI-2).
Principal ComponentsPC 1PC 2PC 3PC 4PC 5PC 6PC 7PC 8
Eigenvalue17.3042.8972.3432.1351.8931.5611.3881.152
% variance50.8958.5216.8916.2805.5684.5914.0833.389
% cumulative variance50.89559.41666.30772.58678.15482.74586.82890.218
Weightage factor0.56410.09440.07640.06960.06170.05090.04530.0376
Factor loadings (Rotated component matrix)
Coarse sand−0.849−0.224−0.203−0.114−0.0460.214−0.0560.068
Fine sand−0.694−0.067−0.449−0.125−0.3280.0400.1320.034
Sand−0.861−0.173−0.338−0.131−0.1810.1540.0260.059
Silt0.640−0.3380.047−0.1060.400−0.129−0.3110.207
Clay0.7420.3810.3930.2140.024−0.1260.123−0.176
BD0.208−0.5120.0570.663−0.3900.1060.0550.241
PD0.3730.0210.372−0.009−0.0200.180−0.4960.452
Porosity−0.0340.5290.136−0.6680.388−0.062−0.2680.056
MWHC0.8650.1680.228−0.0300.0470.0020.221−0.088
VE0.851−0.0710.160−0.0890.0220.0650.2720.138
FC0.8090.269−0.0140.054−0.080−0.1060.070−0.118
PWP0.8030.1560.005−0.0680.005−0.1140.080−0.069
AWC0.7800.408−0.0430.241−0.204−0.0770.043−0.179
pH0.901−0.133−0.066−0.165−0.2120.0240.027−0.072
EC0.695−0.538−0.2030.0060.0560.1950.0040.058
OC−0.4660.2750.5710.067−0.4080.325−0.1330.044
Ex. Ca0.9270.141−0.1630.1530.0330.188−0.0820.064
Ex. Mg0.9280.146−0.0440.134−0.0200.1950.019−0.051
Ex. Na0.570−0.2820.302−0.3820.0760.4630.271−0.067
Ex. K0.4720.215−0.297−0.246−0.0440.146−0.1530.123
Sum of Exchangeable cations0.9610.124−0.1020.1180.0180.215−0.0300.020
Ex. Al+3−0.750−0.0440.0770.1060.4380.2260.1470.099
Ex. H+−0.572−0.0160.2000.4980.4520.1740.014−0.167
Exchangeable Acidity−0.677−0.0260.1790.4180.4860.2050.057−0.098
CEC0.8980.221−0.0920.1500.0790.272−0.127−0.067
BS0.781−0.0610.034−0.061−0.303−0.2050.1630.226
ESP0.254−0.4550.378−0.5270.0490.2800.432−0.058
CEC/clay ratio0.7420.120−0.3750.0630.0730.363−0.231−0.026
CaCO30.877−0.142−0.2480.0740.1660.1430.0330.048
Inorganic carbon0.878−0.142−0.2480.0740.1650.1420.0340.047
Total carbon0.3490.0000.2610.1540.200−0.4120.3110.555
Available N−0.6090.598−0.0140.016−0.1310.3260.2280.239
Total N−0.6170.6000.0040.012−0.1460.3160.2220.237
CN ratio0.076−0.3370.661−0.036−0.3760.096−0.408−0.238
Bold-face factor loadings were considered highly weighted, and the underlined were retained in MDS.
Table 10. Correlation coefficient (Pearson’s) for highly loaded parameters in PC1 (SQI-2).
Table 10. Correlation coefficient (Pearson’s) for highly loaded parameters in PC1 (SQI-2).
SEBCpHCECCaCO3IC
SEBC1
pH0.81 **1
CEC0.96 **0.73 **1
CaCO30.89 **0.80 **0.79 **1
IC0.89 **0.80 **0.79 **1.00 **1
** Correlation is significant at the 0.01 level (2-tailed).
Table 11. Correlation coefficient (Pearson’s) for selected MDS parameters from PCA results (SQI-2).
Table 11. Correlation coefficient (Pearson’s) for selected MDS parameters from PCA results (SQI-2).
CECTotal NCNBDEx. AcidityEx. NaPDTC
CEC1
Total N−0.43 *1
CN−0.08−0.121
BD0.18−0.330.121
Ex. Acidity−0.49 **0.39 *−0.07−0.051
Ex. Na0.47 *−0.350.220.08−0.341
PD0.38−0.140.270.23−0.210.231
TC0.16−0.19−0.110.20−0.160.130.231
** Correlation is significant at the 0.01 level (2-tailed). * Correlation is significant at the 0.05 level (2-tailed).
Table 12. Principal components of soil quality parameters, eigenvalues, and component matrix variables for soil profiles Ap horizons (SQI-3).
Table 12. Principal components of soil quality parameters, eigenvalues, and component matrix variables for soil profiles Ap horizons (SQI-3).
Principal ComponentsPC1PC2PC3PC4PC5PC6PC7PC8
Eigenvalue13.8004.1553.0402.2822.1601.8281.2311.174
% variance40.58912.2208.9406.7116.3545.3773.6223.452
% cumulative variance40.58952.80961.74968.4674.81480.19183.81387.265
Weightage factor0.46510.14000.10240.07690.07280.06160.04150.0396
Factor loadings (Rotated component matrix)
Coarse sand−0.814−0.397−0.0250.0920.0120.1690.135−0.071
Fine sand−0.613−0.157−0.2680.1990.4500.242−0.3190.169
Sand−0.823−0.328−0.1540.1600.2410.230−0.0810.043
Silt0.565−0.2760.203−0.278−0.440−0.0840.0390.115
Clay0.6650.5820.063−0.0230.042−0.2330.076−0.127
BD−0.233−0.1940.524−0.5600.3830.212−0.177−0.157
PD−0.013−0.0030.259−0.186−0.4390.6450.186−0.010
Porosity0.2260.198−0.4330.505−0.5610.0330.2520.155
MWHC0.8480.2800.1200.0500.107−0.197−0.033−0.161
VE0.752−0.0300.0770.152−0.080−0.130−0.420−0.126
FC0.8420.239−0.0850.0400.3530.042−0.1470.140
PWP0.8360.092−0.0850.0090.2700.070−0.1310.128
AWC0.6760.405−0.0670.0780.401−0.008−0.1400.128
pH0.743−0.336−0.251−0.183−0.0910.1740.114−0.059
EC0.129−0.001−0.247−0.199−0.2640.293−0.2900.717
OC−0.1810.8340.119−0.2660.1220.2060.1850.150
Ex. Ca0.921−0.0930.1660.1940.0760.1310.1110.026
Ex. Mg0.9060.0460.2110.0700.044−0.0310.146−0.065
Ex. Na0.3790.066−0.7600.0200.177−0.2470.2100.075
Ex. K0.5510.228−0.042−0.0060.1520.5140.178−0.110
Sum of Ex. cations0.946−0.0370.1700.1510.0710.0740.133−0.008
Ex. Al+3−0.6880.2560.3510.096−0.138−0.2080.013−0.107
Ex. H+−0.5170.1830.4800.3800.126−0.1600.1700.363
Exchangeable Acidity−0.6630.2400.4930.3120.031−0.2030.1270.214
CEC0.851−0.0560.2260.3680.1550.0420.1420.011
BS0.7470.121−0.169−0.402−0.1630.101−0.004−0.079
ESP−0.1540.137−0.789−0.2820.097−0.2680.148−0.009
CEC/clay ratio0.660−0.4150.1410.4350.1110.1800.0680.139
CaCO30.844−0.2950.181−0.039−0.058−0.010−0.0770.085
Inorganic carbon0.844−0.2950.181−0.039−0.058−0.010−0.0770.085
Total carbon0.1440.5270.203−0.118−0.443−0.290−0.4420.106
Available N−0.2020.757−0.1230.261−0.1160.375−0.156−0.191
Total N−0.2060.758−0.1280.265−0.1170.375−0.149−0.189
CN ratio−0.0780.4730.236−0.5580.277−0.0390.3570.300
Bold-face factor loadings were considered highly weighted, and the underlined were retained in MDS.
Table 13. Correlation coefficient (Pearson’s) for highly loaded parameters in PC1 (SQI-3).
Table 13. Correlation coefficient (Pearson’s) for highly loaded parameters in PC1 (SQI-3).
SEBCEx. CaEx. MgCEC
SEBC1
Ex. Ca0.98 **1
Ex. Mg0.95 **0.88 **1
CEC0.94 **0.93 **0.89 **1
** Correlation is significant at the 0.01 level (2-tailed).
Table 14. Correlation coefficient (Pearson’s) for highly loaded parameters in PC2 (SQI-3).
Table 14. Correlation coefficient (Pearson’s) for highly loaded parameters in PC2 (SQI-3).
OCTNAN
OC1
TN0.63 **1
AN0.69 **0.72 **1
** Correlation is significant at the 0.01 level (2-tailed).
Table 15. Correlation coefficient (Pearson’s) for highly loaded parameter in PC4 (SQI-3).
Table 15. Correlation coefficient (Pearson’s) for highly loaded parameter in PC4 (SQI-3).
BDPorosityCN
BD1
Porosity−0.93 **1
CN0.28−0.271
** Correlation is significant at the 0.01 level (2-tailed).
Table 16. Correlation coefficient (Pearson’s) for selected MDS parameters from PCA results (SQI-3).
Table 16. Correlation coefficient (Pearson’s) for selected MDS parameters from PCA results (SQI-3).
CECOCESPBDCNPorosityPDTCEC
CEC1
OC−0.241
ESP−0.39 *0.101
BD−0.200.10−0.211
CN−0.170.78 **0.060.281
Porosity0.21−0.050.11−0.93 **−0.271
PD0.000.14−0.250.250.040.131
TC−0.050.28−0.09−0.140.140.150.041
EC−0.040.030.08−0.110.000.190.200.161
** Correlation is significant at the 0.01 level (2-tailed) * Correlation is significant at the 0.05 level (2-tailed).
Table 17. Soil quality index (SQI) values for 27 soil profiles using different methods.
Table 17. Soil quality index (SQI) values for 27 soil profiles using different methods.
Profile NumberDepth of Soil Profile (cm)SQI Method
Horizon WiseWeighted AverageAp Horizon
Pedon-1280.460.510.55
Pedon-2390.450.550.58
Pedon-3490.520.550.60
Pedon-4790.730.820.89
Pedon-5650.490.540.61
Pedon-61800.700.780.82
Pedon-7670.370.430.63
Pedon-8400.440.530.67
Pedon-9800.600.710.79
Pedon-10350.440.520.57
Pedon-112000.740.850.77
Pedon-12350.300.400.52
Pedon-13260.370.450.50
Pedon-14200.410.480.62
Pedon-15900.420.490.62
Pedon-16460.380.460.58
Pedon-171200.370.440.53
Pedon-181700.410.470.59
Pedon-191300.650.750.71
Pedon-20200.390.480.59
Pedon-21350.430.510.53
Pedon-22300.340.410.48
Pedon-231300.640.730.73
Pedon-24900.450.490.53
Pedon-25550.430.510.63
Pedon-26300.260.320.39
Pedon-271800.570.610.63
Table 18. Pearson’s correlation coefficients (r) of soil quality index (SQI) versus crop yield.
Table 18. Pearson’s correlation coefficients (r) of soil quality index (SQI) versus crop yield.
SoybeanMaizeSorghumGreengram
SQI-10.74 **0.69 **0.69 **0.73 **
SQI-20.63 **0.61 **0.58 **0.67 **
SQI-30.58 **0.58 **0.54 **0.53 **
** Correlation is significant at the 0.01 level (2-tailed).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yadav, M.B.N.; Patil, P.L.; Hebbara, M. Comparative Studies on Soil Quality Index Estimation of a Hilly-Zone Sub-Watershed in Karnataka. Sustainability 2023, 15, 16576. https://doi.org/10.3390/su152416576

AMA Style

Yadav MBN, Patil PL, Hebbara M. Comparative Studies on Soil Quality Index Estimation of a Hilly-Zone Sub-Watershed in Karnataka. Sustainability. 2023; 15(24):16576. https://doi.org/10.3390/su152416576

Chicago/Turabian Style

Yadav, M. Bhargava Narasimha, P. L. Patil, and M. Hebbara. 2023. "Comparative Studies on Soil Quality Index Estimation of a Hilly-Zone Sub-Watershed in Karnataka" Sustainability 15, no. 24: 16576. https://doi.org/10.3390/su152416576

APA Style

Yadav, M. B. N., Patil, P. L., & Hebbara, M. (2023). Comparative Studies on Soil Quality Index Estimation of a Hilly-Zone Sub-Watershed in Karnataka. Sustainability, 15(24), 16576. https://doi.org/10.3390/su152416576

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop