With reference to a more compacted and less conductive upper soil layer overlying a less compacted and more conductive subsoil, a simple three-dimensional (3D) infiltration run is expected to yield more representative results of the upper layer than the subsoil. However, there is the need to quantitatively establish what is meant by more representativeness. At this aim, numerically simulated infiltration was investigated for a theoretically unconfined process under a null ponded head of water (d0H0 setup, with d = depth of ring insertion and H = ponded depth of water) and a practical beerkan run (d1H1 setup, d = H = 1 cm). The considered layered soils differed by both the layering degree (from weak to strong) and the thickness of the upper soil layer (0.5-3 cm). It was confirmed that water infiltration should be expected to be more representative of the upper soil layer when this layer is the less permeable since, for a 2-h experiment, the instantaneous infiltration rates for the layered soil were 1.0-2.1 times greater than those of the homogeneous low permeable soil and 1.3-20.7 smaller than those of the homogeneous coarser soil that constituted the subsoil. Similarity with the homogeneous fine soil increased as expected as the upper layer became thicker. For a weak layering condition, the layered soil yielded an intermediate infiltration as compared with that of the two homogeneous soils forming the layered system. For a strong layering degree, the layered soil was more similar to the homogeneous fine soil than to the homogeneous coarse soil. Using the practical setup instead of the theoretical one should have a small to moderate effect on the instantaneous infiltration rates since all the calculated percentage differences between the d1H1 and d0H0 setups fell into the relatively narrow range of -18.8% to +17.4%. A sequential analysis procedure appeared usable to detect layering conditions but with some modifications as compared with the originally proposed procedure. The practical setup enhanced the possibility to recognize the time at which the characteristics of the subsoil start to influence the infiltration process. In conclusion, this investigation contributed to better interpret both the theoretical and the practically established 3D infiltration process in a soil composed of a less conductive upper soil layer overlying a more conductive subsoil and it also demonstrated that modifying a recently proposed procedure only using infiltration data could be advisable to determine the time when layering starts to influence the process.

In his seminal paper on solution of the infiltration equation, Philip (1957) proposed a gravity time, tgrav, to estimate practical convergence time of his infinite time series expansion, TSE. The parameter tgrav refers to a point in time where infiltration is dominated equally by capillarity and gravity derived from the first two (dominant) terms of the TSE expansion. Evidence that higher order TSE terms describe the infiltration process better for longer times. Since the conceptual definition of tgrav is valid regardless of the infiltration model used, we opted to reformulate tgrav using the analytic approximation proposed by Parlange et al. (1982) valid for all times. In addition to the roles of soil sorptivity (S) and saturated (Ks) and initial (Ki) hydraulic conductivities, we explored effects of a soil specific shape parameter β on the behavior of tgrav. We show that the reformulated tgrav (notably tgrav= F(β) S^2/(Ks - Ki)^2 where F(β) is a β-dependent function) is about 3 times larger than the classical tgrav given by tgrav, Philip= S^2/(Ks - Ki)^2. The differences between original tgrav, Philip and the revised tgrav increase for fine textured soils. Results show that the proposed tgrav is a better indicator for convergence time than tgrav, Philip. For attainment of the steady-state infiltration, both time parameters are suitable for coarse-textured soils, but not for fine-textured soils for which tgrav is too conservative and tgrav, Philip too short. Using tgrav will improve predictions of the soil hydraulic parameters (particularly Ks) from infiltration data as compared to tgrav, Philip.

Lateral saturated soil hydraulic conductivity, Ks,l, is the soil property governing subsurface water transfer in hillslopes, and the key parameter in many numerical models simulating hydrological processes both at the hillslope and catchment scales. Likewise, the hydrological connectivity of lateral flow paths plays a significant role in determining the intensity of the subsurface flow at various spatial scales. The objective of the study is to investigate the relationship between Ks,l and hydraulic connectivity at the hillslope spatial scale. Ks,l was determined by the subsurface flow rates intercepted by drains, and by water table depths observed in a well network. Hydraulic connectivity of the lateral flow paths was evaluated by the synchronicity among piezometric peaks, and between the latter and the peaks of drained flow. Soil moisture and precipitation data were used to investigate the influence of the transient hydrological soil condition on connectivity and Ks,l. It was found that the higher was the synchronicity of the water table response between wells, the lower was the time lag between the peaks of water levels and those of the drained subsurface flow. Moreover, the most synchronic water table rises determined the highest drainage rates. The relationships between Ks,l and water table depths were highly non-linear, with a sharp increase of the values for water table levels close to the soil surface. Estimated Ks,l values for the full saturated soil were in the order of thousands of mm h-1, suggesting the activation of macropores in the root zone. The Ks,l values determined at the peak of the drainage events were correlated with the indicators of synchronicity. The sum of the antecedent soil moisture and of the precipitation was correlated with the indicators of connectivity and with Ks,l. We suggest that, for simulating realistic processes at the hillslope scale, the hydraulic connectivity could be implicitly considered in hydrological modelling through an evaluation of Ks,l at the same spatial scale.