The dominant role of structure for solute transport in soil: Experimental evidence and modelling of structure and transport in a field experiment

Abstract. A classical transport experiment was performed in a field plot of 2.5 m2 using the dye tracer brilliant blue. The measured tracer distribution demonstrates the dominant role of the heterogeneous soil structure for solute transport. As with many other published experiments, this evidences the need of considering the macroscopic structure of soil to predict flow and transport. We combine three different approaches to represent the relevant structure of the specific situation of our experiment: i) direct measurement, ii) statistical description of heterogeneities and iii) a conceptual model of structure formation. The structure of soil layers was directly obtained from serial sections in the field. The sub-scale heterogeneity within the soil horizons was modelled through correlated random fields with estimated correlation lengths and anisotropy. Earthworm burrows played a dominant role at the transition between the upper soil horizon and the subsoil. A model based on percolation theory is introduced that mimics the geometry of earthworm burrow systems. The hydraulic material properties of the different structural units were obtained by direct measurements where available and by a best estimate otherwise. From the hydraulic structure, the 3-dimensional velocity field of water was calculated by solving Richards' Equation and solute transport was simulated. The simulated tracer distribution compares reasonably well with the experimental data. We conclude that a rough representation of the structure and a rough representation of the hydraulic properties might be sufficient to predict flow and transport, but both elements are definitely required.

[1]  Y. Capowiez Differences in burrowing behaviour and spatial interaction between the two earthworm species Aporrectodea nocturna and Allolobophora chlorotica , 2000, Biology and Fertility of Soils.

[2]  S. O. Mohamed,et al.  Morphology and origin of secondary calcite in soils from Beauce, France , 1993 .

[3]  Y. Mualem A New Model for Predicting the Hydraulic Conductivity , 1976 .

[4]  P. Raviart,et al.  A mixed finite element method for 2-nd order elliptic problems , 1977 .

[5]  Hannes Flühler,et al.  Tracer Characteristics of Brilliant Blue FCF , 1995 .

[6]  Y. Capowiez,et al.  A new simulation for modelling the topology of earthworm burrow systems and their effects on macropore flow in experimental soils , 2002, Biology and Fertility of Soils.

[7]  Lynn W. Gelhar,et al.  Stochastic subsurface hydrology from theory to applications , 1986 .

[8]  E. E. Miller,et al.  Physical Theory for Capillary Flow Phenomena , 1956 .

[9]  K. Roth Steady State Flow in an Unsaturated, Two-Dimensional, Macroscopically Homogeneous, Miller-Similar Medium , 1995 .

[10]  Graham Aylmore Tomography of Soil-Water-Root Processes , 1994 .

[11]  K. Beven,et al.  Macropores and water flow in soils , 1982 .

[12]  R. Andricevic,et al.  Geoelectric characterization of the hydraulic conductivity field and its spatial structure at variable scales , 2000 .

[13]  N. J. Jarvis,et al.  The MACRO Model (version 3.1). Technical description and sample simulations. , 1994 .

[14]  Hannes Flühler,et al.  SUSCEPTIBILITY OF SOILS TO PREFERENTIAL FLOW OF WATER : A FIELD STUDY , 1994 .

[15]  V. Hallaire,et al.  Characterization of the burrow system of the earthworms Lumbricus terrestris and Aporrectodea giardi using X-ray computed tomography and image analysis , 1999, Biology and Fertility of Soils.

[16]  S. Redner,et al.  Introduction To Percolation Theory , 2018 .

[17]  S. Finsterle,et al.  Estimating flow parameter distributions using ground-penetrating radar and hydrological measurements , 2004 .

[18]  Mark L. Rivers,et al.  Using X-ray computed tomography in hydrology: systems, resolutions, and limitations , 2002 .

[19]  Muhammad Sahimi,et al.  Dynamics of Fluids in Hierarchical Porous Media. , 1992 .

[20]  T. K. Haraldsen,et al.  Earthworm channels in cultivated clayey and loamy Norwegian soils , 1997 .

[21]  R. J. Millington,et al.  Transport in Porous Media , 1960 .

[22]  Hans-Jörg Vogel,et al.  Validity limits for the van Genuchten–Mualem model and implications for parameter estimation and numerical simulation , 2006 .

[23]  Hannes Flühler,et al.  Quantifying dye tracers in soil profiles by image processing , 2000 .

[24]  M. V. Genuchten,et al.  Mass transfer studies in sorbing porous media. I. Analytical solutions , 1976 .

[25]  B. A. Stewart,et al.  Advances in Soil Science , 1986, Advances in Soil Science.

[26]  M. Th. van Genuchten,et al.  Effect of the shape of the soil hydraulic functions near saturation on variably-saturated flow predictions , 2000 .

[27]  Hans-Jörg Vogel,et al.  Moving through scales of flow and transport in soil , 2003 .

[28]  L. Gelhar Stochastic Subsurface Hydrology , 1992 .

[29]  M. Flury,et al.  Sorption of Brilliant Blue FCF in soils as affected by pH and ionic strength , 2000 .

[30]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[31]  Jean E. Roberts,et al.  Mixed and hybrid finite element methods , 1987 .

[32]  Hans-Jörg Vogel,et al.  Effect of non‐linear adsorption on the transport behaviour of Brilliant Blue in a field soil , 2002 .

[33]  Allan L. Gutjahr,et al.  Cross‐correlated random field generation with the direct Fourier Transform Method , 1993 .

[34]  Johan Bouma,et al.  Soil morphology and preferential flow along macropores , 1981 .

[35]  X‐Ray Tomography of Soil Properties , 1994 .

[36]  M. V. Genuchten,et al.  Mass Transfer Studies in Sorbing Porous Media: II. Experimental Evaluation with Tritium (3H2O)1 , 1977 .

[37]  M. Vanclooster,et al.  Overview of inert tracer experiments in key belgian soil types: Relation between transport and soil morphological and hydraulic properties , 2001 .

[38]  R. LeVeque Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[39]  S. Meyer-Windel,et al.  Adsorption of brilliant blue FCF by soils , 1999 .

[40]  E. Brent The Measured Mobile-Water Content of an Unsaturated Soil as a Function of Hydraulic Regime , 1995 .

[41]  R. Kasteel,et al.  An Improved Dual‐Permeability Model of Water Flow and Solute Transport in the Vadose Zone , 2005 .

[42]  M. V. Genuchten,et al.  A dual-porosity model for simulating the preferential movement of water and solutes in structured porous media , 1993 .

[43]  T. Ligthart,et al.  A method for the three-dimensional mapping of earthworm burrow systems , 1993 .

[44]  Remo Guidieri Res , 1995, RES: Anthropology and Aesthetics.

[45]  N. Jarvis,et al.  The use of tension infiltrometers to assess routes and rates of infiltration in a clay soil , 1987 .

[46]  G. Dagan Statistical Theory of Groundwater Flow and Transport: Pore to Laboratory, Laboratory to Formation, and Formation to Regional Scale , 1986 .

[47]  K. Abbaspour,et al.  Single‐ and dual‐porosity modelling of multiple tracer transport through soil columns: effects of initial moisture and mode of application , 2001 .