Near-surface distributions of soil water and water repellency under three effluent irrigation schemes in a blue gum (Eucalyptus globulus) plantation

Water repellent soils are difficult to irrigate and susceptible to preferential flow, which enhances the potential for accelerated leaching to groundwater of hazardous substances. Over 5 Mha of Australian soil is water repellent, while treated municipal sewage is increasingly used for irrigation. Only if a critical water content is exceeded will repellent soils become wettable. To avoid excessive loss of water from the root zone via preferential flow paths, irrigation schemes should therefore aim to keep the soil wet enough to maintain soil wettability. Our objective was to monitor the near-surface water content and water repellency in a blue gum (Eucalyptus globulus) plantation irrigated with treated sewage. The plantation's sandy soil surface was strongly water repellent when dry. For 4 months, three rows of 15 blue gum trees each received no irrigation, three other rows received 50% of the estimated potential water use minus rainfall, and three more rows received 100%. During this period, 162 soil samples were obtained in three sampling rounds, and their water content (% dry mass) and degree of water repellency determined. Both high and low irrigation effectively wetted up the soil and eliminated water repellency after 2 (high) or 4 (low) months. A single-peaked distribution of water contents was observed in the soil samples, but the water repellency distribution was dichotomous, with 44% extremely water-repellent and 36% wettable. This is consistent with a threshold water content at which a soil sample changes from water repellent to wettable, with spatial variability of this threshold creating a much wider transition zone at the field scale. We characterized this transition zone by expressing the fraction of wettable samples as a function of water content, and demonstrated a way to estimate from this the wettable portion of a field from a number of water content measurements. To keep the plantation soil wettable, the water content must be maintained at a level at which a significant downward flux is likely, with the associated enhanced leaching. At water contents with negligible downward flux, the field is water repellent, and leaching through preferential flow paths is likely. Careful management is needed to resolve these conflicting requirements.

Spatial and temporal distribution of solute leaching in heterogeneous soils: analysis and application to multisampler lysimeter data

Accurate assessment of the fate of salts, nutrients, and pollutants in natural, heterogeneous soils requires a proper quantification of both spatial and temporal solute spreading during solute movement. The number of experiments with multisampler devices that measure solute leaching as a function of space and time is increasing. The breakthrough curve (BTC) can characterize the temporal aspect of solute leaching, and recently the spatial solute distribution curve (SSDC) was introduced to describe the spatial solute distribution. We combined and extended both concepts to develop a tool for the comprehensive analysis of the full spatio-temporal behavior of solute leaching. The sampling locations are ranked in order of descending amount of total leaching (defined as the cumulative leaching from an individual compartment at the end of the experiment), thus collapsing both spatial axes of the sampling plane into one. The leaching process can then be described by a curved surface that is a function of the single spatial coordinate and time. This leaching surface is scaled to integrate to unity, and termed S can efficiently represent data from multisampler solute transport experiments or simulation results from multidimensional solute transport models. The mathematical relationships between the scaled leaching surface S, the BTC, and the SSDC are established. Any desired characteristic of the leaching process can be derived from S. The analysis was applied to a chloride leaching experiment on a lysimeter with 300 drainage compartments of 25 cm2 each. The sandy soil monolith in the lysimeter exhibited fingered flow in the water-repellent top layer. The observed S demonstrated the absence of a sharp separation between fingers and dry areas, owing to diverging flow in the wettable soil below the fingers. Times-to-peak, maximum solute fluxes, and total leaching varied more in high-leaching than in low-leaching compartments. This suggests a stochastic–convective transport process in the high-flow streamtubes, while convection–dispersion is predominant in the low-flow areas. S can be viewed as a bivariate probability density function. Its marginal distributions are the BTC of all sampling locations combined, and the SSDC of cumulative solute leaching at the end of the experiment. The observed S cannot be represented by assuming complete independence between its marginal distributions, indicating that S contains information about the leaching process that cannot be derived from the combination of the BTC and the SSDC.

Quantifying minimum monolith size and solute dilution from multi-compartment percolation sampler data

Preferential flow affects solute transport in natural soils, leading to high spatiotemporal variation of concentration. A multicompartment solute sampler (MCS), yielding multiple breakthrough curves at a given depth, can monitor tracer movement in a heterogeneous soil. We present a technique to estimate from MCS data whether a soil monolith is sufficiently large to capture preferential flow, which is a necessity for tracer breakthrough curves to be representative. For several soils, we estimate that an MCS should be larger than 0.1 to 0.2 m². We also expand dilution theory to analyze the concentration variations of a tracer passing the control plane monitored by the MCS, in addition to the conventional plume spreading analysis. We characterize the set of locally observed breakthrough curves by the entropy-based dilution index. For given first and second-central moment, the spatially uniform log-normal breakthrough curve maximizes the dilution index. The ratio between observed and maximum dilution index is denoted reactor ratio. For a 300-compartment solute sampler, covering an area of 0.75 m², we compute a reactor ratio of 0.665, compared with 0.04 for stochastic-convective and 1 for convective-dispersive transport. With a single, large collector the reactor ratio would be 0.958, severely underestimating concentration variations. Large collector areas are clearly inadequate to estimate dilution. Values of the dilution index and the reactor ratio for individual sampling compartments indicate efficient longitudinal mixing in most but not all cases, and considerable spatial variation of the leaching process.

Investigation of compounds causing water repellency in the rhizosphere of sandy soils from a wide range of locations.

Although soils are generally considered to wet readily, some are actually water repellent at the surface and in the rhizosphere. This phenomenon occurs at low to moderate moisture contents and has been reported from soils under a range of vegetation types and from many regions around the globe. Water repellency in soils can have serious environmental implications including reduced seed germination and plant growth as well as irrigation efficiency, accelerated soil erosion, and enhanced leaching of agrochemicals through preferential flow. it has been proposed that water repellency is caused by the accumulation of hydrophobic organic compounds released as root exudates, microbial byproducts or from decomposing organic matter, which are deposited on mineral or aggregate surfaces, or are present as interstitial matter, Few studies to date have attempted to isolate and characterize these compounds and their structure is therefore only poorly understood, These studies have generally focussed on only a single soil or a small range of samples, have not included non-repellent soils as a control and have not always been able to demonstrate that the substances isolated are indeed responsible for repellency formation.

This study reports on the first part (extraction procedures) of a research programme addressing these gaps in current knowledge by investigating a wide range of severely repellent and wettable ‘control’ samples from different countries, and by including assessments of extraction efficiency and ability of extracts to cause repellency. Analytical methods include DRIFT (Diffuse Reflectance Infrared Fourier Transform Spectroscopy) of soils and IR (Infrared) analysis of extracts.

Key findings are that (i) soil sample heating after extraction is valuable in assessing the effectiveness of the extraction procedure, (ii) Soxhlet extraction using isopropanol/ ammonia (70/30 v/v) was the most effective method in extracting hydrophobic compounds, while leaving the ability of extracted compounds to induce water repellency virtually unaffected, (iii) wettable control soils also contain hydrophobic substances capable of inducing water repellency, (iv) the amount of organic compounds extracted was poorly related to sample repellency, indicating that compounds responsible for repellency may only represent a small fraction of the extract, (v) differences in extraction efficiency between different samples indicate that the compounds responsible may differ generically and/or in terms of their bonding to minerals, and (vi) the combination of repellency assessments with DRIFT on soils and JR on extracts used with internal standards has considerable potential to allow quantification of CH bearing organic matter in the soil, the efficiency of extraction processes for its removal, and its significance in causing water repellency in soils.

Induced polarization, resistivity, and self-potential: a case history of contamination evaluation due to landfill leakage

This article compares the efficiency of induced polarization (IP) and resistivity in characterizing a contamination plume due to landfill leakage in a typical tropical environment. The resistivity survey revealed denser electrical current flow that induced lower resistivity values due to the high ionic content. The increased ionic concentration diminished the distance of the ionic charges close to the membrane, causing a decrease in the IP phenomena. In addition, the self-potential (SP) method was used to characterize the preferential flow direction of the area. The SP method proved to be effective at determining the flow direction; it is also fast and economical. In this study, the resistivity results were better correlated with the presence of contamination (lower resistivity) than the IP (lower chargeability) data.

Dilution of non-reactive tracers in variably saturated sandy structures

Based on a dye tracer experiment in a sand tank we addressed the problem of local dispersion of conservative tracers in the unsaturated zone. The sand bedding was designed to have a defined spatial heterogeneity with a strong anisotropy. We estimated the parameters that characterize the local dispersion and dilution from concentration maps of a high spatial and temporal resolution obtained by image analysis. The plume spreading and mixing behavior was quantified on the basis of the coefficient of variation of the concentration and of the dilution index. The heterogeneous structure modified the flow pattern depending on water saturation. The shape of the tracer plumes revealed the structural signature of the sand bedding at low saturation only. In this case pronounced preferential flow was observed. At higher flow rates the structure remained hidden by a spatially almost homogeneous behavior of the plumes. In this context, we mainly discuss the mechanism of re-distributing a finite mass of inert solutes over a large volume, due to macro- and micro-heterogeneities of the structure. (C) 2001 Elsevier Science Ltd. AU rights reserved.

Transport of iodide in structured clay–loam soil under maize during irrigation experiments analyzed using HYDRUS model

Transport of radioactive iodide 131I− in a structured clay loam soil under maize in a final growing phase was monitored during five consecutive irrigation experiments under ponding. Each time, 27 mm of water were applied. The water of the second experiment was spiked with 200 MBq of 131I− tracer. Its activity was monitored as functions of depth and time with Geiger-Müller (G-M) detectors in 11 vertically installed access tubes. The aim of the study was to widen our current knowledge of water and solute transport in unsaturated soil under different agriculturally cultivated settings. It was supposed that the change in 131I− activity (or counting rate) is proportional to the change in soil water content. Rapid increase followed by a gradual decrease in 131I− activity occurred at all depths and was attributed to preferential flow. The iodide transport through structured soil profile was simulated by the HYDRUS 1D model. The model predicted relatively deep percolation of iodide within a short time, in a good agreement with the observed vertical iodide distribution in soil. We found that the top 30 cm of the soil profile is the most vulnerable layer in terms of water and solute movement, which is the same depth where the root structure of maize can extend.

Tracking water pathways in steep hillslopes by d18O depth profiles of soil water

Assessing temporal variations in soil water flow is important, especially at the hillslope scale, to identify mechanisms of runoff and flood generation and pathways for nutrients and pollutants in soils. While surface processes are well considered and parameterized, the assessment of subsurface processes at the hillslope scale is still challenging since measurement of hydrological pathways is connected to high efforts in time, money and personnel work. The latter might not even be possible in alpine environments with harsh winter processes. Soil water stable isotope profiles may offer a time-integrating fingerprint of subsurface water pathways. In this study, we investigated the suitability of soil water stable isotope (d18O) depth profiles to identify water flow paths along two transects of steep subalpine hillslopes in the Swiss Alps. We applied a one-dimensional advection–dispersion model using d18O values of precipitation (ranging from _24.7 to _2.9‰) as input data to simulate the d18O profiles of soil water. The variability of d18O values with depth within each soil profile and a comparison of the simulated and measured d18O profiles were used to infer information about subsurface hydrological pathways. The temporal pattern of d18O in precipitation was found in several profiles, ranging from _14.5 to _4.0‰. This suggests that vertical percolation plays an important role even at slope angles of up to 46_. Lateral subsurface flow and/or mixing of soil water at lower slope angles might occur in deeper soil layers and at sites near a small stream. The difference between several observed and simulated d18O profiles revealed spatially highly variable infiltration patterns during the snowmelt periods: The d18O value of snow (_17.7 ± 1.9‰) was absent in several measured d18O profiles but present in the respective simulated d18O profiles. This indicated overland flow and/or preferential flow through the soil profile during the melt period. The applied methods proved to be a fast and promising tool to obtain time-integrated information on soil water flow paths at the hillslope scale in steep subalpine slopes.

Tab. 1 Mole % CO2 and 13C values of CO2 at baseline and at the final monitoring event for the eight observation wells

During CO2 storage operations in mature oilfields or saline aquifers it is desirable to trace the movement of injected CO2 for verification and safety purposes. We demonstrate the successful use of carbon isotope abundance ratios for tracing the movement of CO2 injected at the Cardium CO2 Storage Monitoring project in Alberta between 2005 and 2007. Injected CO2 had a d13C value of -4.6±1.1 per mil that was more than 10 per mil higher than the carbon isotope ratios of casing gas CO2 prior to CO2 injection with average d13C values ranging from -15.9 to -23.5 per mil. After commencement of CO2 injection, d13C values of casing gas CO2 increased in all observation wells towards those of the injected CO2 consistent with a two-source end-member mixing model. At four wells located in a NE-SW trend with respect to the injection wells, breakthrough of injected CO2 was registered chemically (>50 mol % CO2) and isotopically 1-6 months after commencement of CO2 injection resulting in cumulative CO2 fluxes exceeding 100000 m**3 during the observation period. At four other wells, casing gas CO2 contents remained below 5 mol % resulting in low cumulative CO2 fluxes (<2000 m**3) throughout the entire observation period, but carbon isotope ratios indicated contributions between <30 and 80% of injected CO2. Therefore, we conclude that monitoring the movement of CO2 in the injection reservoir with geochemical and isotopic techniques is an effective approach to determine plume expansion and to identify potential preferential flow paths provided that the isotopic composition of injected CO2 is constant and distinct from that of baseline CO2.

Application of active heat pulse method with fiber optic temperature sensing for estimation of wetting bulbs and water distribution in drip emitters.

Through the use of the Distributed Fiber Optic Temperature Measurement (DFOT) method, it is possible to measure the temperature in small intervals (on the order of centimeters) for long distances (on the order of kilometers) with a high temporal frequency and great accuracy. The heat pulse method consists of applying a known amount of heat to the soil and monitoring the temperature evolution, which is primarily dependent on the soil moisture content. The use of both methods, which is called the active heat pulse method with fiber optic temperature sensing (AHFO), allows accurate soil moisture content measurements. In order to experimentally study the wetting patterns, i.e. shape, size, and the water distribution, from a drip irrigation emitter, a soil column of 0.5 m of diameter and 0.6 m high was built. Inside the column, a fiber optic cable with a stainless steel sheath was placed forming three concentric helixes of diameters 0.2 m, 0.4 m and 0.6 m, leading to a 148 measurement point network. Before, during, and after the irrigation event, heat pulses were performed supplying electrical power of 20 W/m to the steel. The soil moisture content was measured with a capacitive sensor in one location at depths of 0.1 m, 0.2 m, 0.3 m and 0.4 m during the irrigation. It was also determined by the gravimetric method in several locations and depths before and right after the irrigation. The emitter bulb dimensions and shape evolution was satisfactorily measured during infiltration. Furthermore, some bulb's characteristics difficult to predict (e.g. preferential flow) were detected. The results point out that the AHFO is a useful tool to estimate the wetting pattern of drip irrigation emitters in soil columns and show a high potential for its use in the field.

Indagine modellistica sulla tenuta idraulica delle difese arginali deputate alla protezione idraulica del territorio

La prima parte di questo lavoro di tesi tratta dell’interazione tra un bacino di laminazione e il sottostante acquifero: è in fase di progetto, infatti, la costruzione di una cassa di espansione sul torrente Baganza, a monte della città di Parma. L’obiettivo di tale intervento è di ridurre il rischio di esondazione immagazzinando temporaneamente, in un serbatoio artificiale, la parte più pericolosa del volume di piena che verrebbe rilasciata successivamente con portate che possono essere agevolmente contenute nel tratto cittadino del torrente. L’acquifero è stato preliminarmente indagato e monitorato permettendone la caratterizzazione litostratigrafica. La stratigrafia si può riassumere in una sequenza di strati ghiaioso-sabbiosi con successione di lenti d’argilla più o meno spesse e continue, distinguendo due acquiferi differenti (uno freatico ed uno confinato). Nel presente studio si fa riferimento al solo acquifero superficiale che è stato modellato numericamente, alle differenze finite, per mezzo del software MODFLOW_2005. L'obiettivo del presente lavoro è di rappresentare il sistema acquifero nelle condizioni attuali (in assenza di alcuna opera) e di progetto. La calibrazione è stata condotta in condizioni stazionarie utilizzando i livelli piezometrici raccolti nei punti d’osservazione durante la primavera del 2013. I valori di conducibilità idraulica sono stati stimati per mezzo di un approccio geostatistico Bayesiano. Il codice utilizzato per la stima è il bgaPEST, un software gratuito per la soluzione di problemi inversi fortemente parametrizzati, sviluppato sulla base dei protocolli del software PEST. La metodologia inversa stima il campo di conducibilità idraulica combinando osservazioni sullo stato del sistema (livelli piezometrici nel caso in esame) e informazioni a-priori sulla struttura dei parametri incogniti. La procedura inversa richiede il calcolo della sensitività di ciascuna osservazione a ciascuno dei parametri stimati; questa è stata valutata in maniera efficiente facendo ricorso ad una formulazione agli stati aggiunti del codice in avanti MODFLOW_2005_Adjoint. I risultati della metodologia sono coerenti con la natura alluvionale dell'acquifero indagato e con le informazioni raccolte nei punti di osservazione. Il modello calibrato può quindi essere utilizzato come supporto alla progettazione e gestione dell’opera di laminazione. La seconda parte di questa tesi tratta l'analisi delle sollecitazioni indotte dai percorsi di flusso preferenziali causati da fenomeni di piping all’interno dei rilevati arginali. Tali percorsi preferenziali possono essere dovuti alla presenza di gallerie scavate da animali selvatici. Questo studio è stato ispirato dal crollo del rilevato arginale del Fiume Secchia (Modena), che si è verificato in gennaio 2014 a seguito di un evento alluvionale, durante il quale il livello dell'acqua non ha mai raggiunto la sommità arginale. La commissione scientifica, la cui relazione finale fornisce i dati utilizzati per questo studio, ha attribuito, con molta probabilità, il crollo del rilevato alla presenza di tane di animali. Con lo scopo di analizzare il comportamento del rilevato in condizioni integre e in condizioni modificate dall'esistenza di un tunnel che attraversa il manufatto arginale, è stato realizzato un modello numerico 3D dell’argine mediante i noti software Femwater e Feflow. I modelli descrivono le infiltrazioni all'interno del rilevato considerando il terreno in entrambe le porzioni sature ed insature, adottando la tecnica agli elementi finiti. La tana è stata rappresentata da elementi con elevata permeabilità e porosità, i cui valori sono stati modificati al fine di valutare le diverse influenze sui flussi e sui contenuti idrici. Per valutare se le situazioni analizzate presentino o meno il verificarsi del fenomeno di erosione, sono stati calcolati i valori del fattore di sicurezza. Questo è stato valutato in differenti modi, tra cui quello recentemente proposto da Richards e Reddy (2014), che si riferisce al criterio di energia cinetica critica. In ultima analisi è stato utilizzato il modello di Bonelli (2007) per calcolare il tempo di erosione ed il tempo rimanente al collasso del rilevato.

A method for the in-situ determination of the hydraulic conductivity of gravels as used in constructed wetlands for wastewater treatment

A new instrument and method are described that allow the hydraulic conductivities of highly permeable porous materials, such as gravels in constructed wetlands, to be determined in the field. The instrument consists of a Mariotte siphon and a submersible permeameter cell with manometer take-off tubes, to recreate in-situ the constant head permeameter test typically used with excavated samples. It allows permeability to be measured at different depths and positions over the wetland. Repeatability obtained at fixed positions was good (normalised standard deviation of 1–4%), and results obtained for highly homogenous silica sand compared well when the sand was retested in a lab permeameter (0.32 mm.s–1 and 0.31 mm.s–1 respectively). Practical results have a ±30% associated degree of uncertainty because of the mixed effect of natural variation in gravel core profiles, and interstitial clogging disruption during insertion of the tube into the gravel. This error is small, however, compared to the orders of magnitude spatial variations detected. The technique was used to survey the hydraulic conductivity profile of two constructed wetlands in the UK, aged 1 and 15 years respectively. Measured values were high (up to 900 mm.s –1) and varied by three orders of magnitude, reflecting the immaturity of the wetland. Detailed profiling of the younger system suggested the existence of preferential flow paths at a depth of 200 mm, corresponding to the transition between more coarse and less coarse gravel layers (6–12 mm and 3–6 mm respectively), and transverse drift towards the outlet.

Rainfall infiltration and soil hydrological characteristics below ancient forest, planted forest and grassland in a temperate northern climate

How rainfall infiltration rate and soil hydrological characteristics develop over time under forests of different ages in temperate regions is poorly understood. In this study, infiltration rate and soil hydrological characteristics were investigated under forests of different ages and under grassland. Soil hydraulic characteristics were measured at different scales under a 250-year-old grazed grassland (GL), 6-year-old (6yr) and 48-year-old (48yr) Scots pine (Pinus sylvestris) plantations, remnant 300-year-old individual Scots pine (OT) and a 4000-year-old Caledonian Forest (AF). In situ field-saturated hydraulic conductivity (Kfs) was measured, and visible root:soil area was estimated from soil pits. Macroporosity, pore structure and macropore connectivity were estimated from X-ray tomography of soil cores, and from water-release characteristics. At all scales, the median values for Kfs, root fraction, macroporosity and connectivity values tended to AF>OT>48yr>GL>6yr, indicating that infiltration rates and water storage increased with forest age. The remnant Caledonian Forest had a huge range of Kfs (12 to >4922mmh-1), with maximum Kfs values 7 to 15 times larger than those of 48-year-old Scots pine plantation, suggesting that undisturbed old forests, with high rainfall and minimal evapotranspiration in winter, may act as important areas for water storage and sinks for storm rainfall to infiltrate and transport to deeper soil layers via preferential flow. The importance of the development of soil hydrological characteristics under different aged forests is discussed.

Temporal variation of faecal indicator bacteria in tropical urban storm drains

Human faecal contamination poses a widespread hazard for human health. In urban areas, sewer leakage may be an important cause of faecal pollution to surface water. Faecal indicator bacteria (FIB) are the most widely used indicators to monitor surface water quality. However, assessing whether a water body is meeting water quality criteria is made difficult by the high variability of FIB concentrations over time. In this study, the variation of FIB concentration in surface water from tropical urban catchments is investigated. Eleven urban sub-catchments were sampled hourly over 24-hr and samples analysed for FIB. It was found that FIB show a diurnal pattern that is characterised by daytime FIB concentrations that are significantly higher than nighttime FIB concentrations. This observed diurnal variation of FIB closely follows that of sewer flows and contrasts with observations in rural streams where FIB concentrations are known to be low in the daytime and high during the night. Field tracer tests provide qualitative evidence of sewage exfiltration and transport to drains via preferential flow paths. The diurnal FIB variation and field tracer tests indicate the likelihood of surface water contamination due to leaking sewers. The results further suggest that contamination of surface-water drains is likely a widespread problem in tropical urban areas due to extensive drainage networks and the persistence of FIB under tropical conditions. Because of FIB variation over time, the time at which samples are collected is important in being able to capture the daily maximum and minimum FIB concentrations. The Kruskal-Wallis test shows that hourly sampling from 04:00 to 07:00 and from 12:00 to 15:00 results in significantly different FIB concentration (minimum and maximum, respectively). Furthermore, the Wilcoxon-Mann-Whitney test shows that sampling at 12:00 and 14:00 results in significantly higher FIB concentrations, while sampling at 05:00 and 04:00 or 05:00 and 06:00 results in significantly lower FIB concentrations, than sampling at other hours of the day.

Singularity and multifractal characterization of signals with wavelets. Application to clay soil images

In the field of multiscale analysis of signals, including images, the wavelet transform is one of the most attractive and powerful tool due to its ability to focus on signals structures at different scales. Wavelet Transform at different scales is successfully applied to image characterization (which can be applied to a watermarking scheme) and multiscale singularity detection and processing. In this work we show further research of computation of multifractals properties such as the multifractal spectrum (D(alpha)) applied to dye stained images of natural terrain. This can be useful for statically describing preferential flow path geometry.