African Journal of
Agricultural Research

  • Abbreviation: Afr. J. Agric. Res.
  • Language: English
  • ISSN: 1991-637X
  • DOI: 10.5897/AJAR
  • Start Year: 2006
  • Published Articles: 6863

Full Length Research Paper

Spatial variability of available nutrients in soils of Nainpur tehsil of Mandla district of Madhya Pradesh, India using Geo-statistical approach

G. S. Tagore
  • G. S. Tagore
  • Department of Soil Science and Agricultural Chemistry, Jawaharlal Nehru Krishi Vishwa Vidyalaya, Jabalpur (M. P.), India.
  • Search for this author on:
  • Google Scholar
B. Singh
  • B. Singh
  • Department of Soil Science and Agricultural Chemistry, Rajmata Vijayaraje Scindia Krishi Vishwavidyalaya, Indore (M. P.), India.
  • Search for this author on:
  • Google Scholar
P. S. Kulhare
  • P. S. Kulhare
  • Department of Soil Science and Agricultural Chemistry, Jawaharlal Nehru Krishi Vishwa Vidyalaya, Jabalpur (M. P.), India.
  • Search for this author on:
  • Google Scholar
R. D. Jatav
  • R. D. Jatav
  • Department of Agriculture, Mandla (M.P.) India
  • Search for this author on:
  • Google Scholar


  •  Received: 04 February 2015
  •  Accepted: 16 July 2015
  •  Published: 20 August 2015

 ABSTRACT

Scientific information concerning spatial variability and distribution of soil properties is critical for farmers attempting to increase fertilization efficiency and crop productivity; fertilization based on maps with recommendations related to soil fertility may lead to reduced fertilizer inputs without reducing yield. In the present study, GPS based one hundred fifty surface soil samples (0-15 cm) were collected from dominant cropping system. After processing, the soil samples were analyzed for different soil characteristics in laboratory using standard procedures. The data obtain from laboratory analysis was statistically and geo-statistical interpreted. The results revealed that the 23.6, 28.30, 48.6, 13.9, 25.5 and 54.7% soil samples were found to be deficient in OC, N, P, K, S and Zn, respectively. None of the soil samples were tested low in Cu, Fe, Mn and B. Exponential model was found as the best fit for considered soil parameters whereas, spherical model was found as the best fit for Mn. The best model was used to generate the spatial distribution maps. Spatial maps showed that the soil pH, EC, organic carbon, available N, P, K, S, Zn, Cu, Fe, Mn and B spatially varied and N, P, K, S and Zn were deficient in major areas. Therefore these maps are more useful for guiding site-specific field management for agricultural production and environmental protection. In addition, reduce the losses of nutrients and could be save time and money for fertilizers.

 

Key words: Geo-statistical, soil types, land use, semi-variogram, kriging, nutrient status.


 INTRODUCTION

Soil is the primary source of plant nutrients that control its fertility and thus yield of crops (Qu et al., 2014). Expectation of higher productivity using adequate amount of fertilizer nutrients may lead to become limiting to some micronutrients in the soil and most times due to their over mining by the crops and shortage of which often show the
 
deficiency symptoms and yields are reduced. However, intensive cultivation, indiscriminate use of high analysis chemical fertilizers results deficiency of micronutrients (Singh et al., 2007). The soil must supply all the essential nutrients for desired growth of plants has contributed to a tremendous increase in yields of agricultural crops that feed the world’s population (Hak-Jin et al., 2009). Ideally, application rates should be adjusted based on estimates of the requirements for optimum production at each location because there is high spatial variability of N, P, and K within individual fields (Page et al., 2005; Ruffo et al., 2005). In this context, it is necessary to evaluate the fertility status of the soil and promote the recommendations of soil test for balanced nutrition to maintain soil health. Classical statistics requires the validity of some basic hypotheses, such as the independence between observations, due to the randomness of variations from one place to another. Farmer’s of the tribal areas have very little knowledge of fertility status of their fields; therefore, the desired level of productivity is not being achieved even after utilizing all inputs.
 
In the last few decades, geo-statistics has been used extensively to characterize the spatial variability of soil attributes due to its ability of quantifying and reducing sampling uncertainties and minimizing investigation costs (Cambule et al., 2014). Numerous studies have been conducted based on geo-statistical analysis to characterize the spatial variability of soil physical (Li et al., 2007; Weindorf and Zhu, 2010), chemical (Liu et al., 2013; Lin et al., 2009; Huang et al., 2007), and biochemical (Šnajdr et al., 2008) properties, as well as microbiological processes (Cao et al., 2011). Thus, information on spatial variability of soil nutrients is important for sustainable management of soil fertility (Fraisse et al., 1999). Hence, the present study was conducted. 


 MATERIALS AND METHODS

Description study area

 

Geographical location of study area

 

Mandla district is located in the east-central part of the Madhya Pradesh, India covering an area of 8771 km2 and consists of a rugged high tableland in the eastern part of the Satpura hills.  Geographically, Nainpur is located between 79°56’45.0” to 80°26’15.0” E longitudes and 22°17’31.0” to 22°37’30” N latitudes, having an area of 845.79 km2 of the total geographical area of district.

 

Climate

 

The climate of this district is characterized by hot summer season and general dryness except in the southwest monsoon season. May is the hottest month with the mean daily minimum temperature at 41.3°C and the mean daily minimum at 24°C (Figure 1).

 

 

Land use and cropping pattern

 

Land use map prepared by using Indian remote-sensing satellite-P6, linear imaging self-scanning satellite-III (IRS-P6, LISS-III) satellite imagery dated January-2014, October-2013 and April-2013. The satellite data has the characteristics of 23.5 m spatial resolution, four spectral channels green (0.52 to 0.59 µ), red (0.62 to 0.68 µ), NIR (0.77 to 0.86 µ), and SWIR (1.55 to 1.70 µ) and five days temporal resolution with 141 km swath. The statistics reveals area extent of land mainly under  agricultural (46.37%), followed by dense forest (19.98%), open forest (9.79%), fallow land (9.23%) waste  land- dense (3.22%) and  waste  land open scrub (6.53%) and others etc. Based on ground truth data collected from local agriculture departments and farmers interviewed, the rice and wheat are major food grain crops. Paddy, maize, kodo, kutki, soybean are important crops during Kharif season and wheat, pea, gram, lentil, and mustard crops during Rabi season in tribal areas of Mandla district (Figure 2).

 

 

Soil types

 

The survey of India topographical maps (1:50000) and the map of soil as a secondary data was used from NBSSLUP Nagpur. The highest area occupied under vertisols followed by inceptisols, entisols and alfisols. These soils are fine montmorillonitic, hyperthermic and having high swell shrink potential (Table 1 and Figure 3). 

 

 

 

Soil sampling and analysis

 

Sampling sites were generated using land use and soil association maps. The sites decided randomly distributed over agricultural land of the study area by considering of topography and heterogeneity of the soil type. Field data collection and soil sampling were carried out by using GPS by navigating those points. One hundred fifty soil samples (0 to 15 cm) were collected from farmer’s field during the 2013 off season from the agricultural land. For each main sampling point, 1.0 kg of representative composite soil sample was collected and logged into properly labeled sample bag.

 

 

Laboratory analysis

 

Soil samples collected from the study area were dried and crushed with the help of wooden rod and passed through 2 mm sieve and then used for the determination of soil pH, electrical conductivity, organic carbon, calcium carbonate and macronutrients like N using Subbiah and Asija (1956), P using Olsen et al. (1954) and K content by adopting standard laboratory methods described in Jackson (1973) (Table 2).

 

 

Available micronutrients (Zn, Cu, Fe and Mn) were extracted by DTPA-CaCl2 solution and analyzed using atomic absorption spectrophotometer (Lindsay and Norvell, 1978). Hot water soluble boron in soil was analyzed by azomethine-H method as outlined by Berger and Truog (1939). The available sulphur was extracted by 0.15% CaCl2 solution and the concentration of sulphur was determined by the turbidimetric method using spectrophotometer (Chesnin and Yien, 1951).

 

Nutrient index calculation

 

The nutrient index (NI) values for available nutrients present in the soils were calculated utilizing the formula suggested by Parker et al. (1951) and classified this index as low (<1.67), medium (1.67 to 2.33) and high (>2.33).

 

 

Where: Nl, Nm and Nh are the number of soil samples falling in low, medium and high categories for nutrient status and are given weight age of 1, 2 and 3, respectively. Nt is the total number of sample.

 

Statistical and geo-statistical analysis

 

It is necessary to check whether the available contents of N, P, K and S and micronutrients in soil samples are approximately normally distributed or not because Kriging assumes the normal distribution for each estimated variable. A normal distribution was estimating based on  skewness  values  and  the  variable  datasets

 

having a skewness ranging between -1 and 1 were considered normally distributed. For non-datasets, a logarithmic transformation was performed to achieve a normal distribution for use in the next step of the statistical analysis.

 

Geo-statistical methods were used to analyze the spatial correlation structures of the available contents of N, P, K and S and micronutrients in soil and spatially estimate their values at unsampled locations using geo-statistical tool in GIS 9.3.1 software. The spatial dependency of selected soil parameters was analyzed using semi-variogram analyses with normalized data. Semi variogam analyses have been proven as an excellent approach to exploring the structure of spatial variogram in agricultural soils.

 

The above phenomena is the best accomplished studying the semivariogram (Warrick et al., 1986) which is a plot of semi-variance that characterizes the rate of change of a mapped variable with respect to distance. Semi-variogram  is computed as half the average squared difference between the soil properties of data pairs. The structure of spatial variability was analyzed through semi-variograms (Figure 4).

 

 

 

 

 

 

A semi variogram was calculated for each soil property. The semi variance γ(h) is estimated as:

 

 

where N(h) is the number of data pairs within a given class of distance and direction, z(xi) is the value of the variable at the location xi and z(xi + h) is the value of the variable at a lag of h from the location xi.

 

An experimental semi-variogram was calculated using the measured data. Next, this was generally fitted with a theoretical model, such as Exponential, spherical and Gaussian models (Goovaerts, 1999). Choice of the best-fitted model was based on the lowest residual sum of square (RSS) and the largest coefficient of determination (R2). The model provided information about the spatial structure as well as the input parameters (that is, nugget, sill and range) for the Kriging interpolation. Nugget is the variance at distance zero, sill is the semi- variance value at which the Semi-variogram reaches the upper bound after its initial increase, and range is a value (x axis) at which one variable becomes spatially independent.

 

The nugget to sill ratio was used to define different classes of spatial dependence for the soil properties (Emadi et al., 2008; Zuo et al., 2008). Nugget/sill ratio of 25%, 25 to 75% and >75% were classified as having strong, moderate and weak spatial dependence, respectively, according to Cambardella et al. (1994). Ordinary Kriging was used for the spatial interpolation because it is best at providing an unbiased prediction for specific unsampled locations and minimizing the influence of outliers (Tesfahunegn et al., 2011; Kavianpoor et al., 2012).

 

Evaluation criteria

 

Accuracy of the soil maps was evaluated through cross-validation approach (Davis et al., 1987; Santra et al., 2008). Among three evaluation indices used in this study, mean absolute error (MAE), and mean squared error (MSE) measure the accuracy of pre effectiveness of prediction. MAE is a measure of the sum of the residuals (e.g. predicted minus observed)(Voltz and Webster, 1990).

 

 

where z is the sample mean. If G = 100, it indicates perfect prediction, while negative values indicate that the predictions are less reliable than using sample mean as the predictors.


 RESULTS AND DISCUSSION

Descriptive statistics of soil properties
 
The mean values for pH, EC, OC content were 6.58, 0.21 dSm-1 and 6.30 g kg-1 with a range of 4.8 to 7.90, 0.02 to 0.60 dS m-1 and 3.0 to 13.30 g kg-1, respectively. The available N, P, and K varied from 150 to 500 kg ha-1, 5.88 to 34.89 kg ha-1 and 122.08 to 807.52 kg ha-1 with mean value of 299.43 kg ha-1, 12.31 kg ha-1 and 405.30 kg ha-1, respectively. The S varied from 2.04 to 39.0 mg kg-1 with a mean value of 14.69 mg kg-1. The available micronutrients Zn, Cu, Fe, Mn and B varied from 0.11 to 6.08, 0.41 to 16.35, 5.15 to 101.20, 2.24 to 61.80 mg kg-1 and 1.40 to 3.72 with mean  values  of  0.83,  4.81, 40.81, 35.29 and 2.54 mg kg-1, respectively (Table 3).
 
 
The coefficient of variation, which is the ratio of the standard deviation to mean expressed as a percentage is a useful measure of overall variability. According Hillel (1980) and Yang and Yang (2005), there are three classes about coefficient of variation of soil nutrients (C.V.), weak variation (C.V. <0.1), medium variation (C.V. = 0.1–1.0), strong variation (C.V. > 1.0). The EC had the largest variation (CV = 50.12%) followed by OC (CV = 27.83%) while pH was found least variable (CV = 8.71%). The macronutrients are, NPK, the available P had the highest variability (CV=49.21%) and followed by avail. K (CV = 35.75%) whereas N had the lowest variability (CV= 27.00). The S was found to be moderately variable (CV = 56.84%). Among the micronutrients, the Zn was found to be highly variable (CV = 98.80%), followed by Fe (CV = 86.22%), Cu (CV =60.55 percent) and Mn (CV = 52.57 percent) while B was found least variable (CV= 25.75%). The range of CV for the area suggested different degrees of heterogeneity among the properties studied. All the macronutrients were found to be moderately variable with ranging from 27.00 to 49.21%. All the micronutrients were highly variable with CV ranging from 25.57 to 98.80%. Among the soil fertility parameters, CaCO3, Zn, Cu and Fe were found not normally distributed due to higher value of skewness and kurtosis.
 
Soil nutrient status
           
It is evident from the data presented in Table 4 that soil organic carbon content was found to be  low in 23.6% and only 2% soil samples were medium and rest of soil samples were high (74.5%). It indicates that majority of these soils were low to moderately high in organic carbon content. On the basis of the ratings suggested by Subbiah and Asija (1956), the available N was found to be low (28.30%), medium (1.9%) and high (74.5%) in soil samples. It may be due low organic matter content of soil as well as rapid loss of applied N in soil. On the basis of the limits suggested by Muhr et al. (1963), the available phosphorus content was found to be low 48.6% in soil samples, none of the samples were found in medium and rest of soil samples were fall in high category (51.4%). The low amount of available P may be due to application of lower doses of P fertilizer, fixation of P on clay minerals or CaCO3 surfaces with the time elapsed between fertilizer application and crop uptake. It may also be due low organic matter content of soil (Sanyal et al., 2015). 
 
 
According to Muhr et al. (1963), the available potassium content was found to be low in 13.9 soil samples, high in 86.1% soil samples and non of soil samples were found in medium. For these soils available micronutrients status using critical limit as given in Table 4 for DTPA-extractable Zn content was found to be  deficient in 54.7%, medium in 26.4%  and high in 18.9% soil samples. The data showed the DTPA  extractable  Fe content was found medium in 87.7% soil samples and none of the soil samples were tested low in Cu, Fe, Mn and B content, respectively. All soil samples were found insufficient in case of Cu content. The present study result was supported by Athokpam et al. (2013). Considering soil nutrient index suggested by Parker et al. (1951) (Table 4), soils of Nainpur were found of high status in respect of all fertility status except medium P while low fertility status in case of Zn.
 
Spatial variability maps preparation using geo-statistical approach
                      
It was necessary to normalize the data prior to the geo-statistical analysis because of high skewness (Table 3) and the presence of outliers. Logarithmic transformations were selected to normalize the dataset. The skewness and kurtosis coefficients are zero for a normally distributed random variable. If the data distributions are largely deviated from a normal distribution, data transformations are often  performed  in  order  to  reduce the influence of extreme values on spatial analysis (Webster and Oliver, 2001).
 
where f(x) is the transformed value and x is the value to be transformed. For a given data set (x1, x2,…..xn), the parameter is estimated based on the assumption that the transformed values (y1,y2,…..yn) are normally distributed. When λ = 0, the transformed becomes the logarithmic transformation.
 
The logarithmic transformation resulted in smaller skewness and kurtosis for CaCO3, Zn, Cu and Fe and the transformed data passed the normality test. Log-transformed data were used in the spatial analysis (Table 5).
 
 
Geo-statistical analysis
 
Ordinary Kriging was chosen to create the spatial distribution  maps  of  soil  fertility  parameters,   with   the maximum search radius being set to the autocorrelation range of the corresponding variable. The best model for fitting on experimental variogram was selected based on less RSS value. The exponential model was selected from standard models that are available to fit experimental semi-variograms based on more favorable weighted residual mean squares, and visual fit to the data at short lags. Wasiullah et al. (2010) and Gore et al. (2012) reported the spherical model for Mn. Therefore, the authors recognized the exponential and spherical model to be suitable for estimation of soil properties. The results supported by Ferreira et al. (2015) who reported exponential model is best for most of soil properties. Similar result reported by Eltaib et al. (2002), Nayanaka et al. (2010) and Cao et al. (2011) in case of soil organic carbon (Table 6).
 
 
In this study, no apparent anisotropy was found for any studied variable through experimental variograms. So, all experimental variograms were in isotropic form and ratio of nugget variance to sill expressed in percentages (C0/C+C0) can be regarded as a criterion for classifying the spatial dependence of the soil parameters. If this ratio is less than 25%, then the variable has strong spatial dependence (Shi et al., 2005) 25% and 75%, and higher than 75% correspond to moderate, and weak spatial dependencies, respectively. The OC, Zn, Cu, Fe and  Mn showed strong spatial dependence and the values between 25 and 75% have moderate spatial dependence as shown in table, pH, EC, CaCO3, N, P, K, S and B have moderate spatial structure. According to Cambardella et al. (1994), these spatial auto dependencies may be attributed to both intrinsic factors such as other soil properties and extrinsic factors such as human activities (Table 7). 
 
 
The mean absolute error (MAE) was found to be zero for pH, EC, OC, P, S, Zn, Cu, Fe and Mn. The highest mean squared error (MSE) was observed for K followed by N, Fe, CaCO3, Mn, S, P, Cu, OC, Zn, pH and EC. The performance of the interpolation technique was considered in term of the G value nutrients that showed positive. The goodness of fit (G) values was positive and highest G value was observed for Fe followed by Mn, EC, pH, Cu, Zn. The G values for B and P were negative. The spatial distribution of P could be related to geology, especially parent material of the study area (Wang et al., 2009). Among the micronutrients, available Zn was spatially  correlated for a shorter distance and at a distance less than the range, measured soil property of two samples become more alike with decreasing distance between them (Eltaib  et al., 2002). Spatial dependence of DTPA extractable Zn, Fe, Cu and Mn were also reported by Nayak et al. (2006).


 CONCLUSIONS

On the basis of laboratory data, soil fertility assessment in respect to available macro and micro nutrients have been drawn in terms of nutrient indexing (Low, medium and high). The available N, P, K and S were high, medium, high respectively and micro-nutrients, that is, Zn was low and rest of that were sufficient in their status. It was concluded that it is applicable to investigate the spatial variability of soil properties. Exponential model was found as the best fit of all the parameters except Mn. The spherical model was found as the best fit for Mn. Based on prediction maps, it was realized that all nutrients has no toxic status and that the application of fertilizer will improve crop yields. Hence, study area requires phosphatic and zinc fertilizers for better crop growth and productivity of crops. The spatial variability maps on nutrient status generated geo-statistical techniques will be helpful for the farmers or researchers for the location specific correction of nutrient deficiency. These maps are very much helpful for guiding the farmers to decide the amount and kind of nutrient they can apply for in terms of optimum/economic returns, as the nutrient management will be different for areas having deficiency of one or more nutrients than those having sufficient nutrients. 


 CONFLICT OF INTEREST

The authors have not declared any conflict of interest.



 REFERENCES

Agterberg FP (1984). Trend surface analysis. In Spatial Statistics and Models (eds Gaile, G. L. and Willmott, C. J.), Reidel, Dordrecht, The Netherlands pp. 147–171.
Crossref
 
Athokpam H, Wani SH, Kamei D, Athokpam HS, Nongmaithem J, Kumar D, Singh YK, Naorem BS, Thokchom RD, Devi L (2013). Soil macro- and micronutrient status of Senapati district, Manipur (India). Afr. J. Agric. Res. 8(39):4932-4936.
 
Berger KC, Truog E (1939). Boron determination in soils and plants. Ind. Eng.Chem. Anal. Ed.11: 540-545.
Crossref
 
Cambardella CA, Moorman TB, Novak JM, Parkin TB, Karlen DL, Turco RF, Konopka AE (1994). Field scale variability of soil properties in Central Iowa soils. Soil Sci. Soc. Am. J. 58:1501–1511.
Crossref
 
Cambule AH, Rossiter DG, Stoorvogel JJ, Smaling EMA (2014). Soil organic carbon stocks in the Limpopo National Park, Mozambique: amount, spatial distribution and uncertainty. Geoderma 213:46–56.
Crossref
 
Cao SK, Chena KL, Cao GC, Zhang L, Ma J, Yang L, Lu BL, Chen L, Lu H (2011).The Analysis of Characteristic and Spatial Variability for Soil Organic Matter and Organic Carbon around Qinghai Lake. Procedia Env. Sci. 10:678 –684.
Crossref
 
Chesnin L, Yien CH (1951). TurbidImetric determination of available suiphates. Soil Sci. Soc. Am. Proc. 15:149-151.
Crossref
 
Davis BM (1987). Uses and abuses of cross-validation in geo-statistics. Math. Geol. 19:241–248.
Crossref
 
Eltaib SM, Amin, MSM, Hanfi MM, Sharif, ARM Wayayok A (2002). Spatial variability of N, P and K in rice field in Sawah Sempadan, Malaysia. Songklanakarin J. Sci. Technol. 24:321–328.
 
Emadi M, Baghernejad M, Maftoun M (2008). Assessment of some soil properties by spatial variability in saline and sodic soils in Arsanjan plain. Southern Iran. Pak. J. Biol. Sci. 11(2):238–243.
Crossref
 
Ferreira V, Panagopoulos T, Andrade R, Guerrero C, Loures L (2015). Spatial variability of soil properties and soil erodibility in the Alqueva reservoir watershed. Solid Earth 6:383–392.
Crossref
 
Fraisse CW, Sudduth KA, Kitchen NR, Fridgen JJ (1999). Use of unsupervised clustering algorithms for delineating within Colfield management zones. ASAE Paper 993043, ASAE, Ontario, Canada.
 
Goovaerts P (1999). Geo-statistics in soil science: state-of-the-art and perspectives. Geoderma 89:1-45.
Crossref
 
Gore Yogita D, Nagaraju MSS, Srivastava R, Nasre RA (2012). Spatial variability of soil properties in basaltic terrain for Precision agriculture using geospatial techniques: a case study in savli village of Wardha district of Maharashtra. Proceeding of AIPA.
 
Hak-Jin K, Kenneth K, Sudduth A, Hummel JW (2009). Soil macronutrient sensing for precision agriculture. J. Environ. Monit.11: 1810-1824.
Crossref
 
Hillel D (1980). Fundamentals of Soil Physics. Academic Press, Inc., New York.
 
Huang B, Sun WX, Zhao YC, Zhu J, Yang RQ, Zhou Z, Ding F, Su JP (2007). Temporal and spatial variability of soil organic matter and total nitrogen in an agricultural ecosystem as affected by farming practices. Geoderma 139:336–345.
Crossref
 
Jackson ML (1973). Soil Chemical Analysis. Prentice Hall of India Pvt. Ltd., New Delhi.
 
Kavianpoor H, Esmali Ouri A, Jafarian Jeloudar Z, Kavian A (2012). Spatial variability of some chemical and physical soil properties in nesho mountainous rangelands. Am. J. Environ. Eng. 2(1):34-44.
Crossref
 
Li XG, Li FM, Zed R, Zhan ZY, Singh B (2007). Soil physical properties and their relations to organic carbon pools as affected by land use in an alpine pastureland. Geoderma 15:98–105.
Crossref
 
Lin JS, Shi XZ, Lu XX, Yu DS, Wang HJ, Zhao, YC, Sun, WX (2009). Storage and spatial variation of phosphorus in paddy soils of China. Pedosphere 19:790–798.
Crossref
 
Lindsay WL, Norvell WA (1978). Development of a DTPA soil test for Zn, Fe, Mn and Cu. Soil Sci. Soc. Am. J. 42: 421-428.
Crossref
 
Liu ZP, Shao MA, Wang YQ (2013). Spatial patterns of soil total nitrogen and soil total phosphorus across the entire Loess Plateau region of China. Geoderma 197–198, 67–78.
Crossref
 
Muhr GR, Datta NP, Shankara SN, Dever F, Lecy VK Donahue RR (1963). Soil Testing in India. USAID Mission to India.
 
Nayak AK, Chinchmalatpure AR, Rao G, Khandelwal MK, Tyagi NK (2006). Spatial Variability of DTPA Extractable Micronutrients in the Soils of Bara Tract of Sardar Sarovar Canal Command in Gujarat. J. Indian Soc.Soil Sci.54(2):137-145.
 
Nayanaka VGD, Vitharana WAU, Mapa RB (2010). Geo-statistical Analysis of Soil Properties to Support Spatial Sampling in a Paddy Growing Alfisol. Trop. Agric. Res. 22(1):34–44.
 
Olsen SR, Cole CV, Vatanabe FS, Dean LA (1954). Estimation of available phosphorus in soil by extraction with sodium bicarbonate. Circ. U.S Dept. Agric. 939:1-19.
 
Page T, Haygarth PM, Beven KJ, Jones A, Butler T, Keeler C, Freer J, Owens PN, Wood GA (2005). The spatial variability of soil phosphorus GA in relation to topographic indices and important source areas: samples to assess the risks to water quality. J. Environ. Qual. 34: 2263-2277.
Crossref
 
Parker FW, Nelson WL, Winter Eric, Miller LE (1951). The broad interpretation of soil test informations. Agron. J. 43:105-112.
Crossref
 
Qu M, Li W, Zhang C (2014). Country-Scale Spatial Variability of Macronutrient Availability Ratios in Paddy Soils. Appl. Environ. Soil Sci. 3:1-10.
Crossref
 
Ruffo ML, Bollero GA, Hoeft RG, Bullock DG (2005). Spatial variability of the Illinois soil nitrogen test: Implications for soil sampling. Agron. J. 97:1485-1492.
Crossref
 
Santra P, Chopra UK, Chakraborty D (2008). Spatial variability of soil properties and its application in predicting surface map of hydraulic parameters in an agricultural farm. Curr. Sci. 95(7):937-945.
 
Sanyal SK, Dwivedi BS, Singh VK, Majumdar K, Datta SC, Pattanayak SK, Annapurna K (2015). Phosphorus in relation to dominant cropping sequences in India: chemistry, fertility relations and management options. Curr. Sci. 108:7-10
 
Shi Z, Li Y, Makeschine F, Wang RC (2005). Assessment of temporal and spatial variability of soil salinity in a coastal saline field. Environ Geo. 48:171-178.
 
Singh D, Chhonkar PK, Dwivedi BS (2007). Manual on Soil, Plant and Water Analysis Westville Publishing House, New Delhi.
 
Šnajdr J, Valášková V, Merhautová V, Herinková J, Cajthaml T, Baldrian P (2008). Spatial variability of enzyme activies and microbial biomass in the upper layers of Quercus petraea forest soil. Soil Biol. Biochem. 40:2068–2075.
Crossref
 
Subbiah BV, Asija GL (1956). A rapid procedure for the determination of available nitrogen in soils. Curr. Sci. 25:259-260.
 
Tesfahunegn GB, Tamene L, Vlek Paul LG (2011). Catchment-scale spatial variability of soil properties and implications on site-specific soil management in northern Ethiopia. Soil Till. Res. 117:124–139.
Crossref
 
Voltz M, Webster R (1990). A comparison of kriging, cubic splines and classification for predicting soil properties from sample information. J. Soil Sci. 41:473–490.
Crossref
 
Wang YQ, Zhang XC, Zhang JL, Li SJ (2009). Spatial variability of soil organic carbon in a watershed on the Loess Plateau. Pedosphere 19:486-495.
Crossref
 
Warrick AW, Myers DE, Nielsen DR (1986). Geo-statistical methods applied to soil science. In: Methods of Soil Analysis, Part 1. Physical and Mineralogical Methods. Agronomy Monograph 9(2):53–82.
 
Wasiullah AU, Bhatti KF, Akmal M (2010). Spatial variability and geo-statistics application for mapping of soil properties and nutrients in semi arid district Kohat of Khyber Pakhtunkhwa (Pakistan). Soil Environ. 29(2):159-166.
 
Webster R, Oliver MA (2001). Statistical Methods in Soil Science and Resource Survey. Oxford University Press, New York, NY.
 
Weindorf DC, Zhu Y (2010). Spatial variability of soil properties at Capulin Volcano, Mew Mexico, USA: implications for sampling strategy. Pedosphere 20:185–197.
Crossref
 
Yang Y, Yang J (2005). The trend variability of soil organic matter content in the salinity region of Yucheng city in Shandong Province. Chin. J. Soil Sci. 36(5):647-651.
 
Zuo X, Zhao H, Zhao X, Zhang T, Guo Y, Wang S, Drake S (2008). Spatial pattern and heterogeneity of soil properties in sand dunes under grazing and restoration in Horqin Sandy Land. Northern China. Soil Till. Res. 99:202–212.
Crossref

 




Copyright © 2024 Author(s) retain the copyright of this article.

This article is published under the terms of the Creative Commons Attribution License 4.0

Back to Vol. 10 No. 34
Back to articles
SUBSCRIBE TO RSS
SUBMIT MANUSCRIPT
CONFERENCES
          */?>