Soil and groundwater pollution has become a global issue since the advent of industrialization and mechanized agriculture. Some contaminants such as PAHs may persist in the subsurface for decades and centuries. In a bid to address these issues, protection of groundwater must be based on the quantification of potential threats to pollution at the subsurface which is often inaccessible. Risk assessment of groundwater pollution may however be strongly supported by applying process-based simulation models, which turn out to be particularly helpful with regard to long-term predictions, which cannot be undertaken by experiments. Such reliable predictions, however, can only be achieved if the used modeling tool is known to be applicable. The aim of this work was threefold. First, a source strength function was developed to describe the leaching behavior of point source organic contaminants and thereby acting as a time-dependent upper boundary condition for transport models. For general application of these functions dimensionless numbers known as Damköhler numbers were used to characterize the reaction of the pollutants with the solid matrix. Two functions were derived and have been incorporated into an Excel worksheet to act as a practical aid in the quantification of leaching behavior of organic contaminant in seepage water prognoses. Second, the process based model tool SMART, which is well validated for laboratory scale data, was applied to lysimeter scale data from two research centres, FZJ (Jülich) and GSF (München) for long term predictions. Results from pure forward model runs show a fairly good correlation with the measured data. Finally, the derived source term functions in combination with the SMART model were used to assess groundwater vulnerability beneath a typical landfill at Kwabenya in Ghana. The predicted breakthrough time after leaking from the landfill was more than 200 years considering the operational time of the facility (30 years). Considering contaminant degradation, the landfill would therefore not cause groundwater pollution under the simulated scenarios and the SMART model can be used to establish waste acceptance criteria for organic contaminants in the landfill at Kwabenya / Seit dem Beginn der Industrialisierung und der mechanisierten Landwirtschaft wurde die Boden- und Grundwasserverschmutzung zu einem weltweiten Problem. Einige Schadstoffe wie z. B. PAK können für Jahrzehnte oder Jahrhunderte im Untergrund bestehen. Um diese Probleme behandeln zu können, muss der Schutz des Grundwassers basierend auf der Quantifizierung potentieller Gefährdungen des zumeist unzugänglichen Untergrundes erfolgen. Risikoabschätzungen von Grundwasserverschmutzungen können jedoch durch die Anwendung prozess-basierter Simulationsmodelle erheblich unterstützt werden, die sich besonders im Hinblick auf Langzeitvorhersagen als hilfreich erweisen und nicht experimentell ermittelbar sind. Derart zuverlässige Vorhersagen können jedoch nur erhalten werden, wenn das verwendete Modellierwerkzeug als anwendbar bekannt ist. Das Ziel dieser Arbeit bestand aus drei Teilen. Erstens wurde eine Quellstärke-funktion entwickelt, die das Ausbreitungsverhalten organischer Schadstoffe aus einer Punktquelle beschreibt und dadurch als zeitabhängige obere Randbedingung bei Transportmodellen dienen kann. Im Hinblick auf die allgemeine Anwendbarkeit dieser Funktion werden als Damköhler-Zahlen bekannte, dimensionslose Zahlen verwendet, um die Reaktion von Schadstoffen mit Feststoffen zu charakterisieren. Zwei Funktionen wurden abgeleitet und in ein Excel-Arbeitsblatt eingefügt, das ein praktisches Hilfsmittel bei der Quantifizierung des Freisetzungsverhaltens organischer Schadstoffe im Rahmen der Sickerwasserprognose darstellt. Der zweite Teil dieser Arbeit beinhaltet die Anwendung des prozessbasierten und mittels Laborexperimenten validierten Modellwerkzeugs SMART für Langzeitprognosen auf der Lysimeterskala anhand von Daten zweier Forschungszentren, FZJ (Jülich) und GSF (München). Ergebnisse reiner Vorwärtsmodellierungsläufe zeigten gute Übereinstimmungen mit den gemessenen Daten. Im dritten Teil wurden die erhaltenen Quellstärkefunktionen in Kombination mit dem SMART-Modell eingesetzt, um das Grundwassergefährdungspotential unter einer typischen Deponie in Kwabenya, Ghana, einzuschätzen. Die vorhergesagten Durchbruchszeiten nach einer Leckage in der Deponie betragen über 200 Jahre bei einer Betriebszeit von 30 Jahren. Unter Berücksichtigung des Schadstoffabbaus verursacht die Deponie somit keine Grundwasserverunreinigung im Rahmen der simulierten Szenarien und das SMART-Modell kann verwendet werden, um Schadstoffgrenzwerte für organische Schadstoffe in der Deponie in Kwabenya festzulegen.

Les objectifs de la thèse sont de (i) modéliser le transport de l’eau, d’un traceur et des polluants dans les sols contaminés, (ii) étudier l’estimabilité des paramètres et les corrélations entre les paramètres, (iii) optimiser les paramètres. Les différents modèles implantés dans le logiciel HYDRUS qui permettent de rendre compte ou pas de l’écoulement préférentiel et du transport hors équilibre physique et chimique sont choisis. Les études concernant le transport d’eau dans un lysimètre de terrain ont montré que les données quotidiennes de pressions et de teneurs en eau volumique contiennent plus d’information que les données horaires, que les pressions ont plus d’information que les teneurs en eau volumique, et que les corrélations des paramètres ont fait perturber les résultats de l’optimisation. Sur le même lysimètre, l’étude d’estimabilité des paramètres caractéristiques pour le transport du traceur (bromure) a montré que les concentrations dans les solutions de percolation ne sont pas suffisantes pour estimer le paramètre de transfert de l’eau entre les zones mobile et immobile car ce paramètre est fortement corrélé avec le paramètre de transfert de soluté. Pour le transport des hydrocarbures aromatiques polycycliques (HAP) dans les colonnes de laboratoire sous différentes conditions de saturation en eau, quand le degré de transport hors équilibre chimique des HAP est élevé, les concentrations en HAP dans les solutions de percolation de la colonne non saturée contiennent plus d’information que celles dans la colonne saturée / Preferential flow and nonequilibrium transport are probably the most frustrating in terms of hampering accurate predictions of contaminant transport through the vadose zone. The mathematical description of preferential flow and nonequilibrium transport needs many parameters that are not measurable. Therefore, the inverse method is a promising way to estimate model parameters. The main objectives of this work are to (i) study the water flow using the uniform flow and dual-porosity models, tracer and contaminant transport using the uniform transport model and/or physical and chemical nonequilibrium transport models, (ii) investigate parameter estimability and correlations between different parameters, and (iii) optimize the hydraulic properties and solute transport parameters. The results concerning the water flow in the bare field lysimeter show that daily data contained much more information than hourly data, daily pressure heads contained more information than daily water contents; the correlations between different parameters hamper the optimization results strongly. Basing on the tracer concentrations in the leaching solution of the lysimeter, the first-order rate water transfer coefficient was not estimable since this parameter was highly correlated with the solute transfer coefficient. PAH concentrations in the leaching solution of the contaminated soil column under saturated and nonsaturated flow conditions show that when the degree of chemical nonequilibirum transport is high, the solute leaching of the nonsaturated column contained more information than those of the saturated column. In addition, the fraction of sites with instantaneous sorption and the linear adsorption distribution coefficient always showed a very strong correlation, they were impossible to optimize simultaneously

From the Proceedings of the 1971 Meetings of the Arizona Section - American Water Resources Assn. and the Hydrology Section - Arizona Academy of Science - April 22-23, 1971, Tempe, Arizona / A study was made to determine how the evaporation rate from a bare Adelanto loam soil in Phoenix changes with season and with time since the last irrigation. The evaporation rates were determined by precision lysimeters in a bare field, with measurements being taken in every month of the year for at least a week after irrigation. The data exhibited a cosine-shaped curve, with a maximum evaporation rate of about 5 mm/day in summer and a minimum rate of about 2 mm/day in winter. By the seventh day, seasonal effects virtually disappear, and the evaporation rate is the same in both summer and winter, being about 2 mm/day after the 7th day and about 0.75 mm/day after the 21st day. It is generally accepted that soil dries in 3 stages, and the transition between the 1st and 2nd stages occurs when atmospheric conditions are no longer critical. In previous laboratory studies of soil drying, with constant atmospheric conditions, stage 1 was easily distinguished from stage II, and these results correlated closely with the equations of Gardner and Hillel. The individual drying curves of this field study were qualitatively different from the laboratory studies and did not confirm the predictions of the equations, suggesting that diurnal variations in temperature and other meteorological parameters have caused the difference.

No description available.

L’hétérogénéité de la zone non saturée joue un rôle important dans le transfert d’eau et de soluté car elle accentue à la fois le développement des zones de stockage temporelles et les écoulements préférentiels. Par conséquent, la validation des modèles prédictifs nécessite le développement des outils expérimentaux spécifiques afin d’observer et de quantifier les mécanismes de transport impliqués dans un système non saturé hétérogène. Cette thèse vise à étudier l’effet combiné de la vitesse d’infiltration, de la barrière capillaire et l’angle de la pente d’interface entre deux matériaux sur les processus de l’écoulement de l’eau et du transport de soluté dans un modèle physique, le lysimètre de laboratoire 1x1x1.6 m3, nommé LUGH (Lysimeter for Urban Groundwater Hydrology) et un modèle numérique 3D de ce lysimètre. Le lysimètre LUGH est rempli par un sable fin et un mélange bimodal (50 % sable fin et 50 % gravier) en deux configurations: un profil uniforme de matériau bimodal ou un profil avec deux couches avec une pente de 14o. Ces agencements figurent l’hétérogénéité structurale et texturale observée sur un des sites expérimentaux de l’OTHU (Observatoire de Terrain en Hydrologie Urbaine) : le bassin d’infiltration d’eaux pluviales Django Reinhardt géré par la ville de Lyon (France). Le lysimètre est alimenté en eau et avec un traceur inerte (bromure de potassium, KBr) sur une partie de la surface par un système d’arrosage automatique. Les effluents ont été recueillis dans quinze sorties différentes en bas du lysimètre. La forte hétérogénéité des flux des sorties et des courbes de percée souligne la mise en place des écoulements préférentiels résultant à la fois de l’effet de barrière capillaire et de l’effet de fond du lysimètre. A partir des résultats expérimentaux, la modélisation numérique à l’aide de logiciel COMSOL MultiphysicsTM a permis de mieux comprendre les mécanismes responsables de ces transferts hétérogènes. Lorsque le modèle numérique validé, un test de sensibilité a été conduit pour étudier les effets de la vitesse d’infiltration et de la pente de l’interface sur l’écoulement. Les résultats montrent que la diminution de la vitesse d’infiltration ou l’augmentation de la pente de l’interface favorisent le développement des écoulements préférentiels. Notre étude a donné également des renseignements pertinents sur le couplage entre les processus hydrodynamiques et le transfert des solutés dans les sols non saturés hétérogènes en soulignant le rôle de la géométrie des interfaces ainsi que celui des conditions aux limites comme des facteurs clés pour la quantification des écoulements préférentiels. / The heterogeneity of the unsaturated zone plays an important role in the water and solutes transfer as it accentuates both the development of stagnant zones for water and preferential flow. Therefore, the validation of predictive models requires the development of specific experimental tools to observe and quantify the transport mechanisms involved in a heterogeneous unsaturated system. The aim of this thesis is to describe the combined effect of infiltration, capillary barrier and sloping layered soil on both flow and solute transport processes in a large physical model (1x1x1.6 m3) called LUGH (Lysimeter for Urban Groundwater Hydrology) and a 3D numerical flow model. Sand and a soil composed of a bimodal sand-gravel mixture were placed in the lysimeter represent one of the commun structural and textural elements of the heterogeneity observed in the vadose zone under an infiltration basin of Lyon (France). The soil was compacted in two configurations: a uniform profile and a profile with two layers having a slope of 14°. Water and an inert tracer (KBr) were injected from the top of the lysimeter using a specific water sprinkler system and collected at 15 different outlets at the bottom. The 15 breakthrough curves obtained presented high heterogeneity, emphasising the establishment of a preferential flow resulting from both capillary barrier and soil layer dip effects. Numerical modelling led to better understanding of the mechanisms responsible for these heterogeneous transfers and it was also used to perform a sensitivity analysis of the effects of water velocity (water and solute flux fed by the sprinkler) and the slope interface. The results show that decreasing velocity and increasing the slope of the interface can lead to the development of preferential flows. In addition, the offset of the centre of gravity of the flow distribution at the output increases linearly as a function of the slope angle of the layered soil. This paper allows coupling the hydrodynamic approach with the transfer of pollutants in unsaturated heterogeneous soil and highlighting preferential flow by flow modeling.

De nombreuses friches héritées d’activités industrielles révolues présentent des pollutions organiques persistantes (charbon, goudron…). Alors que la réglementation impose un une évaluation du degré de contamination de ces sites, elle ne prend pas en compte les sous-produits de transformation tels que les composés polaires, peu considérés. Pourtant ils se solubilisent dans la phase aqueuse par percolation des eaux météoriques à travers ces surfaces contaminées. Bien que la littérature ciblant la matière organique dissoute (MOD) d’origine naturelle soit abondante, elle n’est pas directement transposable à la MOD issues des friches contaminées, qui reste à définir pour une meilleure compréhension du devenir de la pollution dans ces sites anthropisés. Une approche multi-techniques a été appliquée pour appréhender le plus précisément possible la MOD issue de terres d’anciennes cokeries et usines à gaz, par le biais d’expériences en laboratoire et des dispositifs de terrain (lysimètres). L’étude de ces derniers montre qu’ils contiennent une forte teneur en MOD aromatique, dont les composés aromatiques polycycliques ne représentent qu’une faible proportion de la MOD totale. Des expériences complémentaires ciblant l’influence de certains paramètres (pH et hydrophobicité) suggèrent un lien fort entre le pH et l’organisation spatiale de la MOD, ainsi qu’une diminution de la masse moléculaire apparente avec l’augmentation de l’hydrophobicité. Des expériences de vieillissement artificiel ont montré un enrichissement en produits polaires condensés entraînant leur forte mobilisation par l’eau / A large amount of wastelands inherited from former industrial activities contains persistent organic contamination (coal, coal tar…). While the regulation requires an evaluation of the contamination degree of these soils, it doesn’t take into account the transformation by-products such as polar compounds, poorly studied. Yet they solubilize in aqueous phase by percolation of meteoric waters through these contaminated sites. Despite the fact that literature targeting the fresh DOM is abundant, it is not directly transposable to the anthropogenic DOM coming from wastelands, which still need to be more precisely defined to improve our knowledge of this specific DOM and its evolution over time. A multi-technical approach was developed to comprehend the anthropogenic DOM coming from former coking and gas plant soils, thanks to a combination of laboratory experiments (under controlled conditions) and on field devices (lysimeters). Their study show that they contained high polycyclic aromatic compound (PAC) contents, whose proportions in polar PACs exceed the ones in soil measurements. However, these PACs only consist of a low proportion of the total DOM. Complementary experiences targeting the influence of some parameters (pH, hydrophobicity) show a major hydrophobic characteristic and suggest a strong link between the pH and the spatial DOM organization. Artificial aging experiences show an enrichment in polar compounds leading to their water mobilization

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.

The disposal of waste from agricultural activities has been recognised as a source of environmental contamination by endocrine disrupting chemicals (EDCs). The New Zealand dairy industry produces a large volume of dairy farm effluent, which contains EDCs in the form of estrogens. Most of this dairy farm effluent is applied onto the land for disposal. Groundwater and soil contamination by estrogens following waste application on the land have been reported overseas, but our understanding of the processes and factors governing the fate of estrogens in the soil is poor. Therefore the main goal of the present study was to better understand the fate and transport of estrogens, in particular 17β-estradiol (E2) and estrone (E1) in soil. In order to quantify E1 and E2 in drainage water and soil samples, chemical analysis by gas-chromatography mass-spectrometry (GC-MS) was carried out. This included sample extraction, sample clean-up through silica gel and gel permeation chromatography, and sample extract derivatisation prior to analysis. In order to develop a reliable method to extract estrogens from soil, research was conducted to optimise E1 and E2 extraction conditions by adjusting the number of sonication and shaking events, as well as the volume and type of solvent. Among five solvents and solvent mixtures tested, the best recovery on spiked and aged soil was obtained using an isopropanol/water (1:1) mix. A microcosm experiment was carried out to determine the dissipation rates of E2 and E1, at 8°C and at field capacity, in the Templeton soil sampled at two different depths (5-10 cm and 30-35 cm). The dissipation rates decreased with time and half-life values of 0.6-0.8 d for E1 and 0.3-0.4 d for E2 were found for the two depths studied. A field transport experiment was also carried out in winter, over three months, by applying dairy farm effluent spiked with estrogens onto undisturbed Templeton soil lysimeters (50 cm in diameter and 70 cm deep). The hormones were applied in dairy farm effluent at 120 mg m⁻² for E2 and 137 mg m⁻² for E1. The results of the transport experiment showed that in the presence of preferential/macropore flow pathways 0.3-0.7% of E2 and 8-13% of E1 was recovered in the leachate at the bottom of the lysimeters after 3 months, and 1-7% of the recovered E2 and 3-54% of the recovered E1 was leached within 2 days of application. These results suggest that leaching of estrogens via preferential/macropore flow pathways is the greatest concern for groundwater contamination. In the absence of preferential/macropore flow pathways, a significant amount (> 99.94%) of both hormones dissipated in the top 70 cm of soil, due to sorption and rapid biodegradation. Surprisingly, in all cases, estrogen breakthrough occurred before that of an inert tracer (bromide). This could not be explained by the advection-dispersion transport of estrogens, nor by their presence as antecedent concentrations in the soil. It was therefore suggested that colloidal enhanced transport of estrogens was responsible for the earlier breakthrough of estrogens and caused the leaching of a fraction of the applied estrogens to a soil depth of 70 cm. A two-phase model, adapted from a state-space mixing cell model, was built to describe the observed estrogen transport processes under transient flow. The model takes into account 3 transport processes namely, advection-dispersion, preferential/macropore flow and colloidal enhanced transport. This model was able to successfully describe the estrogen transport observed from the lysimeters.

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.