ABSTRACT In precision agriculture, determining management zones for soil and plant attributes is a complex process that requires knowledge of several variables, which complicates management and decisionmaking processes. This study evaluated the spatial variability of soybean yield and soil chemical properties using geostatistical and multivariate analyses to define management zones in an Oxisol. The soybean yield and soil chemical properties between 0 to 0.2 and 0.2 to 0.4 m soil depths were sampled at 70 points. Geostatistical and multivariate analyses were then performed on these data. The soil chemical properties showed higher variability at 0.2 to 0.4 m soil depth. The semivariogram parameters of the principal component analysis (PCA) data (PCA 1, PCA 2, and PCA 3) for both depths were more homogeneous than the original data. The maps of soil chemical properties showed high similarity to the soybean yield map. The PCA explained 65.34% (0 to 0.2 m) and 70.50% (0.2 to 0.4 m) of data variability, grouping the soybean yield, organic matter, pH, phosphorous, potassium, calcium, magnesium, and sodium. PCA spatialization allowed for the definition of management zones indicated by PCA 1, PCA 2, and PCA 3 for both depths. The result indicates that the area must be managed using different strategies of soil fertility management to increase soybean yield.