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.     

Temporal variability of solute transport under vadose zone conditions

The heterogeneity of the solute flux field in the horizontal plane at the field scale has been documented in several field studies. On the other hand, little information is available on the persistence of certain solute transport scenarios over consecutive infiltration cycles. This study was initiated to analyse the recurrence of solute leaching behaviour as estimated in two soil column tests emphasizing the preferential flow phenomenon. Twenty-four small-sized soil samples were subjected to two consecutive unsaturated steady-state flow leaching experiments with bromide as tracer. Observed breakthrough curves (BTCs) were analysed by the method of moments and by the advection–dispersion equation (ADE) to classify solute behaviour. Frequency distributions of the parameters indicating the solute velocity were heavily skewed or bimodal, reflecting the broad variability of the leaching scenarios, including some with pronounced preferential solute breakthrough. Exclusion of the preferential flow columns from our calculations revealed an average amount of 37% of immobile water. The large-scale BTCs derived from assembling the individual concentration courses of each run showed similar features, such as an early bromide breakthrough. However, two distinct apices, viz. one preferential and one matrix, were observed only in the first run, whereas the concentration decrease between the peaks was missing from the second run. A change in soil structure with continuous leaching was presumed to modify the interplay of the various flow domains, thereby altering the spreading of the BTCs. Correlation analysis between parameters of both tests suggests that preferential transport conditions are likely to occur at the same locations in the field over several infiltration cycles, whereas the ‘classical’ or expected matrix flow is time variant and therefore seems to be hardly predictable. © 1998 John Wiley & Sons, Ltd.     

Ground-Water Tracing with Injected Helium

Helium has several characteristics that make it attractive for use as a tracer in hydrological studies. Two types of experiments were conducted to investigate applicability of helium as a tracer of ground-water movement. The experiments included studies using laboratory sand and soil columns and field ground-water tracing in a basaltic aquifer. A water helium analyzer comprised of a thin quartz glass membrane and diode ion pump (making use of the preferential permeation of helium through the quartz glass into an evacuated space) was developed and used for the experiments. Results of our studies demonstrated that breakthrough curves of specific conductance and helium were similar under saturated conditions. In the unsaturated sand/soil columns, breakthrough curves of helium were retarded relative to specific conductance reducing the usefulness of helium as a tracer.     

Mechanisms and pathways of lateral flow on aspen‐forested,Luvisolic soils,Western Boreal Plains,Alberta, Canada

Rainfall simulation experiments by Redding and Devito ( 2008 , Hydrological Processes 23: 4287–4300) on two adjacent plots of contrasting antecedent soil moisture storage on an aspen‐forested hillslope on the Boreal Plain showed that lateral flow generation occurred only once large soil storage capacity was saturated combined with a minimum event precipitation of 15–20 mm. This paper extends the results of Redding and Devito ( 2008 , Hydrological Processes 23: 4287–4300) with detailed analysis of pore pressure, soil moisture and tracer data from the rainfall simulation experiments, which is used to identify lateral flow generation mechanisms and flow pathways. Lateral flow was not generated until soils were wet into the fine textured C horizon. Lateral flow occurred dominantly through the clay‐rich Bt horizon by way of root channels. Lateral flow during the largest event was dominated by event water, and precipitation intensity was critical in lateral flow generation. Lateral flow was initiated as preferential flow near the soil surface into root channels, followed by development of a perched water table at depth, which also interacted with preferential flow pathways to move water laterally by the transmissivity feedback mechanism. The results indicate that lateral flow generated by rainfall on these hillslopes is uncommon because of the generally high available soil moisture storage capacity and the low probability of rainfall events of sufficient magnitude and intensity. Copyright © 2010 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.     

Soil moisture dynamics at high temporal resolution in a semiarid Mediterranean watershed with scattered tree cover

Subsurface water flows play a key role in the distribution of water and solutes and thereby in the water availability for plants. However, the characterization of different flow processes (i.e. matrix and preferential flow), the frequency and factors that cause them, is relatively rare. This characterization enables a better understanding of spatio‐temporal variability of water resources and allows for the design of models to be improved. Using a method based on the time derivative of soil moisture variation known as maximum wetting slope, types of soil wetting processes were classified and quantified. For this, capacitance sensors, which registered the volumetric water content at high temporal resolution (30 min) for more than two hydrological years, were installed at different depths and placed in soil moisture stations with different vegetation covers, lithology and topographic position. Results indicated that there is a general behaviour or pattern of soil moisture dynamics in the catchment with a dominant occurrence of slower soil wetting processes (>50%), caused by matrix flows, and a low occurrence of those faster processes (<30%), originated by preferential flows. Nevertheless, when the total volume of water is considered, preferential flow becomes the dominant process, so that the ecological role of both flow types becomes prominent in water‐limited environments. Statistical multivariate analyses based on data‐mining techniques proved that although both flow types depend on variables associated with precipitation and antecedent soil moisture conditions, faster soil wetting processes are mainly related to variables such as rainfall intensity and topography, while slower soil wetting processes are related to flow velocity, soils or vegetation. Copyright © 2015 John Wiley & Sons, Ltd.     

Estimation of water saturation dependence of dispersion in unsaturated porous media: experiments and modelling analysis

In the dispersion theory, a linear relationship has been verified between the coefficient of hydrodynamic dispersion and water velocity, both in saturated and in unsaturated porous media. But for unsaturated soils the variability of flow directions and microscopic velocities can be larger than in saturated soils because of the lower degree of water saturation. This leads to an increased dispersion. Therefore, relationships between water content and relative water velocity fluctuations and water content together with the coefficient of dispersivity in unsaturated porous media respectively have been investigated systematically by displacement experiments in glass beads and coarse-textured sandy soil columns. The breakthrough curves (BTCs) of chloride showed that an increase of solute mixing with a decrease of water content was caused by an increase of flow velocity fluctuations for different pathways. In order to explain the observed tailing effect in unsaturated flow, two mathematical models were used to fit theoretically derived nonlinear functions of water content dependent dispersivities for both porous media. The close agreement between the observed and computed results suggests that the theoretical model of hydrodynamic dispersion can be extended to transport in unsaturated porous media, providing that BTCs of the effluent water are used to estimate representative dispersivity parameters of soils.     

Testing the ‘two water worlds’ hypothesis under variable preferential flow conditions

Widespread observations of ecohydrological separation are interpreted by suggesting that water flowing through highly conductive soil pores resists mixing with matrix storage over periods of days to months (i.e., two ‘water worlds’ exist). These interpretations imply that heterogeneous flow can produce ecohydrological separation in soils, yet little mechanistic evidence exists to explain this phenomenon. We quantified the separation between mobile water moving through preferential flow paths versus less mobile water remaining in the soil matrix after free-drainage to identify the amount of preferential flow necessary to maintain a two water world's scenario. Soil columns of varying macropore structure were subjected to simulated rainfall of increasing rainfall intensity (26 mm h⁻¹, 60 mm h⁻¹, and 110 mm h⁻¹) whose stable isotope signatures oscillated around known baseline values. Prior to rainfall, soil matrix water δ²H nearly matched the known value used to initially wet the pore space whereas soil δ¹⁸O deviated from this value by up to 3.4‰, suggesting that soils may strongly fractionate ¹⁸O. All treatments had up to 100% mixing between rain and matrix water under the lowest (26 mm h⁻¹) and medium (60 mm h⁻¹) rainfall intensities. The highest rainfall intensity (110 mm h⁻¹), however, reduced mixing of rain and matrix water for all treatments and produced significantly different preferential flow estimates between columns with intact soil structure compared to columns with reduced soil structure. Further, artificially limiting exchange between preferential flow paths and matrix water reduced bypass flow under the most intense rainfall. We show that (1) precipitation offset metrics such as lc-excess and d-excess may yield questionable interpretations when used to identify ecohydrological separation, (2) distinct domain separation may require extreme rainfall intensities and (3) domain exchange is an important component of macropore flow.     

FLOW SEPARATION IN UNDISTURBED SOIL USING MULTIPLE ANIONIC TRACERS. PART 3: UNSTEADY CORE-SCALE INFILTRATION EXPERIMENTS

Solute transport through structured, undisturbed soil has been studied in transient, unsaturated experiments using columns from grass and woodland sites on the Lancaster University campus. Three anionic tracers have been used, bromide (Br⁻) and two fluorinated organic acids (pentraflurobenzoic acid and 2,6-diflurobenzoic acid). The process of displacement of stored water from undisturbed columns was investigated using successive inputs of different tracers under similar antecedent conditions. The results indicated that initial breakthrough was rapid, with a relative concentration of 0.8 being reached between 0.4 and 0.5 pore volumes of discharge. It was found that there was an apparent continued discharge of ‘old’ water, stored in the column before any additions of tracer, even after the addition of a total of 4.9 and 5.4 pore volumes of water for the grass and woodland columns, respectively. The implications of the results of these tracer studies for modelling solute transport in structured soils are considered.     

FLOW SEPARATION IN UNDISTURBED SOILS USING MULTIPLE ANIONIC TRACERS. PART 1. ANALYTICAL METHODS AND UNSTEADY RAINFALL AND RETURN FLOW EXPERIMENTS

Tracers provide one of the few ways of obtaining realistic information on the flows of water and solutes in undisturbed structured soils. Three fluorinated organic acids [pentafluorobenzoic acid, 2,6-diflurobenzoic acid and o-(trifluoromethyl) benzoic acid] and bromide were tested as anionic tracers in situations where the separation of different flow components is of interest. The fluorobenzoates were relatively conservative (approximately 90% for loam soil) although, in some instances, co-elution or complexation may lead to apparently non-conservative behaviour. Tracer mixtures applied to soils, were separated and quantified by high-performance liquid chromatography. The relative mobilities of the tracers were studied during unidirectional steady flow in large undisturbed soil columns brought back to the laboratory. The breakthrough curves showed highly dispersive behaviour, with very early breakthrough of solute and a long tail. Contributions to the column outflow of different simultaneous applications of rainfall and return flow, labelled with different tracers, were separated using multicomponent mixing equations. The results show the importance of preferential flow and relatively immobile storage in the transport process.     

Lateral flow thresholds for aspen forested hillslopes on the Western Boreal Plain,Alberta, Canada

To predict the long‐term sustainability of water resources on the Boreal Plain region of northern Alberta, it is critical to understand when hillslopes generate runoff and connect with surface waters. The sub‐humid climate (P ≤ ET) and deep glacial sediments of this region result in large available soil storage capacity relative to moisture surpluses or deficits, leading to threshold‐dependent rainfall‐runoff relationships. Rainfall simulation experiments were conducted using large magnitude and high intensity applications to examine the thresholds in precipitation and soil moisture that are necessary to generate lateral flow from hillslope runoff plots representative of Luvisolic soils and an aspen canopy. Two adjacent plots (areas of 2·95 and 3·4 m²) of contrasting antecedent moisture conditions were examined; one had tree root uptake excluded for two months to increase soil moisture content, while the second plot allowed tree uptake over the growing season resulting in drier soils. Vertical flow as drainage and soil moisture storage dominated the water balances of both plots. Greater lateral flow occurred from the plot with higher antecedent moisture content. Results indicate that a minimum of 15–20 mm of rainfall is required to generate lateral flow, and only after the soils have been wetted to a depth of 0·75 m (C‐horizon). The depth and intensity of rainfall events that generated runoff > 1 mm have return periods of 25 years or greater and, when combined with the need for wet antecendent conditions, indicate that lateral flow generation on these hillslopes will occur infrequently. Copyright © 2008 John Wiley & Sons, Ltd.     

STUDY OF INFILTRATION INTO A HETEROGENEOUS SOIL USING MAGNETIC RESONANCE IMAGING

In this paper, the feasibility of using magnetic resonance imaging (MRI) to study water infiltration into a heterogeneous soil is examined, together with its difficulties and limitations. MRI studies of ponded water infiltration into an undisturbed soil core show that the combination of one- and two-dimensional imaging techniques provides a visual and non-destructive means of monitoring the temporal changes of soil water content and the moisture profile, and the movement of the wetting front. Two-dimensional images show air entrapment in repetitive ponded infiltration experiments. During the early stages of infiltration, one-dimensional images of soil moisture profiles clearly indicate preferential flow phenomena. The observed advance of wetting fronts can be described by a linear relationship between the square root of infiltration time (√t) and the distance of the wetting front from the soil surface. Similarly, the cumulative infiltration is also directly proportional to √t. Furthermore, from the MRI infiltration moisture profiles, it is possible to estimate the parameters that feature in infiltration equations. © 1997 by John Wiley & Sons, Ltd.     

On the controls of preferential flow in soils of different hillslope position and lithological origin

Soils derived from different lithologies and their controls on preferential flow remain underexplored in forested landscapes. In the same lithology, the propensity for preferential flow occurrence at different hillslope positions also remains largely elusive. By utilizing a soil moisture response time method, we compared preferential flow occurrence between a shale site (Shale Hills, silt loam soils) and a sandstone site (Garner Run, sandy loam soils) at four hillslope positions: ridge-top, North- and South-facing mid-slopes and toe slope, for over 2 years. The catchments are neighbouring and covered by temperate forest. For the four hillslope positions, Shale Hills had higher preferential flow frequencies compared to Garner Run. Between these two catchments, the South-facing mid-slope sites showed the highest contrasts in preferential flow frequency (33.5% of events at Shale Hills vs. 8.8% at Garner Run) while the ridge-top sites showed the lowest contrasts (18.7 vs. 13.2%). Additionally, over the unfrozen period, for seven out of eight monitoring sites, drier antecedent conditions tended to be more favourable for preferential flows to occur, with significant (p < .01) relationships at two sites. Except for the South-facing mid-slope sites, both Shale Hills and Garner Run had two preferential flow pathways. The characteristic preferential flow pathways at Shale Hills were the B_w and C horizons, and for Garner Run, preferential flow moved from the E/AE horizon to the B_w horizon. This study shows that shale-derived soils tended to have higher preferential flow occurrence than sandstone soils, but hillslope positions exhibit different levels of contrasts. More effort should be paid to study the impact of lithology on preferential flows in the context of land surface modelling and biogeochemical reactions to improve ecosystem services of headwater catchments.     

Stormflow generation involving pipe flow in a zero‐order basin of Peninsular Malaysia

Hydrological responses in a zero‐order basin (ZOB), a portion of whose discharge emerged via preferential flow through soil pipes, were examined over a 2‐year period in Peninsular Malaysia to elucidate primary stormflow generation processes. Silicon (Si) and specific conductance (EC) in various runoff components were also measured to identify their sources. ZOB flow response was dependent on antecedent precipitation amount; runoff increased linearly with precipitation during events >20 mm in relatively wet antecedent moisture conditions. Runoff derived from direct precipitation falling onto saturated areas accounted for <0·2% of total ZOB flow volume during the study period, indicating the predominance of subsurface pathways in ZOB flow. ZOB flow (high EC and low Si) was distinct from perennial baseflow via bedrock seepage (low EC and high Si) 5 m downstream of the ZOB outlet. Pipe flow responded quickly to ZOB flow rate and was characterized by a threshold flow capacity unique to each pipe. Piezometric data and pipe flow records demonstrated that pipes located deeper in the soil initiated first, followed by those at shallower depths; initiation of pipe flow corresponded to shallow groundwater rise above the saprolite‐soil interface. Chemical signatures of pipe flow were similar to each other and to the ZOB flow, suggesting that the sources were well‐mixed soil‐derived shallow groundwater. Based upon the volume of pipe flow during storms, the combined contribution of the pipes monitored accounted for 48% of total ZOB flow during the study period. Our results suggest that shallow groundwater, possibly facilitated by preferential flow accreted above the saprolite–soil interface, provides dominant stormflow, and that soil pipes play an important role in the rapid delivery of solute‐rich water to the stream system. Copyright © 2006 John Wiley & Sons, Ltd.     

Vertical surface water–groundwater exchange processes within a headwater floodplain induced by experimental floods

Restoring hydrologic connectivity between channels and floodplains is common practice in stream and river restoration. Floodplain hydrology and hydrogeology impact stream hydraulics, ecology, biogeochemical processing, and pollutant removal, yet rigorous field evaluations of surface water–groundwater exchange within floodplains during overbank floods are rare. We conducted five sets of experimental floods to mimic floodplain reconnection by pumping stream water onto an existing floodplain swale. Floods were conducted throughout the year to capture seasonal variation and each involved two replicate floods on successive days to test the effect of varying antecedent moisture. Water levels and specific conductance were measured in surface water, soil, and groundwater within the floodplain, along with surface flow into and out of the floodplain. Vegetation density varied seasonally and controlled the volume of surface water storage on the floodplain. By contrast, antecedent moisture conditions controlled storage of water in floodplain soils, with drier antecedent moisture conditions leading to increased subsurface storage and slower flood wave propagation across the floodplain surface. The site experienced spatial heterogeneity in vertical connectivity between surface water and groundwater across the floodplain surface, where propagation of hydrostatic pressure, preferential flow, and bulk Darcy flow were all mechanisms that may have occurred during the five floods. Vertical connectivity also increased with time, suggesting higher frequency of floodplain inundation may increase surface water–groundwater exchange across the floodplain surface. Understanding the variability of floodplain impacts on water quality noted in the literature likely requires better accounting for seasonal variations in floodplain vegetation and antecedent moisture as well as heterogeneous exchange flow mechanisms. Copyright © 2016 John Wiley & Sons, Ltd.     

Solute transport and water content measurements in clay soils using time domain reflectometry

Abstract

Clayey and saline soils have been shown to be problematic for time domain reflectometry (TDR) measurements. This study presents some of these problems and discusses solutions to them. Thirteen solute transport experiments were carried out in three undisturbed soil columns of swelling clay soil from Tunisia, labelled S1, S2, and S3 respectively. The columns were collected at three different physiographical regions within a catchment. Water fluxes ranged from 1.2 to 7.2 cm day^?1. The large solute transport heterogeneity and large tailing indicated that preferential flow was most pronounced in S1. The preferential flow took place in voids between structural elements and in wormholes. In S3, preferential flow was also evident, but not to the same extent as in S1. In S2, the solute transport was more uniform with little preferential flow. The heterogeneity of the solute transport increased with the water flux in S1 and to a smaller extent in S3, whereas it remained constant in S2. In a previous dye experiment in the field, preferential flow in cracks was observed at those sites where S1 and S3 were collected. In the column experiments, preferential flow in these cracks was less due to the higher initial water content compared to the dye experiments, indicating that the desiccation cracks were closed by the swelling clay.     

Movement of Heavy Metals in Soil through Preferential Flow Paths under Different Rainfall Intensities

Preferential flow in soils deserves attention due to its potential role in accelerating the movement of contaminants to groundwater. This study investigates the movement of Cd, Cu, and Pb through preferential flow paths under different applied rainfall intensities. Artificial acid rain (pH of 4.1) containing CdCl₂, CuCl₂, and Pb(NO₃)₂ was applied to undisturbed soil and repacked sand columns at low and high intensities, and leachate metals and chloride concentrations were measured. Cd was found in the leachate at both low and high rates in all columns, while the increase in Cu concentrations in the leachate was detected only at the high rate of the undisturbed columns. Pb was retained in both columns. For undisturbed columns, the breakthrough curves of Cd and Cu were similar to those of Cl, showing early initial breakthrough by preferential flow and dependency on rainfall intensities. The Cd concentrations were detected in the leachate from repacked columns for high rate rainfall, implying that even homogeneous soil may not be perfectly able to retain metals and the initially wet condition is more harmful for subsurface contaminant transport. In conclusion, the study demonstrated that, despite its highly sorptive nature, the transport of some metals may be as fast as that of a tracer under preferential flow conditions, and the rainfall intensity is a significant factor for the degree of transport.     

FLOW SEPARATION IN UNDISTURBED SOIL USING MULTIPLE ANIONIC TRACERS. PART 2. STEADY-STATE CORE-SCALE RAINFALL AND RETURN FLOWS AND DETERMINATION OF DISPERSION PARAMETERS

A series of experiments designed to study the separation of flow components from two large undisturbed cores under steady-state rainfall (downward) and return (upward) flows under near-saturated conditions is summarized. The experiments were conducted on soil columns collected from Lancaster University and the Slapton Wood catchment, Devon. The use of the relatively conservative tracers, potassium bromide, o-(trifluoromethyl)benzoic acid and 2,6-difluorobenzoic acid and a combination of application rates made it possible to quantify the different sources of water contributing to the discharge hydrographs. There is significant retention of tracer within the cores, despite the application of several pore volumes of water. The use of steady flow conditions allowed the determination of dispersion coefficients, dispersivity and proportion of ‘mobile’ water content parameters of the advection–dispersion equation. It was found that there were significant differences between the dispersivities at different flow-rates under upward and downward flux conditions and that in the undisturbed cores studied here the apparent proportions of mobile pore water ranged between 0.33 and 1.0, with an apparently complex relationship to flux rate. Prediction of transport in undisturbed soil remains problematic and tracer experiments will continue to be needed to provide a fundamental understanding of the complex flow processes involved.     

Tracer and hydrometric study of preferential flow in large undisturbed soil cores from the Georgia Piedmont,USA

We studied the temporal patterns of tracer throughput in the outflow of large (30 cm diameter by 38 cm long) undisturbed cores from the Panola Mountain Research Watershed, Georgia. Tracer breakthrough was affected by soil structure and rainfall intensity. Two rainfall intensities (20 and 40 mm hr⁻¹) for separate Cl⁻ and Br⁻ amended solutions were applied to two cores (one extracted from a hillslope soil and one extracted from a residual clay soil on the ridge). For both low and high rainfall intensity experiments, preferential flow occurred in the clay core, but not in the hillslope core. The preferential flow is attributed to well‐developed interpedal macrochannels that are commonly found in structured clay soils, characteristic of the ridge site. However, each rainfall intensity exceeded the matrix infiltration capacity at the top of the hillslope core, but did not exceed the matrix infiltration capacity at the middle and bottom of the hillslope core and at all levels in the clay core. Localized zones of saturation created when rainfall intensity exceeds the matrix infiltration capacity may cause water and tracer to overflow from the matrix into macrochannels, where preferential flow occurs to depth in otherwise unsaturated soil. Copyright © 1999 John Wiley & Sons, Ltd.