Hydrostructural pedology, culmination of the systemic approach of the natural environment

: 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, validating at the same time the systemic approach which was the framework.


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: 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 (ℎ ( ) = ℎ ( )) and distinguished by their thermodynamic properties. Figure 1: Place of hydrostructural pedology among the earth sciences [10] 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, could be physically redefined and precisely explained because depending on the two levels of organization: molecular and atomic.
In the present article, we will introduce these two levels of organization in the explanation of the terms of the Richards equation (Appendix 1), after having rewritten it in the systemic framework by taking into account the volume variable (volume of the organized system and subsystems). The problem is knowing how to associate in the same equation an extensive variable (volume) and an intensive variable (potential), empirically done today with Darcy's law extended to unsaturated soils. 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 / external energy, internal / external pressure, internal / external chemical potential are in equilibrium relationship, as is well explained in the previous article.

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. Figure 2: Representation of the internal organization of the pedon, hierarchized into its hydro-functional levels of organization: horizons, pedostructure, primary peds (after   [12]).
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. 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].
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, ( ), the water retention curve, ℎ( ), the unsaturated soil water conductivity, k(W) and the swelling curve of primary aggregates as a function of time, ( ) [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 2mm sieved soil from the fractionation of a moderately dry soil sample (<pF3) on a 2mm sieve (can be 4mm when the sample is very clayey with swelling clay). The fine earth is added layer by layer in a cylinder of 5cm 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.

Volume of concern
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).

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 ( ) and the soil water retention curve ℎ( ), 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 ℎ( ) at equilibrium of the two aqueous phases of the pedostructure [8] and, from this, the exact distribution of the two kinds of water content ( and ) in the pedostructure as function of , raises a new and important question about the descriptive variables of the model. This equation ℎ( ) links an extensive variable (water content) to an intensive variable (water suction). Indeed, we use a mini tensiometer (2mm thick) planted laterally in a soil sample at depths z, to simultaneously obtain water suction (ℎ = ℎ = ℎ ), micro and macro water contents ( and ) at thermodynamic equilibrium, locally in the same soil mini-layer. To answer the fundamental question of what exactly is the spatial extension of and corresponding to ℎ and ℎ 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. 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 , 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: 2 ): the "shocks" or meetings of molecules of the fluid phase between them. We can then associate each molecule with an occupancy volume: , 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 ( ) comes in.
is a fractional number equal to Moreover, given that each molecule necessarily has the same chemical potential = ⁄ 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 ( = ⁄ ) and a chemical potential ( = ) that depends on the entropy ( ) 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, , (index void): 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 , an intensive state variable having the same value throughout the phase. At phase equilibrium, since the internal chemical potential of molecules, , is the same everywhere, the overall molecular chemical potential of the phase, = , is therefore an intensive variable characteristic of the phase. It is the same for the external chemical ). 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: where 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 which is equal to the internal molecular pressure ⁄ and the total molecular pressure in the phase: ( )⁄ . However, the suction pressure measured by the tensiometer in soil science, as shown previously [11], has the expression: where ° is the chemical potential of free water under air pressure and standard temperature. According to the equation (3), and because is equal everywhere in all phases of the system at thermodynamic equilibrium, we can substitute the chemical potential of the water ( ) by the intermolecular chemical potential ( ) without changing the value of the pressure ℎ.
This allows to be identified with the pressure potential of the water in the tensiometer, relating ℎ 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 ( ) and not the temperature ( ) 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 of the fluid phases of the system in equilibrium ( = ). The temperature of water molecules, for example, is equal to: 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 below the exhaustive list of hydro-thermodynamic variables qualified as systemic and which cover the 4 levels of organization: macroscopic and microscopic of the aqueous phase then molecular and atomic of the phase. We know that at the macroscopic level of the phases of the pedostructure, the free energies of the two aqueous phases of the pedostructure ( ̅ and ̅ ) are observed constant despite a change in the water content of this phase in the defined system of the pedostructure. Following Sposito [13] ̅ and ̅ was identified before as the free energy ̅ = and ̅ = , and being the water contents micro and macro of the pedostructure. Now, following our previous study [11] which differentiates the intermolecular energy ( ( − 1)) from the total energy ( ), ̅ and ̅ are defined as the intermolecular energy, corresponding to kinetic energy of molecules: This relationship involves the number of water molecules ̅ in the aqueous phase " " (macro or micro), whose molecular mass is . 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: and . The fact that the energy of the phase α (macro or micro), ̅ , is constant despite a change in phase's water content ( ̅ ) or in its chemical potential ( ), 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: or, for each α aqueous phase (macro or micro): 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: This gives the following general equation, since −°= −° and according to (10): Note that this equation is valid only if the saturated state corresponding to h = 0 is set and therefore the potential ° and ° corresponds to °= and °= and that all these equilibrium equations are deduced from the fact that the "free" energies ̅ and ̅ of the two aqueous phases in equilibrium, defined and expressed by the fundamental expressions (7) and (8), are constant and characteristic of the pedostructure.

Definition of the molecular flux 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, automatically an acceleration field derives from the potential and is applied to the molecules placed in this gradient; in the case of the macro phase: 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: ⃗ = ⃗ . This force determines the elementary pressure carried by each molecule of the phase on its environment: ( = ; cf. [11]) which causes a molecular flux through the surface ( ) (occupied by the macro aqueous phase) of the sample section at depth z, such as: where 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 which 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.

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.

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.

Measuring apparatus of the hydric conductivityof 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 2mm), positioned 1cm 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 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 ( ), two aqueous phases ( ) and ( ), and the gas phase ( ) such that the total surface of the section is the sum = + + + . 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 = all contain the same quantity of structural mass. where is a fixed characteristic of the homogeneous sample in terms of structure and its structural mass ; is the actual density of the solid phase, is the total mass of the dry sample and is its length.

T1
The extensive variables, such as water contents and organized volumes and sub-volumes ( , , , ̅ , etc.), are all related to the mass of the local structural phase: = , the mass contained in the same volume as that in which the other variables are defined. We have for example: 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.

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 ( , , etc.) or a discretized subsystem between two depths ( 1− 2 = ∆ 1−2 ) or that considered at the molecular level, that of the horizontal section at depth z. 1 °) The global variables , , 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 ( ) divided by the total mass of the solid phase ( ) constituting the structure of the system, is the mass volume, which we write as ̅ : ̅ = ⁄ . Likewise, the overall water content of the system, , is the mass water content of the system ( ) divided by the structural mass ( ) and written as: = ⁄ . These variables, all related to the same structural mass, are additive: = + . We have the following equalities: 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: ̅ = ( ). 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 ∆ 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 ̅ , using all the descriptive variables of the aqueous and gas phases defined at this 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 = = for a cylindrical sample of height L and structural mass . We define the water content variables in the same way: Note that, since and are constants, any ratio of two occupied surfaces of a section in z ( ⁄ 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: where is the area of the section at z occupied by water ( = + ) and is the total area, at z.

Flux variables at depth z
Just as we have defined the extensive variables at depth z ( ̅ ; ; ; ), we must also define the associated types of fluxes at depth z: where the ratio represents, as in equation (29), the speed of movement of water molecules on the z axis. 3 °) The molecular flux ( ), the real speed of the water molecules of the surface , is the speed (modulus) of the molecules of the macro phase determined by the chemical potential of the phase different from related to the number of micro and macro water molecules leaving the surface during the time dt): Indeed, the speed 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 and therefore directly linked to the state variables of thermodynamic equilibrium of the phase (temperature, chemical potential, entropy etc.)

Concomitant variation of organizational and fluxes variables at z
The two water fluxes, and , are, according to their physical definition (30) and (31), proportional to the molecular flux according to: or, dividing by ̅ : 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 defined by equation (31) gives: The change in sign results from the fact that one passes from a variation of speed, , to a variation of volume calculated with the height taken in the opposite direction of the speed. In the same way, we also have: 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: and we can rewrite the continuity equation in this form: Other relationships between fluxes and water contents are given in Appendix B.

Spatial variation of the product
We saw that is a molecular flux of the aqueous phase macro, and 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 as the speed of each molecule and as the concentration of momentum whose derivative with respect to time is a force.
Consider the molecular expression of the product: , in accordance with equations (24) to (27)  We give for the following physical definition: where is the solid mass element corresponding to and such that: being the element (mass) of the solid phase associated with and such that: Comparing equations (39)

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.

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 ( ) 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 , which drives the water molecules, of molecular mass: , present at the surface at the coordinate z is equal to: They undergo an acceleration of: 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 ( in this case) is always positive.
The products: ( ) and ( ) are therefore respectively: the linear momentum, MLT -1 and the force of inertia, MLT -2 , of the ̅ molecules of mass , both refer to the local structural mass = / as we saw above. Let us derive the linear momentum of molecules of water ( ) with respect to time: By replacing by its equivalent − given by equation (34) and by (50), we obtain: We show in Appendix C that 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 of the macro phase is really the molecular speed of the molecules of this phase when they pass through the section ̅ under the effect of a potential gradient, specific for this phase, determined at z by the equation: The molecules of the micro phase (of potential ) which are found in the clay plasma of the primary aggregates are themselves subjected to the pressure difference (ℎ − ℎ ) which appears between the two phases as soon as there is a change in macro water content at z, i.e. as soon as a variation in the flux of molecules of this phase along the z axis appears ( ⁄ ≠ 0). In this case, ⁄ ≠ 0), the pressure balance between the two phases expressed by ℎ = ℎ is broken and must be reestablished 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, ≠ 0, has the effect of a change in water content at this same depth z, which simultaneously causes a new micro to macro flux totally determined by the equilibrium pressure equation ℎ = ℎ . This equality was studied above, giving rise to relations (10) to (12).
Moreover, starting from the equilibrium condition: ℎ = ℎ , we have every moment The continuity equation (39) becomes: The characteristic retention curve of the sample ℎ( ) is shown in Figure 5, which also shows the two curves (ℎ 1 and ℎ 2 ) of the continuous reading of tensiometers T1 and T2 as a function of the total water content ( ) of the sample. The curves are homothetic: the total water content of the sample ( = ( − )/ ) corresponds to the values of the suction pressures ℎ 1 and ℎ 2 measured by the tensiometers and to the local water contents 1 and 2 that can be read on the retention curve ℎ( ), characteristic of the sample.

Linear relationships between Wz, Wt and Time
. The following relationships are observed: z1 = 1 + 1 and z2 = 2 + 2 75 where and are constants associated with the depth . 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: The 2 graphs in Figure 6 give the values of and of at the two positions of the tensiometers 1 and 2 , which gives us, according to (83), = ⁄ = -8.56 and -8.57.10-5 minutes-1 respectively.

Logarithmic relation between Wz and Wma-Wmaf
We can see in Figure 7 that the relation between and is, for the two cases of z, a simple logarithmic function such as: where C is a dimensionless constant and 1⁄ = is a constant parameter of the exponential of .
What is remarkable is that this logarithmic form of (84) exactly represents the equation ( (1,2) , mean of between z1 and z2. We observe three simple exponentials of , two of which are represented in Figure 8: The constants: , , 2 = − 1 2 ̅ and 2 = 2 1 , are all determined by measurement as we showed above (97).

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 therefore retained 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.. And 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 …). 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 so, 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: = 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 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].

Conclusion
The results of this article confirm the validity of the assertion made ten years earlier that a radical paradigm shift was necessary in soil science to methodologically link the cartographic description of soils in pedology to the physical modeling of the water mobility in the soil. The paradigm shift consisted of recognition of organized volumes as state variables in the same way as a water or air content of a given organized volume. A physical definition of a system follows, involving these nested organizational volumes, and a complete theorisation of the systemic approach of the natural environment, anthropized or not. This opened the door to the development of a new pedology, hydrostructural pedology, which includes a new type of modeling: systemic modeling of the natural environment and a new theory of thermodynamics of soil water. It is mainly thanks to the latter that a long-standing problem in fundamental physics has been solved, experimentally confirmed in this article: that of the dynamics of water in the soil medium.
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.

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: where is the volume concentration of water in the soil (mass of water in the volume V of soil: = ), t is time and the Eulerian velocity of the fluid (water) with density assumed constant. The Eulerian speed can be identified, except for the sign, with a flow that we will call Eulerian flux and