Elsevier

Journal of Hydrology

Volume 310, Issues 1–4, 1 August 2005, Pages 294-315
Journal of Hydrology

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

https://doi.org/10.1016/j.jhydrol.2005.01.010Get rights and content

Abstract

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 (IN3M) to explore the influence of macropore flow variability on infiltration behavior by performing a sensitivity analysis and applying IN3M 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 IN3M 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.

Introduction

Runoff generation at a site determines the flow pathways in catchments and thus the flow hydrograph for rainfall events. Only a realistic simulation and prediction of flow pathways leads to a correct description of the internal hydrological behavior in a catchment (Beven, 2001). Although the prediction of the flow hydrograph with rainfall-runoff models, which do not account for the relevant runoff processes, may be acceptable, we get ‘the right results for the wrong reason’ (Klemes, 1986). A correct description of the internal hydrological behavior is especially important for the simulation of solute transport processes in the catchment (e.g. nitrate or phosphate leaching) and for the assessment of land-use or climate change scenarios. The influence of macropore flow on hillslope runoff generation is well established (Bronstert, 1999, Faeh et al., 1997, McDonnell, 1990, Návar et al., 1995, Smettem et al., 1991, Weiler et al., 1998, Weiler and McDonnell, 2004), particularly for extreme rainfall events with either a very high intensity for a short duration or a lower intensity for a longer duration. (e.g. Faeh et al., 1997, Niehoff et al., 2002). Despite abundant experimental evidence of its significance, macropore flow, a subtype of preferential flow, is a process that is often not considered in rainfall-runoff models. Indirect evidence of macropore flow at the catchment scale can be obtained through model interpretation (Germann and Beven, 1981, Burch et al., 1989, Niehoff et al., 2002, Beckers and Alila, 2004).

Macropore flow causes a rapid downward movement of water in structural pore spaces such as worm channels, shrinkage cracks, and root holes subsequently bypassing portions of the soil profile. For example, earthworm activity, particularly the anecic earthworm species Lumbricus terrestris, causes macropores in soils in humid climate (Langmaack et al., 1999). Macropore flow influences the infiltration of rainfall and therefore the runoff generation and solute transport in natural soils where these structures are common (Larsson, 1999). The impact of macropores is governed to a large extent by water supply to macropores, water flow in macropores, and the water transfer from the macropores into the surrounding soil matrix (Beven and Germann, 1982, Faeh et al., 1997, Buttle and Leigh, 1997). The flow rate in earthworm channels can be very high compared to the rate in the soil matrix. Even for relatively small earthworm channels, the flow rate in macropores seems to be always higher than the rainfall intensity (Bouma et al., 1982, Wang et al., 1994, Shipitalo and Gibbs, 2000, Weiler, 2001). Thus, the flow capacity of the macropore system is usually not the limiting factor during the infiltration process.

Water transfer from macropores into the surrounding soil matrix has been referred to as lateral infiltration from the macropores (Beven and Clarke, 1986). In this study, the term interaction is used to describe this critical process of water flow in structured soils (LogsdonM et al., 1996, Faeh et al., 1997). Despite some experiments to measure interaction in single artificial or natural macropores in the laboratory (Smettem, 1986, Ghodrati et al., 1999) or in the field using dye tracers and soil water measurements (van Stiphout et al., 1987), a consistent description, parameterization and verification of the interaction process remains elusive.

Macropore flow initiation is the process of water supply to macropores. It is a function of initial water content, rainfall intensity, rainfall amount, hydraulic conductivity of the soil matrix, and surface contributing area (Léonard et al., 1999, Trojan and Linden, 1992). Water can flow into a macropore if its water entry pressure is exceeded, either at the soil surface or the nearly saturated subsurface (Ela et al., 1992, Li and Ghodrati, 1997, Weiler and Naef, 2003a). Weiler and Naef (2003b) found a distinct distribution of flow rates in macropores depending on where macropore flow is initiated.

The objective of this study was to develop, evaluate and test a model that combines well-known approaches of simulating initiation and interaction with the yet uncommon initiation-based flow rate distribution. This model was then used as a tool to evaluate the hypothesis that the initiation process is crucial to infiltration in macroporous soils and hence to runoff generations in general. This hypothesis was tested theoretically with a sensitivity analysis evaluating for a synthetic dataset the influence of the initiation process on interaction and on the hydrological response of the soil. The new conceptual framework of the model was tested practically by simulating results from sprinkling experiments at three field sites in Switzerland (Weiler and Naef, 2003a).

Section snippets

Model theory

The INfiltration-INitiation-INteraction Model (IN3M) allows for water flow within discrete macropores, as well as for water flow in the soil matrix based on the dual permeability concept. Analytical solutions of the Green-Ampt equation and a simple accounting scheme serve as the basis for the model. The use of analytical solutions removes any concern regarding the numerical instability observed in most preferential flow codes (Simunek et al., 2003) and results in code that can be more readily

Strategy

Since IN3M can simulate the flow rate distribution in macropores, the influence of initiation on interaction and thus on the vertical flow in the macropores and on the water content change in the soil matrix can be analyzed under the model assumptions. For this purpose the simulation of vertical water flow in the soil matrix was switched off and simulations were run for three different conditions of water flow distribution in the macropores at the upper boundary:

  • (a)

    same flow rate in each macropore

Study sites

At three field sites in Northern Switzerland (Heitersberg, Koblenz, and Niederweningen), we conducted sprinkling experiments with a Brilliant Blue solution at two irrigation rates, 60 mm h−1 (high) and 12 mm h−1 (low) with a total irrigation of ∼75 mm, and two different initial soil moisture conditions referred to as ‘dry’ (no input within 4 weeks) and ‘wet’ (∼75 mm rainfall at the previous day). The dye solution was applied using a sprinkling device consisting of a linearly moving spray bar with 15

Discussion

None of the models that describes flow and/or transport in soils containing macropores, like MACRO (Jarvis et al., 1991), QSOIL (Faeh et al., 1997), RZWQM (Ahuja et al., 1995), tipping bucket model (Emerman, 1995), SAMP (Ewen, 1996), LASOMS (Chen and Wagenet, 1992) and PREFLO (Workman and Skaggs, 1990), have yet implemented the variation of flow rate in macropores due to the initiation process. This study assessed the hypothesis that the flow rate distribution in macropores due to the

Conclusion

Evaluating and testing a model (IN3M) that combines well-known approaches of simulating initiation and interaction with the yet uncommon initiation-based flow rate distribution demonstrated that the initiation process is crucial to correctly simulate infiltration in macroporous soils. Without considering the initiation process, total water flux from the macropores into the soil matrix is overestimated and thus hydrological response is underestimated. Further work is necessary to explore the

Acknowledgements

This work was funded by the Swiss Federal Institute of Technology in Zürich within the project ‘Investigation of the water exchange mechanisms between preferential flow paths and the soil matrix’. I would like to thank Felix Naef for his support during this study and especially Thomy Keller for his help during the field experiments. Kerstin Stahl, Kevin McGuire, and Kellie Vache provided helpful comments on the manuscript.

References (64)

  • M. Weiler et al.

    Inferring flow types from dye patterns in macroporous soils

    Geoderma

    (2004)
  • M. Weiler et al.

    Virtual experiments: a new approach for improving process conceptualization in hillslope hydrology

    J. Hydrol.

    (2004)
  • M. Weiler et al.

    Simulating surface and subsurface initiation of macropore flow

    J. Hydrol.

    (2003)
  • L.R. Ahuja et al.

    Macropore transport of a surface applied bromide tracer: model evaluation and refinement

    Soil Sci. Soc. Am. J.

    (1995)
  • J. Beckers et al.

    A model of rapid preferential hillslope runoff contributions to peakflow generation in a temperate rainforest watershed

    Water Resour. Res.

    (2004)
  • K.J. Beven

    Rainfall-Runoff Modelling: The Primer

    (2001)
  • K. Beven et al.

    Macropores and water flow in soils

    Water Resour. Res.

    (1982)
  • K.J. Beven et al.

    On the variation of infiltration into a homogeneous soil matrix containing a population of macropores

    Water Resour. Res.

    (1986)
  • J. Bouma et al.

    Water infiltration and redistribution in a silt loam subsoil with vertical worm channels

    Soil Sci. Soc. Am. J.

    (1982)
  • A. Bronstert

    Capabilities and limitations of detailed hillslope hydrological modelling

    Hydrologic. Process

    (1999)
  • G.J. Burch et al.

    Soil hydrophobic effects on infiltration and catchment runoff

    Hydrologic. Process

    (1989)
  • W.M. Edwards et al.

    Rainfall intensity affects transport of water and chemicals through macropores in no-till soil

    Soil Sci. Soc. Am. J.

    (1992)
  • S.D. Ela et al.

    Macropore and surface seal interactions affecting water infiltration into soils

    Soil Sci. Soc. Am. J.

    (1992)
  • A.O. Faeh et al.

    A combined field and numerical approach to investigate flow processes in natural macroporous soils under extreme precipitation

    Hydrol. Earth Syst. Sci.

    (1997)
  • M. Flury et al.

    Dyes as tracers for vadose zone hydrology

    Rev. Geophysics

    (2003)
  • H.H. Gerke et al.

    Estimating hydraulic properties of soil aggregate skins from sorptivity and water retention

    Soil Sci. Soc. Am. J.

    (2002)
  • P. Germann et al.

    Water flow in soil macropores III. A statistical approach

    J. Soil Sci.

    (1981)
  • M. Ghodrati et al.

    Characterization of macropore flow mechanisms in soil by means of a split macropore column

    Soil Sci. Soc. Am. J.

    (1999)
  • W.H. Green et al.

    Studies of soil physics. I. The flow of air and water through soils

    J. Agric. Sci.

    (1911)
  • N.J. Jarvis et al.

    Modeling water and solute transport in macroporous soil. I. Model description and sensitivity analysis

    J. Soil Sci.

    (1991)
  • T. Kiang

    Random fragmentation in two and three dimensions

    Zeitschrift für Astrophysik

    (1966)
  • Cited by (88)

    View all citing articles on Scopus
    View full text