Soil Quality and Its Potential Indicators under Di ﬀ erent Land Use Systems in the Shivaliks of Indian Punjab

: The present study assessed the overall state of the soil based on the most sensitive soil attributes under di ﬀ erent land uses—i.e., rainfed agriculture, mixed forest, a ﬀ orestation and non-arable lands—in the lower Shivaliks of Indian Punjab. The soil parameters—i.e., erosion ratio, bulk density and water retention characteristics—and fertility parameters were integrated under di ﬀ erent land uses to identify potential soil quality indicators.The overall state of the soil, based on a weighted average of primary soil functions under di ﬀ erent land uses through fuzzy modeling, was deemed best for agricultural land use (0.515), followed by forests (0.465) and non-arable lands (0.456), and deemed worst under a ﬀ orestation (0.428). Among the di ﬀ erent land use systems, principal component analysis (PCA) clearly separated the agriculture and forest samples from a ﬀ orestation and non-arable lands samples. The contribution of potential indicators such as erosion ratio (ER), phosphorus (P) and potassium (K) toward the soil quality index (SQI) was substantial. The order of contribution of the selected indicators to the SQI was 53.5%, 34.3% and 19.9% for ER, P and K, respectively. These indicators are most inﬂuential for studying real time soil health and ecological processes in the future, under various land use systems in degraded agroecosystems like the Shivaliks.


Introduction
The Shivaliks region, in the lower Himalayas of India, is characterized by its susceptibility to land degradation because of many topographic and land cover characteristics. The young, unconsolidated low Shivaliks hills, which were once dotted with perennial springs and gentle streams, were converted into a tract of ever widening and deepening torrent beds, which bring flash floods from the barren hills [1]. This region receives 800 to 1000 mm of rainfall, which is almost 1.5 times that of the alluvial plain of Punjab, but about 80 percent of the rainfall occurs in three monsoon months, and the major part ends in runoff. In the remaining months, the availability of water for drinking, livestock, domestic use and field crops is extremely limited. The overexploitation and mismanagement of soil resources through deforestation, overgrazing and land clearing for agriculture has resulted in the ecological degradation of the Shivalik hills and contributed immensely toward excessive silt load in rivers of the foothill, reduced recharge of ground water and an overall decline in land productivity, which has also affected the socio-economic conditions of the inhabitants [2,3].This region faces a serious

Study Site Description
The study area lies between 30 • 40 to 32 • 05 N latitude and 75 • 25 to 76 • 58 E longitude, and its elevation varies from 250 to 350 m above mean sea level ( Figure 1). The topography of the area is undulating to rolling, with the angle of the slopes ranging from 5 to 15 degrees. The total area of this belt is 497,520 ha, constituting about 9.88 percent of the geographical area covering the five districts of the Punjab state. The soils occur on the foothill belt of the Shivalik hills known as the kandi area (piedmont area) of Punjab. The benchmark soils of this region have been characterized as deep, imperfectly drained, coarse-loamy and non-calcareous [16]. The soils were classified following the criteria of soil taxonomy [17]. The color of the typifying pedon of the soils is dark brown on the surface horizon and dark yellowish brown on the subsurface horizons. The soils appear to have developed on stratified parent material, as the texture changes with depth ( Table 1). The texture varies from sandy loam through silty loam to loamy sand or sand, from surface to depth. There is no structural development associated with pedogenesis. The organic carbon content also changes irregularly with depth, suggesting the fluventic or stratified nature of soils.
Diagnostics 2020, 10  The climate of the study area is semi-arid, with a mean annual air temperature of 23.7 • C and a mean annual rainfall around 875 mm. More than half of the rainfall just ends in runoff because of the inadequate vegetation cover on the hilly forest and arable sloping land. The soils are mostly erodible and belong to an ustic moisture regime and to a hyperthermic temperature regime. The temperature regime is hyperthermic. These soils are a member of the coarse-loamy, mixed, hyperthermic family of Typic Ustifluvents. Maize-wheat is the main cropping system of the region. Some areas are under green fodder, pulses and fruit plants.

Sample Collection and Analysis
The soil samples (0-0.15 m) were collected from four land use systems, namely rainfed agriculture, mixed forests, afforestation and non-arable lands in the Lower Shivaliks of Punjab. The soil samples were collected from different micro-watersheds depending upon the percentage of area under different land use. A total of 60 composite samples were collected, i.e., 40 from the Agriculture ecosystem (mostly under the Maize-wheat cropping system), 10 from Forests (Mixed forests with moderate to high slopes), 5 from each area prone to afforestation (forest cover reduced to almost zero) and non-arable lands (lands not under use, mostly barren). A higher grid size was kept for Agriculture as compared to other land uses. Composite samples were prepared by taking surface soil samples from random spots (15)(16)(17)(18)(19)(20) using a push probe auger from sites with the different land uses in different parts of the Shivaliks. These composite samples were air-dried at room temperature and ground to pass through a 2 mm sieve after removing all visible residues. They were subsequently stored in a moisture-free environment for the determination of their physico-chemical properties. Dichromate oxidizable SOC was estimated by the rapid titration method, using a diphenyl amine indicator [18]. The pH and EC of the soil suspension-i.e., 1:2 (soil:water)-were measured by a pH meter and a conductivity meter, respectively. The soil samples were analyzed for Olsen P and ammonium acetate (NH 4 OAc)-extractable K by the methods described by Olsen et al. 1954 [19] and Knudsen et al. 1982 [20], respectively. The particle size distribution-i.e., sand, silt and clay contents-was determined by the Pipette method and the bulk density with core method [21]. The soil texture of the study area varied from loamy sand to sandy loam (Agriculture & Afforestation land use), loamy sand (Forest) and sandy loam to loam (Non-arable). The soil moisture characteristic, i.e., field capacity, and the permanent wilting point, i.e., water held at 0.33 and 15 bar of suction, respectively, were determined by the use of a Pressure Plate membrane apparatus [22]. The erosion ratio (ER) was computed by using the relationships suggested by Middleton 1930 [23]. Middleton found that soilshaving anerosion ratio greater than 10 were most frequently susceptible toerosion. It was determined as the ratio of dispersion to the percentage of clay in the moisture equivalent. Dispersion ratio was determined as the ratio of the percentage of silt plus clay in undispersed soil to the percentage of silt plus clay in dispersed soil.

Overall State of Soils
Based on the primary functions of the soils, the relevant indicators were selected for indexing the overall state ofthesoils under different land uses ( Table 2). The observed values of soil attributes (potential soil quality indicators) were converted into a unit less score of 0-1 through fuzzy modeling, using different scoring curves [24]. The detail of this methodology is given elsewhere [2,25].
The notion that "lower is better" was used for the erosion ratio and bulk density, as the low values of these factors indicates a better resistance to erosion. On other hand, inthecase of soil moisture constants and plant growth parameters, that notion that "higher is better" was used. The ranges of indicators and converted scores are presented in Table 2. The relevant indicators were apportioned weights dependent upon their comparative importance ( Table 3). The uppermost weight (0.30) was allocated to physical degradation/erodibility, i.e., erosion ratio. Soil compaction and water retention were assigned a weight of 0.20 (bulk density) and 0.08 (field capacity and permanent wilting point), respectively. A weight of 0.15 was assigned to the organic carbon content, and 0.05 to pH, N, P and K. After this, the individual rating (q) and assigned weights (W) of different soil functions were multiplied to obtain the weighted scores (0 to 1 scale), which were then summed up to get the aggregated score (Q) representing the soil's state or condition (Equation (1)). Q = q er W er + q bd W bd + q wr W wr + q oc W oc + q pg W pg (1) where er represents the erosion ratio; bd the bulk density; wr the water retention; oc the organic carbon; and pg the plant growth functions (pH, N, P and K).

Data Analysis
A principal component analysis (PCA) was employed to define the most significant and important soil properties influenced by the varying state of soils in different land uses. PCA is a mathematical technique that converts a number of correlated variables into a smaller number of uncorrelated or independent variables, based on the eigenvector decomposition of the covariance [26,27]. Each PC explained a certain amount (%) of the variation in the total dataset. The PCs with higher eigen values and variables which had high factor loading were considered as the best representatives of the system's attributes. The correlations among the indicators were assessed by the determination ofthePearson correlation coefficients (r) and probabilities. In all the analyses, significance was accepted at a level of probability (p) of <0.05. The data were statistically analyzed with the analysis of variance (ANOVA) technique using locally developed software [28]. The means for treatment effects were separated based on Fisher's protected least significant difference (LSD) at p < 0.05.

Characteristics of the Soils
The ranges of the selected indicators/soil attributes-i.e., erosion ratio, bulk density, water retention characteristics and fertility parameters-under different land uses are presented in Table 4. The erosion ratio has generally been recognized as the most suitable index to characterize soils according to their erodibility [29]. It is useful for predicting soil's susceptibility to erosion. Particle size distribution and organic matter content are the two most important indicators of erodibility. The erosion ratios in the studied soils ranged from 24.6 to 86.7 in rainfed agriculture lands, 24.5 to 43.0 in forests soils, 51.0 to 51.4 in the afforestation area and 33.8 to 62.6 in non-arable lands. Forest soils have a low value of erosion ratio compared to other land uses [30]. The erosion indices among different land uses [31] in the lower Shivaliks were reported in the following sequence: barren (0.97), cultivated (0.84), grassland (0.74) and forest (0.63). Considering the upper limit oftheerosion ratio as 10 for non-erosive soils [23], all the soils in the Shivaliks are susceptible to erosion. However, the mean value oftheerosion ratio was low in the soils under forests. This may be due to the litter fall and higher root biomass incorporation [30]. Singh et al. 2006 [32] also reported the low value of the erosion ratio for soils under forests and grasslands in southeastern Rajasthan. Due to an increasing demand for timber, pasture and food and to residential dwindling, forests are being degraded or converted to cropland, which results in lower soil quality and, consequently, decreased productivity [33]. and other land uses. This is probably due tothe loss of the soils'organic matter, combined with greater sand content and poorer aggregation; moreover, trampling by grazing livestock on the surface layer resulted in higher bulk density in the natural forests. Higher values of bulk density were observed in torrent-affected areas of the Shivalik-Himalayan region in India, because of a very poor organic carbon content and poor soil aggregation [2]. The soils were foundto bepoor in organic carbon content under all land uses. However, the content was higher in the soils under non-arable (0.30%) and agricultural lands (0.29%), followed by forest soils (0.20%). Variations in organic carbon in soils under various land uses were observed due to varying leaf litter and their rate of decomposition [34]. Forests had a higher basal soil respiration than the other land use systems due to a relatively dense structure of plants and the continuous deposition of organic matter through leaf litter and fine roots [35]. In eroded forest soils of the outer Himalayas, a negative correlation between bulk density with organic carbon and clay content has been reported [36,37]. The higher organic carbon content under cultivated conditions may be ascribed to the application of farm yard manure and the decomposition of crop residues.
The soils under non-arable and cultivated lands exhibited higher water retention than the areas under afforestation and mixed forests. In general, soil water constants-i.e., the field capacity and permanent wilting point of the soils-ranged from 7.0 to 20.8 and from 2.0 to 7.0%, respectively. The arable lands were mostly under the maize-wheat cropping system. The application of farm yard manure and mulching with locally available material helped increase crop productivity, available moisture as well as soil health [38]. For non-arable lands-which are mostly the government's lands, with locally grown wild grasses-water retention was also high due to the fine texture, and slopes were low as compared to forests and other areas under afforestation. In a study on a micro-watershed in the Shivaliks [15], the highest value of moisture retention at field capacity was observed in grassland areas (21.2%). This value was significantly at par with that of cultivated areas (17.7%), followed by horticulture (12.4%), bare areas (12.6%) and forests (11.4%). The soils under study were neutral to alkaline (pH 6.5-9.4). The mean pH was near neutral in forest soils and alkaline in afforestation and non-arable lands, indicating a better nutrient status. The available nitrogen, phosphorus and potassium contents were low in all land uses. However, N, P and K were higher in non-arable lands (185.0 kg ha −1 ), agriculture (10.1 kg ha −1 ) and forests (123.0 kg ha −1 ), respectively. The afforested areas maintained the lower values of N (103.5 kg ha −1 ), P (3.2 kg ha −1 ) and K (108.8 kg ha −1 ). This might be due to the coarser soil texture in the areas prone to afforestation.

Overall State of Soils under Different Land Uses
The aggregated score (Q) representing the overall state or condition of the soils under different land uses varied from 0.359 to 0.900 ( Figure 2, Table 4). Soils under agriculture use had the highest mean aggregated value (Q) of 0.515, followed by forests (0.465), non-arable lands (0.456) and afforestation (0.428). This result aligns with those reported by Bhattacharya et al. 2008 [2], where in an aggregate score was reported to range from 0.41 to 0.77 for different physiographic units of Punjab. The aggregated low scores of the afforestation and non-arable lands indicate variability in the soil's resilience. Regarding the soil state, in terms of the weighted score of different potential indicators, the erosion ratio ( Figure 2) maintained higher values in forests (0.176), followed by agriculture (0.139), afforestation (0.135) and non-arable lands (0.105). The erodibility was found to vary directly with sand and inversely with clay. Forest Soils were coarser in texture in forests and afforestation ecosystems, compared to other ecosystems. Many of the forest sites are located on steep slopes, where the finer soil particles have been selectively removed by erosion, thereby increasing the proportion of coarse particles in the soils [39]. The highest value of bulk density, 0.145, was found in non-arable areas, followed by agriculture (0.141), afforestation (0.130) and forests (0.097). The converted scores of field capacity and permanent wilting point were highest in non-arable areas (0.051 and 0.004) and lowest in forests (0.033 and 0.003). proportion of coarse particles in the soils [39]. The highest value of bulk density, 0.145, was found in non-arable areas, followed by agriculture (0.141), afforestation (0.130) and forests (0.097). The converted scores of field capacity and permanent wilting point were highest in non-arable areas (0.051 and 0.004) and lowest in forests (0.033 and 0.003). Soil organic matter has been reported as the most powerful indicator for assessing the soil quality index in different regions of the world under varied land uses and management interventions [33,40]. The highest score of organic carbon was observed in non-arable and agricultural lands (0.064 and 0.063) and the lowest under afforestation (0.038). It was reported that deforestation and subsequent cultivation decreased organic matter by 48.8% [41]. The converted values were higher in agricultural lands for pH, N, P and K. The aggregated score indicates that the soil's limitations under all the land uses were moderate to low, and that they can be managed with suitable conservation measures to maintain optimum productivity. This conversion of natural forests to other uses, such as cultivation, has created a severe degradation of soil quality, which may lead to a permanent degradation of land productivity [42]. Areaswith these degraded soils require urgent protection and should not be cultivated. The soils under afforestation are more sensitive to soil erosion and require suitable and immediate soil conservation measures. Preferably, these areas should be brought under permanent vegetation or planted to cover crops.

Principal Component Analysis
The amount of variability explained by PC1 was 63.3%, with an eigenvalue of 5.69, which includes the ER with the highest positive factor loading value, i.e., 0.979, and K, i.e., 0.914 (Figure 3, Table 5). The second component (PC2) explained about 28.1% of variance, with an eigenvalue of 2.53 and the highest positive loading value for P with positive factor loading, i.e., 0.828, and pH, i.e., 0.773. The weight of each PC on the basis of the percentage of variance to total variance ranged from 0.31 to 0.69. The weighted factor for the minimum data set (MDS) had the following trend: PC1 (0.9) > PC2 (0.31). Soil organic matter has been reported as the most powerful indicator for assessing the soil quality index in different regions of the world under varied land uses and management interventions [33,40]. The highest score of organic carbon was observed in non-arable and agricultural lands (0.064 and 0.063) and the lowest under afforestation (0.038). It was reported that deforestation and subsequent cultivation decreased organic matter by 48.8% [41]. The converted values were higher in agricultural lands for pH, N, P and K. The aggregated score indicates that the soil's limitations under all the land uses were moderate to low, and that they can be managed with suitable conservation measures to maintain optimum productivity. This conversion of natural forests to other uses, such as cultivation, has created a severe degradation of soil quality, which may lead to a permanent degradation of land productivity [42]. Areas with these degraded soils require urgent protection and should not be cultivated. The soils under afforestation are more sensitive to soil erosion and require suitable and immediate soil conservation measures. Preferably, these areas should be brought under permanent vegetation or planted to cover crops.

Principal Component Analysis
The amount of variability explained by PC1 was 63.3%, with an eigenvalue of 5.69, which includes the ER with the highest positive factor loading value, i.e., 0.979, and K, i.e., 0.914 (Figure 3, Table 5). The second component (PC2) explained about 28.1% of variance, with an eigenvalue of 2.53 and the highest positive loading value for P with positive factor loading, i.e., 0.828, and pH, i.e., 0.773. The weight of each PC on the basis of the percentage of variance to total variance ranged from 0.31 to 0.69. The weighted factor for the minimum data set (MDS) had the following trend: PC1 (0.9) > PC2 (0.31). Sustainability 2020, 12, x FOR PEER REVIEW 9 of 13  In rainfed land use systems, the PCA identified the carbon management index (28.9%), metabolic potential (9.8%) and total organic carbon (5.4%) as the three most dominant and reliable  In rainfed land use systems, the PCA identified the carbon management index (28.9%), metabolic potential (9.8%) and total organic carbon (5.4%) as the three most dominant and reliable indicators of soil health [43]. Under different land uses in submontane Punjab, organic carbon contributed the mostto soil quality (28.5%), followed by the available K (19.4%), electrical conductivity (18.3%), K-factor of universal soil loss equation (USLE) (14.9%), plant available water (10.5%), and clay (8.3%) [15]. In the present study, the contribution of potential indicators like ER, P and K toward the soil quality index (SQI) was substantial. ER plays a dominant role in identifying, based on the soil's characteristics, the areas that may require more careful land management techniques. Particle size distribution and organic matter content are the two most important indicators of erodibility. Nair (1984) [44] reported that agro-forestry, agro-horticulture and grassland systems have the potential to reduce both the runoff and erosion, and to maintain better soil organic matter, which in turn will help improve the fertility status of the soil. The available K is important for osmotic regulations. Potassium provides much of the osmotic pull that draws water into plant roots. Plants that are K-deficient are less able to withstand water deficits, mostly because of their inability to make full use of the available water. The malfunctioning of stomata due to a deficiency of this nutrient has been related to lower rates of photosynthesis and a less efficient use of water, which is a discouraging feature for moisture stress situations in drylands. Further, tillage practices are known to influence K availability by modifying other factors such as oxygen or aeration, temperature, soil moisture and the positional availability of applied K. Despite medium soil test values (<290 kg ha −1 ) and no application of K fertilizer in the present study, K emerged as a key indicator, contributing 34.3% to the SQI.
Among the different land use systems, PCs clearly separated the agriculture and forest samples from the afforestation and non-arable lands samples. The SQI obtained under each land use was maximum under afforestation (1.41), followed by non-arable (1.07), forestry (1.04) and agriculture (0.758), respectively ( Figure 4). This showed that the contribution of ER to the SQI was highest under afforestation (0.633), followed by forestry (0.554), NA (0.534) and agriculture (0.419). For K, the maximum contribution to the SQI was observed under afforestation (0.516), followed by forestry (0.321) and non-arable lands (0.271), and the minimum was observed under agriculture (0.264). P contributed most to the SQI under non-arable lands (NA) (0.292) and least under agriculture (0.075). indicators of soil health [43]. Under different land uses in submontane Punjab, organic carbon contributed the mostto soil quality (28.5%), followed by the available K (19.4%), electrical conductivity (18.3%), K-factor of universal soil loss equation (USLE) (14.9%), plant available water (10.5%), and clay (8.3%) [15]. In the present study, the contribution of potential indicators like ER, P and K toward thesoil quality index (SQI) was substantial. ER plays a dominant role in identifying, based on the soil's characteristics, the areas that may require more careful land management techniques. Particle size distribution and organic matter content are the two most important indicators of erodibility. Nair (1984) [44] reported that agro-forestry, agro-horticulture and grassland systems have the potential to reduce both the runoff and erosion, and to maintain better soil organic matter, which in turn will help improve the fertility status of the soil. The available K is important for osmotic regulations. Potassium provides much of the osmotic pull that draws water into plant roots. Plants that are K-deficient are less able to withstand water deficits, mostly because of their inability to make full use of the available water. The malfunctioning of stomata due to a deficiency of this nutrient has been related to lower rates of photosynthesis and a less efficient use of water, which is a discouraging feature for moisture stress situations in drylands. Further, tillage practices are known to influence K availability by modifying other factors such as oxygen or aeration, temperature, soil moisture and the positional availability of applied K. Despite medium soil test values (<290 kg ha −1 ) and no application of K fertilizer in the present study, K emerged as a key indicator, contributing 34.3% to the SQI. Among the different land use systems, PCs clearly separated the agriculture and forest samples from the afforestation and non-arable lands samples. The SQI obtained under each land use was maximum under afforestation (1.41), followed by non-arable (1.07), forestry (1.04) and agriculture (0.758), respectively (Figure 4). This showed that the contribution of ER to the SQI was highest under afforestation (0.633), followed by forestry (0.554), NA (0.534) and agriculture (0.419). For K, the maximum contribution to the SQI was observed under afforestation (0.516), followed by forestry (0.321) and non-arable lands (0.271), and the minimum was observed under agriculture (0.264). P contributed most to the SQI under non-arable lands (NA) (0.292) and least under agriculture (0.075). The order of contribution of the selected indicators to the SQI was 53.5%, 34.3% and 19.9% for ER, P and K, respectively. This clearly indicates the importance of weighted factor attributes in thePCA. The high weight of ER indicated its highest variability in the data set. The specific contribution of the minimum data set (MDS) to theSQI is presented througharadar plot ( Figure 5). The order of contribution of the selected indicators to the SQI was 53.5%, 34.3% and 19.9% for ER, P and K, respectively. This clearly indicates the importance of weighted factor attributes in the PCA. The high weight of ER indicated its highest variability in the data set. The specific contribution of the minimum data set (MDS) to the SQI is presented through radar plot ( Figure 5). Sustainability 2020, 12, x FOR PEER REVIEW 11 of 13

Conclusions
The overall state of the soil based on the most sensitive soil indicators under different land uses was highestunder agricultural land use, followed by forests and non-arable lands, and lowest under afforestation. Among the different land use systems, PCs clearly separated the agriculture and forest samples from the afforestation and non-arable lands samples.ER, P and K contents were the most sensitive/main indicators to affect the soil quality in the lower Shiwaliks of northwestern India.The usefulness of ER lies mainly in identifying, based on the soil's characteristics, the areas that may require more careful land management techniques. The contribution of the selected indicators to the SQI was highest for ER, followed by P and K. These indicators may be used by researchers for the real-time monitoring of soil health and ecological processes in the future, under various land use systems in degraded agroecosystems. On the basis of the overall state of the soil, a quantitative model can be developed for assessing the site-specific soil loss tolerance limit for sustaining plant growth in the Shivaliks.

Conclusions
The overall state of the soil based on the most sensitive soil indicators under different land uses was highest under agricultural land use, followed by forests and non-arable lands, and lowest under afforestation. Among the different land use systems, PCs clearly separated the agriculture and forest samples from the afforestation and non-arable lands samples.ER, P and K contents were the most sensitive/main indicators to affect the soil quality in the lower Shiwaliks of northwestern India. The usefulness of ER lies mainly in identifying, based on the soil's characteristics, the areas that may require more careful land management techniques. The contribution of the selected indicators to the SQI was highest for ER, followed by P and K. These indicators may be used by researchers for the real-time monitoring of soil health and ecological processes in the future, under various land use systems in degraded agroecosystems. On the basis of the overall state of the soil, a quantitative model can be developed for assessing the site-specific soil loss tolerance limit for sustaining plant growth in the Shivaliks.