Saturated and near-saturated surface soil hydraulic properties influence the partition of rainfall and snowmelt into infiltration and runoff. The goal of this study was to characterize near-saturated surface soil hydraulic properties and water-conducting porosity in sloping landscapes. The specific objectives included exploration of tension and double-ring infiltrometers for estimation of soil hydraulic properties in sloping landscapes, development of an improved method for determining water-conducting porosity, and the application of these methods in characterizing soil hydraulic properties and water-conducting porosity under three land use. Water infiltration from a double-ring infiltrometer and a tension infiltrometer at water pressures between -2.2 and -0.3 kPa was measured in a cultivated field with 0, 7, 15, and 20% slopes at Laura and under three land use (native grass, brome grass and cultivated) at St. Denis in Saskatchewan, Canada. Three-dimensional computer simulation studies were also performed for tension infiltrometer with various disc diameters, water pressures, and surface slopes. Steady infiltration rates and estimated field-saturated hydraulic conductivity (Kfs), hydraulic conductivity-water pressure relationship (K(h)), and inverse capillary length parameter were compared for different slopes and land use. These parameters were not significantly different (p<0.05) among slopes. For specific K(h) functions, a new analytical solution was developed and compared with existing methods for calculating water-conducting porosity. The new method reliably determined water-conducting porosity of surface soils and gave consistent results, regardless of the width of water pressure ranges. At the -0.3 kPa water pressure, hydraulic conductivity of grasslands was two to three times greater than the cultivated lands. Values of inverse capillary length parameter were about two times and values of Kfs about four times greater in grasslands than in cultivated fields. Water-conducting macroporosity of grasslands and cultivated fields were 0.04% and 0.01% of the total soil volume, respectively. Over 40% and 50% of the total water flux at -0.06 kPa water pressure was transmitted through macropores (pores > 1×10-3 m in diameter) of the cultivated land and the grasslands, respectively. Experimental and simulation results of this study indicated that both tension and double-ring infiltrometers are suitable for characterization of saturated and near-saturated surface soil hydraulic properties in landscapes up to 20% slope. The new method can be used to characterize water-conducting porosity from in situ tension and double-ring infiltrometers measurements more adequately and efficiently than the existing methods. Application of these methods for three land use indicated that land use modified surface soil hydraulic properties and consequently may alter the water balance of an area by affecting the partition between, and relative amount of infiltration and surface runoff.

Saturated and near-saturated surface soil hydraulic properties influence the partition of rainfall and snowmelt into infiltration and runoff. The goal of this study was to characterize near-saturated surface soil hydraulic properties and water-conducting porosity in sloping landscapes. The specific objectives included exploration of tension and double-ring infiltrometers for estimation of soil hydraulic properties in sloping landscapes, development of an improved method for determining water-conducting porosity, and the application of these methods in characterizing soil hydraulic properties and water-conducting porosity under three land use. Water infiltration from a double-ring infiltrometer and a tension infiltrometer at water pressures between -2.2 and -0.3 kPa was measured in a cultivated field with 0, 7, 15, and 20% slopes at Laura and under three land use (native grass, brome grass and cultivated) at St. Denis in Saskatchewan, Canada. Three-dimensional computer simulation studies were also performed for tension infiltrometer with various disc diameters, water pressures, and surface slopes. Steady infiltration rates and estimated field-saturated hydraulic conductivity (Kfs), hydraulic conductivity-water pressure relationship (K(h)), and inverse capillary length parameter were compared for different slopes and land use. These parameters were not significantly different (p<0.05) among slopes. For specific K(h) functions, a new analytical solution was developed and compared with existing methods for calculating water-conducting porosity. The new method reliably determined water-conducting porosity of surface soils and gave consistent results, regardless of the width of water pressure ranges. At the -0.3 kPa water pressure, hydraulic conductivity of grasslands was two to three times greater than the cultivated lands. Values of inverse capillary length parameter were about two times and values of Kfs about four times greater in grasslands than in cultivated fields. Water-conducting macroporosity of grasslands and cultivated fields were 0.04% and 0.01% of the total soil volume, respectively. Over 40% and 50% of the total water flux at -0.06 kPa water pressure was transmitted through macropores (pores > 1×10-3 m in diameter) of the cultivated land and the grasslands, respectively. Experimental and simulation results of this study indicated that both tension and double-ring infiltrometers are suitable for characterization of saturated and near-saturated surface soil hydraulic properties in landscapes up to 20% slope. The new method can be used to characterize water-conducting porosity from in situ tension and double-ring infiltrometers measurements more adequately and efficiently than the existing methods. Application of these methods for three land use indicated that land use modified surface soil hydraulic properties and consequently may alter the water balance of an area by affecting the partition between, and relative amount of infiltration and surface runoff.

No description available.

Conduziu-se durante dois anos, em um Latossolo Vermelho Amarelo, em Piracicaba-SP (22o41’00&quot; S; 47o39’00&quot; O; 554 m de altitude), um experimento com a cultura da cana-de-açúcar. Foram utilizados 9 tratamentos com 4 repetições, implantados num delineamento em blocos ao acaso. Os tratamentos foram: tr1 - adubação mineral + calagem; tr2, tr3; tr4 e tr5 - doses de lodo de esgoto (respectivamente, 0; 33; 66 e 99 Mg ha-1, no primeiro ano, e 0; 37; 74 e 112 Mg ha-1, no segundo ano); tr6, tr7; tr8 e tr9 - doses de composto de lixo (respectivamente, 0; 20; 40 e 60 Mg ha-1, no primeiro ano, e 0; 24; 48 e 72 Mg ha-1, no segundo ano). Os resíduos foram aplicados manualmente na superfície do solo e incorporados com enxada rotativa, no primeiro ano em área total e no segundo apenas na entrelinha. Foram determinados, em cada parcela, os seguintes parâmetros do solo: composição granulométrica (amostras coletadas antes da implantação do experimento); o conteúdo de carbono orgânico (amostras coletadas em 15 datas durante os dois anos); conteúdo de argila dispersa em água, densidade e porosidade total (amostras coletadas ao final do primeiro e do segundo ano do ciclo da cultura); curva de retenção de água, condutividade hidráulica do solo saturado e condutividade hidráulica relativa (amostras coletadas ao final do primeiro ano do ciclo da cultura); e condutividade hidráulica do solo saturado e não saturado (determinada no campo ao final de dois anos do ciclo da cultura, utilizando o infiltrômetro de tensão). Os resultados obtidos mostraram que houve redução da densidade do solo e aumento da porosidade total com o aumento das doses dos resíduos, para a camada de 0 a 0,15 m de profundidade, no segundo ano de aplicação. Nas parcelas que receberam aplicações de ambos os resíduos, houve aumento da condutividade hidráulica do solo próximo à saturação e redução à medida que o solo se afasta da saturação, para a faixa de potencial mátrico estudada (0 a 1 kPa). O decaimento da condutividade hidráulica relativa, a partir da saturação, foi inicialmente mais acentuado nas parcelas que receberam aplicações de ambos os resíduos, quando comparadas às que não receberam resíduos. Verificou-se também, através de análises de regressão múltipla, que diferenças na retenção de água e no conteúdo de argila dispersa em água, que pela análise de variância convencional seriam atribuídas exclusivamente aos tratamentos, de fato eram parcial ou totalmente devidas a variações não casuais na composição granulométrica do solo. Pelos resultados pôde-se concluir que a aplicação de resíduos urbanos ao solo leva a modificações de propriedades como a densidade e porosidade total do solo, a condutividade hidráulica do solo saturado e a condutividade hidráulica do solo não saturado. Pôde-se concluir também que a casualização dos tratamentos na área experimental não garante a ausência de covariação entre os tratamentos e outras variáveis independentes, as quais podem interferir nas propriedades de interesse (variáveis dependentes) e, conseqüentemente, na eficiência da análise de variância feita de maneira convencional, sendo recomendado a realização de um &quot;ensaio em branco&quot; para verificar a eficiência desta casualização. / An experiment was carried out during two years on a Red Yellow Latosol cropped with sugar-cane, located in the country of Piracicaba-SP (22o41’00&quot; S; 47o39’00&quot; W; 554 m high). Nine treatments with 4 replications were used in a random-block experimental design. The treatments were: tr1 - mineral fertilization + CaCO3; tr2, tr3; tr4 and tr5 - levels of sewage sludge (respectively, 0; 33; 66 and 99 Mg ha-1, in the first year, and 0; 37; 74 and 112 Mg ha-1, in the second year); tr6, tr7; tr8 and tr9 - levels of composed of garbage (respectively, 0; 20; 40 and 60 Mg ha-1, in the first year, and 0; 24; 48 and 72 Mg ha-1, in the second year). The residues were manually applied on the soil surface and incorporated with a rotative plough. In the first year, the incorporation was done on the total area; in the second year only in inter-rows. The following soil parameters were determinated in each plot: particle size distribution (samples collected before the installation of the experiment); organic carbon content (15 samplings along the two years); natural clay content, soil bulk density and total soil porosity (sampling at the end of the first and second years of the growing cycle); water retention curves, saturated hydraulic conductivity and relative unsaturated hydraulic conductivity (samples collected at the end of the first year of the growing cycle); and saturated and unsaturated hydraulic conductivity (determined in situ at the end of the first and second years of the growing cycle, using tension infiltrometer). Results show a decrease the soil bulk density and increase of total porosity as the levels of residue application increased, for the 0 to 0.15 m soil layer, in the second year of application. On the plots that received applications of both residues, there was increase of the hydraulic conductivity near the soil saturation and decrease for the 0 to 1 kPa matric potential range. The decline of the relative hydraulic conductivity starting from the saturation, was initially more accentuated in the plots that received applications of both residues as compared with the ones that did not. It was also verified, through analyses of multiple regression, that differences in the water retention and in the natural clay content, that they would be attributed exclusively to the treatments by the conventional variance analysis, they were partially or totally owed to non casual variations in the particle size distribution of the soil. It could be concluded that the soil application of urban residues causes modifications of soil physical properties as saturated and unsaturated soil hydraulic conductivity. It could also be concluded that the randomization of the treatments in the experimental area does not guarantee the covariance absence between the treatments and other independent variables, which can interfere in the properties of interest (dependent variables) and, consequently, in the efficiency of the variance analysis done in a conventional way.

Conduziu-se durante dois anos, em um Latossolo Vermelho Amarelo, em Piracicaba-SP (22o41’00&quot; S; 47o39’00&quot; O; 554 m de altitude), um experimento com a cultura da cana-de-açúcar. Foram utilizados 9 tratamentos com 4 repetições, implantados num delineamento em blocos ao acaso. Os tratamentos foram: tr1 - adubação mineral + calagem; tr2, tr3; tr4 e tr5 - doses de lodo de esgoto (respectivamente, 0; 33; 66 e 99 Mg ha-1, no primeiro ano, e 0; 37; 74 e 112 Mg ha-1, no segundo ano); tr6, tr7; tr8 e tr9 - doses de composto de lixo (respectivamente, 0; 20; 40 e 60 Mg ha-1, no primeiro ano, e 0; 24; 48 e 72 Mg ha-1, no segundo ano). Os resíduos foram aplicados manualmente na superfície do solo e incorporados com enxada rotativa, no primeiro ano em área total e no segundo apenas na entrelinha. Foram determinados, em cada parcela, os seguintes parâmetros do solo: composição granulométrica (amostras coletadas antes da implantação do experimento); o conteúdo de carbono orgânico (amostras coletadas em 15 datas durante os dois anos); conteúdo de argila dispersa em água, densidade e porosidade total (amostras coletadas ao final do primeiro e do segundo ano do ciclo da cultura); curva de retenção de água, condutividade hidráulica do solo saturado e condutividade hidráulica relativa (amostras coletadas ao final do primeiro ano do ciclo da cultura); e condutividade hidráulica do solo saturado e não saturado (determinada no campo ao final de dois anos do ciclo da cultura, utilizando o infiltrômetro de tensão). Os resultados obtidos mostraram que houve redução da densidade do solo e aumento da porosidade total com o aumento das doses dos resíduos, para a camada de 0 a 0,15 m de profundidade, no segundo ano de aplicação. Nas parcelas que receberam aplicações de ambos os resíduos, houve aumento da condutividade hidráulica do solo próximo à saturação e redução à medida que o solo se afasta da saturação, para a faixa de potencial mátrico estudada (0 a 1 kPa). O decaimento da condutividade hidráulica relativa, a partir da saturação, foi inicialmente mais acentuado nas parcelas que receberam aplicações de ambos os resíduos, quando comparadas às que não receberam resíduos. Verificou-se também, através de análises de regressão múltipla, que diferenças na retenção de água e no conteúdo de argila dispersa em água, que pela análise de variância convencional seriam atribuídas exclusivamente aos tratamentos, de fato eram parcial ou totalmente devidas a variações não casuais na composição granulométrica do solo. Pelos resultados pôde-se concluir que a aplicação de resíduos urbanos ao solo leva a modificações de propriedades como a densidade e porosidade total do solo, a condutividade hidráulica do solo saturado e a condutividade hidráulica do solo não saturado. Pôde-se concluir também que a casualização dos tratamentos na área experimental não garante a ausência de covariação entre os tratamentos e outras variáveis independentes, as quais podem interferir nas propriedades de interesse (variáveis dependentes) e, conseqüentemente, na eficiência da análise de variância feita de maneira convencional, sendo recomendado a realização de um &quot;ensaio em branco&quot; para verificar a eficiência desta casualização. / An experiment was carried out during two years on a Red Yellow Latosol cropped with sugar-cane, located in the country of Piracicaba-SP (22o41’00&quot; S; 47o39’00&quot; W; 554 m high). Nine treatments with 4 replications were used in a random-block experimental design. The treatments were: tr1 – mineral fertilization + CaCO3; tr2, tr3; tr4 and tr5 - levels of sewage sludge (respectively, 0; 33; 66 and 99 Mg ha-1, in the first year, and 0; 37; 74 and 112 Mg ha-1, in the second year); tr6, tr7; tr8 and tr9 - levels of composed of garbage (respectively, 0; 20; 40 and 60 Mg ha-1, in the first year, and 0; 24; 48 and 72 Mg ha-1, in the second year). The residues were manually applied on the soil surface and incorporated with a rotative plough. In the first year, the incorporation was done on the total area; in the second year only in inter-rows. The following soil parameters were determinated in each plot: particle size distribution (samples collected before the installation of the experiment); organic carbon content (15 samplings along the two years); natural clay content, soil bulk density and total soil porosity (sampling at the end of the first and second years of the growing cycle); water retention curves, saturated hydraulic conductivity and relative unsaturated hydraulic conductivity (samples collected at the end of the first year of the growing cycle); and saturated and unsaturated hydraulic conductivity (determined in situ at the end of the first and second years of the growing cycle, using tension infiltrometer). Results show a decrease the soil bulk density and increase of total porosity as the levels of residue application increased, for the 0 to 0.15 m soil layer, in the second year of application. On the plots that received applications of both residues, there was increase of the hydraulic conductivity near the soil saturation and decrease for the 0 to 1 kPa matric potential range. The decline of the relative hydraulic conductivity starting from the saturation, was initially more accentuated in the plots that received applications of both residues as compared with the ones that did not. It was also verified, through analyses of multiple regression, that differences in the water retention and in the natural clay content, that they would be attributed exclusively to the treatments by the conventional variance analysis, they were partially or totally owed to non casual variations in the particle size distribution of the soil. It could be concluded that the soil application of urban residues causes modifications of soil physical properties as saturated and unsaturated soil hydraulic conductivity. It could also be concluded that the randomization of the treatments in the experimental area does not guarantee the covariance absence between the treatments and other independent variables, which can interfere in the properties of interest (dependent variables) and, consequently, in the efficiency of the variance analysis done in a conventional way.

The increasing commercial and industrial use of ethanol, i.e., in biofuel and gasoline, has generated increased incidents of vadose zone contamination by way of ethanol spills and releases. This has increased the interest in better understanding infiltration behaviours of ethanol in unsaturated porous media and the multiphase interactions in the vadose zone. Solute-dependent capillarity-induced focused flow (SCIFF) is a vertical, highly focused flow infiltration behaviour first reported by Smith et al. (2011) in butanol-contaminated sands. Through the use of highly controlled laboratory experiments, this thesis research investigates focused flow (SCIFF) and related behaviours through water infiltration into ethanol-contaminated unsaturated sand. Focused flow behaviours (SCIFF) were demonstrated through the infiltration of water into an ethanol-contaminated unsaturated sand using both constant flux and constant head methodologies. The observation of focused flow behaviours in ethanol-contaminated sand supported the primary hypothesis of this work. The secondary hypothesis was also supported, as focused flow behaviours were not observed, rather stable semicircular infiltration patterns were observed during ethanol infiltration into water-wet sand. Comparisons between constant flux and constant head application methods under similar flow rates and fluid volumes produced similar results. The zone of lower saturation, or the “halo effect” reported in previous literature, was strongly expressed during water infiltration in ethanol-contaminated sand. This halo effect is affected by the maximum (at 40% to 50%) of aqueous concentration of ethanol. This maximum enhances the zone of lower saturation and stabilizes the solute front. The SCIFF focused flow also overcame the effects of minor heterogeneities in the sand. However, additional laboratory and modelling work is required to further understand the extent of SCIFF behaviour. / Thesis / Master of Science (MSc) / Understanding the behaviour and interaction of water and contaminants in soils is important as environmental contamination and spills can have devastating environmental impacts. In recent decades, ethanol spills and accidental releases onto ground surface have increased as the commercial and industrial use of ethanol has increased. The goals of this work were to qualitatively visualize and quantify the unique nature of water infiltration into the ethanol-contaminated soil and understand the complex mechanisms behind water-ethanol interactions. This research showed that water infiltration creates an uncommon vertical, focused pattern when flowing into sand contaminated by ethanol. However, when ethanol is applied to standard water-wet sand, that behaviour is not observed. This work provided greater insight into the nature of ethanol-contaminated soils. These findings furthered the understanding needed to evaluate impacts that ethanol contamination can have on remedial efforts and the rate of migration of contaminants to groundwater.

Land application of wastes has become increasingly popular, to promote nutrient recycling and environmental protection, with soil functioning as a partial barrier between wastes and groundwater. Dairy shed effluent (DSE), may contain a wide variety of pathogenic micro-organisms, including bacteria (e.g. Salmonella paratyphyi, Escherichia coli. and Campylobacter), protozoa and viruses. Groundwater pathogen contamination resulting from land-applied DSE is drawing more attention with the intensified development of the dairy farm industry in New Zealand. The purpose of this research was to investigate the fate and transport of bacterial indicator-faecal coliform (FC) from land-applied DSE under different irrigation practices via field lysimeter studies, using two water irrigation methods (flood and sprinkler) with contrasting application rates, through the 2005-2006 irrigation season. It was aimed at better understanding, quantifying and modelling of the processes that govern the removal of microbes in intact soil columns, bridging the gap between previous theoretical research and general farm practices, specifically for Templeton soil. This study involved different approaches (leaching experiments, infiltrometer measurements and a dye infiltration study) to understand the processes of transient water flow and bacterial transport; and to extrapolate the relationships between bacterial transport and soil properties (like soil structure, texture), and soil physical status (soil water potential ψ and volumetric water content θ). Factors controlling FC transport are discussed. A contaminant transport model, HYDRUS-1D, was applied to simulate microbial transport through soil on the basis of measured datasets. This study was carried out at Lincoln University’s Centre for Soil and Environmental Quality (CSEQ) lysimeter site. Six lysimeters were employed in two trials. Each trial involved application of DSE, followed by a water irrigation sequence applied in a flux-controlled method. The soil columns were taken from the site of the new Lincoln University Dairy Farm, Lincoln, Canterbury. The soil type is Templeton fine sandy loam (Udic-Ustochrept, coarse loamy, mixed, mesic). Vertical profiles (at four depths) of θ and ψ were measured during leaching experiments. The leaching experiments directly measured concentrations of chemical tracer (Br⁻ or Cl⁻) and FC in drainage. Results showed that bacteria could readily penetrate through 700 mm deep soil columns, when facilitated by water flow. In the first (summer) trial, FC in leachate as high as 1.4×10⁶ cfu 100 mL⁻¹ (similar to the DSE concentration), was detected in one lysimeter that had a higher clay content in the topsoil, immediately after DSE application, and before any water irrigation. This indicates that DSE flowed through preferential flow paths without significant treatment or reduction in concentrations. The highest post-irrigation concentration was 3.4×10³ cfu 100 mL⁻¹ under flood irrigation. Flood irrigation resulted in more bacteria and Br⁻ leaching than spray irrigation. In both trials (summer and autumn) results showed significant differences between irrigation treatments in lysimeters sharing similar drainage class (moderate or moderately rapid). Leaching bacterial concentration was positively correlated with both θ and ψ, and sometimes drainage rate. Greater bacterial leaching was found in the one lysimeter with rapid whole-column effective hydraulic conductivity, Keff, for both flood and spray treatments. Occasionally, the effect of Keff on water movement and bacterial transport overrode the effect of irrigation. The ‘seasonal condition’ of the soil (including variation in initial water content) also influenced bacterial leaching, with less risk of leaching in autumn than in summer. A tension infiltrometer experiment measured hydraulic conductivity of the lysimeters at zero and 40 mm suction. The results showed in most cases a significant correlation between the proportion of bacteria leached and the flow contribution of the macropores. The higher the Ksat, the greater the amount of drainage and bacterial leaching obtained. This research also found that this technique may exclude the activity of some continuous macropores (e.g., cracks) due to the difference of initial wetness which could substantially change the conductivity and result in more serious bacterial leaching in this Templeton soil. A dye infiltration study showed there was great variability in water flow patterns, and most of the flow reaching deeper than 50 cm resulted from macropores, mainly visible cracks. The transient water flow and transport of tracer (Br⁻) and FC were modelled using the HYDRUS-1D software package. The uniform flow van Genuchten model, and the dual-porosity model were used for water flow and the mobile-immobile (MIM) model was used for tracer and FC transport. The hydraulic and solute parameters were optimized during simulation, on the basis of measured datasets from the leaching experiments. There was evidence supporting the presence of macropores, based on the water flow in the post-DSE application stage. The optimised saturated water content (θs) decreased during the post-application process, which could be explained in terms of macropore flow enhanced by irrigation. Moreover, bacterial simulation showed discrepancies in all cases of uniform flow simulations at the very initial stage, indicating that non-equilibrium processes were dominant during those short periods, and suggesting that there were strong dynamic processes involving structure change and subsequently flow paths. It is recommended that management strategies to reduce FC contamination following application of DSE in these soils must aim to decrease preferential flow by adjusting irrigation schemes. Attention needs to be given to a) decreasing irrigation rates at the beginning of each irrigation; b) increasing the number of irrigations, by reducing at the same time the amount of water applied and the irrigation rate at each irrigation; c) applying spray irrigation rather than flood irrigation.

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.

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.

Dynamic soil water repellency is an important soil phenomenon in the vadose zone as it is now recognised that most soils in the world are likely to express some degree of reduced wettability and/or long term hydrophobicity. Fractional wettability and contact angles are, however, rarely discussed or quantified for natural systems. This is particularly the case in the presence of dynamic contact angles. Soil water repellency remains a persistent impediment and challenge to accurate conceptual and numerical models of flow and storage in the vadose zone. This dissertation addresses the opportunity and pressing need for contributions that develop better quantifiable definitions, descriptions, and understanding of soil water repellency. Using materials collected from post wildfire sites, this work employs water and ethanol to identify, isolate, and quantify contact angle dynamics and fractional wettability effects during infiltration. Varied concentrations of water and ethanol solutions were applied to soils and observed through X-ray microtomography, tension infiltration experiments, and moisture content measurements in the laboratory and field. Several analyses from lab and field investigations showed that applications of ethanol and specifically, water-ethanol aqueous solutions provide unique additional insights into proportions of media that remain non-wettable and how those proportions affect overall hydrologic processes, which are not readily observable through water infiltrations alone. Observations include the wetting up of microporous structures, reduced storage, and changes in unsaturated hydraulic conductivities. Challenges which develop as a consequence of variable fluid properties including changes to operational pore assemblages, slow down of wetting fronts, and non-uniqueness relative to infiltration responses are addressed. Important insights and contributions were developed through this approach and water-ethanol mixtures are valuable tools for developing greater quantification and mechanistic data to better inform our models and understanding of dynamic soil water repellency. / Dissertation / Doctor of Philosophy (PhD) / Quantifying fluid behaviours in soils is important for a host of environmental, social, and economic reasons. Over the last 25+ years, one soil phenomenon has garnered increased attention because it interferes with our ability to carry out this work. Soils that are or become water repellent develop all over the world and where hydrophobic or non-wetting substances can enter soil and remain in pore spaces or as coatings on particles. To assist in the tracking and management of its complex effects on water storage and infiltration, the goals of this work were to develop fundamental insights into the manifestation and effects of this variable soil property on key hydrologic properties and processes. This work tests a new conceptual model for understanding these systems through both field and laboratory work and using a number of different technologies. These include X-ray microtomography (μXCT), tension infiltrometry, and more regularly applied techniques which are sensitive to changes in repellency. The works shows how combining fractional wettability and contact angle dynamics generates a stereoscopic conceptual framework which facilitates increased capacity for quantifying and understanding of soil systems expressing dynamic soil water repellency.