Worldwide there is considerable concern over the effects of human activities on the quantity and quality of freshwater. Measurement of infiltration behaviour will be important for improving freshwater management. This study identifies that New Zealand has a sporadic history of measuring soil water movement attributes on a limited number of soil types, although the current practical demand should be large for management of irrigation, dairy farm effluent disposal, as well as municipal / domestic waste- and storm-water disposal. Previous research has demonstrated that infiltration behaviour is governed by the interplay between numerous mechanisms including hydrophobicity and preferential flow, the latter being an important mechanism of contaminant leaching for many NZ soils. Future characterisation will need to recognise the dynamic nature of these interactions, and be able to reliably characterise the key infiltration mechanisms. Since macropores are responsible for preferential flow, it is critical that infiltration studies use a representative sample of the macropore network. The aim of this project was to study the mechanisms governing the infiltration behaviour of a layered soil in large (50 x 70 cm) monolith lysimeters, where the connectivity of the macropore network remains undisturbed. Four lysimeters of the Gorge silt loam were collected, a structured soil with four distinct layers. On each lysimeter there were four separate infiltration experiments, with water applied under suctions of 0, 0.5, 1, and 1.5 kPa by a custom-built tension infiltrometer. Each lysimeter was instrumented with 30 tensiometers, located in arrays at the layer boundaries. There was also a field experiment using ponded dye infiltration to visually define preferential flowpaths. Analysis of dye patterns, temporal variability in soil matric potential (Ψm), and solute breakthrough curves all show that preferential flow is an important infiltration mechanism. Preferential flowpaths were activated when Ψm was above -1.5 kPa. During saturated infiltration, at least 97% of drainage was through the ‘mobile’ pore volume of the lysimeter (θm), estimated among the lysimeters at 5.4 – 8.7 % of the lysimeter volume. Early-time infiltration behaviour did not show the classical square-root of time behaviour, indicating sorptivity was not the governing mechanism. This was consistent across the four lysimeters, and during infiltration under different surface imposed suctions. The most likely mechanism restricting sorptivity is weak hydrophobicity, which appears to restrict infiltration for the first 5 – 10 mm of infiltration. Overall, the Gorge soil’s early-time infiltration behaviour is governed by the dynamic interaction between sorptivity, hydrophobicity, the network of air-filled pores, preferential flow and air encapsulation. Long-time infiltration behaviour was intimately linked to the temporal dynamics of Ψm, which was in turn controlled by preferential flow and soil layer interactions. Preferential flowpaths created strong inter-layer connectivity by allowing an irregular wetting front to reach lower layers within 2 – 15 mm of infiltration. Thereafter, layer interactions dominate infiltration for long-time periods, as Ψm in soil layers with different K(Ψm) relationships self-adjusts to try to maintain a constant Darcy velocity. An important finding was that Ψm rarely attained the value set by the tension infiltrometer during unsaturated infiltration. The results show that ‘true’ steady-state infiltration is unlikely to occur in layered soils. A quasi-steady state was identified once the whole column had fully wet and layer interactions had settled to where Ψm changes occurred in unison through each soil layer. Quasi-steady state was difficult to identify from just the cumulative infiltration curve, but more robustly identified as when infiltration matched drainage, and Ψm measurements showed each layer had a stable hydraulic gradient. I conclude that the in-situ hydraulic conductivity, K(Ψm), of individual soil layers can be accurately and meaningfully determined from lysimeter-scale infiltration experiments. My results show that K(Ψm) is different for each soil layer, and that differences are consistent among the four lysimeters. Under saturated flow the subsoil had the lowest conductivity, and was the restricting layer. Most interestingly this pattern reversed during unsaturated flow. As Ψm decreased below -0.5 to -1 kPa, the subsoil was markedly more conductive, and the topsoil layers became the restricting layers. All four soil layers demonstrate a sharp decline in K(Ψm) as Ψm decreases, with a break in slope at ~ -1 kPa indicating the dual-permeability nature of all layers.

The present work is a contribution to description and understanding of the distribution and movement of water in swelling soils. In order to investigate the moisture distribution in swelling soils a detailed knowledge of volume change properties, flow characteristics and total potential of water in the soil is essential. Therefore, a possible volume change mechanism is first described by dividing the swelling soils into four categories and volume change of a swelling soil is measured under different overburden pressures. The measured and calculated (from volume change data) overburden potential components are used to check the validity of the derivation of a load factor, ∝. Moisture diffusivity in swelling soil under different overburden pressures is measured using Gardner's (1956) outflow method. Behaviour of equilibrium moisture profiles in swelling soils is theoretically explained, solving the differential equation by considering the physical variation of individual soil properties with moisture content and overburden pressure. Using the measured volume change data and moisture potentials under various overburden pressures, the behaviour of possible moisture profiles are described at equilibrium and under steady vertical flows in swelling soils. It is shown that high overburden pressures lead to soil water behaviour quite different from any previously reported.

This research work concentrate on developing digital soil maps to support field based plant phenotyping research. We have developed soil organic matter content (OM), cation exchange capacity (CEC), natural soil drainage class, and tile drainage line maps using topographic indices and aerial imagery. Various prediction models (universal kriging, cubist, random forest, C5.0, artificial neural network, and multinomial logistic regression) were used to estimate the soil properties of interest.