Variation of Physico-Chemical Properties among Different Soil Orders under Different Land Use Systems of the Majha Region in North-Western India

: The impact of different soil orders and land use systems on the distribution of physico-chemical properties is the most critical matter to address in order to maintain sustainable agricultural production. Hence, the present investigation was carried out to study the variation in the physico-chemical characteristics of soil in diverse land use systems (LUSs), i.e., agriculture, horticulture, and forestry, under major soil orders (entisol, inceptisol, and alﬁsol) in the Majha region of Punjab. A total of 225 depth-wise (at 0–20 cm, 20–40 cm, 40–60 cm, 60–80 cm and 80–100


Introduction
The state of Punjab covers an area of 5.03 mha (comprising 1.5% area of the country), of which 4.20 mha of land comes under cultivation.The state has varied climatic conditions, with large variations in nutrients content in soil.The soils of Punjab, in north-western India, developed from alluvium parent material deposited by the Indus system s rivers.The state has been well-known for its agricultural activities and plantations, supplying 40-50% of the rice and 60-65% of the wheat in the central pool.Rising socio-economic demands and an increasing population put pressure on agricultural production systems, resulting in paradigm shifts in patterns of land use [1].As agriculture becomes more advanced and intensive farming systems develop, there have been gradual losses in soil quality that have caused drastic reductions in yield.The rate of degradation of quality is governed by land use systems, the type of soil, topography, and climate factors.Several studies have indicated that inappropriate land use causes physicochemical and biological deterioration of the soil [2,3].Changes in land use are considered to be an important component and a primary cause of global environmental change [4].According to Heluf and Wakene, changes in land use and management practices frequently alter most physical, chemical, and biological characteristics of most soil, which has been reflected in agricultural production [5].Different land use systems viz., agricultural (irrigated and rainfed), horticultural and forestry systems modify the physico-chemical characteristics of the soil as well as the nutrient content [6].Dynamic properties such as soil organic matter (SOM), cation exchange capacity, soil pH, EC and texture are sensitive to land use practices and can provide valuable information about important soil processes such as nutrient cycling, decomposition and formation of SOM and overall productivity potential [7].
Land use systems and their management practices cause qualitative and quantitative changes in soil organic matter, which has a significant impact on the chemical, physical, and biological properties of the soil, which directly influence nutrient availability to plants [8].Anthropogenic activities like input and management of the land use system, as well as the addition of plant litter, can all have a direct impact on soil properties.Land use and management changes can have a major impact on the organic matter properties of soil; conversion from forest to cropland, combined with conventional tillage and a lack of biomass return to soil, has been shown to reduce the degree of humification of the organic matter of soil [9,10].Factors such as vegetation coverage, litter fall, root impact, and disturbance or management system can all contribute to major differences in surface soil organic carbon (SOC) and other soil properties across land use types [11,12].Agricultural land use had significantly more available P and K than agroforestry and grassland systems in different soils of Punjab [13].Datta et al. evaluated several land use types (Guava, Mango, Litchi, Jamun, Prosopis, Eucalyptus, and Rice-wheat cropping system) in alfisols and found that soil pH increased with depth, due to leaching and basic cation accumulation in the lower depth of the profiles in all land uses [14].A few studies on soil orders found a higher OC status in inceptisols and alfisols.This might be due to the accumulation of more organic residues in the soils.The distribution of CaCO 3 is also varied among soil orders and the highest calcium carbonate content was noticed in the aridisols, followed by the inceptisols, entisols, and alfisols [15].Thus, the variability in the distribution of soil properties such as pH, EC, OC, CaCO 3 , and available P and K in different soils is greatly influenced by soil orders under different land use systems in the particular region [16].
The soils of the Majha region of Punjab belong to three soil orders viz., entisols, inceptisols and alfisols, out of twelve soil orders categorized according to the USDA taxonomic classification which are officially recognized as being in Punjab [17].In the Majha region, the dominant area is occupied by inceptisols, followed by entisols and alfisols.Entisols serve as precursors to all other soil orders during their evolutionary history [18].Inceptisols are relatively young soils that have fewer features than mature soils but retain a close resemblance to the parent material.Inceptisols are the most productive soils; they include soils from Ustic to Udic moisture regimes that have altered B horizons due to some chemical weathering processes [19].On the other hand, alfisols are commonly found on stable geomorphic surfaces and developed over a wide range of parent materials.Sharma et al. studied alfisols, which are commonly associated with inceptisols and found on concave mountain slopes (Udic moisture regime), as well as the piedmont and plains (Ustic moisture regime).The major driving processes of alfisols are eluviation and illuviation [20].Physico-chemical properties play an important role in nutrient availability in different group of soils.These properties in soil profiles not only depend on soil formation process, mineralogy, climate and topography of soil but are also influenced by the change in land use system in an area [16].
So far, little attention has been paid to monitoring the distribution of physico-chemical properties in different land use systems in the study area.Soil productivity in various land uses is directly related to the management and conservation of soil fertility over time.A detailed investigation of the processes leading to the delineation of various soil fertility parameters can aid in the breakthrough of ways and means to achieve sustainable high levels of production in various land use systems.Keeping in mind that several distinct land use systems exist in different soil orders with different fertility status, the present study aims to characterize the soil properties in different land use systems under distinct soil orders in the Majha region of Punjab.This investigation fills a knowledge gap regarding the variation of soil properties such as pH, EC, OC, CaCO 3 , available P and K in agriculture, horticulture, and forest land use systems falling in the entisol, inceptisol and alfisol soil orders.The outcome of the research will help students, researchers, scientists, and agricultural administrators in managing various land use systems with the aim of preserving soil fertility and productivity of different land use systems in distinct soil orders in the Majha region of Punjab in north-western India.

Geographical Area and Location
The Punjab state lies between 29 • 33 N and 32 • 32 N latitude and 73 • 53 E to 76 • 56 E longitude (Figure 1).It has an average elevation of 300 m (980 feet) above sea level.An investigation was carried out in order to study the variations of physico-chemical properties under different land use systems in dominant soil orders in the Majha region of Punjab at Punjab Agricultural University, Ludhiana, during 2019-2022.This region consists of four districts, i.e., Amritsar, Gurdaspur, Pathankot, and Tarn-Taran, having three dominant soil orders: entisols, inceptisols, and alfisols.It is characterized by three main land use systems such as agriculture, horticulture, and forestry.Wheat, rice, maize, sugarcane, peas, potato, tomato, oil seeds, and pulses are the most important crops grown in the agricultural system.The main fruit crops in the horticultural system are mango, guava, peach, rock pear and litchi.Mixed forest, poplar, and eucalyptus plantations are the most important forest land use systems in this region.Agricultural systems are characterized by the use of chemical fertilizers and farmyard manure.Crop waste management methods differ depending on the type of agricultural area.Farmers burn crop waste after harvesting their summer crops to clear their fields and prepare for the next sowing.Crop waste, on the other hand, is being incorporated into the soil in some areas.Because of the heavy use of machinery and farm equipment, the subsurface soils under agricultural fields (primarily the rice-wheat cropping system) were extremely hard as depth increased.Horticultural land use systems were characterized by the addition of organic residues through leaf litter and the application of chemical fertilizer and manures.However, some horticultural fields were not properly managed and thus accumulated less organic matter.The regular addition of organic residues in the form of leaves and other debris distinguished the forest land use systems.

Site Selection and Sample Collection
Three dominant land use systems, i.e., agriculture, horticulture and forest, were selected under each soil order (entisol, inceptisol and alfisol) in the Majha regions of Punjab.Depth-wise soil samples (0-20 cm, 20-40 cm, 40-60 cm, 60-80 cm and 80-100 cm) were collected from five locations for each land use system in all soil orders in the region using a stainless-steel auger, and the global positioning system of each profile was recorded.Figure 2 depicts the classification of land use systems within each soil order in the Majha regions of Punjab.In the study area, three soil orders were recognized.Three LUSs were identified in each soil order, five replications were considered in each LUSs, and soil samples were collected at five depths in each replication.Thus, a total of 225 samples (45 locations, 5 depths) were collected from three land use systems in all soil orders for the present study.

Site Selection and Sample Collection
Three dominant land use systems, i.e., agriculture, horticulture and forest, were selected under each soil order (entisol, inceptisol and alfisol) in the Majha regions of Punjab.Depth-wise soil samples (0-20 cm, 20-40 cm, 40-60 cm, 60-80 cm and 80-100 cm) were collected from five locations for each land use system in all soil orders in the region using a stainless-steel auger, and the global positioning system of each profile was recorded.Figure 2 depicts the classification of land use systems within each soil order in the Majha regions of Punjab.In the study area, three soil orders were recognized.Three LUSs were identified in each soil order, five replications were considered in each LUSs, and soil samples were collected at five depths in each replication.Thus, a total of 225 samples (45 locations, 5 depths) were collected from three land use systems in all soil orders for the present study.

Analysis of Soil Samples
The standard methodologies were used to ascertain different physico-chemical parameters for soil.A glass electrode pH meter (Elico, LI-120, Hyderabad, India) was used to determine the pH from soil suspension in water at a ratio of (1:2), after keeping the same supernatant solution at equilibration for 24 h, the EC was determined using a con-−1

Analysis of Soil Samples
The standard methodologies were used to ascertain different physico-chemical parameters for soil.A glass electrode pH meter (Elico, LI-120, Hyderabad, India) was used to determine the pH from soil suspension in water at a ratio of (1:2), after keeping the same supernatant solution at equilibration for 24 h, the EC was determined using a conductivity bridge (Elico, CM-183, Hyderabad, India) and expressed in dS m −1 [21].Soil organic carbon (SOC) levels were determined in soil samples using a potassium dichromate rapid titration method and the results were expressed as a percentage [22].The titration of soil suspension with 0.5 N H 2 SO 4 in the presence of bromothymol blue and bromocresol green indicators yielded the calcium carbonate concentration of soil, which was reported as a percentage [23].The available phosphorus (P) was extracted from the soil using 0.5 M NaHCO 3 and the content of P in the extract was measured by stannous chloride reduced molybdophosphric acid blue color in HCl system by using spectrophotometer measured at 760 nm wavelength and the results were represented in kg ha −1 [24].The available potassium (K) in the soils was extracted using neutral N normal ammonium acetate in a 1:5 soil to extractant ratio and the concentration of K contained in the extract was evaluated using a flame photometer [25].

Statistical Analysis
Using the Agricolae package included in the R programming language, distribution patterns of several physico-chemical characteristics (at various soil orders, land use systems and depth) were subjected to an ANOVA.The Tukeys post-hoc test of significance was used to determine the level of significance (p-value < 0.05).Using R programming, the Pearson correlation coefficient was used to calculate correlations between various soil parameters.Using the Devtools and Factoextra packages included in R programming, principal component analysis (PCA) was carried out on a correlation matrix between the 6 variables of various physico-chemical parameters among different soil orders and land use systems.

Variation of Soil pH under Soil Orders under Different Land Use Systems
The distribution pattern of soil pH at various soil orders under different land use systems is presented in Figure 3. Soils under all three soil orders were neutral to slightly alkaline in nature, and optimal for most of the crops.In the soil profiles of various land use systems, the soil pH is varied from 6.80 to 7.50 in entisols, 7.64 to 8.34 in inceptisols and 6.94 to 7.87 in alfisols.In each soil order, the impact of soil depth and LUSs on soil pH distribution varies greatly.The pH in entisols varies significantly among different land use systems, while depth does not.The highest pH values were found in agricultural and horticultural LUSs, whereas the lowest values were found in FLUS.All the LUSs and depth had a significant impact on soil pH in the inceptisol.When compared to ALUS and HLUS, FLUS exhibited significantly lower soil pH.In all land use systems, soil pH is much lower in surface horizons and increases with depth.Different LUSs significantly affect the soil pH in alfisols, although depth has no significant impact on soil pH.ALUS had the highest soil pH and FLUS had the lowest soil pH in alfisols.A perusal of the data revealed that soil orders and land use systems (LUSs) had a significant impact on soil pH.Among soil order, highest average soil pH (8.09) was observed in inceptisols, followed by alfisols (7.38) and entisols (7.10).Among the examined soils, the highest pH was observed in inceptisols under ALUS.This might be due to a higher concentration of base-forming cations provided by fertilizers and chemical amendments added to the soil.Entisols and alfisols, on the other hand, have a lower pH, possibly owing to agricultural wastes and litter fall from the forest land use system continually adding organic matter to these soils, and the breakdown of the litter also releasing weak organic acids [16].Several studies have found a similar trend in soil pH [7,26].Irrespective of the soil orders, the soil pH was found to increase with depth in the profiles [13,15].The subsurface leaching of basic cations and salts to the deeper layer of soil under all systems results in the gradual rise in the pH of soils with depth [7,27].(7.38) and entisols (7.10).Among the examined soils, the highest pH was observed in inceptisols under ALUS.This might be due to a higher concentration of base-forming cations provided by fertilizers and chemical amendments added to the soil.Entisols and alfisols, on the other hand, have a lower pH, possibly owing to agricultural wastes and litter fall from the forest land use system continually adding organic matter to these soils, and the breakdown of the litter also releasing weak organic acids [16].Several studies have found a similar trend in soil pH [7,26].Irrespective of the soil orders, the soil pH was found to increase with depth in the profiles [13,15].The subsurface leaching of basic cations and salts to the deeper layer of soil under all systems results in the gradual rise in the pH of soils with depth [7,27].

Variation of Soil Electrical Conductivity (EC) under Soil Orders under Different Land Use Systems.
The data provided in Figure 4 demonstrated the variation of soluble salt concentration, i.e., soil EC in profiles of various soil orders under different land use systems.The EC values in the soil profiles of various land use systems ranged from 0.13 to 0.42, 0.19 to 0.54 and 0.19 to 0.46 dS m −1 in entisols, inceptisols and alfisols, respectively.The distribution of soil EC did not differ significantly across soils from different LUSs under entisols; however, soil depth has a significant impact on the variance of soluble salts in the soil

Variation of Soil Electrical Conductivity (EC) under Soil Orders under Different Land Use Systems
The data provided in Figure 4 demonstrated the variation of soluble salt concentration, i.e., soil EC in profiles of various soil orders under different land use systems.The EC values in the soil profiles of various land use systems ranged from 0.13 to 0.42, 0.19 to 0.54 and 0.19 to 0.46 dS m −1 in entisols, inceptisols and alfisols, respectively.The distribution of soil EC did not differ significantly across soils from different LUSs under entisols; however, soil depth has a significant impact on the variance of soluble salts in the soil profile of investigated soils.Significantly higher soil EC was observed under the surface horizon, and soluble salt concentration decreased with an increase in depth in all the land use systems of entisols.The depth of the soil has no influence on the salt concentration in profiles, while the distribution of soil EC varies significantly among distinct LUSs within the inceptisols soil order.ALUS and HLUS had significantly higher soil EC than FLUS.In alfisols, LUSs has a significant influence on the distribution of EC, but there is no significant variation in soil depth.FLUS had significantly higher soil EC when compared to agricultural and horticultural LUSs in alfisols.The distribution of soil EC among LUSs in entisol, inceptisol and alfisols follows the order HLUS > FLUS > HLUS, ALUS > HLUS > FLUS and FLUS > ALUS > HLUS, respectively.The mean values of soil orders are shown to have a nonsignificant influence on the distribution of soil EC.Among different soil orders, the highest value of EC (0.34 dS m −1 ) was recorded in the profile of inceptisols, followed by alfisols (0.29 dS m −1 ) and lowest EC value (0.28 dS m −1 ) was observed in the entisols.The higher EC in inceptisols under ALUS could be attributed to greater evaporation, and excessive usage of fertilizer salts, which results in the deposition of soluble salts [16], resulting in increased salt build-up in soils [28].In general, the lower values of EC in entisols and alfisols could be attributed to sandy loam texture, which allows free leaching of soluble salts, and also due to low amounts of soluble salts in the original deposits [15].The EC in soils changed significantly as the soil depth increased.Irrespective of soil orders and LUSs, there was no specific pattern in the distribution of soil EC among soil profiles, but it showed a subsequent decrease as depth increased.The EC of the surface layer was higher, which might mean more nutrient ions at the surface; however, as depth increased, the EC in the profile decreased as nutrient ions declined [29].Similar patterns of decreasing soil EC with an increase in depth of all the land uses have been observed [7,13].
followed by alfisols (0.29 dS m −1 ) and lowest EC value (0.28 dS m −1 ) was observed in the entisols.The higher EC in inceptisols under ALUS could be attributed to greater evaporation, and excessive usage of fertilizer salts, which results in the deposition of soluble salts [16], resulting in increased salt build-up in soils [28].In general, the lower values of EC in entisols and alfisols could be attributed to sandy loam texture, which allows free leaching of soluble salts, and also due to low amounts of soluble salts in the original deposits [15].The EC in soils changed significantly as the soil depth increased.Irrespective of soil orders and LUSs, there was no specific pattern in the distribution of soil EC among soil profiles, but it showed a subsequent decrease as depth increased.The EC of the surface layer was higher, which might mean more nutrient ions at the surface; however, as depth increased, the EC in the profile decreased as nutrient ions declined [29].Similar patterns of decreasing soil EC with an increase in depth of all the land uses have been observed [7,13].

Variation of Soil Organic Carbon (SOC) under Soil Orders under Different Land Use Systems
The status of soil organic carbon in profiles of different soil orders is shown in Figure 5.The results revealed that in each soil order, depth had a significant impact on the distribution of SOC content, whereas all the LUSs in the entisol and alfisol soil orders have a

Variation of Soil Organic Carbon (SOC) under Soil Orders under Different Land Use Systems
The status of soil organic carbon in profiles of different soil orders is shown in Figure 5.The results revealed that in each soil order, depth had a significant impact on the distribution of SOC content, whereas all the LUSs in the entisol and alfisol soil orders have a substantial impact on SOC except for the inceptisols.In all the soil orders, the mean values of SOC were highest in FLUS, followed by HLUS and comparatively lowest in ALUS.Irrespective of soil orders and LUSs, the surface soils had substantially larger SOC concentrations than subsurface layers in all soil profiles investigated.The SOC varied substantially between 0 to 20 and 20 to 40 cm soil depth, but not between 40 to 60, 60 to 80 and 80 to 100 cm soil depth.In all soil orders, the 80 to 100 cm depth had the lowest soil OC levels.The organic carbon content in the Majha region of Punjab typically declines due to the hyperthermic soil temperature regime.The soil OC in soil profiles of various land use systems ranged from 0.14 to 0.99% in entisols, 0.21 to 0.69% in inceptisols and 0.15 to 0.72% in alfisols.The distribution of SOC is not significantly affected among different soil orders in the region.Among the soil orders studied, entisols under FLUS had a higher SOC content than other soil orders due to there being more organic residues via leaf litter and other plant residues (roots, stubbles), whereas alfisols had the lowest concentration due to their limited vegetation and widespread erosion impact, resulting in a lower OC status in these soils [16].Furthermore, the prolonged and intensive cultivation and other agricultural activities that increased the mineralization of organic residues, along with the activity of soil microorganisms, may be the cause of the low level of SOC in alfisols and entisols [15,30].Nevertheless, higher OC in entisols attributed to adequate moisture availability, which led to better vegetation growth.As a result, these soils have higher levels of organic residue and consequently, more OC in the soil [30].In addition, lower temperatures and prolonged rainfall intensity and duration favor higher OC in soil [20].FLUS are more productive than ALUS in terms of resource capture and biomass growth, resulting in increased carbon inputs to the soil.Soil carbon in FLUS is facilitated by inputs from litter production and turnover, as well as soil stabilization processes [3,4].The maximum OC level was found at the surface, owing to litter and organic residue deposition, regardless of soil orders or land use systems and this decreased with depth [29].A similar pattern of decreasing OC content with depth has been documented [28].In general, the organic carbon content of arid and semi-arid soils was low, which might be attributed to the poor biological activity or the quick degradation of organic biomass under frequent torrid conditions [31].
lower OC status in these soils [16].Furthermore, the prolonged and intensive cultivation and other agricultural activities that increased the mineralization of organic residues, along with the activity of soil microorganisms, may be the cause of the low level of SOC in alfisols and entisols [15,30].Nevertheless, higher OC in entisols attributed to adequate moisture availability, which led to better vegetation growth.As a result, these soils have higher levels of organic residue and consequently, more OC in the soil [30].In addition, lower temperatures and prolonged rainfall intensity and duration favor higher OC in soil [20].FLUS are more productive than ALUS in terms of resource capture and biomass growth, resulting in increased carbon inputs to the soil.Soil carbon in FLUS is facilitated by inputs from litter production and turnover, as well as soil stabilization processes [3,4].The maximum OC level was found at the surface, owing to litter and organic residue deposition, regardless of soil orders or land use systems and this decreased with depth [29].A similar pattern of decreasing OC content with depth has been documented [28].In general, the organic carbon content of arid and semi-arid soils was low, which might be attributed to the poor biological activity or the quick degradation of organic biomass under frequent torrid conditions [31].

Variation of Soil Calcium Carbonate (CaCO 3 ) under Soil Orders under Different Land Use Systems
The variation of CaCO 3 in different soil orders are given in Figure 6.The CaCO 3 in soils of different land use systems varied from 0.75 to 2.07% in entisols, 1.07 to 3.32% in inceptisols and 0.93 to 2.29% in alfisols.The data illustrated that there is a significant variation in CaCO 3 content among the different soil orders at various depths.The interesting feature of the study was the presence of calcium carbonate in all the profiles of the soil orders.In each soil order, the impact of land use systems and depth on calcium carbonate concentration is highly variable.Land use system and depth did not significantly affect the distribution of CaCO 3 under inceptisols.ALUS and FLUS both had the highest concentrations; however, FLUS had a lower concentration.The CaCO 3 composition of the entisol and alfisol profiles is significantly influenced by land use systems and depth.Among the three land uses in entisol soil order, FLUS had a significantly higher CaCO 3 concentration, followed by HLUS and the lowest was recorded in ALUS.In the alfisols, however, significantly higher values were detected in FLUS and ALUS than to HLUS.Thus, there was no specific trend in the distribution of calcium carbonate in soils of distinct LUSs.Based on the results of the study, the distribution of CaCO 3 differs significantly among land use systems in all the soil orders except inceptisols.Thus, the results indicate the impact of land use systems on the distribution of CaCO 3 in each soil order in the Majha region.The distribution of CaCO 3 is also affected by soil formation process, mineralogy and climate of a particular region [28].The highest percentage of CaCO 3 (1.75%) was recorded in inceptisols, followed by alfisols (1.43%) and entisols (1.20%).When compared to other soil orders, the higher concentration of inceptisols may be the result of a greater accumulation of CaCO 3 throughout the soil profiles.However, there was no trend in the depth-wise distribution of calcium carbonate in soils of different LUSs.The literature revealed that the irregular distribution of CaCO 3 with depth was caused by the soil s alluvial nature.Regardless of the soil orders and land uses, the values of CaCO 3 varied considerably from one another and increased with depth [32].The leaching of calcium, which then precipitated as carbonate at a lower depth, suggested that a larger concentration of CaCO 3 was associated with depth.Similar findings were also reported in the literature [33].The CaCO 3 leaching might be caused by high permeability and sand content in all land uses.As a consequence, the drainage system, water supply and soil depth, all played a role in the variation of CaCO 3 in different LUSs among soil orders [34].In addition to rainfall and temperature, irrigation water quality, soil texture and parent material of an area all affect the distribution of calcium carbonate in different soil orders under various LUSs [35].

Variation of Available Phosphorus (P) under Soil Orders under Different Land Use Systems
The distribution of available P in studied soils was found to be significantly influenced by soil orders and land use systems at various depths (Figure 7).The available P in soil profiles of various land use systems ranged from 7.77 to 41.84, 10.56 to 40.23 and 7.24 to 39.51 kg ha −1 in entisols, inceptisols and alfisols, respectively.The distribution of P in land use systems of all soil orders did not follow any trend.The pattern of available P in entisol, inceptisol and alfisols is as follows: HLUS > ALUS > FLUS, FLUS > HLUS > ALUS and ALUS > FLUS > HLUS.ALUS had higher available phosphorus content in alfisols than in forest and horticulture, which could be attributed to the use of P fertilizers in agricultural fields [36].Forest and horticultural LUSs had the highest P levels in inceptisols and entisols, which could be attributed to the frequent addition of organic residues, during decomposition which releases more available P in soils.This variability effectively shows the impact of land uses and their management strategies on soil available P in different soil orders in an area.Soil P is much greater in the surface layers of all soil profiles studied, and it declines as depth rises.When compared to the land use systems of entisols and alfisols, agricultural, horticultural and forest systems in inceptisols had a comparatively higher quantity of available P, resulting in significantly higher available P in inceptisols, followed by entisols (20.21 kg ha −1 ) and alfisols (19.72 kg ha −1 ).The higher amount of P in inceptisols might be due to the addition of organic matter through leaf litter (from HLUS and FLUS) and the continuous application of fertilizers in ALUS in the soils [37].Several factors, including soil amendments, farming methods and weather pat-

Variation of Available Phosphorus (P) under Soil Orders under Different Land Use Systems
The distribution of available P in studied soils was found to be significantly influenced by soil orders and land use systems at various depths (Figure 7).The available P in soil profiles of various land use systems ranged from 7.77 to 41.84, 10.56 to 40.23 and 7.24 to 39.51 kg ha −1 in entisols, inceptisols and alfisols, respectively.The distribution of P in land use systems of all soil orders did not follow any trend.The pattern of available P in entisol, inceptisol and alfisols is as follows: HLUS > ALUS > FLUS, FLUS > HLUS > ALUS and ALUS > FLUS > HLUS.ALUS had higher available phosphorus content in alfisols than in forest and horticulture, which could be attributed to the use of P fertilizers in agricultural fields [36].Forest and horticultural LUSs had the highest P levels in inceptisols and entisols, which could be attributed to the frequent addition of organic residues, during decomposition which releases more available P in soils.This variability effectively shows the impact of land uses and their management strategies on soil available P in different soil orders in an area.Soil P is much greater in the surface layers of all soil profiles studied, and it declines as depth rises.When compared to the land use systems of entisols and alfisols, agricultural, horticultural and forest systems in inceptisols had a comparatively higher quantity of available P, resulting in significantly higher available P in inceptisols, followed by entisols (20.21 kg ha −1 ) and alfisols (19.72 kg ha −1 ).The higher amount of P in inceptisols might be due to the addition of organic matter through leaf litter (from HLUS and FLUS) and the continuous application of fertilizers in ALUS in the soils [37].Several factors, including soil amendments, farming methods and weather patterns, have a substantial impact on the widely varying P concentration in distinct soil orders under various land use systems.The amount of P was higher on the surface soils but decreased with an increase in depth of all soil orders, which might be attributable to its higher release proportion to the larger accumulation of organic matter on the surface of the soil [7].The decrease in the concentration with depth might be due to the constant decline in content and organic matter addition in soils [31].Intensive farming, inadequate application of organic manures and P fertilizers, ongoing exclusion of crop residues and resulting nutrient decline in various soils may all be contributing factors to the decline in available P in the soil profiles of several LUSs [38].

Variation of Available Potassium (K) under Soil Orders under Different Land Use Systems
The results presented in Figure 8 indicated the variability of available K in soil orders of the Majha region of Punjab.The available K in all of the examined soils was greater in soil profiles across all land use systems and ranged from 98.37 to 334.68, 94.51 to 230.18 and 93.01 to 367.39 kg ha −1 in entisols, inceptisols and alfisols, respectively.In each soil order, the variation of soil K is differed among land use systems and depth.The study reported that land use systems and depth strongly affected the available K in alfisols, but there was no significant variation was observed in entisols and inceptisols.In alfisols and inceptisols, the highest available K was found in FLUS, followed by HLUS and the lowest amount was found in ALUS; however, in entisols, ALUS recorded the highest amount of available K, followed by FLUS and HLUS.The high available K values in FLUS might be attributable to nutrient recycling caused by dropping leaves of different tree species [13].Soil orders were found to have a non-significant effect on available K distribution in various soils.Among the soil orders investigated, higher available K was observed in entisols (200.59kg ha −1 ), followed by alfisols (178.58 kg ha −1 ) and inceptisols (157.76 kg ha −1 ).The distribution of available K in soil follows a distinct geomorphic pattern and is heavily influenced by parent material, weathering, particle size and management practices.The majority of potassium in soil is in the form of minerals such as feldspars and micas [39].Entisols had the highest amount of available K when compared to inceptisols and alfisols, which correlates strongly to the governance of minerals rich in K as well as the addition of a higher proportion of OM, which aided in the restoration of soil nutrient status [15].

Variation of Available Potassium (K) under Soil Orders under Different Land Use Systems
The results presented in Figure 8 indicated the variability of available K in soil orders of the Majha region of Punjab.The available K in all of the examined soils was greater in soil profiles across all land use systems and ranged from 98.37 to 334.68, 94.51 to 230.18 and 93.01 to 367.39 kg ha −1 in entisols, inceptisols and alfisols, respectively.In each soil order, the variation of soil K is differed among land use systems and depth.The study reported that land use systems and depth strongly affected the available K in alfisols, but there was no significant variation was observed in entisols and inceptisols.In alfisols and inceptisols, the highest available K was found in FLUS, followed by HLUS and the lowest amount was found in ALUS; however, in entisols, ALUS recorded the highest amount of available K, followed by FLUS and HLUS.The high available K values in FLUS might be attributable to nutrient recycling caused by dropping leaves of different tree species [13].Soil orders were found to have a non-significant effect on available K distribution in various soils.Among the soil orders investigated, higher available K was observed in entisols (200.59kg ha −1 ), followed by alfisols (178.58 kg ha −1 ) and inceptisols (157.76 kg ha −1 ).The distribution of available K in soil follows a distinct geomorphic pattern and is heavily influenced by parent material, weathering, particle size and management practices.The majority of potassium in soil is in the form of minerals such as feldspars and micas [39].Entisols had the highest amount of available K when compared to inceptisols and alfisols, which correlates strongly to the governance of minerals rich in K as well as the addition of a higher proportion of OM, which aided in the restoration of soil nutrient status [15].Regardless of soil orders, the surface layer recorded the highest quantity of available K, but it decreased as depth increased.The replenishing of soil with organic matter through crop residues and the addition of animal manures were both responsible for the greatest amount of K in the surface horizons [21].The decrease in the concentration with depth might be due to the constant decline in content and organic matter addition in soils [40].Our findings are in agreement with many other studies [41,42].The soils are very rich in K.It is therefore advisable that farmers apply the requisite amount of potassium to replenish its removal by the respective crops of the system, so that pressure on soil K reserve may be avoided, which may otherwise result in a negative balance in soil in near future.Thus, multiple factors such as depth, soil amendments and climatic conditions under different land-use systems all contributed to increased soil available K content.

The Correlation Studies between Various Physico-Chemical Parameters
The correlation analysis was carried out between soil physico-chemical parameters of the soils examined (Table 1).The analysis revealed that soil available P showed a significantly positive correlation with EC (p < 0.01 and r = 0.47) and OC (p < 0.01 and r = 0.66).However, it showed a negative correlation between with pH (r = −0.

The Correlation Studies between Various Physico-Chemical Parameters
The correlation analysis was carried out between soil physico-chemical parameters of the soils examined (Table 1).The analysis revealed that soil available P showed a significantly positive correlation with EC (p < 0.01 and r = 0.47) and OC (p < 0.01 and r = 0.66).However, it showed a negative correlation between with pH (r = −0.03)and CaCO 3 (p < 0.01 and r = −0.38).The available soil K was positively associated with EC (r = 0.   9).High eigenvalue principal components (PCs) were considered to represent the maximum variations among different soil properties [43].The findings show that the first two principal components (PC1 and PC2) with >1 eigenvalues explained 69.83 percent of the total variances in the original dataset.The amount of variability expressed by PC1 47.19%, with an eigenvalue of 2.83, which indicates that soil pH had the highest positive factor loading value, i.e., 0.55, and that it contributes 30.44% in the total variations of soil properties.The PC1 was specified to the pH factor of the soil.The second principal component (PC2) explained about 22.63% of the variance, with an eigenvalue of 1.35 and the highest positive loading value was for OC (0.68), which contributes 46.92% in the total variation among other soil variables.The PC2 was specified to the SOC factor.In accordance with the PCA analysis of the various soil parameters, soil pH, followed by SOC, were the most significant factor that contributes to and represents the most variability in determining the soil quality across different LUSs and soil orders in the Majha region.The soil pH and OC are two of the most important basic soil properties which influence the availability of macro and micronutrients under different LUSs.The SOC is a measure of the quality and quantity of soil organic matter under agricultural, horticultural and forest land uses in all the soil orders.The frequent accumulation of soil organic matter under forest and horticulture land use systems increases SOC content, thereby increasing the availability of nutrients and maintaining soil fertility and productivity.In the first principal component of PCA-biplot, among various land uses and soil orders, the ALUS and HLUS under inceptisols had the highest positive loadings and represent the greatest variability in the examined soils.Thus, soils of inceptisols have greater soil pH, which reflects the maximum variability among various soil parameters and governs the availability of most of the essential nutrients to the crops in the Majha region.In the second component, the FLUS under alfisols and entisols had the maximum positive loadings, while ALUS had the highest negative loadings in all soil orders.Poor land use and management strategies made significant contributions to soil degradation, which lowers soil fertility and productivity [44,45].Hence, FLUS has the potential to improve soil fertility and reduce soil degradation, by improving the physico-chemical properties of soil in horticultural and agricultural systems.

Conclusions
The distribution of physico-chemical properties varies under different soil order land use systems at different depths in profiles of the soil samples examined.The r of the research indicated that soil pH, CaCO3 and available P content in soil signific varied among different soil orders, whereas EC, SOC and available K were not si cantly influenced by soil orders.The highest soil pH, EC and CaCO3 values were rep in inceptisols, followed by alfisols and the lowest values were found in entisols.Th tern of SOC variation in soil orders is as follows: entisol > inceptisols > alfisols.Incep had the highest available P, followed by entisols and alfisols.Entisols had the hi amount of available K, followed by alfisols and inceptisols.In all soil orders, the dis tion of the physico-chemical properties of soil under different land use systems di follow any specific trend.When compared to other land uses, the agriculture system the highest values of soil pH, EC, CaCO3 and the lowest amount of SOC in the soils ied.The forest land use system reported maximum soil OC and available K status region.In general, the values of pH and CaCO3 increased with depth, attributed t leaching and accumulation of basic cations in the sub-surface layers of the soil profi

Conclusions
The distribution of physico-chemical properties varies under different soil orders and land use systems at different depths in profiles of the soil samples examined.The results of the research indicated that soil pH, CaCO 3 and available P content in soil significantly varied among different soil orders, whereas EC, SOC and available K were not significantly influenced by soil orders.The highest soil pH, EC and CaCO 3 values were reported in inceptisols, followed by alfisols and the lowest values were found in entisols.The pattern of SOC variation in soil orders is as follows: entisol > inceptisols > alfisols.Inceptisols had the highest available P, followed by entisols and alfisols.Entisols had the highest amount of available K, followed by alfisols and inceptisols.In all soil orders, the distribution of the physico-chemical properties of soil under different land use systems did not follow any specific trend.When compared to other land uses, the agriculture system had the highest values of soil pH, EC, CaCO 3 and the lowest amount of SOC in the soils studied.The forest land use system reported maximum soil OC and available K status in the region.In general, the values of pH and CaCO 3 increased with depth, attributed to the leaching and accumulation of basic cations in the sub-surface layers of the soil profile.As the depth of the soil increases, soil properties such as OC, P and K levels gradually decrease in all the land uses.The constant decline in content and addition of organic matter in soils may be the cause of the decrease in concentration with depth.The principal component analysis (PCA) reported that soil pH and OC were the most variable soil parameters, which influence the availability of other physico-chemical properties under different soil orders and land use systems.The study s findings indicated that in agriculture systems, SOC levels are lower because soil organic matter is in relatively poor condition, and that in order to preserve soil fertility and productivity, proper management strategies are considered necessary.The soils contained a large amount of K in all the soils examined.Therefore, it is advised that farmers apply the necessary amount of potassium to replace that which is lost by the various crops in the system.The physico-chemical properties of soil play an important role in the availability of nutrients in the soil, which is a major concern for efficient crop production.Hence, it may be concluded that the variation in the distribution of soil physico-chemical properties in different soil orders is highly affected by the land use systems of the regions.However, there is a need to investigate the other physical, chemical and biological properties of soil profiles for a better understanding of the effects of existing land use systems in different soil orders of the Majha region.The results of the current study might be of great importance to researchers, state agricultural officers and farmers in implementing nutrient management strategies in different land use systems for sustainable crop production.

Sustainability 2023 , 17 Figure 1 .
Figure 1.Locations for soil sampling in the Majha region of Punjab.

Figure 1 . 17 Figure 2 .
Figure 1.Locations for soil sampling in the Majha region of Punjab.Sustainability 2023, 15, x FOR PEER REVIEW 5 of 17

Figure 2 .
Figure 2. Classification of land use systems and soil orders in the Majha region of Punjab, India.

Figure 3 .
Figure 3. Variation in soil pH under different soil orders and land use systems in the Majha region of Punjab.

Figure 3 .
Figure 3. Variation in soil pH under different soil orders and land use systems in the Majha region of Punjab.

Figure 4 .
Figure 4. Variation in soil electrical conductivity (dS m −1 ) under different soil orders and land use systems in the Majha region of Punjab.

Figure 4 .
Figure 4. Variation in soil electrical conductivity (dS m −1 ) under different soil orders and land use systems in the Majha region of Punjab.

Figure 5 .
Figure 5. Variation of soil organic carbon (%) content under different soil orders and land use systems in the Majha region of Punjab.

Figure 5 .
Figure 5. Variation of soil organic carbon (%) content under different soil orders and land use systems in the Majha region of Punjab.

Sustainability 2023 , 17 Figure 6 .
Figure 6.Variation in calcium carbonate content (%) under different soil orders and land use systems in the Majha region of Punjab.

Figure 6 .
Figure 6.Variation in calcium carbonate content (%) under different soil orders and land use systems in the Majha region of Punjab.

Sustainability 2023 , 17 Figure 7 .
Figure 7. Variation in soil available phosphorus (kg ha −1 ) under different soil orders and land use systems in the Majha region of Punjab.

Figure 7 .
Figure 7. Variation in soil available phosphorus (kg ha −1 ) under different soil orders and land use systems in the Majha region of Punjab.

Sustainability 2023 ,
15, x FOR PEER 12 of 17Thus, multiple factors such as depth, soil amendments and climatic conditions under different land-use systems all contributed to increased soil available K content.

Figure 8 .
Figure 8. Variation in soil available potassium (kg ha −1 ) under different soil orders and land use systems in the Majha region of Punjab.
03) and CaCO3 (p < 0.01 and r = −0.38).The available soil K was positively associated with EC (r = 0.26), OC (p < 0.01 and 0.46), but significantly negatively correlated with soil pH (p < 0.01 and r = −0.46)and CaCO3 (p < 0.05 and r = −0.35)content in soil.Depth-wise distribution of CaCO3 was significantly negatively associated to OC (p < 0.05 and r = −0.35)but positively correlated with soil pH (p < 0.01 and r = 0.52) and EC (r = 0.01).Soil OC was found to be significantly negatively correlated with pH (p < 0.05 and r = −0.31),but showed a positive correlation with EC (p < 0.05 and r = 0.38).The distribution of EC with pH showed a non-significant positive correlation (r =0.15) in the soils studied.

Figure 8 .
Figure 8. Variation in soil available potassium (kg ha −1 ) under different soil orders and land use systems in the Majha region of Punjab.
26), OC (p < 0.01 and 0.46), but significantly negatively correlated with soil pH (p < 0.01 and r = −0.46)and CaCO 3 (p < 0.05 and r = −0.35)content in soil.Depth-wise distribution of CaCO 3 was significantly negatively associated to OC (p < 0.05 and r = −0.35)but positively correlated with soil pH (p < 0.01 and r = 0.52) and EC (r = 0.01).Soil OC was found to be significantly negatively correlated with pH (p < 0.05 and r = −0.31),but showed a positive correlation with EC (p < 0.05 and r = 0.38).The distribution of EC with pH showed a non-significant positive correlation (r =0.15) in the soils studied.

Figure 9 .
Figure 9. Principal component analysis for soil physico-chemical properties from different lan systems under distinct soil orders in the Majha region of Punjab.

Figure 9 .
Figure 9. Principal component analysis for soil physico-chemical properties from different land use systems under distinct soil orders in the Majha region of Punjab.

Table 1 .
Linear correlation (Pearson) coefficient between soil physico-chemical properties of different land use systems under major soil orders in the Majha region of Punjab.

Table 1 .
Linear correlation (Pearson) coefficient between soil physico-chemical properties of different land use systems under major soil orders in the Majha region of Punjab.*,*indicatecorrelation is significant at 1 and 5% significance level.3.8.Principal Component Analysis for the Determination of Most Variable Factor among Several Soil Parameters under Different Soil Orders and Land Use Systems of Majha RegionPrincipal component analysis (PCA) exhibited the general sensitivity pattern of soil parameters and most variable parameters based on the factor loadings from each principal component (PC) (Table2 and Figure *