No Brasil pouco se estuda sobre as estratégias de amostragem em áreas com mais de um tipo de solo e várias feições de paisagem, estratégias estas, importantes para mapear a variabilidade espacial dos atributos físicos: frações granulométricas, como também a correlação destas com os parâmetros geomorfológicos: altitude e inclinação máxima da área amostrada da extensão multivariada da krigagem, a cokrigagem. Essa extensão multivariada permite a associação de duas variáveis dependentes gerarem novas informações correlacionadas através do modelo de variância cruzada. O objetivo deste trabalho foi avaliar a correlação espacial entre a distribuição da fração granulométrica com os parâmetros de altitude e a inclinação máxima através de uma densidade amostral agregada (100x100 e 50x50m) com dois níveis de densidade amostral, 20x20 e 40x40m. O estudo se deu numa área de 180 ha, onde se coletou frações de solo numa grade amostral agregada, ou seja, um arranjo amostral 100x100m e três sub-grades de 50x50m distribuídas aleatoriamente na área. Foi realizada análise estatística descritiva dos dados processados no laboratório e a análise espacial através de modelos de variância. De posse dos modelos ajustados na análise dos variogramas, realizou a interpolação pela krigagem para as frações granulométricas e através dos variogramas cruzados pela cokrigagem, as correlações dessas frações com os dados, obtidos pelo modelo digital de elevação, da altitude e a inclinação máxima do terreno da área e comparada, através dos mapas, essas mesmas correlações em densidades amostrais diferentes. A partir dos modelos de mapas gerados, chegou aos seguintes resultados: houve uma moderada variabilidade das frações granulométricas e uma moderada dependência espacial. De acordo com os parâmetros da variância cruzada a distribuição espacial da fração areia correlacionada com a altitude e com a inclinação máxima não diferenciou significativamente em relação aos níveis de densidade amostral 20x20m e 40x40m, praticamente não houve correlação entre a fração areia e a altitude, pois não apresentou um padrão definido de distribuição espacial. Pode-se concluir que a declividade influenciou na distribuição dos atributos físicos do solo. Como também foi avaliada a disponibilidade hídrica na área do pivô da Fazenda, foi constatada a forte correlação do atributo argila+silte com a capacidade de água disponível (CAD), enfatizando a dependência da disponibilidade hídrica em relação aos atributos físicos do solo.

In Brazil, few studied are related about soil sampling strategies in areas with more than one type soil and various features of the landscape, these strategies are important to map the spatial variability of physical attributes: size fractions, as well as the correlation of these with geomorphologic parameters: altitude and maximum slope of the sampled area of the multivariate extension of kriging, cokriging. This extension allows multivariate combination of two dependent variables to generate new information correlated with the model cross variance. The objective of this study was to evaluate the spatial correlation between the distribution of particle size fraction with the parameters of altitude and terrain slope through a maximum aggregate sample density (100x100 and 50x50m) with two levels of sample density, 20x20 and 40x40m. The study was an area of 180 ha, where fractions of collected soil at a sampling grid aggregate, ie an arrangement 100x100m sample and three sub-grids 50x50m randomly distributed in the area. We performed a descriptive statistical analysis of data processed in the laboratory and spatial analysis through models of variance. In possession of the adjusted models in the analysis of variogram for kriging interpolation performed for the fractions and through the variogram crossed by cokriging, the correlations of these fractions with the data obtained by digital elevation model, altitude and the maximum slope of the land area and compared, through maps, these same correlations in different sampling densities. From the models generated maps, it was reached the following results: there was a moderate variability in size fractions and a moderate spatial dependence. In accordance with parameters of cross (semi)variance the spatial distribution of the sand fraction correlated with altitude and the maximum slope did not differ significantly from the levels of sample density 20x20m and 40x40m, virtually no correlation between the sand fraction and altitude because it did not show a definite pattern of spatial distribution. It can be concluded that the slope influenced the distribution of soil physical attributes studieds. It was also quantified the water availability in the area under the center pivot equipment, it was found a strong correlation of the attribute clay + silt with TAW (total available water), emphasizing the dependence on water availability in relation to soil physical properties.