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.

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.

Phosphorus retention in forest soils and the functioning of buffer zones used in forestry

Phosphorus (P) retention properties of soils typical for boreal forest, i.e. podzolic soil and peat soils, vary significantly, but the range of this variation has not been sufficiently documented. To assess the usefulness of buffer zones used in forestry in removing P from the discharge by chemical sorption in soil, and to estimate the risk of P leaching after forestry operations, more data is needed on soil P retention properties. P retention properties of soils were studied at clear-cut areas, unharvested buffer zones adjoining the clear-cut and at peatland buffer zone areas. Desorption-sorption isotherms were determined for the humus layer, the mineral soil horizons E, B and C of the Podzol profile and for the surface layer peat (0-15 cm) and the subsurface layer peat (15-30 cm). The efficiency of buffer zones in retaining P was studied at six peatland buffer zone areas by adding P-containing solute in the inflow. A tracer study was conducted at one of the buffer zone areas to determine the allocation of the added P in soil and vegetation. Measured sorption or desorption rather than parameter values of fitted sorption equations described P desorption and sorption behaviour in soil. The highest P retention efficiency was in the B horizon and consequently, if contact occurred or was established between the soluble P in the water and the soil B horizon, the risk of P leaching was low. Humus layer was completely incapable of retaining P after clear-cutting. In the buffer zones, the decrease in P retention properties in the humus layer and the low amount of P sorbed by it indicated that the importance of the layer in the functioning of buffer zones is low. The peatland buffer zone areas were efficient in retaining soluble P from inflow. P sorption properties of the peat soil at the buffer zone areas varied largely but the contribution of P sorption in the peat was particularly important during high flow in spring, when the vegetation was not fully developed. Factors contributing to efficient P retention were large buffer size and low hydrological load whereas high hydrological load combined with the formation of preferential flow paths, especially during early spring or late autumn was disadvantageous. However, small buffer zone areas, too, may be efficient in reducing P load.

Physical properties of a boreal clay soil under differently managed perennial vegetation

The physical properties of surface soil horizons, essentially pore size, shape, continuity and affinity for water, regulate water entry into the soil. These properties are prone to changes caused by natural forces and human activity. The hydraulic properties of the surface soil greatly impact the generation of surface runoff and accompanied erosion, the major concern of agricultural water protection. The general target of this thesis was to improve our understanding of the structural and hydraulic properties of boreal clay soils. Physical properties of a clayey surface soil (0 - 10 cm, clay content 51%), with a micaceous/illitic mineralogy subjected to three different management practices of perennial vegetation, were studied. The study sites were vegetated buffer zones located side by side in SW Finland: 1) natural vegetation with no management, 2) harvested once a year, and 3) grazed by cattle. The soil structure, hydraulic properties, shrinkage properties and soil water repellency were determined at all sites. Two distinct flow domains were evident. The surface soil was characterized by subangular blocky, angular blocky and platy aggregates. Hence, large, partially accommodated, irregular elongated pores dominated the macropore domain at all sites. The intra-aggregate pore system was mostly comprised of pores smaller than 30 μm, which are responsible for water storage. Macropores at the grazed site, compacted by hoof pressure, were horizontally oriented and pore connectivity was poorest, which decreased water and air flux compared with other sites. Drying of the soil greatly altered its structure. The decrease in soil volume between wet and dry soil was 7 - 10%, most of which occurred in the moisture range of field conditions. Structural changes, including irreversible collapse of interaggregate pores, began at matric potentials around -6 kPa indicating, instability of soil structure against increasing hydraulic stress. Water saturation and several freezethaw cycles between autumn and spring likely weakened the soil structure. Soil water repellency was observed at all sites at the time of sampling and when soil was dryer than about 40 vol.%. (matric potential < -6 kPa). Therefore, water repellency contributes to water flow over a wide moisture range. Water repellency was also observed in soils with low organic carbon content (< 2%), which suggests that this phenomenon is common in agricultural soils of Finland due to their relatively high organic carbon content. Aggregate-related pedofeatures of dense infillings described as clay intrusions were found at all sites. The formation of these intrusions was attributed to clay dispersion and/or translocation during spring thaw and drying of the suspension in situ. These processes generate very new aggregates whose physical properties are most probably different from those of the bulk soil aggregates. Formation of the clay infillings suggested that prolonged wetness in autumn and spring impairs soil structure due to clay dispersion, while on the other hand it contributes to the pedogenesis of the soil. The results emphasize the dynamic nature of the physical properties of clay soils, essentially driven by their moisture state. In a dry soil, fast preferential flow is favoured by abundant macropores including shrinkage cracks and is further enhanced by water repellency. Increase in soil moisture reduces water repellency, and swelling of accommodated pores lowers the saturated hydraulic conductivity. Moisture- and temperature-related processes significantly alter soil structure over a time span of 1 yr. Thus, the pore characteristics as well as the hydraulic properties of soil are time-dependent.