An Assessment of Soil Quality and Agricultural Production Status in the Alluvial Soil Region: A Case Study in Koch Bihar District, West Bengal, India Current World

Soil refers to the upper layer of the Earth’s surface, which is made up of a mixture of organic residues, clay, and rock particles, and that's where plants grow. The soil quality is the environmental aspect that is most significant in agricultural activities, as well as for the concern of the safety of agricultural produces. At present-day, soil quality assessment becomes the most important issue because of the raising food security awareness. This study was assessed in Koch Bihar district, West Bengal, India, to quantify soil quality using the weighted Soil Quality Index (SQI) approach. For this assessment, the soil-related data were collected from the Soil Health Card (SHC) of the Agricultural Development Offices in all blocks of Koch Bihar and then analyzed using the principal component analysis (PCA) and expert opinion (EO) method. The weights of the selected soil quality indicators were determined using the integrated Fuzzy-AHP model. According to this method, 12 indicators, i.e., soil pH, soil electrical conductivity (EC), soil organic carbon content, nitrogen content, potassium content, phosphors content, soil texture, the groundwater level in the winter season, cropping intensity, and drainage frequency were considered based on literature review. The SQI assessment was done accordingly for each spatial unit (Block), and the spatial variability of the soil quality map was produced by GIS spatial analysis module. The findings may aid in promoting non-harmful produce production, the provision of scientific data for agricultural structure adjustment, and the maintenance of agricultural sustainability.


Introduction
The dynamic quality of soil properties affected by good management is known to as soil quality. The soil quality concept was familiarizing for correct grouping and allocation of agricultural inputs materials'. 1,2 The ability of the soil to work within a regulated or native ecosystem limits in order to preserve or improve the productivity of the ecosystem's plants and animals is referred to as soil quality. 3,4,5 As a result, the most significant environmental parameter in the ecosystem is soil quality, which helps the yields and safety of agricultural production. 6 The worldwide environmental excellence and food supply were maintained by the soil, and this awareness excites the evaluation and measurement of soil quality. 7,8,9 The soil quality assessment is very complex as the soil properties may be composed of several features, and there is no set process for selecting soil quality indicators. 8 It is frequently difficult to choose how many and which are considered the soil quality indicators for the assessment. 8,10,11 To assess the Soil Quality Index (SQI), it is necessary to secure the indicators and determine their respective weight value. Indicators such as soil pH, soil electrical conductivity (EC), soil organic carbon content, nitrogen content, potassium content, phosphors content, sand, silt, clay, groundwater level in the winter season, cropping intensity and drainage frequency were used to determining the SQI. Some statistical processes such as Principal Component Analysis (PCA), Fuzzy set Theory, Analytic Hierarchy Process (AHP), and R (Regression) were expanding to find these indicators selection and weight assignment process. 12,13,14 Different types of soils can be found in India, and their creation is influenced by various elements, including elevation, climate, excessive rainfall, and so on. Different types of soil are found in multiple parts of the country. Previously, the soil was classed according to its fertility. The soil was either 'Urvara', which meant fertile, or 'Usara', which meant non-fertile or sterile; however, today, many features are taken into account, and the soil type is classed based on texture, colour, or moisture content, and chemical composition. The Government of India established the Soil Survey of India in 1956 as an entity to study soil and its features. Mountain soil, Laterite soil, alluvial soil, black soil, desert soil, saline soil, red soil and peat soil are some of the soil types found in the country. 39,40 On the basis of chemical composition (mainly p H ) there are three types of soil such as, Neutral soil (p H = 7.0), Acidic soil (p H < 7.0) and Alkaline soil (p H > 7.0). 39,40 As alluvial soil is most important for crop production and maintaining food supply over the country, hence the assessment of this type of soil is very much necessary in the present-day context. Alluvial soils are generated by the sediments carried downstream by the rivers. In the Terai and Teesta alluvial zone, most of the rivers originated in the Himalayas and transported a large amount of silt. Particles such as silt, sand, and clay make up the soil, and it has enough phosphorus, potassium, and lime. There are two types of alluvial soil: (a) old alluvium known as 'bhangar', and (b) young alluvium known as 'khadar'. It's (khadar) also the most important type of soil in the country, accounting for roughly 40% of the total land area. It stretches from Punjab to West Bengal and Assam in the northern plains physiographic unit. In peninsular India, it can also be found in the deltas of rivers like the Krishna, Godavari, Kaveri, and Mahanadi. Alluvial soil is light grey in colour and quite fertile. Wheat, rice, maize, sugarcane, legumes, and oilseeds are the most often grown crops in this soil region.
The Terai and Teesta alluvial zones cover all twelve blocks in the Koch Bihar district, according to the categorization. The creation of effective land-use management systems requires an understanding of the distribution patterns of soil property in spatial aspects as well as in temporal aspects. 14,15,16 Several investigations are needed to create an acceptable and practical SQI protocol for the wide and harsh tropical monsoon area soils; the current work is an attempt in that regard. It was carried out in a tropical monsoon region soil to assess the soil quality index. Hence, for this study, the authors took Koch Bihar district in India as an example, using Soil Health Card (SHC) data for the assessment of the spatial variation of the SQI. The SHC indices, those uses for calculation SQI were p H , EC, soil organic carbon content, nitrogen content, potassium content, phosphors content. This evaluation of the SQI was carried out by the integrated Fuzzy-AHP soil quality index model for this study. This SQI assessment can make fertilization more efficient and assist underdeveloped regions and backward rural areas, heavily reliant on the soil for crop production. 14,17,18,19

Materials And Methods Study Area
Koch Bihar district belongs to the Terai-Tista Alluvial Zone. Physiographically, this district has two parts; one is an active alluvial plain, another is a recent alluvial plain, and both are formed by the Tista River. In order, the soil of this study area falls under Entisol, and there is the presence of 4 soil series, i.e., Lotafela, Matiarkuthi, Balarampur, and Rajpur series. The elevation of the study area varying between 30-50 meters above mean sea level and the main rivers run across a level terrain with a small southeastern slope. (Tista, Torsa, Kaljani, Gadadhar, Raidak, Sankosh) are flow. The climatic condition of this district is humid to peri-humid with an average of 250-300 mm annual rainfall. The maximum temperature in the summer and winter is roughly 33°C in May, with the lowest temperature being 7°C in January. Geographically, this district is located (Fig. 1) in the northeastern part of West Bengal. The latitudinal and longitudinal extension of the study is 25°57ʹ57ʹʹN to 26°32ʹ58ʹʹN and 88°45ʹ28ʹʹE to 89°51ʹ50ʹʹE, respectively. The district's overall geographical area is 338700 ha, with 1146 ha of net sown area and 246491 ha of gross cropped land. The district is economically agricultural-based, and most inhabitants depend on agricultural and allied activities. According to Census 2011, the district has 34.74% agricultural labourers and 32.34% cultivators (Koch Bihar District profile, http://www.coochbehar. nic.in/HTMfiles/dist_profile.html).

Data sources
The calculation of soil quality index was conducted utilizing a variety of indicators such as soil electrical conductivity (EC), pH, soil organic carbon content, nitrogen content, potassium content, phosphors content, the textural classification (amount of sand, silt, clay), the groundwater level in the winter season, cropping intensity and drainage frequency. The data about the parameters such as 'EC', soil pH, 'OC', 'N' content, 'P' content, and 'K' content were collected from the Soil Health Card (SHC) of existing Agricultural Development Offices in all blocks of Koch Bihar district.
The textural classification (amount of sand, silt, clay) of soil were collected from SATSA report Koch Bihar Unit; data about the groundwater level in the winter season were collected from CGWB, Government of India; cropping intensity data were collected from Annual agricultural plan report of Koch Bihar district (2017), and drainage frequency was calculated from Topographical map (No. 125). The available N, P, and K were measured in kilograms per hectare (kg/ha). In contrast, the soil textural elements such as sand, silt, clay content and organic carbon content were measured in percentages (Table 1).

Evaluation of Soil Quality Index
Four steps (Fig. 2) were taken to calculate SQI: (i) identify the management aim, (ii) selecting indicators, (iii) evaluating the chosen indicators, and (iv) computing SQI. Soil functions include biological diversity preservation, water and solute flow regulation and partitioning, buffering and detoxification of organic and inorganic contaminants, nutrient storage, and cycling the plant support. 4,20,21 However, it is widely acknowledged that crop output is stagnating in any region due to inadequate management, largely rain-fed conditions, and a scarcity of resources. 22 The production activity is given top priority among the soil functions in this study, despite the fact that crop cultivation in any region is influenced by the edaphic and non-edaphic parameters of the soil. 23

Defining Management Goal
The main goal is to create a Soil Quality Index (SQI). that spatially varies in the study area and manages crop cultivation. The SQI indicators should be selected based on the important soil properties for cultivation. Individual management goals are common, but these should be socioeconomic based. The farm management is very crucial because land degradation, soil pollution are increase by excessive use of pesticides.

Indicator Selection
The principal component analysis (PCA) and expert opinion (EO) were chosen to determine soil quality parameters.

Principal Component Analysis
Using SPSS (version 20.3), PCA was used to reduce the number of soil quality parameters and derive indices for the soil quality index. The target of PCA was to decrease the number of variables in data while avoiding information loss. 24 PCs (Principal component) with higher eigenvalue were determined to be the prime choices for understanding the variability. 25 PCs which have the eigenvalue of greater than one or equal to one were adopted because, the eigen value less than one accounts for less variation than a single variable performs. 26 Thus for soil quality index indicators, highly valued variables were chosen from each PCA analysis. The descriptive statistics and PCA analysis for these soil quality indicators are shown in Table 2 and Table 3, respectively.

Expert Opinion
On the other hand, experts believe that characteristics of the study area location, such as climate, rainfall, and associated pedogenic processes that affect soil health, should be taken into account when selecting variables as indicators. Furthermore, it is critical that the indicators were chosen accurately reflect the soil's complexity and function. 21 In conclusion, soil quality parameters were chosen depending on available data and literature about the soil of the study area.

Scoring of Indicators
Using the linear scoring approach, selected indicators from a minimum data set of soil were evaluated into dimensionless values from 0 to 1. 21,27 The positive or negative relation of a soil parameter with soil quality index was evaluated in ascending or descending order. For this, in the case of positive relation indicators, each value was divided by the highest value, and in the case of negative relation indicators, the lowest value is divided by each particular indicators value. In this study, the Fuzzyset technique was employed to score the indicators, and AHP was employed to weight the indicators.

Fuzzy Set Theory
Fuzzy mathematics, fuzzy measurements, fuzzy integrals, and other concepts are all part of the Fuzzy Set Theory. 28 Fuzzy logic is one feature of fuzzy mathematics, and it is in this study for scoring the indicators. In classical set theory, a set's identity indicates a right or wrong statement represented as 0 or 1. On the other hand, a Fuzzy set theory allows participants to express themselves on a continuous scale spanning from 0 to 1. 29,30 In this context, the question arise about which value is better? It is depends on the relationship between the particular indicator and the SQI. If the relationship is positive then the value of 1 is better and if the relationship is negative then the value of 0 is better. The fuzzified or normalized value of soil quality indicators is described in Table 4.  31 On each basis, a nine-point scale is utilised to score each criterion's relative preferences. 32 The weightage value of each criterion has been assigned according to Saaty's relative importance scale which has been shown in Table-5. The calculation part of the AHP method has been evaluated sequentially with four respective processes to verify the result is consistent. The processes which are involved in the AHP method are the pairwise comparison matrix process, weight normalization tabulation, weights estimation, and the last stage of the method to check consistency. 31,32 The indicators (i.e., soil electrical conductivity (EC), pH, organic carbon (OC) content, nitrogen (N), phosphorus (P) and potassium (K) content, the textural classification (amount of sand, silt, clay), the groundwater level in the winter season, cropping intensity, and drainage frequency), and the factor weights are summarised in the pair-wise comparison matrix (Table 6). Moderate to strong importance 5 Equally preferred 6 Strong to very strong importance 7 Very strong importance 8 Very strong importance 9 Extreme importance  The consistency ratio (CR) is the proportion of the consistency index to the resulting consistency index. 14,32 If the CR is less than or equal to 0.1, then it is allowed; however, if the CR is greater than 0.1, the subjective judgment must be changed. Where, RI = the resultant consistency index's average, which depends on CI = index of consistency, is expressed as: CI=(λmax-n)/(n-1) ... (2) Where, λmax = the largest individual value in the matrix is max, and the order of the matrix is n -1.
The CR is a scale from 0 to 1, with 1 representing the likelihood of a randomly generated matrix and CR less than 0.10 suggesting a good amount of homogeneity. 32 The weights that resulted have an acceptable CR, as indicated in Table 6. The weights should sum up to 1.0 according to the linear weighted combination computation rule. 14 ∑ n j=1 wj=1 ...(3)

SQI Calculation
To calculate SQI, we used data about soil profiles from the Soil Health Card (SHC). Soil parameters at the surface (0-15 cm) are more dynamic, indicating that the soil is in a state of flux. 25 SQI for surface depths was calculated through a weighted additive index by the Fuzzy-AHP model.

Fuzzy-AHP Model
The Fuzzy-AHP method is made up of the AHP weight of the factor and the Fuzzy values for each parameter. The final soil quality map was made using the sum of the individual soil quality parameters' values generated by the Fuzzy-AHP method. The following formula is described the method: ... (4) Where, μ= the participant function for each parameter, W = the weight of each parameter. Weights were assigned to the indicators chosen by PCA and EO methods depending on the relative value of each indication determined by the Fuzzy AHP approach. The weighting factor was chosen so that the sum of all elements equals one of the soil functions. Then the weighted variables were added to derive the Soil Quality Index (SQI) value for all the blocks of the entire district. The weighted layers were combined in GIS to create the final soil quality zonation maps for the study area. Table 7 shows the statistical values of the soil quality index indicators developed using the Fuzzy-AHP model.

Validation of Soil Quality Index (SQI)
By computing correlation coefficients, the SQIs predicted from the Fizzy-AHP model were validated against Paddy, Jute, Potato, Maize, and Mustard yield data. The SQI values were also compared and found to be related to the above-mentioned crop's block-wise yield rate.

Assessment of Spatial Variation of Soil Quality
The assessment of the spatial variation of soil quality was done by using the spatial autocorrelation among the soil quality index value of 12 blocks in the study area. According to Salima and de Bellefon, the positive or negative correlation of a variable with itself due to the spatial location of the observation, is known as spatial autocorrelation. 34 The most familiar Moran's Index is used to check the spatial autocorrelation between any variable and spatial location. The Moran's index value ranged from positive 1 to negative 1. The negative 1 value indicates perfectly spatial heteroscedasticity, i.e., the value is spatially dispersed, and the positive 1 value indicates perfectly homoscedasticity, i.e., the value is spatially clustered. The Moran's index calculation equation is present below: ... (5) Where, N= number of the spatial unit, Σi Σj wij= summation of neighborhood relation through the matrix, Σ i (X i −X̅ )2 = summation of deviation, Σ i Σ j wij (X i −X̅ )(X j −X̅ ) = summation of neighborhood deviation through the matrix.

Spatial Distribution of Soil Quality Indicators
The indicators used to calculate the soil quality index are p H , EC, soil organic carbon content, nitrogen content, potassium content, phosphors content are shown in Fig. 3 and the textural classification (amount of sand, silt, clay), the groundwater level in the winter season, cropping intensity, drainage frequency is showed in Fig. 4. Sometime fertilizers are used widely, resulting in an unbalanced distribution of these inputs, and it frequently occurs due to a lack of information on the variability of soil quality indicators in different blocks of the study area. Some locations receive more than what is required, while others receive less. So, for the spatial variation of soil quality index, the spatial variation of its' indicators assessment is necessary for the study area, and it is discussed as follows:

Fig. 4: Soil quality indicators i.e. (g) Soil sand content (h) Soil silt content (i) Soil clay content (j) Ground water level (k) Cropping intensity (l) Drainage frequency
The soil electrical conductivity (EC) values were found to be distributed in a range of 0.16 ds/m to 0.39 ds/m. The EC distribution zones in the study area are classified into four classes based on natural break classification. The results of the classification revealed that the soil EC values of 0.16 to 0.22 ds/m belong in the north eastern part of the district, whereas the highest class, i.e., 0.31 to 0.39 ds/m, has belonged in the north western part of the district. The pH of the soil was generally basic and within the crucial ranges ranging from 5.33 in larger parts of the north western and south western part to 6.55 in distribution northeastern of the plain of the district. The organic carbon (OC) content of the soil was observed to be distributed in a range of 0.49 to 0.84 percent. In some areas, such as the northeast ward and lower in all other parts of the study area the organic carbon was higher. In the majority of the study area, total nitrogen content in soil was low, ranging from 181.08 kg/h to 307.95 kg/h. N total had a similar spatial distribution pattern as organic carbon, with more concentrations in the northeastern section and less concentrations in the rest of the area. The phosphorous content of the soil ranged between 60 kg/h to 174.30 kg/h in the study area, and this was beyond the critical value. The spatial distribution pattern of P was highly concentrated in the northeastern part, which is 135kg/h to 174 kg/h, and the lowest concentration of P was distributed in the middle and western part of the study area, which is approximately 60 kg/h to 91 kg/h. The potassium content of the soil ranged between 65.21kg/h to 619.88 kg/h in the study area. Higher concentrations of accessible K were detected in the northwestern and southern parts of the study area, while lower concentrations were found mostly in the western, northeastern, and southeastern patches in that region. Soil texture refers to the proportion of sand, silt, and clay-sized particles that make up the mineral component of the soil. Available sand ranged between 16.80% to 55%, and higher concentration sand are distributed in the western and northwestern part, whereas lower concentration sand is distributed in the southern and northeastern part of the study area. Available silt ranged between 25.10 % to 59.99 %, the distribution pattern of silt is approximately the same as sand distribution. Available clay content is ranged between 13.70 % to 23.20 %, and higher concentration clay is distributed in the southern and northeastern parts, whereas lower concentration clay is distributed in the central portion of the study area. The relation between soil quality and groundwater level in winter is reciprocal, i.e., if the groundwater DAS & PAUL, Curr. World Environ., Vol. 17(1) 268-283 (2022) level (mbsl) increases, then the soil quality decrease 35 The range of groundwater levels in the winter season in the study area is 2.37 mbsl to 4.25. mbsl in the study area. Groundwater level below sea level is higher in the northeastern portion and lower through the western and southwestern parts of the study area. Generally, where the cropping intensity is high, the quality of soil is better.
T h e c r o p p i n g i n t e n s i t y r a n g e d b e t w e e n 149.78 to 256.09, whereas it is high in the western and southeastern parts and low in the northwestern part of this district. Drainage water management is expected to result in a minor increase in soil organic matter content, which would be beneficial to soil quality. 36 Drainage frequency is high in northern and northeastern parts and low in the southwestern part of the study area. Throughout the district, the data about soil quality indicators of 12 blocks are revealed in Table 1.

Spatial Variability of Soil Quality Status
The spatial variability of soil quality in the district was determined with the help of Fuzzy-AHP model and Moran's index. Fuzzy-AHP model was used to determine soil quality indicators weights and Moran's index was use to analysis the spatial variation of soil quality. According to this method, 12 indicators such as electrical conductivity (EC), soil pH, organic carbon (OC) content, available nitrogen (N) content, phosphorus (P) content, potassium (K) content, amount of sand, silt, clay, the groundwater level in the winter season (GWL), cropping intensity (CI) and drainage frequency (Df) were selected. 8    The soil quality indices predicted by the Fuzzy-AHP in the study region are shown on a scale of 0 to 1, with soil quality indices ranging from 0.36 to 0.59 in the study area it was covered all the blocks of the district, as displayed in Fig. 5 (a). The soil quality classification was revealed in a few isolated places, and poor soil quality was observed in the western, middle, and southeasterly parts of the study area and moderate towards the southwest, northern, and eastern portion of the plain. Among the 12 blocks, the soil quality index is spatially varied from one block to another such as, the SQI of Sitai  Table 9. The Moran's index of 12 spatial units (block) is -0.22 (Fig. 5 b), which means the distribution pattern of SQI values is heterogeneous. The soil quality is spatially varied from one spatial unit to another through the study area. . The agricultural productions are varied with the variation of soil quality index in the study area. The cultivation preference orders from the selected major produced crops such as paddy, jute, maize, potato, wheat, and mustard are analyzed by the correlation between crop yield and soil quality index in each spatial unit in this district.
The correlation values between crop yield and SQI ranged from 0.0009 to 0.2028, according to the analysis (Fig. 6). In the study area, paddies (0.2028) have the highest preference for cultivation, and it is followed by jute (0.1172), wheat (0.0222), maize (0.0133), potato (0.0035), and mustard (0.0009) which is shown in Table 10.

Conclusion
There are two reasons why soil quality is crucial. First, improper soil use can harm both the soil and the ecosystem; as a result, land management must be tailored to the soil's capabilities. Second, we must have a baseline understanding of soil quality in order to detect changes as they occur. So, the main objective of soil quality assessment is to safeguard and increase agricultural production, water quality, and shelters for all animals, including humans, over the long run. Soil quality research has recently concentrated on the connections between management techniques and systems, visible soil properties, soil activities, and soil utility efficiency. When selecting the relevant soil qualities to include in an index, site-specific and user oriented management goals must be considered, as well as focusing on crop productivity, consider sustainability. These measures would be beneficial in determining the soil's weakness and identifying how better management could improve its resilience.
This work made a methodological addition by developing an integrated Fuzzy-AHP with GIS technique based on the degree of certainty for all the soil quality criteria. The foundations of the PCA, AHP, and Fuzzy-AHP approaches were introduced for this aim. The soil quality indices for paddy, jute, maize, potato wheat, and mustard production ranged from low to moderate, according to the Fuzzy-AHP analysis in this district which belongs to the alluvial soil zone. AHP and Fuzzy-AHP both proposed identical decision-making alternatives, which is that the Koch Bihar district is somewhat moderately suitable for cash crops. The pH and electrical conductivity (EC) distributions were within the threshold value for the below and slightly over the required accordingly in this study region. Furthermore, the most relevant parameters on soil quality were discovered to be pH, EC, OC, and N content concentration, and soil quality indices predicted by Fuzzy-AHP models were correlated with the selected crops' yield. The crops sequence of paddy, jute, wheat, maize, potato, and mustard is most desirable in the studied district; according to the Fuzzy-AHP analysis, when compared to these different types of crops production, the validation employing linear regression analyses was found to be very simple and transparent for decisionmakers to understand. The results showed that the Fuzzy-AHP decision problem solution was more accurate and less sensitive to the criteria weights than the solution achieved using other methods. So, a prospective crop for cultivation was identified in this research, and the Smallholder farm productivity in the area is gradually increased. It also showed how soil quality indexing can be used to detect changes in soil and crop management that are needed for more effective and long-term use of land resources.