Effects of Soil Compaction and Organic Carbon Content on Preferential Flow in Loamy Field Soils

Abstract: Preferential flow and transport through macropores affect plant water use efficiency and enhance leaching of agrochemicals and the transport of colloids, thereby increasing the risk for contamination of groundwater resources. The effects of soil compaction, expressed in terms of bulk density (BD), and organic carbon (OC) content on preferential flow and transport were investigated using 150 undisturbed soil cores sampled from 15 × 15–m grids on two field sites. Both fields had loamy textures, but one site had significantly higher OC content. Leaching experiments were conducted in each core by applying a constant irrigation rate of 10 mm h−1 with a pulse application of tritium tracer. Five percent tritium mass arrival times and apparent dispersivities were derived from each of the tracer breakthrough curves and correlated with texture, OC content, and BD to assess the spatial distribution of preferential flow and transport across the investigated fields. Soils from both fields showed strong positive correlations between BD and preferential flow. Interestingly, the relationships between BD and tracer transport characteristics were markedly different for the two fields, although the relationship between BD and macroporosity was nearly identical. The difference was likely caused by the higher contents of fines and OC at one of the fields leading to stronger aggregation, smaller matrix permeability, and a more pronounced pipe-like pore system with well-aligned macropores.

A new modelling approach to simulate preferential flow and transport in water repellent media: parameter sensitivity, and effects on crop growth and solute leaching

A modified version of the popular agrohydrological model SWAP has been used to evaluate modelling of soil water flow and crop growth at field situations in which water repellency causes preferential flow. The parameter sensitivity in such situations has been studied. Three options to model soil water flow within SWAP are described and compared: uniform flow, the classical mobile-immobile concept, and a recent concept accounting for the dynamics of finger development resulting from unstable infiltration. Data collected from a severely water-repellent affected soil located in Australia were used to compare and evaluate the usefulness of the modelling options for the agricultural management of such soils.

The study shows that an assumption of uniform flow in a water-repellent soil profile leads to an underestimation of groundwater recharge and an overestimation of plant transpiration and crop production. The new concept of modelling taking finger dynamics into account provides greater flexibility and can more accurately model the observed effects of preferential flow compared with the classical mobile–immobile concept. The parameter analysis indicates that the most important factor defining the presence and extremity of preferential flow is the critical soil water content.

Comparison of the modelling results with the Australian field data showed that without the use of a preferential flow module, the effects of the clay amendments to the soil were insufficiently reproduced in the dry matter production results. This means that the physical characteristics of the soil alone are not sufficient to explain the measured increase in production on clay amended soils. However, modelling with the module accounting for finger dynamics indicated that the preferential flow in water repellent soils that had not been treated with clay caused water stress for the crops, which would explain the decrease in production.

Concepts and dimensionality in modeling unsaturated water flow and solute transport

Many environmental studies require accurate simulation of water and solute fluxes in the unsaturated zone. This paper evaluates one- and multi-dimensional approaches for soil water flow as well as different spreading mechanisms to model solute behavior at different scales. For quantification of soil water fluxes,Richards equation has become the standard. Although current numerical codes show perfect water balances, the calculated soil water fluxes in case of head boundary conditions may depend largely on the method used for spatial averaging of the hydraulic conductivity. Atmospheric boundary conditions, especially in the case of phreatic groundwater levels fluctuating above and below a soil surface, require sophisticated solutions to ensure convergence. Concepts for flow in soils with macro pores and unstable wetting fronts are still in development. One-dimensional flow models are formulated to work with lumped parameters in order to account for the soil heterogeneity and preferential flow. They can be used at temporal and spatial scales that are of interest to water managers and policymakers. Multi-dimensional flow models are hampered by data and computation requirements.Their main strength is detailed analysis of typical multi-dimensional flow problems, including soil heterogeneity and preferential flow. Three physically based solute-transport concepts have been proposed to describe solute spreading during unsaturated flow: The stochastic-convective model (SCM), the convection-dispersion equation (CDE), and the fraction aladvection-dispersion equation (FADE). A less physical concept is the continuous-time random-walk process (CTRW). Of these, the SCM and the CDE are well established, and their strengths and weaknesses are identified. The FADE and the CTRW are more recent,and only a tentative strength weakness opportunity threat (SWOT)analysis can be presented at this time. We discuss the effect of the number of dimensions in a numerical model and the spacing between model nodes on solute spreading and the values of the solute-spreading parameters. In order to meet the increasing complexity of environmental problems, two approaches of model combination are used: Model integration and model coupling. Amain drawback of model integration is the complexity of there sulting code. Model coupling requires a systematic physical domain and model communication analysis. The setup and maintenance of a hydrologic framework for model coupling requires substantial resources, but on the other hand, contributions can be made by many research groups.

Studio dei metodi di valutazione della pericolosità associata ai "debris flow" in ambiente alpino

The present work concerns with the study of debris flows and, in particular, with the related hazard in the Alpine Environment. During the last years several methodologies have been developed to evaluate hazard associated to such a complex phenomenon, whose velocity, impacting force and inappropriate temporal prediction are responsible of the related high hazard level. This research focuses its attention on the depositional phase of debris flows through the application of a numerical model (DFlowz), and on hazard evaluation related to watersheds morphometric, morphological and geological characterization. The main aims are to test the validity of DFlowz simulations and assess sources of errors in order to understand how the empirical uncertainties influence the predictions; on the other side the research concerns with the possibility of performing hazard analysis starting from the identification of susceptible debris flow catchments and definition of their activity level. 25 well documented debris flow events have been back analyzed with the model DFlowz (Berti and Simoni, 2007): derived form the implementation of the empirical relations between event volume and planimetric and cross section inundated areas, the code allows to delineate areas affected by an event by taking into account information about volume, preferential flow path and digital elevation model (DEM) of fan area. The analysis uses an objective methodology for evaluating the accuracy of the prediction and involve the calibration of the model based on factors describing the uncertainty associated to the semi empirical relationships. The general assumptions on which the model is based have been verified although the predictive capabilities are influenced by the uncertainties of the empirical scaling relationships, which have to be necessarily taken into account and depend mostly on errors concerning deposited volume estimation. In addition, in order to test prediction capabilities of physical-based models, some events have been simulated through the use of RAMMS (RApid Mass MovementS). The model, which has been developed by the Swiss Federal Institute for Forest, Snow and Landscape Research (WSL) in Birmensdorf and the Swiss Federal Institute for Snow and Avalanche Research (SLF) takes into account a one-phase approach based on Voellmy rheology (Voellmy, 1955; Salm et al., 1990). The input file combines the total volume of the debris flow located in a release area with a mean depth. The model predicts the affected area, the maximum depth and the flow velocity in each cell of the input DTM. Relatively to hazard analysis related to watersheds characterization, the database collected by the Alto Adige Province represents an opportunity to examine debris-flow sediment dynamics at the regional scale and analyze lithologic controls. With the aim of advancing current understandings about debris flow, this study focuses on 82 events in order to characterize the topographic conditions associated with their initiation , transportation and deposition, seasonal patterns of occurrence and examine the role played by bedrock geology on sediment transfer.

Combined effects of heterogeneity, anisotropy, and saturation on steady state flow and transport: A laboratory sand tank experiment

Field soils show rather different spreading behavior at different water saturations, frequently caused by layering of the soil material. We performed tracer experiments in a laboratory sand tank. Such experiments complement and help comprehension of field investigations. We estimated, by image analysis, the first two moments of small plumes traveling through a two-dimensional, heterogeneous medium with strongly anisotropic correlation structure. Three steady state regimes were analyzed. Two main conclusions were drawn. First, low saturation led to very large heterogeneity and to strong preferential flow. Thus the description of the flow paths and the prediction of the solute arrival times require, in this case, more accurate knowledge about the topological structure. Second, saturation-dependent macroscopic anisotropy is an essential element of transport in unsaturated media. For this reason, small structural soil features should be properly upscaled to give appropriate effective soil parameters to be input in transport models.

Modelling susceptibility of grassland soil to macropore flow

Investigating preferential flow, including macropore flow, is crucial to predicting and preventing point sources of contamination in soil, for example in the vicinity of pumping wells. With a view to advancing groundwater protection, this study aimed (i) to quantify the strength of macropore flow in four representative natural grassland soils on the Swiss plateau, and (ii) to define the parameters that significantly control macropore flow in grassland soil. For each soil type we selected three measurement points on which three successive irrigation experiments were carried out, resulting in a total of 36 irrigations. The strength of macropore flow, parameterized as the cumulated water volume flowing from macropores at a depth of 1 m in response to an irrigation of 60 mm h−1 intensity and 1 h duration, was simulated using the dual-permeability MACRO model. The model calibration was based on the key soil parameters and fine measurements of water content at different depths. Modelling results indicate high performance of macropore flow in all investigated soil types except in gleysols. The volume of water that flowed from macropores and was hence expected to reach groundwater varied between 81% and 94% in brown soils, 59% and 67% in para-brown soils, 43% and 56% in acid brown soils, and 22% and 35% in gleysols. These results show that spreading pesticides and herbicides in pumping well protection zones poses a high risk of contamination and must be strictly prohibited. We also found that organic carbon content was not correlated with the strength of macropore flow, probably due to its very weak variation in our study, while saturated water content showed a negative correlation with macropore flow. The correlation between saturated hydraulic conductivity (Ks) and macropore flow was negative as well, but weak. Macropore flow appears to be controlled by the interaction between the bulk density of the uppermost topsoil layer (0–0.10 m) and the macroporosity of the soil below. This interaction also affects the variations in Ks and saturated water content. Further investigations are needed to better understand the combined effect of all these processes including the exchange between micropore and macropore domains.

In situ investigation of rapid subsurface flow : Temporal dynamics and catchment-scale implication

Non peer reviewed

A three dimensional numerical investigation of hillslope flow processes

Physically based distributed models of catchment hydrology are likely to be made available as engineering tools in the near future. Although these models are based on theoretically acceptable equations of continuity, there are still limitations in the present modelling strategy. Of interest to this thesis are the current modelling assumptions made concerning the effects of soil spatial variability, including formations producing distinct zones of preferential flow. The thesis contains a review of current physically based modelling strategies and a field based assessment of soil spatial variability. In order to investigate the effects of soil nonuniformity a fully three dimensional model of variability saturated flow in porous media is developed. The model is based on a Galerkin finite element approximation to Richards equation. Accessibility to a vector processor permits numerical solutions on grids containing several thousand node points. The model is applied to a single hillslope segment under various degrees of soil spatial variability. Such variability is introduced by generating random fields of saturated hydraulic conductivity using the turning bands method. Similar experiments are performed under conditions of preferred soil moisture movement. The results show that the influence of soil variability on subsurface flow may be less significant than suggested in the literature, due to the integrating effects of three dimensional flow. Under conditions of widespread infiltration excess runoff, the results indicate a greater significance of soil nonuniformity. The recognition of zones of preferential flow is also shown to be an important factor in accurate rainfall-runoff modelling. Using the results of various fields of soil variability, experiments are carried out to assess the validity of the commonly used concept of `effective parameters'. The results of these experiments suggest that such a concept may be valid in modelling subsurface flow. However, the effective parameter is observed to be event dependent when the dominating mechanism is infiltration excess runoff.

Complementary methods to investigate the development of clogging within a horizontal sub-surface flow tertiary treatment wetland

A combination of experimental methods was applied at a clogged, horizontal subsurface flow (HSSF) municipal wastewater tertiary treatment wetland (TW) in the UK, to quantify the extent of surface and subsurface clogging which had resulted in undesirable surface flow. The three dimensional hydraulic conductivity profile was determined, using a purpose made device which recreates the constant head permeameter test in-situ. The hydrodynamic pathways were investigated by performing dye tracing tests with Rhodamine WT and a novel multi-channel, data-logging, flow through Fluorimeter which allows synchronous measurements to be taken from a matrix of sampling points. Hydraulic conductivity varied in all planes, with the lowest measurement of 0.1 md1 corresponding to the surface layer at the inlet, and the maximum measurement of 1550 md1 located at a 0.4m depth at the outlet. According to dye tracing results, the region where the overland flow ceased received five times the average flow, which then vertically short-circuited below the rhizosphere. The tracer break-through curve obtained from the outlet showed that this preferential flow-path accounted for approximately 80% of the flow overall and arrived 8 h before a distinctly separate secondary flow-path. The overall volumetric efficiencyof the clogged system was 71% and the hydrology was simulated using a dual-path, dead-zone storage model. It is concluded that uneven inlet distribution, continuous surface loading and high rhizosphere resistance is responsible for the clog formation observed in this system. The average inlet hydraulic conductivity was 2 md1, suggesting that current European design guidelines, which predict that the system will reach an equilibrium hydraulic conductivity of 86 md1, do not adequately describe the hydrology of mature systems.

A comparison of in situ constant and falling head permeameter tests to assess the distribution of clogging within horizontal subsurface flow constructed wetlands

Clogging is the main operational problem associated with horizontal subsurface flow constructed wetlands (HSSF CWs). The measurement of saturated hydraulic conductivity has proven to be a suitable technique to assess clogging within HSSF CWs. The vertical and horizontal distribution of hydraulic conductivity was assessed in two full-scale HSSF CWs by using two different in situ permeameter methods (falling head (FH) and constant head (CH) methods). Horizontal hydraulic conductivity profiles showed that both methods are correlated by a power function (FH= CH 0.7821, r 2=0.76) within the recorded range of hydraulic conductivities (0-70 m/day). However, the FH method provided lower values of hydraulic conductivity than the CH method (one to three times lower). Despite discrepancies between the magnitudes of reported readings, the relative distribution of clogging obtained via both methods was similar. Therefore, both methods are useful when exploring the general distribution of clogging and, specially, the assessment of clogged areas originated from preferential flow paths within full-scale HSSF CWs. Discrepancy between methods (either in magnitude and pattern) aroused from the vertical hydraulic conductivity profiles under highly clogged conditions. It is believed this can be attributed to procedural differences between the methods, such as the method of permeameter insertion (twisting versus hammering). Results from both methods suggest that clogging develops along the shortest distance between water input and output. Results also evidence that the design and maintenance of inlet distributors and outlet collectors appear to have a great influence on the pattern of clogging, and hence the asset lifetime of HSSF CWs. © Springer Science+Business Media B.V. 2011.

Use of tritium and helium to define groundwater flow conditions in Everglades National Park

The concentrations of tritium (3H) and helium isotopes (3He and4He) were used as tracers of groundwater flow in the surficial aquifer system (SAS) beneath Everglades National Park (ENP), south Florida. From ages determined by 3H/3He dating techniques, groundwater within the upper 28 m originated within the last 30 years. Below 28 m, waters originated prior to 30 years before present with evidence of mixing at the interface. Interannual variation of the 3H/3He ages within the upper 28 m was significant throughout the 3 year investigation, corresponding with varying hydrologic conditions. In the region of Taylor Slough Bridge, younger groundwater was consistently detected below older groundwater in the Biscayne Aquifer, suggesting preferential flow to the lower part of the aquifer. An increase in 4He with depth in the SAS indicated that radiogenic 4He produced in the underlying Hawthorn Group migrates into the SAS by diffusion. Higher Δ4He values in brackish groundwaters compared to fresh waters from similar depths suggested a possible enhanced vertical transport of4He in the seawater mixing zone. Groundwater salinity measurements indicated the presence of a wide (6–28 km) seawater mixing zone. Comparison of groundwater levels with surface water levels in this zone indicated the potential for brackish groundwater discharge to the overlying Everglades surface water.