An infiltration model based on flow variability in macropores: development, sensitivity analysis and applications

Simulating infiltration in soils containing macropores still provides unsatisfactory results, as existing models seem not to capture all relevant processes. Recent studies of macropore flow initiation in natural soils containing earthworm channels revealed a distinct flow rate variability in the macropores depending on the initiation process. When macropore flow was initiated at the soil surface, most of the macropores received very little water while a few macropores received a large proportion of the total inflow. In contrast, when macropore flow was initiated from a saturated or nearly saturated soil layer, macropore flow rate variation was much lower. The objective of this study was to develop, evaluate, and test a model, which combines macropore flow variability with several established approaches to model dual permeability soils. We then evaluate the INfiltration–INitiation–INteraction Model (IN³M) to explore the influence of macropore flow variability on infiltration behavior by performing a sensitivity analysis and applying IN³M to sprinkling and dye tracer experiments at three field sites with different macropore and soil matrix properties. The sensitivity analysis showed that the flow variability in macropores reduces interaction between the macropores and the surrounding soil matrix and thus increases bypass flow, especially for surface initiation of macropore flow and at higher rainfall intensities. The model application shows reasonable agreement between IN³M simulations and field data in terms of water balance, water content change, and dye patterns. The influence of macropore flow variability on the hydrological response of the soil was considerable and especially pronounced for soils where initiation occurs at the soil surface. In future, the model could be applied to explore other types of preferential flow and hence to get a generally better understanding of macropore flow.     

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.     

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.     

ABSENCE OF PREFERENTIAL FLOW IN THE PERCOLATING WATERS OF A CONIFEROUS FOREST SOIL

Evidence for the functioning of macropores and the presence of preferential flow in forest soils is equivocal. This is partly because many workers use only one diagnostic technique to indicate whether or not macropore flow occurs. In this paper three lines of evidence are used to suggest that preferential flow does not occur in the percolating waters of a coniferous forest soil under the range of hydrological conditions that prevail in the field. To simulate field conditions, realistic rainfall intensities were used in conservative solute transport experiments on four undisturbed soil columns. A method is described in which breakthrough data can be used to calculate the percentage of antecedent water displaced from a soil column during frontal-type breakthrough experiments. Calculations based on this method using the experimental data show that as little as five percent of the antecedent water was immobile. The simple form of the functional advection–dispersion equation, based on a single value for linear velocity and the dispersion coefficient was fitted to two of the breakthrough curves with reasonable accuracy, further suggesting that preferential flow did not occur in the experiments. Finally, soil moisture characteristic curves were determined for replicate soil samples from the forest soil. The operational water contents of the columns during the breakthrough experiments were compared with the soil moisture characteristics and it was found that pores exerting pressure heads greater than −0·5 kPa did not appear to contribute to flow through the columns, again suggesting an absence of preferential flow. © 1997 by John Wiley & Sons, Ltd.     

Parametric study of a physically-based,plot-scale hillslope hydrological model through virtual experiments

Abstract

A physically-based hillslope hydrological model with shallow overland flow and rapid subsurface stormflow components was developed and calibrated using field experiments conducted on a preferential path nested hillslope in northeast India. Virtual experiments were carried out to perform sensitivity analysis of the model using the automated parameter estimation (PEST) algorithm. Different physical parameters of the model were varied to study the resulting effects on overland flow and subsurface stormflow responses from the theoretical hillslopes. It was observed that topographical shapes had significant effects on overland flow hydrographs. The slope profiles, surface storage, relief, rainfall intensity and infiltration rates primarily controlled the overland flow response of the hillslopes. Prompt subsurface stormflow responses were mainly dominated by lateral preferential flow, as soil matrix flow rates were very slow. Rainfall intensity and soil macropore structures were the most influential parameters on subsurface stormflow. The number of connected soil macropores was a more sensitive parameter than the size of macropores. In hillslopes with highly active vertical and lateral preferential pathways, saturation excess overland flow was not evident. However, saturation excess overland flow was generated if the lateral macropores were disconnected. Under such conditions, rainfall intensity, duration and preferential flow rate governed the process of saturation excess overland flow generation from hillslopes.

Editor D. Koutsoyiannis; Associate editor C. Perrin     

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.     

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.     

The relevance of preferential flow in catchment scale simulations: Calibrating a 3D dual-permeability model using DREAM

The occurrence of preferential flow in the subsurface has often been shown in field experiments. However, preferential flow is rarely included in models simulating the hydrological response at the catchment scale. If it is considered, preferential flow parameters are typically determined at the plot scale and then transferred to larger-scale simulations. Here, we successfully used the optimization algorithm DiffeRential Evolution Adaptive Metropolis (DREAM) to calibrate a 3D physics-based dual-permeability model directly at the catchment scale. In order to keep computational costs of the optimization routine at a reasonable level, we limited the number of parameters to be calibrated to the ones that had been shown before to be most influential for the simulation of discharge. We also calibrated parameters of the matrix domain and the macropore domain with a fixed parameter ratio between soil layers instead of calibrating every layer separately. These ratios reflected observed depth profiles of soil hydraulic properties at our study site. The dual-permeability parameter sets identified during calibration were able to simulate observed discharge time series satisfactorily but did not outperform a calibrated single-domain reference model scenario. Saturated hydraulic conductivities of the macropore domain were calibrated such that they became very similar to matrix saturated hydraulic conductivities, thereby effectively removing the effect of macropores. This suggests that the incorporation of vertical preferential flow as represented by the dual-permeability approach was not relevant for reproducing the hydrometric response reasonably well in the studied catchment. We also tested the scale-invariance of the calibrated dual-permeability parameter sets by using the parameter sets performing best at catchment scale to simulate plot-scale bromide depth profiles obtained from tracer irrigation experiments. This parameter transfer proved to be not successful, indicating that soil hydraulic parameters are scale-variant, independent of the direction of parameter transfer.     

Estimating macropore and matrix flow using the hydrograph separation procedure in an experimental forest plot

This study was conducted to estimate macropore space, macropore flow and matrix flow in an experimental forest plot in the Ouachita Mountains of Arkansas. Lateral soil water fluxes and soil capillary potentials were observed in the isolated plot during applied rainfall experiments. Rainfalls were applied 17 times during the period 17 July to 10 October 1991. The subsurface hydrograph separation technique was used to estimate macropore space, macropore flux and matrix flux. The boundary between macropore and matrix flow was statistically determined by covariance analysis. The maximum estimated lateral macropore space was approximately 0.006 (cm³ cm^?3). The maximum estimated lateral macropore and matrix flow were 0.042 and 0.00066 cm s^?1, respectively. This report also emphasizes the need for further research on the hydrograph separation procedure for estimating macropores and macropore flow.     

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.     

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.     

A physically based distributed subsurface–surface flow dynamics model for forested mountainous catchments

This study was designed to develop a physically based hydrological model to describe the hydrological processes within forested mountainous river basins. The model describes the relationships between hydrological fluxes and catchment characteristics that are influenced by topography and land cover. Hydrological processes representative of temperate basins in steep terrain that are incorporated in the model include intercepted rainfall, evaporation, transpiration, infiltration into macropores, partitioning between preferential flow and soil matrix flow, percolation, capillary rise, surface flow (saturation‐excess and return flow), subsurface flow (preferential subsurface flow and baseflow) and spatial water‐table dynamics. The soil–vegetation–atmosphere transfer scheme used was the single‐layer Penman–Monteith model, although a two‐layer model was also provided. The catchment characteristics include topography (elevation, topographic indices), slope and contributing area, where a digital elevation model provided flow direction on the steepest gradient flow path. The hydrological fluxes and catchment characteristics are modelled based on the variable source‐area concept, which defines the dynamics of the watershed response. Flow generated on land for each sub‐basin is routed to the river channel by a kinematic wave model. In the river channel, the combined flows from sub‐basins are routed by the Muskingum–Cunge model to the river outlet; these comprise inputs to the river downstream. The model was applied to the Hikimi river basin in Japan. Spatial decadal values of the normalized difference vegetation index and leaf area index were used for the yearly simulations. Results were satisfactory, as indicated by model efficiency criteria, and analysis showed that the rainfall input is not representative of the orographic lifting induced rainfall in the mountainous Hikimi river basin. Also, a simple representation of the effects of preferential flow within the soil matrix flow has a slight significance for soil moisture status, but is insignificant for river flow estimations. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

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.     

The role of overland flow and subsurface flow on the spatial distribution of soil moisture in the topsoil

Many investigations show relationships between topographical factors and the spatial distribution of soil moisture in catchments. However, few quantitative analyses have been carried out to elucidate the role of different hydrological processes in the spatial distribution of topsoil moisture in catchments. A spatially distributed rainfall—runoff model was used to investigate contributions of subsurface matric flow, macropore flow and surface runoff to the spatial distribution of soil moisture in a cultivated catchment. The model results show that lateral subsurface flow in the soil matrix or in macropores has a minor effect on the spatial distribution of soil moisture. Only when a perched groundwater table is maintained long enough, which is only possible if the subsurface is completely impermeable, may a spatial distribution in moisture content occur along the slope. Surface runoff, producing accumulations of soil moisture in flat flow paths of agricultural origin (field boundaries), was demonstrated to cause significant spatial variations in soil moisture within a short period after rainfall (<2 days). When significant amounts of surface runoff are produced, wetter moisture conditions will be generated at locations with larger upstream contributing areas. Copyright © 2001 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.     

Using time‐domain reflectometry to characterize shallow solute transport in an oak woodland hillslope in northern California,USA

The natural heterogeneity of water and solute movement in hillslope soils makes it difficult to accurately characterize the transport of surface‐applied pollutants without first gathering spatially distributed hydrological data. This study examined the application of time‐domain reflectometry (TDR) to measure solute transport in hillslopes. Three different plot designs were used to examine the transport of a conservative tracer in the first 50 cm of a moderately sloping soil. In the first plot, which was designed to examine spatial variability in vertical transport in a 1·2 m² plot, a single probe per meter was found to adequately characterize vertical solute travel times. In addition, a dye and excavation study in this plot revealed lateral preferential flow in small macropores and a transport pattern where solute is focused vertically into preferential flow pathways. The bypass flow delivers solute deeper in the soil, where lateral flow occurs. The second plot, designed to capture both vertical and lateral flow, provided additional evidence confirming the flow patterns identified in the excavation of the first plot. The third plot was designed to examine lateral flow and once again preferential flow of the tracer was observed. In one instance rapid solute transport in this plot was estimated to occur in as little as 3% of the available pore space. Finally, it was demonstrated that the soil anisotropy, although partially responsible for lateral subsurface transport, may also homogenize the transport response across the hillslope by decreasing vertical solute spreading. Copyright © 2002 John Wiley & Sons, Ltd.     

Dual-porosity and kinematic wave approaches to assess the degree of preferential flow in an unsaturated soil

Abstract

The purpose of this study was to assess the degree of preferential flow in an unsaturated soil column using two different models: the dual-porosity model, MACRO, and the kinematic wave approach (KWA) based on boundary-layer flow theory. The soil column experiments consisted of six infiltrations with intensities varying from 15 to 101 mm h^?1. Bromide solution was also infiltrated at an intensity of 79 mm h^?1 and a concentration of 80 mg l^?1. Both MACRO and the KWA indicated the absence of pure preferential flow. The KWA indicated intermediate flow with dispersion of the wetting front with depth, whereas MACRO indicated flow dominated by the diffusion of capillary potential. These results shed light on the transition between flows dominated by momentum dissipation and by diffusion of capillary potential. The absence of pure macropore flow in the structured sandy soil is mainly due to efficient lateral mass exchange in this material.