This work aimed to study the spatial variability of clay and total sand size fractions of Oxisol cultivated with eucalyptus. Soil samples were collected at depths of 0-0.20 m and 0.20-0.40 m in a sampling grid of 33 x 33 m, totaling 94 georeferenced points. Data were subjected to descriptive statistics, geostatistics and kriging and cokriging analyzes. In cokriging depth values of 0-0,20 m were used to estimate depth values of 0.20-0.40 m. Highest clay fraction mean value was found at a depth of 0.20-0.40 m and for total sand at a depth of 0-0.20 m. The spherical model fitted to the data in the methods of kriging and cokriging. The higher ranges of spatial dependence were found at a depth of 0.20-0.40 m. The mean errors found for the two methods showed similar and negative values indicating superstimation of sample values.
Geostatistic; Spatial Variability; Multivariate
