O objetivo do presente trabalho foi estudar a variabilidade espacial dos parâmetros da equação da condutividade hidráulica determinada no campo em função do conteúdo de água no solo, pelo método do perfil instantâneo, e avaliar a estabilidade temporal da armazenagem da água no solo, com vistas a identificação de um local de amostragem ótimo para a determinação da densidade de fluxo da água pela equação de Darcy-Buckingham. O estudo foi desenvolvido num Latossolo Vermelho Amarelo, textura areno-argiloso, em área experimental do Campus "Luiz de Queiroz", Universidade de São Paulo, município de Piracicaba, Estado de São Paulo, Brasil. Suas coordenadas geográficas são: 22° 42’ 43,3’’ de latitude sul, 47° 37’ 10,4’’ de longitude oeste e 456 m de altitude. A parcela experimental apresentava um comprimento de 45 m e uma largura de 15 m, na qual foram instalados 40 tubos de alumínio para acesso de uma sonda de nêutrons para medida do conteúdo da água no solo nas profundidades 0,20 m; 0,40 m; 0,60 m; 0,80 m e 1,00 m e cálculo da armazenagem da água no perfil 0- 1,0 m. A distribuição desses tubos foi feita na forma de uma grade de quatro colunas por dez linhas, com cada tubo distando de seu vizinho de 5 m. As funções K(θ) nos 40 pontos foram determinadas a partir das análises de regressão de θ em função de lnt e h_z em função de lnt, durante o processo de redistribuição da água no solo. De posse das 40 funções K(θ) nas cinco profundidades estudadas, foram aplicadas as técnicas geoestatitisticas para o estudo da variabilidade espacial dos parâmetros θ₀, K₀ e γ da função. Foi realizada também a análise da estabilidade temporal da armazenagem da água no solo, ao longo do período de redistribuição da água no solo para verificar quais e quantos são os locais adequados para o monitoramento da água precisão aceitável e reduzido esforço amostral. Pelos resultados obtidos, pôde-se concluir que: a) os métodos geoestatísticos utilizados foram adequados para descrever a estrutura de dependência espacial dos parâmetros da relação da condutividade hidráulica em função da umidade do solo: θ₀ , K₀ e γ; b) para as três variáveis estudadas ( θ₀ , K₀ e γ) houve estrutura de dependência espacial dos parâmetros na escala utilizada com alcances máximos de 8,33, 19,09 e 10,35 m; c) a técnica da estabilidade temporal que possibilita identificar, no campo, o ponto, ou os pontos, que representam a média e os pontos que superestimam ou subestimam a média real de determinada variável, identificou, para os valores de armazenagem da água no solo do presente estudo, os pontos 19, 20 e 21 como os mais representativos da média geral do campo; d) a metodologia com base nos coeficientes de correlação de Spearman também permitiu concluir que os valores de armazenagem da água no solo, foram estáveis no tempo, para os 40 pontos amostrados. Palavras-chave: perfil instantâneo, estabilidade temporal, sonda de nêutrons, armazenagem

The objective of this work was to study the spatial variability of the parameters of the hydraulic conductivity equation determined in the field as a function of soil water content, by the instantaneous profile method, and to evaluate the soil water storage time stability, in order to identify a sampling location for the determination of soil water flux by means of Darcy-Buckingham equation. The study was carried out in a clay sandy Oxisol in the country of Piracicaba, State of São Paulo, Brazil (22º42’ 43,3’’ S, 47° 37’ 10,4’’ W, 546 m). The dimensions of the experimental plot were 45 m x 15 m in wich 40 aluminium tubes were installed to acess a neutron probe to measure the soil water content at the depths of 0.2, 0.4, 0.6, 0.8 and 1.0 m and, then, calculate the soil water storage of the 0 - 1.0 m soil layer. The distribution of these tubes were made in grid of four columns by ten lines in spacing of 5m x 5m. The Kθ functions were determined in the 40 points from regression analyses of θ as function lnt and h_z as a function of lnt, being K the hydraulic conductivity, θ the volumetric soil water content, z h the soil water storage in the 0 - Z m layer and t the soil water redistribution time. With the 40 Kθ functions in the five studied soil depths, geostatistical techniques were used to evaluate the spatial variability of the parameters 0 θ , 0 K and γ of the function ( ) [ ] 0 0 exp θ θ γ − = K K , where the subindex o means lnt = 0. Soil water storage time stability analysis were also carried out along the soil water redistribution period in order to verify which and how many are the adequate locations for water monitoring with acceptable precisi on and reduced sample effort. From the results, the follwing could conclude: (i) the used geostatistical methods were adequate to described the spatial dependence structure of the parameters θ₀ , Kθ₀ and γ; (ii) there were a spatial dependence structure of the parameters in the used scale with maximum ranges of 8.33, 19.09 and 10.35 m; (iii) from the time stability technique, it was possible to identify the points 19,20 and 21 as the more representative of the field overall mean; (iv) the rank Spearmam coefficients also showed that the soil water storage of the 40 sampled points were in the time.