Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Pre-monsoon spatial distribution of available micronutrients and sulphur in surface soils and their management zones in Indian Indo-Gangetic Plain

  • A. K. Shukla,

    Roles Conceptualization, Investigation, Writing – review & editing

    Affiliation ICAR- Indian Institute of Soil Science, Bhopal, Madhya Pradesh, India

  • S. K. Behera ,

    Roles Conceptualization, Data curation, Investigation, Writing – original draft, Writing – review & editing

    sanjibkumarbehera123@gmail.com

    Affiliation ICAR- Indian Institute of Soil Science, Bhopal, Madhya Pradesh, India

  • V. K. Singh,

    Roles Data curation, Formal analysis, Investigation, Validation

    Affiliation ICAR- Indian Agricultural Research Institute, New Delhi, India

  • C. Prakash,

    Roles Data curation, Formal analysis, Software, Visualization

    Affiliation ICAR- Indian Institute of Soil Science, Bhopal, Madhya Pradesh, India

  • A. K. Sachan,

    Roles Formal analysis, Investigation, Supervision

    Affiliation Chandra Sekhar Azad University of Agriculture & Technology, Kanpur, Uttar Pradesh, India

  • S. S. Dhaliwal,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Validation

    Affiliation Punjab Agricultural University, Ludhiana, Punjab, India

  • P. C. Srivastava,

    Roles Investigation, Methodology, Visualization, Writing – review & editing

    Affiliation Govind Ballabh Pant University of Agriculture and Technology, Pantnagar, Udham Singh Nagar, Uttarakhand, India

  • S. P. Pachauri,

    Roles Data curation, Investigation, Methodology

    Affiliation Govind Ballabh Pant University of Agriculture and Technology, Pantnagar, Udham Singh Nagar, Uttarakhand, India

  • A. Tripathi,

    Roles Methodology, Software, Validation

    Affiliation ICAR- Indian Institute of Soil Science, Bhopal, Madhya Pradesh, India

  • J. Pathak,

    Roles Formal analysis, Investigation, Supervision

    Affiliation Banda University of Agriculture and Technology, Banda, Uttar Pradesh, India

  • A. K. Nayak,

    Roles Conceptualization, Supervision, Validation, Visualization, Writing – review & editing

    Affiliation ICAR-National Rice Research Institute, Cuttack, Odisha, India

  • A. Kumar,

    Roles Data curation, Formal analysis, Investigation

    Affiliation ICAR-Indian Agricultural Statistics Research Institute, New Delhi, India

  • R. Tripathi,

    Roles Formal analysis, Software, Writing – original draft

    Affiliation ICAR-National Rice Research Institute, Cuttack, Odisha, India

  • B. S. Dwivedi,

    Roles Conceptualization, Project administration, Writing – review & editing

    Affiliation ICAR- Indian Agricultural Research Institute, New Delhi, India

  • S. P. Datta,

    Roles Formal analysis, Investigation, Resources, Validation, Writing – review & editing

    Affiliation ICAR- Indian Agricultural Research Institute, New Delhi, India

  • M. C. Meena,

    Roles Data curation, Formal analysis, Investigation, Resources

    Affiliation ICAR- Indian Agricultural Research Institute, New Delhi, India

  • S. Das,

    Roles Investigation, Resources, Validation, Visualization

    Affiliation International Zinc Association, New Delhi, India

  •  [ ... ],
  • V. Trivedi

    Roles Formal analysis, Investigation, Methodology

    Affiliation ICAR- Indian Agricultural Research Institute, New Delhi, India

  • [ view all ]
  • [ view less ]

Abstract

The efficient (site-specific) management of soil nutrients is possible by understanding the spatial variability in distribution of phyto-available nutrients (here after called available nutrients) and identifying the soil management zones (MZs) of agricultural landscapes. There is need for delineating soil MZs of agricultural landscapes of the world for efficient management of soil nutrients in order to obtain sustainability in crop yield. The present study was, therefore, undertaken to understand the spatial distribution pattern of available micronutrients (zinc (Zn), boron (B), iron (Fe), manganese (Mn) and copper (Cu)), available sulphur (S), and soil properties (soil acidity (pH), electrical conductivity (EC) and organic carbon (SOC) content) in soils of intensively cultivated Indo-Gangetic Plain (IGP) of India and to delineate soil MZs for efficient management of soil nutrients. Totally, 55101 soil samples from 0–15 cm depth were obtained from 167 districts of IGP during 2014 to 2017 and were analysed for different soil parameters. Soil pH, EC and SOC content varied from 4.44 to 9.80, 0.02 to 2.13 dS m-1 and 0.10 to 1.99%, respectively. The concentration of available Zn, B, Fe, Mn, Cu and S varied from 0.01 to 3.27, 0.01 to 3.51, 0.19 to 55.7, 0.05 to 49.0, 0.01 to 5.29 and 1.01 to 108 mg kg-1, respectively. Geostatistical analysis resulted in varied distribution pattern of studied soil parameters with moderate to strong spatial dependence. The extent (% area) of nutrient deficiencies in IGP followed the order: S > Zn > B > Mn > Cu > Fe. Principal component analysis and fuzzy c-means clustering produced six distinctly different soil MZs of IGP for implementation of zone-specific soil nutrient management strategies for attaining sustainability in crop yield. The developed MZ maps could also be utilized for prioritization and rationalization of nutrients supply in IGP of India.

Introduction

The Indo-Gangetic Plain (IGP) of India, one of the intensively-cultivated agricultural landscapes of the world, plays a pivotal role in food production of the country. It occupies 52.01 m ha of land area and produces nearly 50% of total food grain production of the country [1]. It is characterized by availability of deep and fertile soils, favourable climatic conditions and sufficient water supply which sustain better agricultural productivity. The IGP played significant role in enhancing food grain production of the country during green revolution era due to growing of nutrient responsive and high yielding crop varieties and adoption of better crop management practices [2]. However, post-green revolution scenario witnessed the decline in factor productivity in IGP, predominantly because of receding ground water table and soil degradation especially due to secondary and micronutrient deficiencies [3]. The emerging deficiency of phyto-available (here after called available) secondary and micronutrients in different parts of the country including IGP is due to less or nil application of organics, over-dependence on straight fertilizers and imbalanced application of nutrients ignoring the replenishment requirement of mined nutrients [45].

There are reports of deficiencies of available micronutrients viz., zinc (Zn), boron (B), iron (Fe), copper (Cu) and manganese (Mn) [6] and available sulphur (S) [7] in various crops and soils of world. A recent analysis of Indian soils revealed an average deficiency level of 36.5% for Zn, 23.2% for B, 12.8% for Fe, 7.1% for Mn, 4.2% for Cu and 28.5% for S [8]. Availability of micronutrients and S in soils the result of combined influence of native soil nutrient status and important soil properties (pH, electrical conductivity (EC), and soil organic carbon (SOC)). The reduction in crop yield owing to deficit concentration of available micronutrients and S and crop responses to micronutrients and S application in soils of IGP of India have been reported [910].

Imbalanced addition of micronutrients and S to soil without the knowledge of their spatial distribution leads to unsustainable crop production. Therefore, proper understanding of spatial distribution of these nutrients and associated soil properties is needed for adoption of efficient (site-specific) soil nutrients management options through variable rate application in order to obtain sustainable crop production [1112]. Several researchers have investigated spatial distribution variability of available nutrients and associated soil properties in different soils of the world [1318] at field to regional scale using geostatistics [1920]. The spatial variations of micronutrients and S availability in soils of IGP is expected to be high primarily due to varied soil types, climatic conditions, crops and crop husbandries. The knowledge pertaining to spatial distribution of available S and micronutrients in soils of IGP is limited.

Geostatistics helps in effective evaluation of the spatial distribution variability of soil properties and available nutrients [21]. Geostatistical estimation predicts the values at unsampled location by establishing spatial correlation between sampled and estimated points and by reducing the estimation error and cost of investigation [22]. The effective way to address the spatial distribution variability of available nutrients in soil for site-specific nutrient management is by delineating soil management zones (MZs) of a particular area [23]. The delineation of MZs involves the techniques like principal component analysis (PCA) and fuzzy clustering, especially fuzzy c-means algorithm [24]. Soil MZs have been delineated for enhancing rice crop production in eastern part of India [25], for augmenting oil palm productivity in southern India [18], for improved nitrogen management in wheat in Argentina [26], and for enhancing corn productivity in Chile [27]. However, there is lack in information related to the zone-wise nutrients management in IGP of India. Keeping this in view, we conducted the present study (i) to understand the distribution variability of soil pH, EC, SOC, available S, and micronutrients (Zn, B, Fe, Mn, and Cu) and (2) to delineate potential soil MZs of IGP for site-specific S and micronutrient management, using geostatistical tools.

Materials and methods

Study area

Indo-Gangetic Plain (located at 21.583° to 32.467° N, 73.833° to 89.817° E) is an extensive fluvial plain spreading in West Bengal, Bihar, Uttar Pradesh, Haryana, and Punjab states of India. The eastern part of IGP is at lower elevation compared to western part. The area is having many rivers which are the sources of alluvium deposited in the plain. For the study purpose, surface soil samples from 0 to 15 cm depth were obtained from the farm lands of IGP (Fig 1). The study area experiences arid (western part), semi-arid (south-western part) and sub-humid (northern, southern and eastern part) climate with average annual rainfall of 300 to 600 mm in arid and semi-arid part and of 600 to 1000 mm in sub-humid part. Majority portion of IGP of India is having hyperthermic temperature regime and experiences mean average temperature of 40°C in summer months and of 10°C in winter months. Soils predominantly fit into Inceptisols, Entisols, Mollisols, Alfisols, and Aridisols orders [28] and having sandy to sandy loam texture [29]. Rice (Oryza sativa L.) and wheat (Triticum aestivum L.) are the prime crops grown in the area. Other crops include maize (Zea mays L.), cotton (Gossypium spp.), sugarcane (Saccharum officinarum L.), pulses, oil seeds and vegetables.

Soil sampling and analysis

Under the aegis of “All India Coordinated Research Project on Micro and Secondary Nutrients and Pollutant Elements in Soils and Plants” (AICRP-MSPE), a total of 55101 soil samples were collected by stratified random sampling from both rain-fed and irrigated agricultural soils of small, medium and large land holdings in 167 districts of IGP of India, during April to June months (pre-monsoon period) of 2014 to 2017. Prior to sample collection, necessary permissions were obtained from the owners of the land holdings. It was presumed that there were non-significant changes in available micronutrients and S in soil samples during sampling period. A hand-held stainless-steel auger was used for collection of soil samples. The geographical coordinates (longitude, latitude, and altitude) of each sampling point were recorded using a global positioning system. Composite samples were obtained to reduce the effect of sampling and to enhance the prediction accuracy [30]. Two to 3 subsamples for small holding, 5 to 6 subsamples for medium and 9 to 10 subsamples for large land holdings were collected to make a composite sample. The collected samples were air-dried under shade in a dust-free environment. Stones and debrises were removed before grinding of samples to pass through a sieve of 2 mm size. Processed samples were stored in polythene bottles for analysis.

For conducting analysis of soil samples in laboratory, prior approval was obtained from the Director, ICAR-Indian Institute of Soil Science, Bhopal, India. Soil pH and EC were determined in soil-water suspension (1: 2.5 weight/volume) [31] using pH meter (Make: Eutech, Model: pH 510) and conductivity meter (Make: Hanna, Model: HI 2300), respectively. Estimation of SOC was done by wet oxidation method [32]. Available S concentration was determined by extraction of soil samples with 0.15% calcium chloride (CaCl2) [33] and estimation using spectrophotometer (Make: Shimadzu, Model: UV-1800). Available B concentration was estimated through spectrophotometer after extraction with hot water [34]. Available Zn, Fe, Mn and Cu concentration in soil samples were estimated by extracting soils with diethylene-triamine-penta-acetic-acid (DTPA) extractant [35]. Measurement of the micronutrients concentration in the extract was carried out using atomic absorption spectrophotometer (Make: Varian, Model: AA240FS).

Descriptive statistics

The descriptive statistics viz., mean, minimum, maximum, standard deviation (SD), coefficient of variation (CV), skewness and kurtosis of studied soil parameters were obtained using SAS 9.2 software package [36]. Relationships among the soil parameters were visualized by carrying out Pearson’s correlation analysis.

Geostatistical analysis

Before geostatistical analysis, the data set was checked for normal distribution by Kolmogorov–Smirnov test. The data set followed the normal distribution. The trend analysis revealed no trend of the data. The analysis for semi-variogram (Eq 1) was obtained using ArcGIS software 10.4.1 for understanding the spatial structure of the soil parameters and to understand the interpolation function [22]. A semi-variogram plots the variance of spatially separated points of data and the separating distance (lag).

(1)

Where, γ (h) is semi-variogram at h distance interval; m(h) is the sample pair value at h distance interval; Z(Xi), Z(Xi+h) are the sample points separated h distance. Ordinary kriging (OK) interpolation technique was used for developing distribution maps [37]. In the kriging process, several semi-variogram models were tested for best-fitting. The exponential, stable, K-Bessel and circular models were found best-fitted. These models were selected through cross validation technique which measures the accuracy of the prediction. The root mean square error (RMSE) (Eq 2) which compares the estimated values from semi-variogram and the observed values was used for cross validation.

(2)

Where, z(xi, yi), z*(xi, yi) and n denote observed value, predicted value and number of observations, respectively.

Other semi-variogram parameters viz., nugget, sill and range were obtained for available S and micronutrients and associated soil properties. The ratio of nugget: sill is the criterion to describe the nature of spatial dependence [38]. The nugget: sill ratio values of ≤0.25, >0.25 to ≤0.75 and >0.75 denote strong, moderate, and weak nature of spatial dependence, respectively. The distances in which the values of soil parameters are inter-correlated are called ranges.

Principal component analysis was carried out using correlation analysis values as input using SPSS software (version 26.0). Principal components (PCs) with eigenvalues >0.90 were considered for developing the MZs in the present study. A bi-plot using altitude and studied soil parameters was drawn to indicate the effect of altitude on soil nutrients content and to examine relations among altitude and available S and micronutrients. Two to 8 clusters were obtained from the dataset through fuzzy c-means clustering using FUZME software [39]. The membership in each cluster was determined through an iterative process beginning with a random set of cluster means. Each of the observation was provided to the nearest of cluster means. The new mean for each cluster was re-estimated depending upon the distance of the observation from the cluster mean. The distance of the data points to the cluster centre was calculated using the euclidean distance. The optimum cluster number was determined by deriving fuzzy performance index (FPI) (extent of fuzziness) and normalized classification entropy (NCE) (degree of disorganization of specific class) (Eqs 3 and 4). The parameters c, n, μik, loga represents cluster number, observation number, fuzzy membership and natural logarithm, respectively. The differences in the mean values of soil parameters in different MZs were evaluated by variance analysis procedure.

(3)(4)

Result and discussion

Overall variability of soil properties and available nutrients

The IGP of India had acidic (4.04) to alkaline (9.80) soil pH (Table 1) with mean value of 7.41±0.94. About 0.8, 5.1, 8.9, 33.4, 45.7 and 6.1% samples had soil pH ranging ≤4.5, >4.5 to ≤5.5, >5.5 to ≤6.5, >6.5 to ≤7.5, >7.5 to ≤8.5 and >8.5, respectively (Fig 2). The variation in soil pH of the region is attributed to the various soil forming factors especially types of parent materials [40], and prevailing climatic parameters such as mean annual precipitation, temperature and evapotranspiration [41]. The soils of IGP region consisted of both younger and older alluvium [42] and had different units of soil-geomorphic units surrounded by different rivers and landforms. Acid soils were from the parent materials like granite, sandstone and shale. Whereas, limestone parent material gave rise to neutral to basic soils. Moreover, the proportion of mean annual precipitation and evaporation affects soil pH changes. Excess of evaporation compared to precipitation in IGP region resulted in accumulation of Ca2+ in soil surface resulting in higher soil pH. Whereas higher precipitation compared to evaporation results in leaching of Ca2+ and accumulation of Al3+ in surface soil resulting in lower soil pH. Soils were non-saline (EC 0.02 dS m-1) to slightly saline (EC 2.13 dSm-1) in nature with mean EC value of 0.35±0.26 dS m-1. About 37.4 and 45.2% soil samples had EC value of ≤0.25 and >0.25 to ≤0.50 dS m-1, respectively. Some parts of IGP had higher salt concentration in soil because of low precipitation, higher evaporation and irrigation of crops with saline water and poor soil-crop management. Soil organic carbon is an integral part of soil organic matter which influences soil physical, biological and chemical properties. The level of SOC in IGP region varied widely from 0.10% to 1.99% with mean value of 0.58±0.24% (Table 1). SOC content in 36.2, 42.7 and 14.1% samples were >0.25 to ≤0.50, >0.50 to ≤0.75 and >0.75 to ≤1.00%, respectively. Most parts of the IGP had low SOC levels because of imbalanced use of fertilizers, nil or lees crop residue addition, and adoption of excessive tillage practices [43].

thumbnail
Fig 2. Frequency distribution of soil samples (%) and area (%) under different ranges of available nutrients and associated soil properties of IGP, India.

https://doi.org/10.1371/journal.pone.0234053.g002

thumbnail
Table 1. Descriptive statistics of available nutrients and associated soil properties of IGP, India (n = 55101).

https://doi.org/10.1371/journal.pone.0234053.t001

The soils of IGP had mean concentration of 20.85±17.87 mg kg-1 for available S (varied from 1.01 to 107.52 mg kg-1), 1.06±0.61 mg kg-1 for available Zn (varied from 0.01 to 3.27 mg kg-1) and 0.97±0.68 mg kg-1 for available B (varied from 0.01 to 3.51 mg kg-1) (Table 1). The mean concentrations were 17.94±12.61 mg kg-1 for available Fe (varied from 0.19 to 55.71 mg kg-1), 12.47±10.36 mg kg-1 for available Mn (varied from 0.05 to 49.02 mg kg-1) and 1.74±1.19 mg kg-1 for available Cu (varied from 0.01 to 5.29 mg kg-1). The mean concentrations of available nutrients followed the order: available S > available Fe > available Mn > available Cu > available Zn > available B. According to the critical concentration ranges proposed by Shukla and Tiwari [44], the concentration of available S in 16.9% sample of IGP region was ≤7.5 mg kg-1, in 31.7% sample >7.5 to ≤15.0 and in 20.8% sample >15.0 to ≤22.5 mg kg-1 (Fig 2). The concentration of available micronutrients was as follows: ≤0.6 mg Zn kg-1 in 26.2% sample, ≤0.5 mg B kg-1 in 25.99% sample, ≤4.5 mg Fe kg-1 in 11.8% sample, ≤3.0 mg Mn kg-1 in 11.4% sample and ≤0.2 mg Cu kg-1 in 2.7% sample. In line with our observations, Singh et al. [45] reported available S concentration of 1.4 to 90.3 mg kg-1 in north-western IGP soils and west-Himalayan soils of India. Shukla et al [17] also recorded the concentration of available S ranging 0.55 to 130 mg kg-1 in soils of Shiwalik Himalayan region (SHR), India. The concentrations of available Zn, Fe, Mn and Cu were recorded ranging 0.10 to 8.00, 0.12 to 48.8, 0.53 to 26.6 and 0.10 to 7.97 mg kg-1 respectively, in soils of Trans-Gangetic Plain (TGP), India [12] and ranging 0.14 to 2.35, 0.90 to 28.5, 0.81 to 24.4 and 0.09 to 2.34 mg kg-1 respectively, in the soils of a Deccan Plateau Region (DPR), India [47]. The variations in the concentration of available B (hot water-soluble B) ranging 0.01 to 2.92 mg kg-1 in some acid soils [46] and ranging 0.19 to 6.11 mg kg-1 in some Indian Alfisols and Vertisols [47] were also reported. This variation in S and micronutrient concentration in IGP region is primarily because of nature of parent material, types of crops grown and amount and type of fertilizers applied in the region. The native supply of S and micronutrients is determined by weathering of the parent material [48] The organic acids released from the decomposition of crop residue as well as by the microbes facilitate the weathering of soil minerals and thus nutrient release. Several researchers have also reported variations in S and micronutrient concentration in soils of IGP due to diversity in crops grown and adoption of different of soil-crop management practices in the region [4951]. The CV values of studied soil parameters varied from 12.64% (pH) to 85.72% (available S) revealing their moderate variability [52]. This is in parallel to the observations of Bogunovic et al. [16] who recorded low, moderate, and high variability for pH, EC and SOC, respectively, in Rasa river valley soils of Croatia. Similarly, low variability for soil pH, and moderate variability for SOC and available Fe were reported by Tesfahunegn et al. [13] in northern Ethiopian soils. Wang et al. [53] recorded moderate variability for available Zn, Fe, Cu, and Mn in paddy soils of China.

Relationship among soil properties and available nutrients

There were significant correlations among the studied soil parameters (Table 2) although some correlations had low correlation coefficient values. Lower but significant correlation coefficient values are because of larger sample population. There was positive correlation (p ≤ 0.01) of soil pH with EC, available S, and B and negative correlation (p ≤ 0.01) with SOC, available Zn, Fe, Mn and Cu (Table 2) in the IGP region. Soil pH is the key variable influencing soil chemical, biological and physical properties [54]. It modifies soil chemical processes and thereby nutrient forms and their phyto-availability. The formation of free metallic cations and protonated anions is favoured by low pH condition [55, 56]. Whereas higher pH supports carbonate and hydroxyl complexes formation. Hence, the availability of cations increases with increase in soil acidity and availability of anions decreases with decrease in soil pH. The negative correlation of soil pH with available Zn, Fe, Mn and Cu explained the reduction in availability of these nutrients in IGP region with increase in soil pH. The negative relationship of soil pH with available cationic micronutrients was recorded by Wei et al., [57] in loess plateau soils of China and by Katyal et al. [58] and Shukla et al., [12] in some Indian soils. Soil EC was positively correlated (p ≤ 0.01) with all the studied soil parameters except SOC and available Mn. Corwin & Lesch, [59] reported EC as an indicator of phyto-available nutrients and soil salinity. The SOC content was negatively correlated (p ≤ 0.01) with available S and B and positively correlated (p ≤ 0.01) with available Zn, Fe, Mn and Cu (Table 2). The negative correlation of SOC with available S is in contradiction with the observation that soil organic fractions often serve as an important source of plant available S in soil [60, 61]. However, availability of S in soil depends upon the mineralization and immobilization process which in turn influenced by types of soil, crop managements and microbial activities [62]. The positive correlation of SOC with available cationic micronutrients indicates higher availability of these nutrients with increase in SOC content. Soil organic matter (in which SOC is a main component) releases chelated which enhances phyto-availability of these nutrients. Available S was positively correlated (p ≤ 0.01) with available Zn and B and negatively correlated (p ≤ 0.01) with available Fe, Mn and Cu. The correlations of available Zn with available B, Fe, Mn and Cu (p ≤ 0.01) were positive. There was negative correlation of available B with available Fe (p ≤ 0.01), Mn (p ≤ 0.05) and Cu (p ≤ 0.01). Available Fe was positively correlated (p ≤ 0.01) with available Mn and Cu. There was positive correlation of available Mn with available Cu (p ≤ 0.01). The positive correlations among the cationic micronutrients indicate that similar sets of factors influence distribution of these nutrients in IGP region. Behera and Shukla [63] also recorded positive correlations among the phyto-available cationic micronutrients in Indian acid soils.

thumbnail
Table 2. Correlation (Pearson’s) coefficients revealing relationship among available nutrients and associated soil properties of IGP, India (n = 55101).

https://doi.org/10.1371/journal.pone.0234053.t002

Spatial structure of soil properties and available nutrients

Geostatistical analysis revealed the best fitted exponential model for soil pH, available Fe, Mn and Cu, stable model for EC, SOC and available Zn, K-Bessel for available S and circular for available B with lower RMSE values (Table 3, Fig 3). The nugget value, indicating micro-variability, was higher for available S, Fe and Mn. This is attributed to inability of the sampling distance to capture spatial dependence. Nugget: sill ratio varied from 0.25 (SOC) to 0.61 (available Fe) with strong (SOC) to moderate (rest of soil parameters) spatial dependence (Table 3). The strong spatial dependence is because of intrinsic characters like soil mineralogy. On the other hand, moderate spatial dependence is due to combined influence of both intrinsic soil characters and extrinsic factors like effect of crops grown and fertilization.

thumbnail
Fig 3. Semi-variograms of available nutrients and associated soil properties of IGP, India.

https://doi.org/10.1371/journal.pone.0234053.g003

thumbnail
Table 3. Parameters of semi-variogram models for available nutrients and associated soil properties of IGP, India.

https://doi.org/10.1371/journal.pone.0234053.t003

The range value indicates the distance within which the samples are related spatially. The range values were 13115 m for pH, 36000 m for EC, available Fe and Mn, 30000 m for SOC, 33109 m for available S, 48000 m for available Zn, 54000 m for available B and 60000 m for available Cu. Soil parameters with higher range values are affected by natural as well as anthropogenic factors to a greater distance than the soil parameter with lower range value [64]. The variations in range values of pH, EC, SOC, available S, Zn, B, Fe, Mn and Cu is attributed to joint action of parent material, climatic condition and various land management practices. In line with the current study, other authors recorded ranges values of 132,000 m for pH, 65,000 m for EC, 59,000 m for SOC, 82,000 m for available S, 66 000 m for available Zn and 82,000 m for available Fe in SHR [17] and of 32, 490 m for available Zn, 61, 400 m for available Cu, 5, 370 m for available Mn and 140, 000 m for available Fe in intensively cultivated TGP, India [12]. The range values of the studied soil parameters in this study could be used as guide for designing future sampling strategies in similar regions. Though the ideal sampling interval needs to be less than half of semi-variogram range [30], it is hereby recommended to have shorter sampling distance than the range value obtained in this investigation for future studies characterizing spatial variability of pH, EC, SOC, available S, Zn, B, Fe, Mn and Cu in similar regions.

The distribution maps (generated by OK) displayed varied distribution patterns of soil properties, available S and micronutrients in IGP region of India (Fig 4). About 8.7, 6.8, 39.8 and 43.0% area were having soil pH ranging >4.5 to ≤5.5, >5.5 to ≤6.5, >6.5 to ≤7.5 and >7.5 to ≤8.5, respectively. The area having soil pH >4.5 to ≤6.5 is predominantly found in eastern part of IGP. Soils with pH values between >6.5 to ≤7.5 is optimum for plant growth and development. However, the areas having soil pH ≤6.5 (15.5%) and >7.5 (44.7%) need soil management practices for better crop production. There is need to adopt befitting liming technologies along with suitable crop cultivars in acid soils areas. Majority portion of IGP had EC ranging ≤0.25 (34.0% area) and >0.25 to ≤0.50 (52.0% area) dS m-1. Soil salinity was recorded in some pockets of IGP. However, soil sodicity (viz. presence of higher percentage of exchangeable sodium) is prevalent in northern and north-western part of IGP. The soil forming process of clay illuviation leads to development of soil sodicity [65]. The repeated cycles of drying and wetting of soils cause hydrolysis of feldspar causing release of alkalis. This results in calcium carbonate precipitation at high soil pH and generation of sub-soil acidity. Proper soil-crop management practices such as irrigation with good quality water, addition of soil amendments, and growing of salt tolerant crop cultivars with matching crop-nutrient management options need to be adopted in IGP areas having higher soil pH (> 7.5) [66]. The status of SOC was low (≤0.50%) in 34.3%, medium (>0.5 to ≤0.75%) in 47.3% and high (>0.75%) in 18.4% area of IGP. A large portion of area in western and southern IGP had low SOC status. The status of SOC was medium in most area of northern, central and eastern IGP. The variations in SOC status in different parts of IGP is the combined influence of soil types, climatic conditions, types of crops grown and adoption of different soil-crop management practices. Bhattacharyya et al. [67] recorded variations in SOC status in different parts of IGP. Efforts must be done to enhance SOC content in the areas of IGP having low and medium SOC status. It could be achieved by addition of organics, adoption of crop diversification practices by including pulse crops in cereal based cropping systems [68] and by following conservation agriculture practices [69].

thumbnail
Fig 4. Spatial distribution maps of available nutrients and associated soil properties of IGP, India.

https://doi.org/10.1371/journal.pone.0234053.g004

Available S was acutely deficient in 14.0% area, deficient in 30.2% area and latently deficient in 23.5% area. More area in eastern IGP had S deficiency. The variations in available S concentrations in different parts of IGP is because of differences in inherent S content in soils, soil properties, climatic conditions, soil-crop management practices and mismatch between S addition and S removal under various crops and/or cropping systems [70]. The management of soil S needs to be prioritized in the areas having S deficiency. About 2.5, 24.6 and 19.8% area had available Zn concentration ranging ≤0.3 (acute deficient), >0.3 to ≤0.6 (deficient) and >0.6 to ≤0.9 mg kg-1 (latent deficient), respectively. The deficiency Zn was prevalent in almost all parts of IGP except eastern and western parts. The concentration of available B was deficient (including acute deficient) (≤0.5 mg kg-1) in 22.5% area. About 22.2% area was latently deficient (>0.5 to ≤0.7 mg kg-1)) in available B. Boron deficiency was prevalent in eastern and some pockets of northern, southern and western part of IGP. Available Fe was deficient in 6.5% area and latently deficient in 7.2% area of IGP. Some pockets of southern and central part of IGP had Fe deficiency. About 5.1 and 13.2% area had available Mn concentration ranging >1.0 to ≤3.0 (deficient) and >3.0 to ≤5.0 mg kg-1 (latent deficient), respectively. The deficiency of Mn was found in some parts of western IGP. The concentration of available Cu was ≤0.2 mg kg-1 (acute deficient) in 1.1% area, >0.2 to ≤0.4 mg kg-1 (deficient) in 5.7% area and >0.4 to ≤0.6 mg kg-1 (latent deficient) in 8.5% area. Some pockets of western, northern and southern IGP had Cu deficiency. The spatial diversity in micronutrients concentration in different parts of IGP is because of physiographic variations. Moreover, cultivation of different crops, adoption of different soil-crop management practices and non-addition of micronutrient fertilizers add to the different spatial distribution scenario of micronutrients. The crops responses to S and micronutrient application varied to a greater extent in different parts of IGP (Table 4) [71]. This is due to varied concentration of phytoavailable S and micronutrients, nature crop grown and amount of nutrients applied. For example, the response of cereal crops ranged from 0.05 to 1.06 t ha-1 for S, 0.01 to 5.47 t ha-1 for Zn, 0.01 to 4.40 t ha-1 for Fe, 0.01 to 3.78 t ha-1 for Mn, 0.01 to 1.78 t ha-1 for Cu and 0.01 to 1.67 t ha-1 for B. However, the response rates of other group of crops were different. This warrants site-specific S and micronutrients application based on their phyto-availability and crop requirement for sustainable crop production. Therefore, the developed spatial distribution maps could be used by the farmers, farm managers, policy makers, fertilizer industries, planners and extension agencies for understanding the deficiencies of S and micronutrients in different parts of IGP. Accordingly, right kind and quantity of S and micronutrients fertilizer could be produced and distributed for their rational and site-specific application.

thumbnail
Table 4. Crop responses to S and micronutrient application in different states in IGP, India.

https://doi.org/10.1371/journal.pone.0234053.t004

Soil management zones

Principal component analysis, aggregating and summarizing the variability in the nine studied soil parameters and altitude, resulted in ten PCs (Table 5). Out of which, first five PCs having eigenvalue > 0.90 and accounting for 69.62% of total variance were considered (Table 5, Fig 5). The PC1 explained 25.26% of total variance and dominated with pH, available Mn, Fe and Cu and altitude. Additionally, PC2, PC3, PC4 and PC5 explained 15.13, 10.35, 9.86 and 9.01%, respectively, of total variance. Principal component 2 was dominated with available S, Zn and B, PC3 by EC, PC 4 by SOC and PC5 by available Zn. PC1 and PC2 bi-plot of revealed 3 groupings (Fig 6). Soil pH, altitude, EC, available S and B formed one group and available Zn, Fe, Mn and Cu another. Contrastingly, Shukla et al. [17] recorded the grouping of altitude with SOC, available Zn, Mo, Cu, Mn, and Fe in SHR, India. Soil pH, EC, available S and B formed another group. This is mainly attributed to the differences in soil types and prevailing climatic conditions of both the regions.

thumbnail
Fig 6. PC1 and PC2 bi-plot of soil parameters and altitude of soil samples.

https://doi.org/10.1371/journal.pone.0234053.g006

thumbnail
Table 5. Principal component analysis of soil parameters and loading coefficient for the first five PCs.

https://doi.org/10.1371/journal.pone.0234053.t005

The clustering was carried out by considering first 5 PCs. This resulted in 6 clusters considering minimum FPI and NCE values (Table 6) and six nutrient MZs (Fig 7). This result is in line with the findings of Fu et al. [72] (in Sanjiang plain of China), Davatgar et al. [73] (in paddy cultivated areas of Iran), Tripathi et al. [25] (in paddy growing soils of India), Shukla et al. [17] (in SHR, India) and Shukla et. al. [46] (in a DPR, India). The % of area in different MZs was in the order: MZ5 (24.4%) > MZ3 (21.7%) > MZ4 (16.2%) > MZ6 (14.8%) > MZ2 (14.1%) >MZ1 (8.8%). The 6 MZs were different form one another and they had different soil properties and available S and micronutrient concentration (Table 7). This is attributed to the variations in soils, agro-ecological conditions and adoption of soil-crop management practices in different MZs. Altitude is another factor in IGP affecting soil parameters. Hence, appropriate soil management strategies need to be devised in different MZs. For example, the level of SOC content needs to be enhanced in MZ3 and MZ4 through adoption of different management practices for better soil function and nutrient availability. Though mean concentrations of available S and micronutrients in different MZs were higher than the critical range of deficiency for respective nutrients, different levels of deficiency exist in the MZs (Table 7). The MZ1 had the highest area (54.6%) with S deficiency followed by MZ5 (53.6%), MZ2 (49.8%), MZ4 (41.3%), MZ6 (35.8%), and MZ3 (32.8%). Similarly, Zn deficiency area in different MZs followed the order: MZ1 (35.8%) < MZ5 (27.5%) < MZ6 (26.7%) < MZ2 (26.5%) < MZ3 (25.3%) < MZ4 (25.1%). The MZ3 had the highest area with Fe (9.80%) and Mn (8.9%) deficiency whereas; the MZ1 had the highest area with Cu (12.3%) deficiency. The highest area under B deficiency was recorded in MZ5 (32.6%) followed by MZ1 (23.0%), MZ6 (22.7%), MZ4 (22.1%), MZ2 (21.6%) and MZ3 (11.7%). The states in IGP are having different soil problems [29] and have varied average consumption of S and micronutrients fertilizers [74] (Table 8). All the states of IGP had all 6 MZs (Table 8). Punjab had the highest % of area in MZ3 (59.4%) whereas Haryana had the largest area in MZ3 (29.3%) and MZ4 (28.1%). The state of Uttar Pradesh had higher area in MZ3 (21.8%) and MZ (21.7%). Bihar had MZ6 in 40.9% of area and West Bengal had 62.0% of area in MZ5.

thumbnail
Fig 7. Sulphur and micronutrient management zones of IGP, India.

https://doi.org/10.1371/journal.pone.0234053.g007

thumbnail
Table 6. FPI and NCE values deciding optimum cluster numbers for the study area.

https://doi.org/10.1371/journal.pone.0234053.t006

thumbnail
Table 7. Mean values of soil properties and available nutrients in different management zones.

https://doi.org/10.1371/journal.pone.0234053.t007

thumbnail
Table 8. States in IGP with different predominant soil problems, average consumption of S and micronutrients fertilizer and MZs.

https://doi.org/10.1371/journal.pone.0234053.t008

Based on this information nutrient management decisions and supply of nutrients to different zones could be prioritized (Table 9). The MZs having higher deficiency level of a particular nutrient need to receive first attention followed by other MZs having subsequently low levels of deficiency. The MZ3 needs to be paid more attention compared to other MZs for Fe and Mn management. Similarly, B management needs to be prioritized in MZ5, MZ1 and MZ6 compared to other MZs. This will help in optimum utilization of resources. Considering the available S and micronutrient status in soils and per cent area of deficiency, different quantities of customized multi-nutrient mixture fertilizers in different grades could be provided to different MZs for efficient nutrient management. In addition, the farmers and farm mangers could suitable be advised to grow efficient and inefficient crop cultivars in different MZs based on soil nutrient status and resource availability. Nutrient-efficient and nutrient-inefficient cultivars of the crops behave differently (in terms of crop growth and yield) under different soil nutrient status and management practices. For example, resource poor farmers could grow nutrient efficient crop cultivars in nutrient deficient soil and could obtain good crop yield without application of that nutrient and vice versa. It is, therefore, pertinent for the farmers and farm managers of IGP to take cognizance of soil parameters in different MZs to devise simple, easy, cost efficient soil-crop manipulation ways for higher and sustainable crop production.

thumbnail
Table 9. Nutrient management priority for different management zones.

https://doi.org/10.1371/journal.pone.0234053.t009

Conclusion

The study revealed wide spatial variability with moderate (except strong for SOC) spatial dependence for phyto-available S, Zn, B, Fe, Mn and Cu and associated soil properties (pH, EC and SOC) in IGP, India. The concentration of S and micronutrients variedly widely with deficiency level of 67.7% for S, 46.9% for Zn, 44.7% for B, 137% for Fe, 18.3% for Mn and 15.3% for Cu. The range values of the semi-variograms of different soil parameters could be considered for designing future soil sampling strategies in IGP of India. Principal component analysis and fuzzy c-means clustering resulted in 6 MZs of IGP having significantly different values of soil parameters. The generated MZ maps could be used for zone-specific manipulation of available S and micronutrients for sustainable crop production. Further, this study also revealed that soil MZs could be delineated in other cultivated regions of the world for zone specific nutrient supply and management.

Acknowledgments

The authors thank the academic editor and reviewers for their constructive suggestions which helped in improving the quality of the manuscript.

References

  1. 1. Aggarwal PK, Joshi PK, Ingram JSI, Gupta RK. Adapting food systems of the Indo-Gangetic plains to global environmental change: key information needs to improve policy formulation. Environ Sci Policy. 2004; 7: 487–498.
  2. 2. Kumar P, Jha D, Kumar A, Chaudhary MK, Grover RK, Singh RK, et al. Economic analysis of total factor productivity of crop sector in Indo-Gangetic Plain of India by district and region. Agricultural Economics Research Report 2. 2002; IARI, New Delhi, India.
  3. 3. Pingali PL. Green revolution: impacts, limits, and the path ahead. Proc Nat Acad Sci USA. 2012; 109(31): 12302–12308. pmid:22826253
  4. 4. Takkar PN. Micronutrient research and sustainable agricultural productivity in India. J Indian Soc Soil Sci. 1996; 44, 562–581.
  5. 5. Shukla AK, Behera SK. Micronutrients Research in India: Retrospect and Prospects. Preprint, FAI Annual Seminar. 2017; The Fertilizer Association of India, New Delhi. SII-4/1-SII-4/17.
  6. 6. Alloway BJ. Micronutrient deficiencies in global crop production (1st ed.). 2008; Dordrecht: Springer. https://doi.org/10.1007/978‐1‐4020‐6860‐7
  7. 7. Wilhelm SH. Sulfur in soils. J Plant Nutri Soil Sci. 2009; 172: 326–335.
  8. 8. Shukla AK, Behera SK. All India research project on micro- and secondary nutrients and pollutant elements in soils and plants: Research achievements and future thrusts. Indian J Fert. 2019; 15(5): 522–543.
  9. 9. Khurana MPS, Sadana US, Bijay Singh. Sulfur Nutrition of Crops in the Indo-Gangetic Plains of South Asia. In: Jez J, editor, Sulfur: A Missing Link between Soils, Crops, and Nutrition, Agron. Monogr. 2008; 50. ASA, CSSA, SSSA, Madison, WI. p. 11–24. https://doi.org/10.2134/agronmonogr50.c2
  10. 10. Singh VK, Dwivedi BS, Mishra RP, Shukla AK, Timsina J, Upadhyay PK, et al. Yields, soil health and farm profits under a rice-wheat system: Long-term effect of fertilizers and organic manures applied alone and in combination. Agronomy. 2019; 9(1): 1.
  11. 11. Adhikari P, Shukla MK, Mexal JG. Spatial variability of soil properties in an arid ecosystem irrigated with treated municipal and industrial waste water. Soil Sci. 2012; 177(7): 458–469.
  12. 12. Shukla AK, Behera SK, Lenka NK, Tiwari PK, Prakash C, Malik RS. Spatial variability of soil micronutrients in the intensively cultivated Trans-Gangetic Plains of India. Soil Till Res. 2016; 163: 282–289.
  13. 13. Tesfahunegn GB, Tamene L, Vlek PLG. Catchment-scale spatial variability of soil properties and implications on site-specific soil management in northern Ethiopia. Soil Till Res. 2011; 117: 124–139.
  14. 14. Behera SK, Shukla AK. Spatial distribution of surface soil acidity, electrical conductivity, soil organic carbon content and exchangeable potassium, calcium and magnesium in some cropped acid soils of India. Land Deg Dev. 2015; 26: 71–79.
  15. 15. Ferreira V, Panagopoulos T, Andrade R, Guerrero C, Loures L. Spatial variability of soil properties and soil erodibility in the Alqueva reservoir watershed. Solid Earth. 2015; 6: 383–392.
  16. 16. Bogunovic I, Pereira P, Brevik EC. Spatial distribution of soil chemical properties in an organic farm in Croatia. Sci Total Env. 2017; 584‐585: 535–545.
  17. 17. Shukla AK, Sinha NK, Tiwari PK, Prakash C, Behera SK, Lenka NK, et al. Spatial distribution and management zones for sulfur and micronutrients in Shiwalik Himalayan region of India. Land Deg Dev. 2017; 28 (3): 959–969.
  18. 18. Behera SK, Mathur RK, Shukla AK, Suresh K, Prakash C. Spatial variability of soil properties and delineation of soil management zones of oil palm plantations grown in a hot and humid tropical region of southern India. Catena. 2018; 165: 251–259.
  19. 19. Warrick AW, Meyers DE, Nielsen DR. Geostatistical methods applied to soil science. pp. 53–82. In Klute A. (Ed.) Methods of Soil Analysis. Part 1. 2nd ed. 1986; Agronomy Monograph 9. ASA and SSSA, Madison, WI.
  20. 20. White JG, Welch RM, Norvell WA. Soil zinc map of the USA using geostatistics and geographical information system. Soil Sci Soc Am J. 1997; 61: 185–194.
  21. 21. Mueller TG, Hartsock NJ, Stombaugh TS, Shearer SA, Cornelius PL, Barnhise RI. Soil electrical conductivity map variability in limestone soil overlain by loess. Agron J. 2003; 95: 496–507.
  22. 22. Goovaerts P. Geostatistical tools for characterizing the spatial variability of microbiological and physio-chemical soil properties. Biol Fert Soils. 1998; 27: 315–334.
  23. 23. Xin-Zhong W, Guo-Shun L, Hong-Chao H, Zhen-Hai W, Qing-Hua L, Xu-Feng L, et al. Determination of management zones for a tobacco field based on soil fertility. Comp Electr Agri. 2009; 65: 168–175.
  24. 24. Fu W, Tunney H, Zhang C. Spatial variation of soil nutrients in a dairy farm and its implications for site-specific fertilizer application. Soil Till Res. 2010; 106: 185–193.
  25. 25. Tripathi R, Nayak AK, Shahid M, Lal B, Gautam P, Raja R., et al. Delineation of soil management zones for a rice cultivated area in eastern India using fuzzy clustering. Catena. 2015; 133: 128–136.
  26. 26. Peralta NR, Costa JL, Balzarini M, Franco MC, Cordoba M, Bullock D. Delineation of management zones to improve nitrogen management of wheat. Comp Electr Agri. 2015; 110: 103–113.
  27. 27. Ortega RA, Santibanez OA. Determination of management zones in corn (Zea mays L.) based on soil fertility. Comp Electr Agri. 2007; 58(1): 49–59.
  28. 28. Soil Survey Staff. Keys to Soil Taxonomy (12th ed.). 2014; Washington, DC: USDA-Natural Resources Conservation Service.
  29. 29. Pal DK, Bhattacharyya T, Srivastava P, Chandran P, Ray SK. Soils of the Indo-Gangetic Plains: their historical perspective and management. Curr Sci. 2009; 96: 1193–1202.
  30. 30. Kerry R, Oliver MA, Frogbrook ZL. Sampling in precision agriculture. In Oliver M. A. (Ed.), Geostatistical Applications for Precision Agriculture (pp. 35–63). 2010; Dordrecht: Springer. https://doi.org/10.1007/978‐90‐481‐9133‐8‐2
  31. 31. Jackson ML. Soil Chemical Analysis. 1973; New Delhi: Prentice Hall of India Pvt. Ltd.
  32. 32. Walkley AJ, Black IA. An examination of the Degtjareff method for determining soil organic matter and a proposed modification of the chromic acid titration method. Soil Sci. 1934; 37: 29–38.
  33. 33. Williams CH, Steinbergs A. Soil sulphur fractions as chemical indices of available sulphur in some Australian soils. Aust J Agri Res. 195; 10: 340–352.
  34. 34. Gupta UC. A simplified method for determining hot-water soluble boron in podzol soils. Soil Sci. 1967; 103: 424–428.
  35. 35. Lindsay WL, Norvell WA. Development of a DTPA soil test for zinc, iron, manganese and copper. Soil Sci Soc Am J. 1978; 42: 421–448.
  36. 36. SAS Institute. The SAS system for Windows. 2011; Release 9.2. Cary, NC: SAS Inst.
  37. 37. Goovaerts P. Geostatistics for Natural Resources Evaluation. 1997; New York: Oxford Univ. Press.
  38. 38. Cambardella CA, Moorman TB, Novak JM, Parkin TB, Karlen DL, Turco RF, et al. Field-scale variability of soil properties in central Iowa soils. Soil Sci Soc Am J. 1994; 58: 1501–1511.
  39. 39. Minasny B, McBratney AB. A conditioned Latin hypercube method for sampling in the presence of ancillary information. Comp Geosci. 2006; 32(9): 1378–1388.
  40. 40. Gruba P, Socha J. Effect of parent material on soil acidity and carbon content in soils under silver fir (Abies alba Mill.) stands in Poland. Catena. 2016;140: 90–95.
  41. 41. Zhang YY, Wu W, Liu H. Factors affecting variations of soil pH in different horizons in hilly regions. PLoS ONE. 2019; 14(6): e0218563. pmid:31216328
  42. 42. Srivastava P, Pal DK, Aruche KM, Wani SP, Sahrawat KL. Soils of the Indo-Gangetic Plains: a pedogenic response to landscape stability, climatic variability and anthropogenic activity during the Holocene. Earth-Sci Rev. 2015; 140: 54–71.
  43. 43. Lal R. Soil carbon sequestration in India. Climatic Change. 2004; 65: 277–296.
  44. 44. Shukla AK. Tiwari PK. Micro and secondary nutrients and pollutant elements research in India. Coordinators Report, pp. 1–196. 2016; AICRP on Micro and Secondary Nutrients and Pollutant Elements in Soils and Plants, ICAR-IISS, Bhopal.
  45. 45. Singh VK, Dwivedi BS, Shukla AK, Kumar V, Gangwar B, Rani M, et al. Status of available sulfur in soils of north-western Indo-Gangetic plain and western Himalayan region and responses of rice and wheat to applied sulfur in farmer’s fields. Agri Res. 2015; 4: 76–92.
  46. 46. Shukla AK, Sinha NK, Tiwari PK, Prakash C, Behera SK, Babu PS, et al. Evaluation of spatial distribution and regional zone delineation for micronutrients in a semi-arid deccan plateau region of India. Land Deg Dev. 2018; 29(8): 2449–2459.
  47. 47. Sahrawat KL, Shirisha K, Murthy KVS, Gajbhiye PN, Kundu S, Wani SP, et al. Comparative evaluation of inductively coupled plasma-atomic emission spectroscopy and colorimetric methods for determining hot-water-extractable boron in soils. Commun Soil Sci Plant Anal. 2012; 43: 948–957.
  48. 48. Anderson DW. The effect of parent material and soil development on nutrient cycling in temperate ecosystems. Biogeochemistry. 1988; 5(1): 71–97.
  49. 49. Chandi KS, Takkar PN. Effects of agricultural cropping systems in micronutrient transformation. I. Zinc. Plant Soil. 1982; 69: 423–436.
  50. 50. Katyal JC, Rao NH, Reddy MN. Critical aspects of organic matter management in the Tropics: the example of India. Nutr Cycl Agroeco. 2001; 61:77–86.
  51. 51. Singh B, Singh V. Productivity and fertility of soils in the Indo-Gangetic Plains of South Asia. Arch Agron Soil Sci. 2012; 58: S33–S40.
  52. 52. Nielsen DR, Bouma J. Soil spatial variability. Proceedings of a Workshop of the ISSS and the SSSA, Las Vegas, USA. 1985; 30 November to 1 December 1984. Pudoc: Wageningen, The Netherlands; 243.
  53. 53. Wang L, Wu JP, Liu YX, Huang HQ, Fang QF. Spatial variability of micronutrients in rice grain and paddy soil. Pedosphere. 2009; 19(6): 748–755.
  54. 54. Brady NC, Weil RR. The Nature and Properties of Soil, 13th ed. 2002; Springer Netherlands, p.249.
  55. 55. Neilsen D, Hoyt PB, Mackenzie AF. Distribution of soil Zn fractions in British Columbia interior orchard soils. Canadian J Soil Sci. 1986; 66: 445–454.
  56. 56. Rengel Z. Availability of Mn, Zn and Fe in the rhizosphere. J Soil Sci Plant Nutri. 2015; 15(2): 397–409.
  57. 57. Wei X, Hao M, Shao M, Gal WJ. Changes in soil properties and the availability of soil micronutrients after 18 years of cropping and fertilization. Soil Till Res. 2006; 91: 120–130.
  58. 58. Katyal JC, Vlek PLG. Micronutrient problems in tropical Asia. Fert Res. 1985; 7: 131–150.
  59. 59. Corwin DL, Lesch SM. Apparent soil electrical conductivity measurements in agriculture. Comp Electr Agri. 2005; 46: 11–43.
  60. 60. Freney JR, Melville GE, Williams CH. Soil organic matter fractions as sources of plant-available sulphur. Soil Biol Biochem.1975; 7(3): 217–221.
  61. 61. De Bona FD, Monteiro FA. Nitrogen and sulfur fertilization and dynamics in a Brazilian Entisol under pasture. Soil Sci Soc Am J. 2010; 74(4): 1248–1258.
  62. 62. Kertesz MA, Mirleau P. The role of soil microbes in plant sulphur nutrition. J Ex Bot. 2004; 55(404): 1939–1945. pmid:15181108
  63. 63. Behera SK, Shukla AK. Depth-wise distribution of zinc, copper, manganese and iron in acid soils of India and their relationship with some soil properties. J Indian Soc Soil Sci. 2013; 61(3): 244–252.
  64. 64. Lopez‐Granados F, Jurado‐Exposito M, Atenciano S, Garcia‐Ferrer A, De la Orden MS, Garcia‐Torres L. Spatial variability of agricultural soil properties in southern Spain. Plant Soil. 2002; 246: 97–105.
  65. 65. Pal DK, Srivastava P, Durge SL, Bhhatacharya T. Role of microtopography in the formation of sodic soils in the semi-arid part of the Indo-Gangetic Plains, India. Catena. 2003; 51: 3–31.
  66. 66. Singh YP, Mishra VK, Singh S, Sharma DK, Singh D, Singh US, et al. Productivity of sodic soils can be enhanced through the use of salt tolerant rice varieties and proper agronomic practices. Field Crops Res. 2016; 190: 82–90. pmid:27212787
  67. 67. Bhattacharyya T, Pal DK, Easter M, Batjes NH, Milne E, Gajbhiye KS, et al. Modelled soil organic carbon stocks and changes in the Indo-Gangetic Plains, India from 1980 to 2030. Agri Ecosyst Env. 2007;122 (1): 84–94.
  68. 68. Ghosh PK, Hazra KK, Venkatesh MS, Nath CP, Singh J, Nadarajan N. Increasing soil organic carbon through crop diversification in cereal–cereal rotations of Indo-Gangetic plain. Proc National Acad Sci India Section B: Biol Sci. 2019; 89: 429–440.
  69. 69. Jat HS, Datta A, Sharma PC, Kumar V, Yadav AK, Choudhary M, et al. Assessing soil properties and nutrient availability under conservation agriculture practices in a reclaimed sodic soil in cereal-based systems of North-West India. Arch Agron Soil Sci. 2018; 64(4): 531–545. pmid:30363929
  70. 70. Tandon HLS. Sulphur in soils, crops and fertilizers. 2011; Fertilizer Development and Consultation Organization, New Delhi, India, pp 204 + X.
  71. 71. Takkar PN, Chhibba IM, Mehta SK. Twenty Years of Coordinated Research on Micronutrients in Soils and Plants. 1989; Bhopal, India: Indian Institute of Soil Science.
  72. 72. Behera SK, Shukla AK, Singh M, Dwivedi BS. Extractable boron in some acid soils of India: status, spatial variability and relationship with soil properties. J Indian Soc Soil Sci. 2016; 64(2): 183–192.
  73. 73. Davatgar N, Neishabouri MR, Sepaskhah AR. Delineation of site-specific nutrient management zones for a paddy cultivated area based on soil fertility using fuzzy clustering. Geoderma. 2012; 173–174: 111–118.
  74. 74. Fertilizer Association of India. Speciality fertilizer and micronutrient statistics 2016–2017. 2017; The Fertilizer Association of India, New Delhi, India, p 207.