O estudo da variação espacial da fertilidade do solo cresceu em importância, em discussões sobre manejo agrícola, uma vez que o tratamento uniforme das áreas cultivadas gera resultados inconsistentes, capazes de alterar o rendimento das culturas. Utilizando-se métodos de geoestatística multivariada, objetivou-se avaliar a relação espacial entre o estoque de nutrientes e a densidade de um Latossolo Vermelho-Amarelo húmico. O estudo foi realizado em área cultivada com Coffea arabica L. (cv. Catucaí), em uma malha amostral com 100 pontos. A amostragem do solo foi realizada na camada de 0-0,20 m, para determinar a densidade e os teores de P, K, Ca e Mg, os quais foram transformados em valores de estoque. Após esta transformação, os dados foram analisados por meio de estatística descritiva e exploratória, determinando-se a correlação entre a densidade e os estoques. Posteriormente, realizou-se análise geoestatística uni e multivariada, seguidas de interpolação por krigagem ordinária e cokrigagem. A distribuição dos estoques de nutrientes mostrou-se inversamente proporcional à da densidade do solo, com exceção do estoque de P. A cokrigagem permitiu caracterizar o solo, favorecendo a interpretação do comportamento de distribuição espacial da fertilidade.
The study of the soil fertility spatial variation has grown in importance, in discussions on agricultural management, once the uniform approach of cultivated areas provides inconsistent results, capable of affecting the crops yield. By using multivariate geostatistics methods, this study aimed at evaluating the spatial relationship between the stock of nutrients and density of a humic Rhodic Hapludox. The study was conducted in an area cultivated with Coffea arabica L. (Catucaí cultivar), in a sample grid with 100 points. The soil sampling was performed in the 0-0.20 m layer, in order to determine the soil density and the P, K, Ca and Mg contents, which were converted into stock values. After this transformation, the data were analyzed by using descriptive and exploratory statistics, determining the correlation between density and stocks. Then, the univariate and multivariate geostatistics were used, followed by interpolation with ordinary kriging and cokriging. The distribution of nutrients stocks was inversely proportional to the soil density, except for the P stock. The cokriging allowed to characterize the soil, favoring the interpretation of the fertility spatial distribution behavior.