Effective conductivity of inertial flows through porous media.

We study two-dimensional incompressible inertial flows through porous media. At core (small) scale, we prove that the constitutive, nonlinear model can be rewritten into a linear one by means of a new parameter K^{★} which encompasses all the inertial effects. In natural (large-scale) formations, K^{★} is erratically changing, and we analytically compute its counterpart, which is coined generalized effective conductivity, by the self-consistent approach (SCA). In spite of its approximate nature, the SCA leads to simple results that are in good agreement with Monte Carlo simulations.

Effects of Hydraulic Soil Properties on Vegetation Pattern Formation in Sloping Landscapes.

Current models of vegetation pattern formation rely on a system of weakly nonlinear reaction-diffusion equations that are coupled by their source terms. While these equations, which are used to describe a spatiotemporal planar evolution of biomass and soil water, qualitatively capture the emergence of various types of vegetation patterns in arid environments, they are phenomenological and have a limited predictive power. We ameliorate these limitations by deriving the vertically averaged Richards' equation to describe flow (as opposed to "diffusion") of water in partially saturated soils. This establishes conditions under which this nonlinear equation reduces to its weakly nonlinear reaction-diffusion counterpart used in the previous models, thus relating their unphysical parameters (e.g., diffusion coefficient) to the measurable soil properties (e.g., hydraulic conductivity) used to parameterize the Richards equation. Our model is valid for both flat and sloping landscapes and can handle arbitrary topography and boundary conditions. The result is a model that relates the environmental conditions (e.g., precipitation rate, runoff and soil properties) to formation of multiple patterns observed in nature (such as stripes, labyrinth and spots).

A boundary-layer solution for flow at the soil-root interface.

Transpiration, a process by which plants extract water from soil and transmit it to the atmosphere, is a vital (yet least quantified) component of the hydrological cycle. We propose a root-scale model of water uptake, which is based on first principles, i.e. employs the generally accepted Richards equation to describe water flow in partially saturated porous media (both in a root and the ambient soil) and makes no assumptions about the kinematic structure of flow in a root-soil continuum. Using the Gardner (exponential) constitutive relation to represent the relative hydraulic conductivities in the Richards equations and treating the root as a cylinder, we use a matched asymptotic expansion technique to derive approximate solutions for transpiration rate and the size of a plant capture zone. These solutions are valid for roots whose size is larger than the macroscopic capillary length of a host soil. For given hydraulic properties, the perturbation parameter used in our analysis relates a root's size to the macroscopic capillary length of the ambient soil. This parameter determines the width of a boundary layer surrounding the soil-root interface, within which flow is strictly horizontal (perpendicular to the root). Our analysis provides a theoretical justification for the standard root-scale cylindrical flow model of plant transpiration that imposes a number of kinematic constraints on water flow in a root-soil continuum.

Macrodispersion by diverging radial flows in randomly heterogeneous porous media.

Radial flow takes place in a heterogeneous porous formation where the transmissivity T is modelled as a stationary random space function (RSF). The steady flow is driven by a given rate, and the mean velocity is radial. A pulse-like of a tracer is injected in the porous formation, and the thin plume spreads due to the fluctuations of the velocity which results a RSF as well. Transport is characterized by the mean front, and by the second spatial moment of the plume. We are primarily interested in tracer macrodispersion modelling. With the neglect of pore-scale dispersion, macrodispersion coefficients are computed at the second order of approximation, without neglecting the head-gradient fluctuations. Although transport is non-ergodic at the source, it is shown that ergodicity is achieved at small distances from the source. This is due to the fact that close to the source local velocities are quite large, and therefore solute particles become uncorrelated very soon. Under ergodic conditions, we compare macrodispersion mechanism in radial flows with that occurring in mean uniform flows. At short distances the spreading effect is highly enhanced by the large variability of the flow field, whereas at large distances transport exhibits a lesser dispersion due to the reduction of velocities. This supports the explanation provided by Indelman and Dagan (1999) to justify why the macrodispersivity is found smaller than that pertaining to mean uniform flows. The model is tested against a tracer transport experiment (Fernàndez-Garcia et al., 2004) by comparing the theoretical and experimental breakthrough curves. The accordance with real data, that is achieved without any fitting to concentration values, strengthens the capability of the proposed model to grasp the main features of such an experiment, the approximations as well as experimental uncertainties notwithstanding.

Darcian preferential water flow and solute transport through bimodal porous systems: experiments and modelling.

Soils often exhibit a variety of small-scale heterogeneities such as inter-aggregate pores and voids which partition flow into separate regions. In this paper a methodological approach is discussed for characterizing the hydrological behaviour of a heterogeneous clayey-sandy soil in the presence of structural inter-aggregate pores. For the clay soil examined, it was demonstrated that, coupling the transfer function approach for analyzing BTCs and water retention data obtained with different methods from laboratory studies captures the bimodal geometry of the porous system along with the related existence of fast and slow flow paths. To be effectively and reliably applied this approach requires that the predominant effects of the soil hydrological behaviour near saturation be supported by accurate experimental data of both breakthrough curves (BTCs) and hydraulic functions for high water content values. This would allow the separation of flow phases and hence accurate identification of the processes and related parameters.

Analytical solutions for reactive transport under an infiltration-redistribution cycle.

Transport of reactive solute in unsaturated soils under an infiltration-redistribution cycle is investigated. The study is based on the model of vertical flow and transport in the unsaturated zone proposed by Indelman et al. [J. Contam. Hydrol. 32 (1998) 77], and generalizes it by accounting for linear nonequilibrium kinetics. An exact analytical solution is derived for an irreversible desorption reaction. The transport of solute obeying linear kinetics is modeled by assuming equilibrium during the redistribution stage. The model which accounts for nonequilibrium during the infiltration and assumes equilibrium at the redistribution stage is termed partial equilibrium infiltration-redistribution model (PEIRM). It allows to derive approximate closed form solutions for transport in one-dimensional homogeneous soils. These solutions are further applied to computing the field-scale concentration by adopting the Dagan and Bresler [Soil Sci. Soc. Am. J. 43 (1979) 461] column model. The effect of soil heterogeneity on the solute spread is investigated by modeling the hydraulic saturated conductivity as a random function of horizontal coordinates. The quality of the PEIRM is illustrated by calculating the critical values of the Damköhler number which provide the achievable accuracy in estimating the solute mass in the mobile phase. The distinguishing feature of transport during the infiltration-redistribution cycle as compared to that of infiltration only is the finite depth of solute penetration. For irreversible desorption, the maximum solute penetration W/theta(r) is determined by the amount of applied water W and the residual water content theta(r). For sorption-desorption kinetics, the maximum depth of penetration z(r)(e, infinity ) also depends on the ratio between the rate of application and the column-saturated conductivity. It is shown that z(r)(e, infinity ) is bounded between the depths W/(theta(r)+K(d)) and W/theta(r) corresponding to the maximum solute penetration for equilibrium transport and for irreversible desorption, respectively. This feature of solute penetration explains the unusual phenomena of plume contraction after an initial period of spreading [Lessoff, S.C., Indelman, P., Dagan, G., 2002. Solute transport in infiltration-redistribution cycles in heterogeneous soils. In Raats, P.A.C., Smiles, D.,Warrick, A.W. (Eds), Environmental Mechanics: Water, Mass and Energy Transport in the Biosphere. American Geophysical Union, pp. 133-144]. Unlike transport under equilibrium conditions, when the solute is completely concentrated at the front, the solute under nonequilibrium conditions is spread out behind the front. Heterogeneity leads to additional spreading of the plume.