Next Article in Journal
Building System Capacity with a Modeling-Based Inquiry Program for Elementary Students: A Case Study
Previous Article in Journal
Application of Systems-Approach in Modelling Complex City-Scale Transdisciplinary Knowledge Co-Production Process and Learning Patterns for Climate Resilience
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrostructural Pedology, Culmination of the Systemic Approach of the Natural Environment

1
Institut de Recherche pour le Développement, (IRD), F13572 Marseille, France
2
Valorhiz—Parc Scientifique Agropolis, F34980 Montferrier sur Lez, France
3
Faculty of Agricultural and Food Sciences, American University of Beirut, Beirut 11-0236, Lebanon
4
Department Biological and Agricultural Engineering (BAEN), Zachry Department Civil and Environmental Engineering (CVEN), Texas A&M University, College Station, TX 77843, USA
*
Author to whom correspondence should be addressed.
Submission received: 1 December 2020 / Revised: 13 January 2021 / Accepted: 19 January 2021 / Published: 22 January 2021

Abstract

:
The subject of this article is the dynamics of water in a soil pedostructure sample whose internal environment is subjected to a potential gradient created by the departure of water through surface evaporation. This work refers entirely to the results and conclusions of a fundamental theoretical study focused on the molecular thermodynamic equilibrium of the two aqueous phases of the soil pedostructure. The new concepts and descriptive variables of the hydro-thermodynamic equilibrium state of the soil medium, which have been established at the molecular level of the fluid phases of the pedostructure (water and air) in a previous article, are recalled here in the systemic paradigm of hydrostructural pedology. They allow access to the molecular description of water migration in the soil and go beyond the classical mono-scale description of soil water dynamics. We obtain a hydro-thermodynamic description of the soil′s pedostructure at different hydro-functional scale levels including those relating to the water molecule and its atoms. The experimental results show a perfect agreement with the theory, at the same time validating the systemic approach that was the framework.

1. Introduction

The problem of water transfer equations in soil dates back to the beginning of soil science. The best-known equation and the basis of all models of water circulation in the soil, is the “Richards” equation, which associates Euler′s law of continuity and Darcy′s law extended to unsaturated soils:
d θ d t = d f d z   and   f = k d h d z
where θ is the volumetric water content (unit less), t is time, f the flow, z the vertical coordinate, k the water conductivity and h the soil water retention pressure.
We resume here the study of the water transfer equation in the soil with a completely new approach: that of the systemic approach we recently theorized [1,2,3] from the work of Bertalanffy, initiator of the general theory of systems [4] and Le Moigne [5], author of the General System model. The application of this systemic approach applied to pedology has created a new paradigm of characterization, water modeling, and representation of the natural environment (multi-scale mapping). It is named hydrostructural pedology [1,2,3] and is presented schematically in Figure 1.
Hydrostructural pedology allowed the development of a new physics of soil water, qualified as systemic, based on the recognition of the pedostructure and the two types of water associated with it [6,7,8,9]. We demonstrated that these two types of water in the pedostructure are two aqueous phases in pressure equilibrium ( h m i W m i = h m a W ) and distinguished by their thermodynamic properties.
Figure 1. Place of hydrostructural pedology among the earth sciences [10].
Figure 1. Place of hydrostructural pedology among the earth sciences [10].
Systems 09 00008 g001
This new physics of water in the soil has led to reconsideration the fundamental bases of classical thermodynamics; in particular, those of Gibbs free energy [11]. This study, associated with the new concept of “system organized into organized subsystems that are molecules and their atoms”, has made it possible to develop a new vision of the thermodynamic equilibrium of the soil medium. The liquid and gas phases are all recognized as subsystems organized in molecules and themselves in atoms, relative to the solid phase that makes up the soil structure. Classical variables such as temperature, entropy, pressure, internal energy, and Gibbs free energy, can then be physically redefined and precisely explained accounting the two levels of organization: molecular and atomic.
In the present article, we will introduce these two levels of organization into the current description of the pedostructure organization for re-examining the terms of the Richards equation (Appendix A) after having rewritten it in the systemic framework of the hydrostructural pedology [1,2,3]. The new formulation will allow for taking into account all functional volumes as thermodynamics variables (volume of the whole system and volumes of hydro-functional subsystems). We will see then that the challenge becomes how to associate in the same equation an extensive variable (volume) and an intensive variable (potential). This is currently done empirically using Darcy′s law extended to unsaturated soils. In this article, we will see that the organization level, in which the extensive and intensive variables meet, is the molecular level of organization. At this level, the terms internal or external energy, internal or external pressure and internal or external chemical potential are all in equilibrium relationship, as explained in [11].

2. Updated Theoretical Background of Hydrostructural Pedology

2.1. The Pedostructure, Test Body of Hydrostructural Pedology

2.1.1. Preparation of a Standard Sample of Pedostructure

The pedostructure is the fundamental concept at the basis of hydrostructural pedology. Materially, it constitutes the first two levels of organization of the soil horizon: that of the clay plasma and that of the assembly of primary aggregates between them and possibly with other mineralogical or biological grains of sand size. Pedostructure is present in all soil horizons (Figure 2); its volume percentage in the horizon and its specific hydrostructural properties, due to the clay plasma that makes up the primary aggregates, characterize the hydrostructural behavior of a soil horizon.
The variables that are used for the systemic description of the pedostructure of any soil horizon are listed in Table 1. They all are reported to the same mass of solids: the pedostructural mass of the considered total soil volume.
It is, therefore, necessary to define the representative sample of pedostructure in the laboratory: that sample upon which all the measurements of the hydro-functional curves of the soil will relate as well as the determination of their parameters. These characteristic curves are the shrinkage curve, V W , the water retention curve, h W , the unsaturated soil water conductivity, k(W) and the swelling curve of primary aggregates as a function of time, W m i t [2,9].
In our study, a standard laboratory pedostructure sample is a soil sample that is reconstituted with what is traditionally called “fine earth”, the 2 mm sieved soil from the fractionation of a moderately dry soil sample (<pF3) on a 2 mm sieve (can be 4 mm when the sample is very clayey with swelling clay). The fine earth is added layer by layer in a cylinder of 5 cm diameter and 5 cm height, placed on a damp terry cloth; each layer added wets along with the filling. The objective is to obtain a homogeneous sample in terms of structure and hydrostructural behavior. The soil cylinder is then subjected to 2 cycles of desiccation-humidification, the desiccation being carried out using either the Richards press at pF3 (15 bar) or evaporation in ambient air, the sample being positioned, in this case, so that the evaporation occurs simultaneously on both sides of the cylinder.
These preparation standards for the pedostructure sample are at the same level of importance as the oven temperature standard of 105° for the definition of dry soil. The term “pedostructural mass” is the mass of the solid phase that constitutes the pedostructure of the sample: it constitutes the universal benchmark for the extensive variables of a soil horizon (water content, salts etc. referred to the pedostructural mass).

2.1.2. Characterization and Modelling of the Hydrostructural Properties of the Soil

The parameters of the two equilibrium equations of the hydrostructural state of the pedostructure, the shrinkage curve V W and the soil water retention curve h W , are determined from the curves measured on the standard sample using the TypoSoil® device which can simultaneously process up to 8 samples [12]. This characterization is totally accepted by the soil water model Kamel [9,13] and was fully established within the systemic paradigm of hydrostructural pedology with constant reference to the notions of nested organizations, hydro-functional levels of organizations (primary aggregate, pedostructure, soil horizon, pedon, primary soil unit, etc.), and using only variables, functions and parameters, said to be systemic because they are defined in the systemic description network of the hydrostructural pedology [1,2,3]. All the extensive variables of the studied homogeneous organized system, in particular the cut volume of the sample taken, are related to the fixed mass of the solid phase comprising the structure cut out in this volume.
However, the exact thermodynamic formulation of the water retention curve h W at equilibrium of the two aqueous phases of the pedostructure [8] and, from this, the exact distribution of the two kinds of water content ( W m a and   W m i ) in the pedostructure as function of W , raises a new and important question about the descriptive variables of the model. This equation h W links an extensive variable (water content) to an intensive variable (water suction). Indeed, we use a mini tensiometer (2 mm thick) planted laterally in a soil sample at depths z, to simultaneously obtain water suction ( h = h m i = h m a ), micro and macro water contents ( W m i and W m a ) at thermodynamic equilibrium, locally in the same soil mini-layer. To answer the fundamental question of what exactly the spatial extension of is W m i and W m a corresponding to h m i and h m a , we must recall the principal results of the previous article [11] on the hydro-thermodynamic equilibrium of the pedostructure at the molecular and atomic levels.

2.2. Molecular Thermodynamic Equilibrium of the Fluid Phases of the Pedostructure

2.2.1. Internal Molecular Organization of the Aqueous Phases at Equilibrium

The water molecule and simple gas molecules constitute the material point of the fluid phases of the natural environment. These molecules have a specific energy volume V m , the sum of the volumes occupied by their constituent atoms, and a mass, the sum of the masses of their atoms. The free energy that the atom develops in the parent molecule is of the oscillatory type: 1 2 m λ 2 ν 2 . The sum of the atomic free energies of the molecule constitutes the oscillatory energy of the molecule contained in its volume V m   =   V a i . It is this internal “free” energy of the molecule that has been identified as the temperature of a molecule. Thus, because of the mathematical property of fractions: a b = c d = a + c b + d , not only is the internal pressure equal everywhere in the molecule ( P = T V = T a i V a i ), but also the chemical potential defined by μ m = T m = T a i m a i = T a i m a i .
In the aqueous or gaseous fluid phases in which the molecules, while being optionally ionized, maintain their chemical compositions in atoms, the internal oscillatory energies T = 1 2 m λ ν 2 of the molecules are balanced with the kinetic energies acquired and maintained by what can be called thermal agitation ( E m v = T S 1 = 1 2 m v 2 ): the “shocks” or meetings of molecules of the fluid phase between them. We can then associate each molecule with an occupancy volume: V m t , that contains the two types of energy: oscillatory and linear kinetics (Figure 3).
At equilibrium, the volume concentrations of the internal and external energies (pressures) of the molecules are the same and this is where the molecular entropy of the phase S a comes in. S a is a fractional number equal to V m t V m , which makes thermodynamic equilibrium possible where all molecules have the same internal and external pressure.
Moreover, given that each molecule necessarily has the same chemical potential μ m i = T i / m i as that of its atoms and the same mass volume, equilibrium is achieved if the chemical potential of the molecules is the same in all the phases. Molecules of chemical masses have different temperatures in the fluid phase in equilibrium, but have the same molecular chemical potential μ m = T i / m i and a chemical potential μ m α = μ m S α that depends on the entropy S α of the α phase. We can define the chemical potential of the phase that concerns only the external kinetic energy of the molecules, which we will call the intermolecular chemical potential, μ v α , (index v void):
μ v α = μ m S α 1
The big difference with statistical thermodynamics is in the definition of temperature and the understanding of the free energy of the thermodynamic system (homogeneous liquid and gaseous fluid phases in equilibrium). The temperature of the phase is not the statistical average of a variable temperature around a value but, rather, the exact average of the temperatures of a finite number of molecules of different chemical species, the molecules of the same species having the same temperature (internal energy).
A phase is characterized by its entropy S α , an intensive state variable having the same value throughout the phase. At phase equilibrium, since the internal chemical potential of molecules, μ m , is the same everywhere, the overall molecular chemical potential of the phase, μ α = μ m S α , is, therefore, an intensive variable characteristic of the phase. It is the same for the external chemical potential of molecules of the phase ( μ v α ), more specifically the characteristic kinetic potential of the phase ( 1 2 v 2 ). This means that all molecules, regardless of their mass, have the same linear speed in the phase.
We, therefore, have a fundamental relationship between the chemical potential of the molecule, the entropy of the phase and the speed of the molecules in the phase:
μ v α = μ m S α 1 = μ m α μ m = v α 2 2
where v α is the linear speed of molecules of the phase at the thermodynamic equilibrium state. The pressure of the water molecules in this intermolecular space is T w S α 1 / V m t w V m w , which is equal to the internal molecular pressure T w / V m w and the total molecular pressure in the phase: T w S α / V m t w .
However, the suction pressure measured by the tensiometer in soil science, as shown previously [11], has the expression:
h = ρ w μ w μ w °
where μ w ° is the chemical potential of free water under air pressure and standard temperature. According to the Equation (3), and because μ m is equal everywhere in all phases of the system at thermodynamic equilibrium, we can substitute the chemical potential of the water ( μ w ) by the intermolecular chemical potential ( μ v w ) without changing the value of the pressure h .
h = ρ w μ v w μ v w °
This allows μ v α to be identified with the pressure potential of the water in the tensiometer, relating h to the speed squared of molecules in the phase (Equations (3) and (5)).
The particularity of this speed is being the same for all molecules of the phase, whatever their mass, at thermodynamic equilibrium state since the criterion of intra and inter phase thermodynamic equilibrium is the molecular chemical potential ( μ m ) and not the temperature ( T i ) as was shown in [11]. The temperature, identified with molecular energy, is in fact different for each chemical species of the phase since it is obtained from the molecular chemical potential μ m of the fluid phases of the system in equilibrium ( T i = m i μ m ). The temperature of water molecules, for example, is equal to:
T w = m w μ m = m w μ m α S α = m w μ v α S α 1
With these state variables of the thermodynamic equilibrium of the fluid phases (liquid and gaseous) of the pedostructure defined at the two organization levels, the molecule and the phase, we are able to search for the existing relationships between these newly defined variables and the usual ones (flow, water content, pressures, suction, etc.). It should be remembered that this link is only possible between variables defined according to the systemic approach, whether intensive or extensive.
We give in Table 2 the exhaustive list of hydro-thermodynamic variables qualified as systemic and which cover the four levels of organization: macroscopic and microscopic of the aqueous phase then molecular and atomic of the phase.

2.2.2. Identification of Constants E ¯ m a and E ¯ m i as Intermolecular Free Energies of the Pedostructure

We know that at the macroscopic level of the phases of the pedostructure, the free energies of the two aqueous phases of the pedostructure ( E ¯ m i and E ¯ m a ) are observed constant despite a change in the water content of this phase in the defined system of the pedostructure. Following Sposito [13], E ¯ m i and E ¯ m a was identified before as the free energy G ¯ m i = μ w m i W m i and G ¯ m a = μ w m a W m a , W m i and W m a being the water contents micro and macro of the pedostructure. Now, following our previous study [11], which differentiates the intermolecular energy ( T S α 1 ) from the total energy ( T S α ), E ¯ m i and E ¯ m a are defined as the intermolecular energy, corresponding to kinetic energy of molecules:
E ¯ w m a = 1 2 m ¯ w m a v m a 2 = μ v m a m ¯ w m a = n ¯ w m a m w μ v m a
E ¯ w m i = 1 2 m ¯ w m i v m i 2 = μ v m i m ¯ w m i = n ¯ w m i m w μ v m i
This relationship involves the number of water molecules n ¯ w α in the aqueous phase “ w α ” (macro or micro), whose molecular mass is m w . The horizontal line above the extensive variables signifies the ratio to the solid mass present in the same elementary volume.
Recall that these two aqueous phases coexist in the section (elementary layer) at z, one surrounding the other, in the pedostructure sample. They are indexed: W m a and W m i . The fact that the energy of the phase α (macro or micro), E ¯ w α , is constant despite a change in phase’s water content m ¯ w α or in its chemical potential μ w α , appears as a displacement constraint for the molecules since a potential gradient is created as soon as a deficit of water molecules appears in the system. This constraint is written:
d E ¯ w α = 0 = m ¯ w α d μ v α + μ v α d m ¯ w α = 0
or, for each α aqueous phase (macro or micro):
d μ v α μ v α = d m ¯ w α m ¯ w α   similar   to   d log μ v α = d log m ¯ w α
As the soil medium of the pedostructure is in thermodynamic equilibrium, we have, at every point of the medium, equality of the retention pressures between the two phases:
h z = h m a = h m i = ρ w μ v m a μ v m a ° z = ρ w μ v m i μ v m i ° z
This gives the following general equation, since μ v m a μ v m a ° = μ v m i μ v m i ° and according to (10):
μ v m a d log m ¯ w m a = μ v m i d log m ¯ w m i = d μ v m a = d μ v m i
Note that this equation is valid only if the saturated state corresponding to h = 0 is set and therefore the potential μ v m a ° and μ v m i ° corresponds to W m a ° = W m a s a t and W m i ° = W m i s a t and that all these equilibrium equations are deduced from the fact that the “free” energies E ¯ w m a and E ¯ w m i of the two aqueous phases in equilibrium, defined and expressed by the fundamental expressions (7) and (8), are constant and characteristic of the pedostructure.

2.2.3. Definition of the Molecular Flux f m a in the Pedostructure in Thermodynamic Equilibrium

Consider the two aqueous phases of the pedostructure, one in the inter-aggregate space (macro phase), surrounding the other (micro phase) in the clay plasma of the primary peds. Only the macro phase is in contact with air and in capillary continuity throughout the sample. It has the possibility of moving according to a potential gradient created by the departure of water molecules from the surface. The micro phase has no capillary continuity and is in contact only with the macro phase with which it locally comes into pressure equilibrium by exchanging water molecules.
In the absence of a potential gradient, the sum of the “speed” vectors of all molecules is zero (zero divergence). When a potential gradient is set up, on the z axis for example, an acceleration field automatically derives from the potential and is applied to the molecules placed in this gradient; in the case of the macro phase:
γ ma z = d μ v m a d z
In the case of our standard experiment of drying a soil cylinder by evaporation of water on the surface, the force exerted on all molecules of the phase is, in accordance with Newton′s second law: F i = m i γ w m a . This force determines the elementary pressure carried by each molecule of the phase on its environment: ( p = ρ w γ w m a ; cf. [11]), which causes a molecular flux f m a through the surface s m a (occupied by the macro aqueous phase) of the sample section at depth z, such as:
γ ma z = d f m a d t = d μ v m a d z
where f m a is the molecular flux of the macro phase, identified to the molecular speed of phase molecules along a gradient line if it exists (non-zero divergence). It is, therefore, an intensive phase variable that has the same value for each molecule in the phase. This is what conceptually distinguishes this type of (molecular) flux from the flows usually considered as flow in the Euler equation, defined as the speed of passage of a volume of n molecules through a chosen surface and not the actual passage of molecules. Moreover, as said above, the argument of velocity vectors of molecules of a phase is equal throughout the phase in thermodynamic equilibrium, whatever their molecular masses.

3. Materials and Methods

3.1. Soils

3.1.1. Provenance

All tested soils in this study come from Martinique; they were collected and characterized as part of a IRD project to establish a SIG of Soils of the Martinique [14]. The goal of the project was to physically characterize the hydrostructural properties of the soils described in the notice of the very detailed existing soil map of Martinique. Soils are clayey of volcanic origin, differentiated by pedogenesis according to the geomorphological situation and geographical position they occupy around the ancient volcano. These lead, over small distances, to well-differentiated pedohydric regimes and different microclimates and plant cover on the surface.

3.1.2. Hydrostructural Characterization

Hydrostructural characterization consists of measuring the shrinkage curve and the potential curve performed at the same time on the same sample [7]. At the time of the project, TypoSoil did not exist yet and the hydrostructural characterization was accomplished by measuring the shrinkage curve V (W) and the suction curve h (W) (or water retention) of the soil on two separate samples. The porous ceramic of the tensiometer had to be placed in the center of the sample, and the water had to evaporate uniformly over the entire surface so that the curve was representative of the sample.

3.2. Measuring Apparatus of the Hydric Conductivity of the Pedostructure

The apparatus used was manufactured to measure the water conductivity of the soil pedostructure (Figure 4). It is composed of a balance on which rests a metal cylinder of 5 cm in diameter and 5 cm in height containing the soil sample, collected at a moisture state close to the field capacity. The cylinder is provided with two holes that allow the introduction of two mini-tensiometers T1 and T2 (diameter 2 mm), positioned 1 cm and 2 cm, respectively, from the surface. The sample, reshuffled or not, is first brought to saturation with its cylinder on a sandbox. The upper surface is then leveled at the edge of the cylinder while the lower surface is covered with plastic film to prevent any evaporation on this side. Finally a flat ring, with an outer diameter equal to the diameter of sample and forming a strip a few millimeters in width, is laid on its upper face to limit the lateral evaporation that occurs after reduction in the diameter by shrinking.
The tensiometers are connected to a box of pressure sensors, in turn connected to a computer. The assembly is placed in a thermostatic chamber at 34 °C. Tensions and weights are recorded every 5 min. We deduce the values of the overall water content W t , the water potential μ m i ,   μ m a of the micro and macro phases at two positions of the tensiometer as well as the corresponding local water contents W m i ,   W m a by using the following relationships:
W t = M t M S M S ;   h = h m a = h m i = ρ w μ m a μ m a ° = ρ w μ m i μ m i °
h = ρ w E ¯ m a 1 W m a 1 W m a S a t = ρ w E ¯ m i 1 W m i 1 W m i S a t
E ¯ m a = W m a μ W m a   and   E ¯ m i = W m i μ W m i
where E ¯ m i , E ¯ m a , W m a S a t and W m i S a t are the characteristic parameters of the retention curve h (W) given by the previous hydrostructural characterization of a soil horizon.

3.3. Systemic Variables Used for Modelling the Water Movement at the Different Organization Levels of the Pedostructure

3.3.1. Hierarchical Arrangement of the Solid, Liquid and Air Phases at a Given Depth

The standard object of the study being a sample of pedostructure of cylindrical shape by which we study the variation of the water state along the z axis, all the descriptive variables used must relatable to the same level of scale: that of the horizontal section of the sample at depth z, then allowing us to relate the two microscopic and macroscopic aspects of the sample. We can imagine this section, over an elementary height dz, surfaces occupied by the well identified phases: the solid phase ( s s ), two aqueous phases ( s m i ) and ( s m a ), and the gas phase ( s a i r ) such that the total surface of the section s t is the sum s t = s s + s m i + s m a + s a i r . The entire surface is homogeneously filled with these 4 phases with the essential constraint that the arrangement of the 4 phases between them is respected: the solid phase is surrounded by micro water, which is surrounded by the macro phase, and the macro phase surrounded by the air phase.
Recall that the basic assumption of the systems approach is to consider the solid phase of the structure as invariant in mass in the discretization of space: the elementary horizontal layers of volume δ V = s t d z all contain the same quantity of structural mass.
δ m s d z = ρ s s s d z d z = ρ s s s = M s L
where ρ s s s is a fixed characteristic of the homogeneous sample in terms of structure and its structural mass M s ; ρ s is the actual density of the solid phase, M s is the total mass of the dry sample and L is its length.
The extensive variables, such as water contents and organized volumes and sub-volumes ( W ,   W m i ,   W m a ,   V ¯ , etc.), are all related to the mass of the local structural phase: δ m s = ρ s s s d z , the mass contained in the same volume δ V as that in which the other variables are defined. We have for example:
W z = ρ w s w d z ρ s s s d z = ρ w s w z ρ s s s z
In the systemic approach, this structural mass is the fixed reference to which are attached all the variables of the same volume and which, in turn, is variable with the water content. This ensures consistency in the creation and definition of descriptive variables as well as their correct use.

3.3.2. Extensive Variables such as Volumes and Water Contents

The crucial problem is indeed that defining the extensive variables of water content and of volumes that depend on the organizational scale at which the variable is considered: the macroscopic level of the entire system ( V t ,   W t ,   etc . ) or a discretized subsystem between two depths ( V z 1 z 2 = Δ V 1 2 ) or that considered at the molecular level, that of the horizontal section at depth z.
(1) The global variables W t , V t , etc. marked with the index t, are considered homogeneous and in thermodynamic equilibrium over the entire pedostructural system, the subject of the study. The volume ( V t ) divided by the total mass of the solid phase ( M s t ) constituting the structure of the system, is the mass volume, which we write as V ¯ t : V ¯ t = V t / M s t . Likewise, the overall water content of the system, W t , is the mass water content of the system M w divided by the structural mass M s and written as: W t = M w / M s t . These variables, all related to the same structural mass, are additive: W t = W m a t + W m i t . We have the following equalities:
W t = ρ w V w t M s t = ρ w V ¯ w t = M w M s t
The mass volume of the pedostructure (mass pedostructural volume) is variable with its water content in accordance with a characteristic property of the soil: its shrinkage curve: V ¯ t = f W t .
(2) A second type of variable is the local variable, defined for a delimited part of the homogeneous medium of the total system. In our case of a standard sample of pedostructure, this is a horizontal layer of the sample between the dimensions z1 and z2. This layer is a subsystem of the overall system, and has the same structural characteristics but the extensive variables of volume and water content of the various mobile phases of the system (aqueous and gaseous) are not related to the total structural mass of the system. They are related to the local structural mass of the medium, between the depths z1 and z2. This type of variable with an extensive character is indexed with Δ z
W Δ z = ρ w V w Δ z M s Δ z = ρ w V ¯ w Δ z = M w Δ z M s Δ z
(3) The third category of variables of an extensive nature is the “molecular” variable, a function of z. These variables are attached to the horizontal section of the sample identified by the corresponding z score. Let us redefine the mass volumes, like V ¯ z , using all the descriptive variables of the aqueous and gas phases defined at this depth z:
V ¯ z = lim d z 0 s t ρ s   s s d z d z = s t C M s 1 L 3
where
C = ρ s   s s   d z d z = d m s d z = ρ s   s s   at   depth   z .
Subsequently, the z index, indicating that the variable is molecular, will generally be omitted.
Assuming that the medium is homogeneous from the perspective of its structure, we can consider c as a constant that can be estimated at C = d m s d z = M s L for a cylindrical sample of height L and structural mass M s . We define the water content variables in the same way:
W m a = ρ w ρ s s m a s s = ρ w   C s m a = ρ w V ¯ m a
W m i = ρ w ρ s s m i s s d z d z = ρ w   C s m i = ρ w V ¯ m i
W = ρ w ρ s s w s s d z d z = ρ w C s w = ρ w V ¯ w
Note that, since ρ w and C are constants, any ratio of two occupied surfaces of a section in z ( s m a / s m i at depth z, for example) is equal to the ratio of the volumes based on these surfaces and for height dz. This ratio of two differentiated surfaces at the molecular level of a horizontal section of the sample, at z, can be considered equivalent to the ratio of the corresponding extensive variables at the same z, for example:
s w s t = V w V = V ¯ w V ¯ = θ w = W ρ w   V ¯
where s w is the area of the section at z occupied by water ( s w = s m a + s m i ) and s t is the total area, at z.

3.3.3. Concomitant Variation of Organizational and Fluxes Variables at z

The two water fluxes, f e w and f w , are, according to their physical definition (30) and (31), proportional to the molecular flux according to:
ρ w   V ¯ f e w = W m a f w = W f m a
or, dividing by ρ w   V ¯ :
f e w = θ w m a f w = θ w f m a
However, these Equations (32) and (33) do not provide any information on their reciprocal variations in time and space: it must be the same dz for the same dt that makes up the equations defining the three types of variable containing the equation (flux, volume and mass). It is possible to resolve the uncertainty about dz by considering the derivative of these variables with respect to z. The derivative of f m a defined by Equation (31) gives:
d f m a d z = d d z s m a s m a d z d t = d d t s m a d z s m a d z = d W m a d t W m a = d ln W m a d t
The change in sign results from the fact that one passes from a variation of speed, d z d t , to a variation of volume calculated with the height d z taken in the opposite direction of the speed. In the same way, we also have:
d f w d z = d d z s w s m a d z d t = s w s m a d d z s w s w d z d t = W W m a d ln W d t = d W W m a d t
d f e w d z = d s w d z d z   s t   d t = d ln W d t s w s t = d ln W d t W ρ w V ¯ = d W ρ w V ¯ d t
Note that the relations (34) to (36) participate in the definition of local extensive variables (at depth z); we have in fact:
The indeterminacy having been lifted, it is allowed to relate these 3 equations to each other, which gives:
d f w d z / d f m a d z = d W d t / d W m a d t = 1 + d W m i / d t d W m a / d t
d f e w d z / d f w d z = W m a ρ w V ¯ = θ w m a
and we can rewrite the continuity equation in this form:
d W d t = ρ w V ¯ d f e w d z = W m a d f w d z = W m a d f m a d z 1 + d W m i / d t d W m a / d t
Other relationships between fluxes and water contents are given in Appendix B.

3.3.4. Spatial Variation of the Product f m a W m a

We saw that f m a is a molecular flux of the aqueous phase macro, and W m a is the water content of this phase at depth z, given by Equation (24). The problem is the constant C that makes reference to the solid phase. By taking the correct expression for the solid phase, we can then consider f m a as the speed of each molecule and f m a W m a as the concentration of momentum whose derivative with respect to time is a force.
Consider the molecular expression of the product: f m a W m a , in accordance with Equations (24)–(27) and (31) of the physical definition of the two variables and their derivatives with respect to z:
  f m a W m a = d z d t ρ w s m a C = ρ w s m a d z ρ s   s s d z d z d t = ρ w δ V m a δ m s d z d t = ρ w d V m a α z m s z m a d t = ρ w d V ¯ m a α z d t = d W m a α z d t
where C is the constant of the material C = ρ s   s s defined by Equations (22) and (23), δ V m a the element of volume equal to δ V m a = s m a d z , δ m s = ρ s   s s d z , the mass of the solid phase concomitant, and m s z m a the mass of the solid phase at the level of the section at z associated with the volume variation d V m a = d s m a α z ; α z = C / m s z m a L 1 .
We give for α z the following physical definition:
ρ w s m a d z   C d t = ρ w δ V m a   δ m s d z d t = ρ w d V m a   α z m s d t = ρ w d V ¯ m a   α z d t = d W m a α z d t
where δ m s is the solid mass element corresponding to δ V m a and such that:
δ m s d z = α z m s m a
m s m a being the element (mass) of the solid phase associated with d V m a and such that:
d V m a m s m a = d V ¯ m a
Comparing Equations (39) and (40), we deduce that:
α z f m a W m a = d   f m a d z W m a = d W m a d t
thus,
α z f m a = d   f m a d z
Since f m a = d z d t , we also have:
f m a d W m a d z = d z d t d W m a d z = d W m a d t
and from (44), the general equation:
α z W m a f m a = d f m a d z W m a = d W m a d t = f m a d W m a d z
The consequence of Equation (47) is that:
d ln f m a d z = d ln W m a d z = α z   and   d ln W m a f m a d z = 0

3.3.5. Flux Variables at Depth z

Just as we have defined the extensive variables at depth z V ¯ ; W ; W m a ; W m i , we must also define the associated types of fluxes at depth z:
(1) The global or Eulerian flux ( f e w ): the flux of water crossing the entire horizontal section ( s t ), of the sample, without distinction of the surface actually crossed by this section:
f e w = m w ρ w d n w s t d t = d V w s t d t = s w s t d z d t = s t s t d l w d t = d l w d t
where m w is the molecular mass of water and n w is the number of water molecules in the elementary water volume d V w = s w d z , which is in the elementary soil volume d V t = s t d z . This defines the elementary high of water ( d l w ) such that s t d l w = d V w = s w d z . It follows that d z d t is the rate of transfer of water molecules through the surface ( s t ), while d l w / d t is the rate of drainage of the height of water in the volume element ( d V = s t d z ).
This implies the equivalence: s w d z = s t d l w and the different forms of writing of the Eulerian flux that we obtain by using Equation (28) of equivalences with the ratio of molecular surfaces:
  f e w = d l w d t = s w s t d z d t = θ w d z d t = W ρ w   V ¯ d z d t .
(2) The real flux,   f w , transfer speed of water molecules through their real surface of passage ( s m a ), the surface occupied by the macro water molecules, the aqueous phase external to the primary aggregates:
f w = m w ρ w d n w s m a d t = d V w s m a d t = s t s m a d l w d t = ρ w V ¯ W m a f e w = f e w θ m a = W W m a d z d t
where the ratio d z d t represents, as in Equation (29), the speed of movement of water molecules on the z axis.
(3) The molecular flux ( f m a ), the real speed of the water molecules of the surface s m a , is the speed (modulus) of the molecules of the macro phase determined by the chemical potential of the phase different from f w related to the number of micro and macro water molecules leaving the surface during the time (dt):
f m a = m w ρ w d n w m a s m a d t = s m a d z s m a d t = d z d t
Indeed, the speed d z d t defined by this Equation (31) is that of the water molecules of the phase, namely, the speed of agitation of the molecules of the macro phase.
Thus, the Eulerian and real fluxes of the water in the pedostructure are linked to the speed of agitation of the molecules of the macro aqueous phase through the intermediary of the molecular flux f m a and therefore directly linked to the state variables of thermodynamic equilibrium of the phase (temperature, chemical potential, entropy etc.)

3.4. Writing of the Physical Process at a z-Section Level of Scale

Having defined the descriptive variables of the organization of the internal environment of the pedostructure at depth z, we can now introduce the physical processes that govern the movement of water (fundamental equation of dynamics) and the regulation of liquid phases by relative to the solid phase (thermodynamic equilibrium) due to evaporation of surface water.

3.4.1. Application of Newton’s Law

The relation of flux with time, when it comes to a speed of movement, goes through the fundamental law of mechanics and Newton’s 2nd law, mentioned above (14). These laws apply to the molecular flux ( f m a ) which proceeds from the chemical potential gradient of the macro aqueous phase, the relationship of which is known with the water content of the phase at depth z:
The upward force F m a , which drives the water molecules, of molecular mass: m w , present at the s m a surface at the coordinate z is equal to:
F m a = m w γ m a = m w d f m a d t = + m w d μ d z
They undergo an acceleration of:
d f m a d t = d μ v a d z
The + sign of Equation (49) is negative in the literature but must be corrected as positive. In fact, the negative sign arises from the fact that the potential, μ, is taken negative in a standard way, in accordance with current thinking about potentials. However, we showed [11] that the chemical potential of the thermodynamic phases ( μ v a in this case) is always positive.
The products: f m a W m a and W m a d f m a d t are, therefore, respectively: the linear momentum, MLT−1 and the force of inertia, MLT−2, of the n ¯ m a molecules of mass m w , both refer to the local structural mass m s z = C / α z as we saw above.
Let us derive the linear momentum of molecules of water f m a W m a with respect to time:
d d t f m a W m a = W m a d f m a d t + f m a d W m a d t
By replacing d W m a d t by its equivalent W m a d f m a d z given by Equation (34) and d f m a d t by d μ m a d z (50), we obtain:
d d t f m a W m a = W m a d μ m a d z f m a W m a d f m a d z = W m a d μ m a d z d f m a 2 2 d z
We show in Appendix C that
d μ m a d z d f m a 2 2 d z = α t f m a
where α t is the constant of dimension t-1 of the time exponential in the relation:
W m a = A t = 0 exp α t t + B
Thus, according to (52), we get the important relationship:
d d t f m a W m a = α t f m a W m a
which leads to:
d ln W m a   f m a d t = d ln f m a d t + d ln W m a d t = α t = d d t ln W m a B
Furthermore, since d f m a W m a d z = 0 according to (48), the derivative of (52) with respect to z is null:
d 2 f m a W m a d z d t = d α t d z W m a f m a + α t W m a d f m a d z + α t f m a d W m a d z = 0 ,  
meaning that α t is constant with depth.

3.4.2. Equations of W m a and f m a and Their Derivatives

We deduce from the above the equation of W m a depending on t and z:
W m a W m a f = W ma z = 0 exp α z z W m a f                                                            =   W m a ° W m a f exp α t t exp α z z
where B = W m a f and A t = 0 = W m a ° W m a f in Equation (54), W m a f being the macro-water content of the soil surface at equilibrium with atmosphere, at the end of the experiment, and W m a ° being the macro-water content at saturation at the beginning of the experiment with the condition that W m a f W m a W m a ° .
Thus, we have:
d W m a d z = W ma z = 0 exp α z z = α z W m a W m a f
d W m a d t = α t W m a ° W m a f exp α t t exp α z z = α t W m a W m a f
d 2 W m a d z d t = α t α z W m a ° W m a f exp α t t exp α z z = d 2 W m a d t d z
From (56) we can calculate the derivatives of f m a as functions of W m a :
d ln f m a d t = α t d W m a d t W m a = α t α t W m a B W m a = α t B W m a
and, since d ln W m a d t = d f m a d z = α z f m a ,
d f m a d z = α z f m a = d ln W m a d t = d ln f m a d t α t = α t B W m a 1 ,
which gives:
f m a = α t α z B W m a W m a = d W m a / d t d W m a / d z = d z d t
Finally, retaking (62) and using (64), we obtain:
d f m a d t = α t B W m a f m a = α t 2 α z B W m a B W m a W m a

3.4.3. Application of the Equilibrium Equations between the Two Pedostructure Aqueous Phases

The pressure balance between the two aqueous phases which is established simultaneously with the migration of macro water to the surface can be seen as follows.
The molecular flux f m a of the macro phase is really the molecular speed of the molecules of this phase when they pass through the section s ¯ m a under the effect of a potential gradient, specific for this phase, determined at z by the equation:
d H d z = d h m a d z ρ w g = ρ w d μ m a d z g
The molecules of the micro phase (of potential μ m i ) that are found in the clay plasma of the primary aggregates are themselves subjected to the pressure difference h m a h m i , which appears between the two phases as soon as there is a change in macro water content W m a at z, i.e., as soon as a variation in the flux of molecules of this phase along the z axis appears ( d f m a / d z 0 ).
In this case, ( d f m a / d z 0 ) , the pressure balance between the two phases expressed by h m a = h m i is broken and must be re-established by a lateral flux of molecules from the micro phase to the macro phase.
We can then describe the process of water migration in the sample following the evaporation of water at the surface as follows: the variation of the flux of inter-aggregate water at depth z, d f m a d z 0 , has the effect of a change in water content W z at this same depth z, which simultaneously causes a new micro to macro flux totally determined by the equilibrium pressure equation h m a = h m i . This equality was studied above, giving rise to relations (10) to (12).
Moreover, starting from the equilibrium condition: h m i = h m a , we have every moment
μ m i e q μ m i S a t = μ m a e q μ m a S a t
which is written, according to the local values of W m i and W m a :
E ¯ m i W m i E ¯ m i W m i S a t = E ¯ m a W m a E ¯ m a W m a S a t
By setting the constant parameters of the shrinkage curve [6]:
A = μ m a S a t μ m i S a t = E ¯ m a W m a S a t E ¯ m i W m i S a t   and   E ¯ = E ¯ m i   +   E ¯ m a
we get the following equalities:
W m i E ¯ m i = W m a E ¯ m a A W m a = W E ¯ A W m a ;
showing that the ratios W m i W m a , W W m i and W W m a are all functions of W m a alone. Therefore, we have:
W = W m a E ¯ A W m a E ¯ m a A W m a   et   W m i = E ¯ m i W m a E ¯ m a A W m a  
Having the distribution of W m a in space and its variation with time, we automatically have the values and variations of W and W m i in any point of the medium.
The following equations can be verified:
d W m i d W m a = E ¯ m a E ¯ m i E ¯ m i 2 E ¯ m a A W m a 2 = E ¯ m a E ¯ m i W m i W m a 2
leading to:
d W d t = d W m a d t 1 + E ¯ m a E ¯ m i W m i W m a 2 = α t W m a W m a f 1 + E ¯ m a E ¯ m i W m i W m a 2
and
d W d z = d W m a d z 1 + E ¯ m a E ¯ m i W m i W m a 2 = α z W m a 1 + E ¯ m a E ¯ m i W m i W m a 2
The continuity Equation (39) becomes:
d W d t = ρ w V ¯ d f e w d z = W m a d f w d z = W m a d f m a d z 1 + E ¯ m a E ¯ m i W m i W m a 2

4. Results

4.1. Linear Relationships between Wz, Wt and Time

The characteristic retention curve of the sample h W is shown in Figure 5, which also shows the two curves ( h 1 and h 2 ) of the continuous reading of tensiometers T1 and T2 as a function of the total water content ( W t ) of the sample. The curves are homothetic: the total water content of the sample W t = M M s / M s corresponds to the values of the suction pressures h 1 and h 2 measured by the tensiometers and to the local water contents W 1 and W 2 that can be read on the retention curve h W , characteristic of the sample.
The following relationships are observed:
W z 1 = a 1 W t + b 1   and   W z 2 = a 2 W t + b 2
where a i and b i are constants associated with the depth z i . At a given depth, the local water content is in constant proportion to the overall water content of the sample. The same applies to the difference in water content between two depths:
W z 2 W z 1 = W t a 2 a 1 b 2 b 1
Let us find the relationship between a i and b i . At water saturation W 0 of the sample,
W z s a t = a i W t s a t + b i   and   W t s a t = W z s a t = W 0
thus,
b i = W 0 1 a i
and
a i = W z W 0 W t W 0
For a i to be time independent, whatever z, it is necessary that:
d a i d t = d W t d t W z W 0 W t W 0 2 + d W z d t W t W 0 W t W 0 2 = 0
d W t d t W z W 0 = d W z d t W t W 0
or else
d ln W t W 0 d t = d ln W z W 0 d t
The relation between the water content local, W Z , and total, W t , is such that:
d W z d t / d W t d t = W z W 0 W t W 0 = a i
W z being defined as the ratio of the areas s w = s m a + s m i and s s at z: W z = s m a + s m i s s z .
The 2 graphs in Figure 6 give the values of a i and of d W z d t at the two positions of the tensiometers z 1 and z 2 , which gives us, according to (83), d W t d t = d W z / d t a i = −8.56 × 10−5 min−1 and −8.57 × 10−5 min−1, respectively.

4.2. Logarithmic Relation between Wz and Wma-Wmaf

We can see in Figure 7 that the relation between W z and W m a is, for the two cases of z, a simple logarithmic function such as:
W z W c = ln W m a z W m a f + C
where C is a dimensionless constant and 1 / W c = α w is a constant parameter of the exponential of W z .
What is remarkable is that this logarithmic form of W z (84) exactly represents the Equation (71) of W z function of W m a z :
W z W c = ln W ma z W m a f + C = 1 W c W m a E ¯ A W m a E ¯ m a A W m a
By differentiating (84) with respect to time and using the relation (60) giving d W m a d t , we obtain:
W m a W m a f W c d W z d t = d W m a d t = α t W m a W m a f
We, therefore, have whatever z in the unsaturated zone:
d W z d t = α t W c = c t e
and according to the relation (73) that exists between d W z d t and d W m a d t :
W c = W m a W m a f 1 + E ¯ m a E ¯ m i W m i W m a 2 = W m a W m a f R z

4.3. The Fundamental Relationships between Flux, Water Potential and Water Content at Macroscopic Scale

4.3.1. Central Role of W ma

Figure 8 shows the experimental result of the relationship between the macro water content gradients Δ W m a Δ z and the pressure Δ h Δ z of a thin soil layer and the average water content of this layer.
We are in the case of a systemic discretization of the soil medium to apply transfer equations [9]. The values of W m a 1 and W m a 2 are calculated from the data of tensiometers T1 and T2 using the characteristic soil retention curve. They are supposed to represent the average water contents in the 1 cm thick slices around tensiometers 1 and 2. The products W m a Δ h ρ w Δ z , Δ h ρ w Δ z Δ W m a Δ z and Δ ln W m a Δ z named F1, F2 and F3 are presented as functions of W m a 1 , 2 , mean of W m a between z1 and z2.
We observe three simple exponentials of W m a , two of which are represented in Figure 8:
F 1 W m a Δ H ρ w Δ z = k 1 exp α 1 W m a ;   k 1   =   3.123   Jm 1   and   α 1   =   54.04   kg _ soil / kg _ water  
F 2 Δ W m a ρ w Δ z Δ H Δ z = k 2 exp α 2 W m a ;   k 2 =   14.629   Jm 2   and   α 2 =   107.4   kg _ soil / kg _ water
F 3 Δ ln W m a Δ z = k 3 exp α 3 W m a ;   k 3 =   4.721   m 1   and   α 3 =   53.61   kg _ soil / kg _ water
To interpret these results, recall that the basic variables, suction pressure (h), chemical potential ( μ m a ), and molecular flux of the mobile phase ( f m a ), are linked by the relation:
d h ρ w d z = d μ m a d z = d μ m i d z = d f m a d t
and knowing that
d ln μ m a = d ln W m a   et   W m a μ m a = E ¯ m a
we then have the following equalities:
W m a d h ρ w d z = W m a d μ m a d z = W m a d f m a d t = E ¯ m a d ln μ m a d z = E ¯ m a d ln W m a d z
Assuming that the discretization is fine enough to maintain at the macroscopic scale the relationship observed at the molecular scale between extensive and intensive variables, we should observe, after (89)–(91) and (94): F 1 = E ¯ m a F 3 and F 1 F 3 = F 2
F 1 F 3 = F 2 Δ W m a Δ z Δ h ρ w Δ z = k 1 2 E ¯ m a exp 2 α 1 W m a
F 3 Δ ln W m a Δ z = k 1 E ¯ m a exp α 1 W m a
We can notice the good accordance between the measured values of parameters and the theoretical relationships between them
k 2 = k 1 2 E ¯ m a = k 1 k 3   and   α 1 = α 3 = α 2 / 2

4.3.2. Pedostructure Water Conductivity k ps

Thus, we have all the physical equations determining the space–time relationship of variation of the three variables describing the dynamics of the medium: the fluxes, water contents and chemical potentials of the two aqueous phases.
Recall the equation of continuity (74) that takes account of the thermodynamic equilibrium. Using the relation f m a = d z d t , we can write:
f m a d W d z = d z d t d W d z = d W d t
so the equation of continuity can be written such as:
d W d t = d W m a d t 1 + E ¯ m i E ¯ m a μ m a μ m i 2 = W m a d f w d z = f m a d W d z = ρ w V ¯ d f e w d z
Furthermore, experience has shown that
d W d t = α t W c = c t e
W m a being an exponential of time and space, we deduce from the fact that W m a d f w d z = d W d t = c t e (Equations (98) and (99) that d f w d z and therefore also f w are simple exponentials with the same coefficients as W m a . Thus, as the soil water conductivity by definition is written:
k p s = f e w d h / d z = f w d h / d z f e w f w = f w d h / d z θ m a ,
by multiplying (101) up and down by α z W m a = d W m a d z (47) we get:
k p s = α z W m a f w α z W m a d h / d z θ m a = W m a d f w / d z W m a / d z d h / d z θ m a = d W / d t F 2 θ m a
and using Equation (100):
k p s = α t W c k 2 θ m a exp α 2 W m a
The constants: α t ,   W c ,   k 2 = k 1 2 E ¯ m a and α 2 = 2 α 1 , are all determined by measurement as we showed above (97).

5. Discussion

The systemic modeling of the hydrostructural soil water properties by the model Kamel [9,12] already accounted for the levels of internal organization of the “soil factory”. It precisely identified the pedostructure as assembly of primary peds containing two thermodynamically distinct aqueous phases, intra-aggregate (macro) phase and inter-aggregate (micro) phase. However, this modeling still ignored the lower levels of organization (molecular and atomic) of the fluid phases of the pedostructure. It retained, therefore, a semi-empirical character because it is at these two levels of organization that the variables of temperature, pressure, entropy and chemical potential have their basis of definition, as our previous study [11] showed. Furthermore, all these variables intervene in the hydro-thermodynamic equilibria of the soil at higher levels of organization (soil suction gradient, soil water retention curve, hydric conductivity, etc.). By understanding the internal organization of the aqueous phases and their roles in the evaporation process, we were able to relate the variables at each level to each the others in a comprehensive and orderly manner.
The present study focused on the notions of flux: molecular flux (fma), water flux (fw, fma and fmi,) and Eulerian flux (few). The fma, flux is said to be molecular because it is equal to the celerity of the molecules of the phase and, therefore, linked to the chemical potential of the phase as we have shown it. Application of the Newton’s second law makes it possible to identify the gradient of the chemical potential of the inter-aggregates aqueous phase (macro) to the time derivative of its molecular flux, leading to the relation: d μ m a d z = d f m a d t , which has the dimensions of an acceleration (LT−2). This important relationship could not be obtained without the acknowledgement of these molecular and atomic organization levels. This allows us to say that we have solved the Navier–Stokes equation for the particular case of the water flux in pedostructure during its drying by evaporation at its upper surface.
Above this molecular level, there are the nested levels of organization that we have already dealt with exhaustively in hydrostructural pedology [1,2]. The present study has defined and highlighted the junction point of both worlds by studying the molecular and non-molecular descriptive variables (intensive and extensive) attached to the z-depth where each descriptive variable can be considered punctual and continuous.
Finally, the equation of water conductivity (103) is given here theoretically for the first time from a combination of hydro-thermodynamic variables. It has the same form as the one we found semi-empirically and published recently in [13]. This result solves one of the many points of soil science questioned [14] and which was not yet fully resolved by hydrostructural pedology; in particular, the physical equation of the water conductivity of unsaturated soils and its measurement [15]. The resolution of this last point confirms the validity of the “systemic approach” such as it was redefined by the authors [1,2,3] to face the “black-box” modelling in soil hydrodynamics and thermodynamics. It also justifies the new methods of characterization of the water functions of the soil [16,17,18].

6. Conclusions

A long-standing problem in fundamental physics has been solved and experimentally confirmed in this article: that of the dynamics of water in unsaturated soil. The exact physical equation for non-saturated hydric conductivity has been found; this by reconsidering the equations of Euler and Richards with a new approach, systemic and molecular, of thermodynamics. The methodology for obtaining parameters of the pedostructure hydric conductivity curve of soil is also given in the article, valid for all types of soil, swelling or not. At the same time, we discovered that the link between the two domains of the soil science: soil hydrodynamics and soil water thermodynamics is provided by Newton′s 2nd law, which links the acceleration of molecules of the macro aqueous phase to the gradient of the chemical potential of this phase.
These results confirm the validity of two recent theories in environmental science thanks to which they were obtained: the systemic approach which gave birth to hydrostructural pedology and the systemic and molecular thermodynamics of soil water.

Author Contributions

Conceptualization, writing—original draft preparation, E.B.; writing—review and editing, R.H.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. Field sampling and laboratory data were acquired in 2003–2005 in Martinique by the first author in the frame of its research work at IRD.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to non-prepared presentation.

Acknowledgments

The authors warmly thank Hassan Boukcim (agro-biologist, CEO of Valorhiz) and Amjad T. Assi, research assistant professor at Texas A&M University for their long and constant support for the applied aspects of this work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. The Richards Equation, an Eulerian Point of View

The Richards equation is the combination of two physical equations of different conceptual origins:
(a) conservation law of the matter which is expressed mathematically by the Euler equation
d θ w d t = d f e w d z ,
where θ w is the volume water content of the soil and where f e w is the Eulerian flow, understood as the volume flow velocity on z (upward positive);
(b) Darcy′s law extended to unsaturated porous media, the classical expression of which is
f e w = k Δ H / Δ z
where Δ H Δ z = Δ h Δ z ρ w g Δ z Δ z is the pressure gradient of the soil water with respect to z, sum of the water retention pressure gradients and the pressure gravity.
This relationship is called “Darcy′s law extended to unsaturated soils” and its application to the Euler equation constitutes the Richards equation used in all soil-water models of water dynamics in soil:
d θ w d t = z K d H d z
The question that arises when one goes from the verified Darcy law for saturated soils to the extension of the law for unsaturated soils concerns the surface area of the flux at depth z. We generally take s t the total surface of the sample but would it not be s w , the surface occupied by the water molecules, which decreases with the water content, or else only s m a , the surface corresponding to the inter-aggregate of water, which should be considered? The systems approach helps answer this question as we will see below. In fact, the variable θ w , which is the ratio of the volume of water to the volume of soil θ w = V w / V t , is a non-systemic variable defined according to the notion of REV (Representative Elementary Volume), which implies that the variable f e w cannot be defined by Equation (8): we do not know what it is spatially represents. The notions of surfaces occupied by water molecules and solids at the depth z of the cylinder and of common thickness of the soil layer do not appear in the concept of REV.
The so-called Richards Equation (A3) is therefore empirical and cannot be related to fundamental physics. To write it correctly, we must first understand the exact physical meaning of the flux variables above, f w et f e w , but also that of the thermodynamic variables that describe the energy equilibria that are established at the molecular scale between the molecules of each fluid phase as well as at the global scale between the thermodynamic phases of the pedostructure. Indeed, these thermodynamic variables, like chemical potential, temperature and entropy, take on their full physical significance at the molecular scale as we have shown in [11]. We have to take this level of molecular scale into account to describe the process of migration of water molecules in the soil (its pedostructure) submitted to a potential gradient at thermodynamic equilibrium.

Appendix A.2. Euler′s Equation for Conservation of the Mass

As mentioned above, the Richards equation comes from the introduction of Darcy′s law “extended to unsaturated media” into the Euler continuity equation. We need to study the terms, including the meaning of flow, to relate it exactly to the speed of molecules and the thermodynamic variables that were defined at the molecular level in the first part.
Euler′s equation for the conservation of the mass of water moving in the soil is written:
ρ t + ρ w v = 0
where ρ is the volume concentration of water in the soil (mass of water in the volume V of soil: ρ = M w V ), t is time and v the Eulerian velocity of the fluid (water) with density ρ w assumed constant. The Eulerian speed can be identified, except for the sign, with a flow that we will call Eulerian flux and write: f e w ( ρ w = M w V w ,   v = f e w = ± d l w d t ).
Let us rewrite this Equation (A4) using the systemic variables ( W and V ¯ = V / m s ) rather than the non-systemic variables ( ρ   and   θ = V w / V ) in order to discern the internal process variables involved in the migration of water in the sample subjected to evaporation. Let us first recall the existing relationships between these types of variables, m w and m s being the mass of water and solids contained in the volume V of soil (pedostructure):
ρ = m w V = ρ w V w V = ρ w θ = m w / m s V / m s = W V ¯ = ρ w V ¯ w V ¯
d ρ d t = d W /   V ¯ d t = ρ w d θ d t
since ρ w = c t e at constant temperature and pressure.
The second term is written such as:
ρ w v = ρ w v = ρ w d f e w d z + d f e w d y + d f e w d x
where f e w = v is the Eulerian flow which has the dimensions of a velocity, in LT-1.
Thus, the Euler equation that is known in hydrology is:
d θ d t = d f e w d z

Appendix B

Appendix B.1. Equations of W , f w and Their Derivatives According to W m a and f m a

Consider now the products f w W and f e w ρ w V ¯ ; we have, according to the definition of the derivatives of fluxes (39–41):
d f w d z = 1 W m a d W d t
and
d f e w d z = 1 ρ w V ¯ d W d t
We then obtain the general equation, similar to the Euler equation written with the systemic variables:
d W d t = ρ w V ¯ d f e w d z = W m a d f w d z
Moreover, using the relation f m a = d z d t we can write:
f m a d W d z = d z d t d W d z = d W d t
so we have
W m a d f w d z = f m a d W d z = d W d t
that we can compare to
W m a d f m a d z = f m a d W m a d z = d W m a d t
By subtracting the two equations term by term, we obtain:
W m a d f w d z d f m a d z = f m a d W d z d W m a d z = d W d t d W m a d t
equivalent to:
W m a d f w d z d f m a d z = f m a d W m i d z = d W m i d t
Defining f m i = f w f m a as the virtual speed of molecules of the micro phase at z, we have the following fundamental relationships:
W m a d f m i d z = f m a d W m i d z = d W m i d t
We find here the central role of f m a and W m a :
f m a = d W m i / d t d W m i / d z = d W m a / d t d W m a / d z = d W / d t d W / d z
and
W m a = d W m i / d t d f m i / d z = d W m a / d t d f m a / d z = d W / d t d f w / d z

Appendix C

Appendix C.1. Application of Newton’s Law, d f m a W m a d t = α t f m a W m a Demonstration

Let us try to determine the relation between the two terms d μ m a d z and d f m a 2 2 d z of relation (52). The product f m a W m a written with its fundamental variables is equal to d W m a α z d t from Equation (40). The derivative with respect to time is, therefore, the second derivative of W m a :
d f m a W m a d t = d 2 W m a α z d t 2
From Equations (52) and (A19) we therefore have:
d 2 W m a α z dt 2 = d f m a W m a d t = d f m a d t W m a d W m a d t f m a = W m a d μ m a d z d f m a 2 2 d z
Suppose that d W m a d t is an exponential function of time, in its most general form, as will be verified experimentally:
W m a = A t = 0 exp α t t + B   and   d 2 W m a d t 2 = α t 2 A t = 0 exp α t t = α t d W m a d t = α t 2 W m a B
Equation (A20) is then written:
d 2 W m a α z d t 2 + d W m a d t f m a = d f m a d t W m a
which, according to (55) and (50), can be put in the form:
d W m a d t α t α z + f m a = d f m a d t W m a = d μ m a d z W m a
By dividing all the members of Equation (A23) by W m a and using Equation (34), we obtain the relation which links together d μ m a d z , d f m a d t and d f m a d z :
d f m a d t = d μ m a d z = d ln W m a d t α t α z + f m a = d f m a d z α t α z + f m a
Using the relation d f m a d z = α z f m a (45), (A24) becomes:
d f m a d t = d μ m a d z = α t f m a + 1 2 d f m a 2 d z
From where we get
d μ m a d z 1 2 d f m a 2 d z = α t f m a
which, reported in Equation (A20), gives:
d d t f m a W m a = α t W m a f m a

References

  1. Braudeau, E. Modélisation systémique du couple eau-sol: La pédologie hydrostructurale. Res-Systemica 2015, 13, 94. [Google Scholar]
  2. Braudeau, E.; Assi, A.T.; Mohtar, R.H. Emergence of a new scientific discipline: Hydrostructural pedology. In Hydrostructural Pedology; ISTE: London, UK; Wiley: Hoboken, NJ, USA, 2016; ISBN 978-1-84821-994-6. [Google Scholar]
  3. Braudeau, E.; Mohtar, R.H. Modeling the soil system: Bridging the gap between pedology and soil-water physics. Glob. Planet. Chang. J. 2009, 67, 51–61. [Google Scholar] [CrossRef]
  4. Pouvreau, D. Une Histoire De La “Systémologie Générale” De Ludwig Von Bertalanffy—Généalogie, Genèse, Actualisation Et Postérité D’un Projet Herméneutique, Thèse, Spécialité “Histoire Des Sciences”, Ecole Des Hautes Etudes En Sciences Sociales (Ehess), 2013, Paris. Available online: http://tel.archives-ouvertes.fr/tel-00804157 (accessed on 22 January 2021).
  5. Le Moigne, J.-L. La Théorie du Système Général, Théorie de la Modélisation, 4th ed.; PUF: Paris, France, 1994. [Google Scholar]
  6. Braudeau, E.; Assi, A.T.; Boukcim, H.; Mohtar, R.H. Thermodynamic formulation of the pedostructure water retention and shrinkage curves. Front. Environ. Sci. 2014, 2. [Google Scholar] [CrossRef] [Green Version]
  7. Assi, A.T.; Accola, J.; Hovhannissian, G.; Mohtar, R.H.; Braudeau, E. Pedostructure characterization through measurement and modeling of the soil moisture characteristic curves. Front. Environ. Sci. 2014, 2, 5. [Google Scholar] [CrossRef] [Green Version]
  8. Braudeau, E.; Hovahnnissian, G.; Assi, A.T.; Mohtar, R.H. Soil water thermodynamic to unify water retention curve by pressure plates and tensiometer. Front. Environ. Sci. 2014, 2. [Google Scholar] [CrossRef] [Green Version]
  9. Braudeau, E.; Mohtar, R.H. A framework for soil-water modeling using the pedostructure and structural representative elementary volume (SREV) concepts. Front. Environ. Sci. 2014, 2, 24. [Google Scholar] [CrossRef] [Green Version]
  10. Braudeau, E.; Boukcim, H.; Assi, A.T.; Mohtar, R.H. Hydrostructural pedology, new scientific discipline allowing for physical modelling of ‘Green Water’ dynamics in the soil-plant-atmosphere system. J. Agric. Ecol. Res. Int. 2018, 15, 1–8. [Google Scholar] [CrossRef] [Green Version]
  11. Braudeau, E. Systemic modelling of soil water thermodynamics under natural conditions of air temperature and pressure. Eur. J. Appl. Sci. 2020, 8, 46–72. [Google Scholar] [CrossRef]
  12. Braudeau, E.; Mohtar, R.H.; El Ghezal, N.; Crayol, M.; Salahat, M.; Martin, P. A multi-scale “soil water structure” model based on the pédostructure concept. Hydrol. Earth Syst. Sci. Discuss. 2009, 6, 1111–1163. [Google Scholar]
  13. Braudeau, E.; Assi, A.T.; Mohtar, R.H. Methods for determining the characteristic parameters. In Hydrostructural Pedology; ISTE: London, UK; Wiley: Hoboken, NJ, USA, 2016; ISBN 978-1-84821-994-6. [Google Scholar]
  14. Hunt, A.G.; Ewing, R.P.; Horton, R. What’s wrong with soil physics? Soil Sci. Soc. Am. J. 2013, 77, 1877–1887. [Google Scholar] [CrossRef]
  15. Schelle, H.; Heise, L.; Janicke, K.; Durner, W. Water retention characteristics of soils over the whole moisture range: A comparison of laboratory methods. Eur. J. Soil Sci. 2013, 64, 814–821. [Google Scholar] [CrossRef]
  16. Blake, J.; Assi, A.T.; Mohtar, R.H.; Braudeau, E.; Morgan, C.L.S. Quantification of available water capacity comparing standard methods and a pedostructure method on a weekly structured soil. Trans. ASABE 2019, 62, 289–301. [Google Scholar] [CrossRef]
  17. Assi, A.T.; Mohtar, R.H.; Braudeau, E. Soil pedostructure-based method for calculating the soil water-holding properties. MethodsX 2018, 5, 950–958. [Google Scholar] [CrossRef] [PubMed]
  18. Assi, A.T.; Blake, J.; Mohtar, R.H.; Braudeau, E. Soil aggregates structure-based approach for quantifying the field capacity, permanent wilting point and available water capacity. Irrig. Sci. 2019, 37, 511–522. [Google Scholar] [CrossRef]
Figure 2. Representation of the internal organization of the pedon, hierarchized into its hydro-functional levels of organization: horizons, pedostructure, primary peds (after Braudeau et al. 2009 [12]).
Figure 2. Representation of the internal organization of the pedon, hierarchized into its hydro-functional levels of organization: horizons, pedostructure, primary peds (after Braudeau et al. 2009 [12]).
Systems 09 00008 g002
Figure 3. Conceptual model of the thermodynamic equilibrium at molecule scale; v is the celerity of the molecule and ν 2 is the number of shocks per unit of time [11]. The figure represents one molecule (in red) of volume V m , in its occupational volume, V m t , of diameter D 2 (red + blue). The external volume of the molecule, in blue, is V m v = V m t V m , the volume of its external energy: E m v =   T S 1   = 1 2 m v 2 . In a fluid phase, the volume concentration of energy, the pressure, is the same everywhere, in atoms, molecule and the intermolecular space.
Figure 3. Conceptual model of the thermodynamic equilibrium at molecule scale; v is the celerity of the molecule and ν 2 is the number of shocks per unit of time [11]. The figure represents one molecule (in red) of volume V m , in its occupational volume, V m t , of diameter D 2 (red + blue). The external volume of the molecule, in blue, is V m v = V m t V m , the volume of its external energy: E m v =   T S 1   = 1 2 m v 2 . In a fluid phase, the volume concentration of energy, the pressure, is the same everywhere, in atoms, molecule and the intermolecular space.
Systems 09 00008 g003
Figure 4. Device for measuring the hydric conductivity of the pedostructure.
Figure 4. Device for measuring the hydric conductivity of the pedostructure.
Systems 09 00008 g004
Figure 5. Representation on the same graph of the water retention characteristic curve h W , measured and modelled, and of the tensiometric reading h 1 and h 2 according to the total water content W t . Modelled curves of h W , h 1 W t and h 2 W t used the same pedostructural parameters characteristic of the sample: W m a S a t ,   W m i S a t ,   E ¯ m a   and   E ¯ m i .
Figure 5. Representation on the same graph of the water retention characteristic curve h W , measured and modelled, and of the tensiometric reading h 1 and h 2 according to the total water content W t . Modelled curves of h W , h 1 W t and h 2 W t used the same pedostructural parameters characteristic of the sample: W m a S a t ,   W m i S a t ,   E ¯ m a   and   E ¯ m i .
Systems 09 00008 g005
Figure 6. Linear relationships between Wz, Wt and Time. (a) represents the local water content at depths z1 and z2 according to the global water content of the sample; (b) represents the time dependence of the two local water contents (at depths z1 and z2) and the global water content. The z-area ratio a i can be read on the figure (a): 2.29 for z = z1 and 1.82 for z = z2.
Figure 6. Linear relationships between Wz, Wt and Time. (a) represents the local water content at depths z1 and z2 according to the global water content of the sample; (b) represents the time dependence of the two local water contents (at depths z1 and z2) and the global water content. The z-area ratio a i can be read on the figure (a): 2.29 for z = z1 and 1.82 for z = z2.
Systems 09 00008 g006
Figure 7. Relationship between W t , W z and W m a z W m a f . Rivière Lézarde Ponterre (halloysite soil) 60–65 cm. (a) represents the two local water contents ( W z and W m a z ) relationship at two given depths (b) represents the global water content of the sample as a function of the local macro water content at a given depth z. The constant W c in Equation (85) is read on the left figure is 0.04995 and 0.0500 kg of water/kg of soil.
Figure 7. Relationship between W t , W z and W m a z W m a f . Rivière Lézarde Ponterre (halloysite soil) 60–65 cm. (a) represents the two local water contents ( W z and W m a z ) relationship at two given depths (b) represents the global water content of the sample as a function of the local macro water content at a given depth z. The constant W c in Equation (85) is read on the left figure is 0.04995 and 0.0500 kg of water/kg of soil.
Systems 09 00008 g007
Figure 8. Measurement result of functions F1 and F2 defined hereafter in the text, which are simple exponentials of W m a : (a) represents F1 as function of W m a and (b) F2 as function of W m a .
Figure 8. Measurement result of functions F1 and F2 defined hereafter in the text, which are simple exponentials of W m a : (a) represents F1 as function of W m a and (b) F2 as function of W m a .
Systems 09 00008 g008
Table 1. Pedostructure state variables. Subscripts mi, ma and s; refer to micro, macro and solids; ip, st, bs and re, refer to the name of the corresponding shrinkage phase of the shrinkage curve: interpedal, structural, basic and residual [12].
Table 1. Pedostructure state variables. Subscripts mi, ma and s; refer to micro, macro and solids; ip, st, bs and re, refer to the name of the corresponding shrinkage phase of the shrinkage curve: interpedal, structural, basic and residual [12].
Volume of ConcernSpecific Volume
[dm3/kg]
Specific Pore Volume
[dm3/kg]
Specific Water Content
[kgwater/kgsoil]
Non Saturating Water
[kgwater/kgsoil]
Saturating Water
[kgwater/kgsoil]
Suction
[kPa]
Pedostructure V ¯ V p ¯ W h
Interpedal porosity V p ¯ m a W m a w s t w i p h m a ,   h i p
Primary peds V ¯ m i V p ¯ m i W m i w r e w b s h m i
Primary particles V ¯ s
Table 2. Primary thermodynamic variables and their units.
Table 2. Primary thermodynamic variables and their units.
VariablesSymbolsEquationTypeUnits
Total molecular energy E m t 1 2 m λ ν 2 + m v 2 Energyjoule
Temperature T = E m 1 2 m λ ν 2 Energyjoule
Entropy SS V m v / V m inter-mol. Vol./molecular Vol.number
Pressure P P = P m = P m v = P m t T V m ,   T S 1 V m v ,   T S V m t Energy/Volumejoule/m3; Pa
Chemical potential µ μ m ,   μ m v , μ α T m ,   T S 1 m ,   T S m Energy/Massjoule/kg
Temperature is pseudo-intensive, the others are true intensive variables and they represent each point of the medium and stay unchanged across scales [11].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Braudeau, E.; Mohtar, R.H. Hydrostructural Pedology, Culmination of the Systemic Approach of the Natural Environment. Systems 2021, 9, 8. https://doi.org/10.3390/systems9010008

AMA Style

Braudeau E, Mohtar RH. Hydrostructural Pedology, Culmination of the Systemic Approach of the Natural Environment. Systems. 2021; 9(1):8. https://doi.org/10.3390/systems9010008

Chicago/Turabian Style

Braudeau, Erik, and Rabi H. Mohtar. 2021. "Hydrostructural Pedology, Culmination of the Systemic Approach of the Natural Environment" Systems 9, no. 1: 8. https://doi.org/10.3390/systems9010008

APA Style

Braudeau, E., & Mohtar, R. H. (2021). Hydrostructural Pedology, Culmination of the Systemic Approach of the Natural Environment. Systems, 9(1), 8. https://doi.org/10.3390/systems9010008

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop