Field study of macropore flow processes using tension infiltration of a dye tracer in partially saturated soils

Macropores are important preferential pathways for the migration of water and contaminants through the vadose zone. The objective of this study was to examine small‐scale preferential flow processes during infiltration in macroporous, low permeability soils. A series of tension infiltration tests were conducted using Brilliant Blue dye tracer at two field sites in southwestern Ontario, Canada. The maximum applied pressure head was varied for each test and the resulting dye stain patterns and macropore networks were characterized by excavation, mapping, photography, and image analysis. Worm burrows were the dominant macropore type, with average macropore densities exceeding 400 m^?2 and peak densities of more than 750 m^?2 at 30 cm depth at both sites. Flow in macropores became significant at infiltration pressures > ? 3 cm, with corresponding increases in infiltration rate, soil water content variability (spatially and temporally), and depth of dye staining. The results demonstrated clear evidence for partially saturated macropore flow under porewater tension conditions and the associated importance of macropore–matrix interaction in controlling this flow. Field observations of transient infiltration showed that film and rivulet flow along macropores yielded vertical flow velocities exceeding 40 m d^?1. Simple calculations showed that film flow along the walls and corners of irregularly shaped macropores could explain the observed results. Copyright © 2009 John Wiley & Sons, Ltd.     

Effects of antecedent moisture and macroporosity on infiltration and water flow in frozen soil

Infiltration into frozen soil plays an important role in soil freeze–thaw and snowmelt-driven hydrological processes. To better understand the complex thermal energy and water transport mechanisms involved, the influence of antecedent moisture content and macroporosity on infiltration into frozen soil was investigated. Ponded infiltration experiments on frozen macroporous and non-macroporous soil columns revealed that dry macroporous soil produced infiltration rates reaching 10³ to 10⁴ mm day⁻¹, two to three orders of magnitude larger than dry non-macroporous soil. Results suggest that rapid infiltration and drainage were a result of preferential flow through initially air-filled macropores. Using recorded flow rates and measured macropore characteristics, calculations indicated that a combination of both saturated flow and unsaturated film flow likely occurred within macropores. Under wet conditions, regardless of the presence of macropores, infiltration was restricted by the slow thawing rate of pore ice, producing infiltration rates of 2.8 to 5.0 mm day⁻¹. Reduced preferential flow under wet conditions was attributed to a combination of soil swelling, due to smectite-rich clay (that reduced macropore volume), and pore ice blockage within macropores. In comparison, dry soil column experiments demonstrated that macropores provided conduits for water and thermal energy to bypass the frozen matrix during infiltration, reducing thaw rates compared with non-macroporous soils. Overall, results showed the dominant control of antecedent moisture content on the initiation, timing, and magnitude of infiltration and flow in frozen macroporous soils, as well as the important role of macropore connectivity. The study provides an important data set that can aid the development of hydrological models that consider the interacting effects of soil freeze–thaw and preferential flow on snowmelt partitioning in cold regions.     

Effect of macropores on soil freezing and thawing with infiltration

An understanding of heat transport and water flow in unsaturated soils experiencing freezing and thawing is important when considering hydrological and thermal processes in cold regions. Macropores, such as cracks, roots, and animal holes, provide efficient conduits for enhanced infiltration, resulting in a unique distribution of water content. However, the effects of macropores on soil freezing and thawing with infiltration have not been well studied. A one‐directional soil‐column freezing and thawing experiment was conducted using unsaturated sandy and silt loams with different sizes and numbers of macropores. During freezing, macropores were found to retard the formation of the frozen layer, depending on their size and number. During thawing, water flowed through macropores in the frozen layer and reached the underlying unfrozen soil. However, infiltrated water sometimes refroze in a macropore. The ice started to form at near inner wall of the macropore, grew to the centre, and blocked flow through the macropore. The blockage ice in the macropore could not melt until the frozen layer disappeared. Improving a soil freezing model to consider these macropore effects is required. Copyright © 2016 John Wiley & Sons, Ltd.     

Sensitivity analysis on processes affecting bypass flow

Bypass flow in structured soils is dominated by soil hydrological processes, such as rain intensity, initial pressure head of the soil, surface storage of rain, horizontal contact area and absorption rate, and hydraulic conductivity of the soil matrix. This study was conducted to determine the relative impact of these processes in different soil types. A quasi 3-dimensional simulation model was used to calculate the effects of these soil hydrological input parameters on surface infiltration, macropore flow (with related horizontal absorption) and drainage. For light textured soils, surface infiltration was the most important term in the water balance. Heavy textured soils, in contrast, had drainage as the main term. In the latter soils bypass flow, when occurring, was almost equal to the amount of rain applied, indicating that absorption processes were strongly reduced. Lateral absorption on macropore walls was a minor fraction in the total mass balances, due to limited contact area and relatively weak diffusivity forces. Surface infiltration is a crucial parameter in bypass flow and is mainly dependent on rain intensity, initial pressure head and conductivity of the soil matrix. This requires measurement methods for hydraulic conductivity that specifically consider the effect of macropores.     

Synergistic effects of water repellency and macropore flow on the hydraulic conductivity of a burned forest soil,south‐east Australia

Research shows that water repellency is a key hydraulic property that results in reduced infiltration rates in burned soils. However, more work is required in order to link the hydrological behaviour of water repellent soils to observed runoff responses at the plot and hillslope scale. This study used 5 M ethanol and water in disc infiltrometers to quantify the role of macropore flow and water repellency on spatial and temporal infiltration patterns in a burned soil at plot (<10 m²) scale in a wet eucalypt forest in south‐east Australia. In the first summer and winter after wildfire, an average of 70% and 60%, respectively, of the plot area was water repellent and did not contribute to infiltration. Macropores (r > 0·5 mm), comprising just 5·5% of the soil volume, contributed to 70% and 95%, respectively, of the field‐saturated and ponded hydraulic conductivity (K_p). Because flow occurred almost entirely via macropores in non‐repellent areas, this meant that less than 2·5% of the soil surface effectively contributed to infiltration. The hydraulic conductivity increased by a factor of up to 2·5 as the hydraulic head increased from 0 to 5 mm. Due to the synergistic effect of macropore flow and water repellency, the coefficient of variation (CV) in K_p was three times higher in the water‐repellent soil (CV = 175%) than under the simulated non‐repellent conditions (CV = 66%). The high spatial variability in K_p would act to reduce the effective infiltration rate during runoff generation at plot scale. Ponding, which tend to increase with increasing scale, activates flow through macropores and would raise the effective infiltration rates at larger scales. Field experiments designed to provide representative measurements of infiltration after fire in these systems must therefore consider both the inherent variability in hydraulic conductivity and the variability in infiltration caused by interactions between surface runoff and hydraulic conductivity. Copyright © 2010 John Wiley & Sons, Ltd.     

Topographic controls upon soil macropore flow

Macropores are important components of soil hydrology. The spatial distribution of macropore flow as a proportion of saturated hydraulic conductivity was tested on six humid–temperate slopes using transects of tension infiltrometer measurements. Automated water table and overland flow monitoring allowed the timing of, and differentiation between, saturation‐excess overland flow and infiltration‐excess overland flow occurrence on the slopes to be determined and related to tension‐infiltrometer measurements. Two slopes were covered with blanket peat, two with stagnohumic gleys and two with brown earth soils. None of the slopes had been disturbed by agricultural activity within the last 20 years. This controlled the potential for tillage impacts on macropores. The proportion of near‐surface macropore flow to saturated hydraulic conductivity was found to vary according to slope position. The spatial patterns were not the same for all hillslopes. On the four non‐peat slopes there was a relationship between locations of overland flow occurrence and reduced macroporosity. This relationship did not exist for the peat slopes investigated because they experienced overland flow across their whole slope surfaces. Nevertheless, they still had a distinctive spatial pattern of macropore flow according to slope position. For the other soils tested, parts of slopes that were susceptible to saturation‐excess overland flow (e.g. hilltoes or flat hilltops) tended to have least macropore flow. To a lesser extent, for the parts of slopes susceptible to infiltration‐excess overland flow, the proportion of macropore flow as a component of infiltration was also smaller compared with the rest of the slope. The roles of macropore creation and macropore infilling by sheet wash are discussed, and it is noted that the combination of these may result in distinctive topographically controlled spatial patterns of macropore flow. Copyright © 2008 John Wiley & Sons, Ltd.     

A model to couple overland flow and infiltration into macroporous vadose zone

Most vegetated land surfaces contain macropores that may have a significant effect on the rate of infiltration of water under ponded conditions on the ground surface. Owing to the small-scale variations of the land topography (microtopography), only portions of the land area may get ponded during the process of overland flow. As the macropores transmit water at much higher rates than the primary soil matrix, higher macropore activation in ponded areas produces larger effective infiltration rates into the soil. Therefore, overland flow and infiltration into the macroporous vadose zone are interrelated. Representing the microtopographic variation of the land surface by a simple sine wave function, a method was developed to relate the ponding area to the average ponding depth which was determined by overland flow. A numerical model coupling overland flow and infiltration into the macroporous vadose zone was developed. Overland flow was simulated using the St. Venant equations with the inertia terms neglected. A single macropore model was used to simulate the infiltration into the macroporous vadose zone. The interaction between overland flow and the infiltration into the macroporous vadose zone was analyzed for a hypothetical watershed. The sensitivity analysis revealed that the interaction of macropore flow and overland flow is significant. For the conditions tested, the macropore flow and the overland flow were found to be more sensitive to the macroporosity and less sensitive to the microtopographic surface variation.     

Snowmelt‐driven macropore flow and soil saturation in a semiarid forest

Lateral subsurface flow is generally assumed to occur as a result of the development of a saturated zone above a low‐permeability interface such as at the soil–bedrock contact, and it is often augmented by macropore flow. Our objective was to evaluate the development of lateral subsurface flow and soil saturation at a semiarid ponderosa pine forest in New Mexico with respect to the conceptual model of saturation building above the soil–bedrock contact. At this site, we have long‐term observations of the water budget components, including lateral flow. A 1·5 m deep by 7 m long trench was constructed to observe lateral subsurface flow and development of saturation directly. Our observations are based on flow resulting from a melting snowdrift. The edge of the drift was about 7 m upslope from the trench. Lateral subsurface flow only occurred from root macropores in the Bt soil horizon. Saturation developed and grew outward from flowing root macropores, rather than growing upward from the soil–bedrock interface. This macropore‐centred saturation resulted in a highly heterogeneous distribution of water content until enough macropores began flowing and individual macropore saturated zones grew large enough to coalesce and saturate large volumes of the soil. Our observations are based on one snowmelt event and a relatively short hillslope flow path, and thus do not represent a full range of hydrologic conditions. Nevertheless, the observed behaviour did not conform to the traditional model of soil–bedrock control of saturation and lateral flow. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulation of water and chemicals in macropore soils Part 1. Representation of the equivalent macropore influence and its effect on soilwater flow

Macropores are a relatively small proportion of the soil volume, but they play an important role in the movement of water and chemicals owing to occasional rapid fluxes through them. The occurrence of macropore flow does not depend on the water content (or potential) of the bluk matrix unless the soil is close to saturation, but depends instead principally upon surface boundary conditions. Accordingly, three control situations of infiltration are recognized: macropore control, application control, and matrix control. These three situations indicate that the two-domain system may be a proper approach for the simulation of macropore soil. In this conceptualization, macropores are defined as channeling pores of different radii in which the flux density (with unit hydraulic gradient) occurring in the minimum sizes of such pores is greater than or equal to the saturated matrix hydraulic conductivity. Recognizing the two structural domains of the macropore and matrix, and possible water flow situations, three flow regions are suggested: matrix, macropore, and transaction. The matrix and the macropore are the two domains, and the transaction represents the exchange of water between the matrix and the macropore. The classic approach of the Richards equation is applicable to describe the flow in the matrix domain. The Hagen-Poiseuille and the Chezy-Manning equations for tube flow can be applied to represent the relationship between the hydraulic conductivity of the macroporosity and the total macroporosity, where the total macroporosity is defined as the ratio of the summed macropore cross-sectional area and the total soil cross-sectional area. An equation describing water flow in the macropore domain is then obtained.     

Modelling investigation of water partitioning at a semiarid ponderosa pine hillslope

The effects of vegetation root distribution on near‐surface water partitioning can be two‐fold. On the one hand, the roots facilitate deep percolation by root‐induced macropore flow; on the other hand, they reduce the potential for deep percolation by root‐water‐uptake processes. Whether the roots impede or facilitate deep percolation depends on various conditions, including climate, soil, and vegetation characteristics. This paper examines the effects of root distribution on deep percolation into the underlying permeable bedrock for a given soil profile and climate condition using HYDRUS modelling. The simulations were based on previously field experiments on a semiarid ponderosa pine (Pinus ponderosa) hillslope. An equivalent single continuum model for simulating root macropore flow on hillslopes is presented, with root macropore hydraulic parameterization estimated based on observed root distribution. The sensitivity analysis results indicate that the root macropore effect dominates saturated soil water flow in low conductivity soils (K_matrix below 10^?7 m/s), while it is insignificant in soils with a K_matrix larger than 10^?5 m/s, consistent with observations in this and other studies. At the ponderosa pine site, the model with simple root‐macropore parameterization reasonably well reproduces soil moisture distribution and some major runoff events. The results indicate that the clay‐rich soil layer without root‐induced macropores acts as an impeding layer for potential groundwater recharge. This impeding layer results in a bedrock percolation of less than 1% of the annual precipitation. Without this impeding layer, percolation into the underlying permeable bedrock could be as much as 20% of the annual precipitation. This suggests that at a surface with low‐permeability soil overlying permeable bedrock, the root penetration depth in the soil is critical condition for whether or not significant percolation occurs. Copyright © 2010 John Wiley & Sons, Ltd.     

Soil macroporosity and infiltration characteristics of a forest podzol

Spatial distribution of soil macroporosity was determined for a forest podzol from tension infiltrometer measurements at the soil surface. Surface‐derived macroporosity values were compared with point infiltration characteristics obtained from soil water content and soil water chemistry measurements during an experimental irrigation, and with parameters of a kinematic wave model applied to soil water content data. Macroporosity estimated by the tension infiltrometer ranged from 0·00087 to 0·0219% of soil volume, and infiltration at these two sites was dominated by propagation of a well‐defined wetting front through the soil profile and bypass flow via soil macropores, respectively. Infiltration at sites with intermediate macroporosities reflected a combination of these two processes, although results were inconclusive at one site owing to lateral flow at the base of the soil profile. There was no agreement between macroporosities estimated by the tension infiltrometer and the kinematic wave model. The maximum soil conductance parameter within the profile at a site, however, was related directly to the surface‐derived macroporosity. The partial agreement between surface‐derived macroporosity estimates and point infiltration characteristics shown here supports the use of tension infiltrometry as a rapid, non‐destructive method of assessing spatial variations in the relative contribution of macropore flow to the infiltration process. Copyright © 2000 John Wiley & Sons, Ltd.     

The influence of preferential flow on hillslope hydrology in a semi‐arid watershed (in the Spanish Dehesas)

Preferential flow is known to influence hillslope hydrology in many areas around the world. Most research on preferential flow has been performed in temperate regions. Preferential infiltration has also been found in semi‐arid regions, but its impact on the hydrology of these regions is poorly known. The aim of this study is to describe and quantify the influence of preferential flow on the hillslope hydrology from small scale (infiltration) to large scale (subsurface stormflow) in a semi‐arid Dehesa landscape. Precipitation, soil moisture content, piezometric water level and discharge data were used to analyse the hydrological functioning of a catchment in Spain. Variability of soil moisture content during the transition from dry to wet season (September to November) within horizontal soil layers leads to the conclusion that there is preferential infiltration into the soils. When the rainfall intensity is high, a water level rapidly builds up in the piezometer pipes in the area, sometimes even reaching soil surface. This water level also drops back to bedrock within a few hours (under dry catchment conditions) to days (under wet catchment conditions). As the soil matrix is not necessarily wet while this water layer is built up, it is thought to be a transient water table in large connected pores which drain partly to the matrix, partly fill up bedrock irregularities and partly drain through subsurface flow to the channels. When the soil matrix becomes wetter the loss of water from macropores to the matrix and bedrock decreases and subsurface stormflow increases. It may be concluded that the hillslope hydrological system consists of a fine matrix domain and a macropore domain, which have their own flow characteristics but which also interact, depending on the soil matrix and macropore moisture contents. The macropore flow can result in subsurface flow, ranging from 13% contribution to total discharge for a large event of high intensity rainfall or high discharge to 80% of total discharge for a small event with low intensity rainfall or low discharge. During large events the fraction of subsurface stormflow in the discharge is suppressed by the large amount of surface runoff. Copyright © 2008 John Wiley & Sons, Ltd.     

Differences in preferential flow with antecedent moisture conditions and soil texture: Implications for subsurface P transport

Preferential flowpaths transport phosphorus (P) to agricultural tile drains. However, if and to what extent this may vary with soil texture, moisture conditions, and P placement is poorly understood. This study investigated (a) interactions between soil texture, antecedent moisture conditions, and the relative contributions of matrix and preferential flow and (b) associated P distributions through the soil profile when fertilizers were applied to the surface or subsurface. Brilliant blue dye was used to stain subsurface flowpaths in clay and silt loam plots during simulated rainfall events under wet and dry conditions. Fertilizer P was applied to the surface or via subsurface placement to plots of different soil texture and moisture condition. Photographs of dye stains were analysed to classify the flow patterns as matrix dominated or macropore dominated, and soils within plots were analysed for their water‐extractable P (WEP) content. Preferential flow occurred under all soil texture and moisture conditions. Dye penetrated deeper into clay soils via macropores and had lower interaction with the soil matrix, compared with silt loam soil. Moisture conditions influenced preferential flowpaths in clay, with dry clay having deeper infiltration (92 ± 7.6 cm) and less dye–matrix interaction than wet clay (77 ± 4.7 cm). Depth of staining did not differ between wet (56 ± 7.2 cm) and dry (50 ± 6.6 cm) silt loam, nor did dominant flowpaths. WEP distribution in the top 10 cm of the soil profile differed with fertilizer placement, but no differences in soil WEP were observed at depth. These results demonstrate that large rainfall events following drought conditions in clay soil may be prone to rapid P transport to tile drains due to increased preferential flow, whereas flow in silt loams is less affected by antecedent moisture. Subsurface placement of fertilizer may minimize the risk of subsurface P transport, particularily in clay.     

Capabilities and limitations of detailed hillslope hydrological modelling

Hillslope hydrological modelling is considered to be of great importance for the understanding and quantification of hydrological processes in hilly or mountainous landscapes. In recent years a few comprehensive hydrological models have been developed at the hillslope scale which have resulted in an advanced representation of hillslope hydrological processes (including their interactions), and in some operational applications, such as in runoff and erosion studies at the field scale or lateral flow simulation in environmental and geotechnical engineering. An overview of the objectives of hillslope hydrological modelling is given, followed by a brief introduction of an exemplary comprehensive hillslope model, which stimulates a series of hydrological processes such as interception, evapotranspiration, infiltration into the soil matrix and into macropores, lateral and vertical subsurface soil water flow both in the matrix and preferential flow paths, surface runoff and channel discharge. Several examples of this model are presented and discussed in order to determine the model's capabilities and limitations. Finally, conclusions about the limitations of detailed hillslope modelling are drawn and an outlook on the future prospects of hydrological models on the hillslope scale is given.The model presented performed reasonable calculations of Hortonian surface runoff and subsequent erosion processes, given detailed information of initial soil water content and soil hydraulic conditions. The vertical and lateral soil moisture dynamics were also represented quite well. However, the given examples of model applications show that quite detailed climatic and soil data are required to obtain satisfactory results. The limitations of detailed hillslope hydrological modelling arise from different points: difficulties in the representations of certain processes (e.g. surface crusting, unsaturated–saturated soil moisture flow, macropore flow), problems of small‐scale variability, a general scarcity of detailed soil data, incomplete process parametrization and problems with the interdependent linkage of several hillslopes and channel–hillslope interactions. Copyright © 1999 John Wiley & Sons, Ltd.     

MACROPORE FLOW AND MASS WASTING OF GULLIES IN THE LOESS PLATEAU, CHINA

Due to the extensive gullying from historically excessive erosion in the loess plateau of China, much of this region is being converted to native grass and shrub vegetation. Tunnel scour and mass wasting are important gully erosion processes resulting from preferential flow through macropores （ pores 〉 1 mm diameter）. The objective of this study is to assess the changes with time in macropore flow characteristics of soils on the Loess Plateau following conversion to grass vegetation and the associated degree of mass wasting of gully faces. Ridge areas that had been revegetated for 1 year, 6 years, and 〉 15 years following tilling, and for 6 years following contour-ditching and the adjacent gully faces were characterized for their macropore and soil matrix properties on a 50 cm by 50 cm area. The total number of macropores increased from 11.6/m^2 to 39.6/m2 from 1 to 6 years and to 51.6/m2 after 15 years of revegetation following tillage. The macroporosity increased from 0.0008 m^3/m^3 to 0.0018 m^3/m^3 from 1 to 6 years of revegetation following tillage but the lowest macroporosity （0.0005 m3/m3） was 6 years of revegetation following contour-ditching. The contourditched area had the lowest infiltration rate （95 m/d） through the soil matrix （areas without macropores） with the tilled areas having similar infiltration rates regardless of the number of years of revegetation （averaged 146 m/d）. Due to tunnel scour erosion of macropores during infiltration into the area revegetated for 1 year, pore diameters enlarged by more than 200% resulting in this condition having the highest individual macropore infiltration rates （7967 m/d）. Macropores in all other areas were stable with no tunnel scour erosion of macropores. The total capacity for infiltration through macropores increased significantly with time following revegetation. The number of macropores on the gully faces was triple （92.8/m2） and the macroporosity quadruple （0.004 m3/m3） that of the ridge surfaces. The upper gully faces exhibited 1.1 slumps m^-1 for a total soil loss of 48622 kg per ha.     

Infiltration pattern in a regolith–fractured bedrock profile: field observation of a dye stain pattern

We examined the infiltration pattern of water in a regolith–bedrock profile consisting of two overburdens (OB1 and OB2), a buried rice paddy soil (PS), two texturally distinctive weathered materials (WM1 and WM2) and a fractured sedimentary rock (BR), using a Brilliant Blue FCF dye tracer. A black‐coloured coating in conducting fractures in WM1, WM2 and BR was analysed by X‐ray diffraction and scanning electron microscopy. The dye tracer penetrated to greater than 2 m depth in the profile. The macropore flow and saturated interflow were the major infiltration patterns in the profile. Macropore flow and saturated interflow were observed along fractures in WM1, WM2 and BR and at the dipping interfaces of PS–WM1, PS–WM2 and PS–BR respectively. Heterogeneous matrix flow occurred in upper overburden (OB1) and PS. Compared with OB1, the coarser textured OB2 acted as a physical barrier for vertical flow of water. The PS with low bulk density and many fine roots was another major conducting route of water in the profile. Manganese oxide and iron oxide were positively identified in the black coating material and had low crystallinity and high surface area, indicating their high reactivity with conducting contaminants. Copyright © 2006 John Wiley & Sons, Ltd.     

Soil hydrophobic effects on infiltration and catchment runoff

After dry summers or drought, eucalypt forest soils at two sites in southeastern Australia developed hydrophobic or non-wetting surface characteristics that reduced infiltration, measured using a sprinkling infiltrometer. At one site the development of hydrophobic conditions caused the rainfall to runoff conversion efficiency of a forested catchment to increase from 5 per cent to 15 per cent. Under non-hydrophobic conditions at this site, grassland always generated more runoff than forest. However, one major rainfall-runoff was recorded at a time of highly hydrophobic forest soil conditions and this storm generated greater runoff on the forested catchment than the grassland catchment. At the second site forest soils have naturally highly conductive surface layers because of a dense network of macropores and pathways for preferential flow. Hydrophobic conditions produced by drought caused soil water movement to be confined to only a few of the larger macropores exposed to surface ponded water. Even so, infiltration rates remained relatively high so that the impacts of hydrophobic soils were not translated into increased catchment runoff as at the first site.     

Microtopography,water storage and flow patterns in a fine‐textured soil under agricultural use

Agricultural use of soils implies tillage and often compaction and therefore influences processes on soil surface and affects infiltration of water into the subsoil. Although many studies on soil surface processes or flow patterns in soils exist, works relating both are rare in literature. We did two tracer experiments with Brilliant Blue FCF on a tilled and compacted plot and a non‐tilled one to investigate water storage on the soil surface during simulated rainfall and changes of soil microtopography, to analyse the associated flow patterns in the soil and to relate both to tillage and compaction. Our results show that storage was larger on the tilled and compacted plot than on the non‐tilled one. After tillage, transport processes above the plough pan were partly disconnected from those underneath because macropores were disrupted and buried by the tillage operation. However, preferential flow along cracks occurred on both plots and the macropores buried below the tillage pan still functioned as preferential flow paths. Therefore, we conclude that the studied soil is susceptible to deep vertical solute propagation at dry conditions when cracks are open, irrespective of tillage and compaction. Copyright © 2012 John Wiley & Sons, Ltd.